北京市城市边界提取流程.pptxVIP

北京市城市边界提取流程

目录CONTENTS1数据获取及预处理2北京城市边界提取3提取结果分析

一、数据获取及预处理1.研究区遥感影像获取：将北京市行政区划边界导入GoogleEarthEngine平台，编写代码进行去云等操作获取2019年及2023年北京市哨兵二号卫星遥感影像数据Level-2A（15m分辨率）。2.研究区遥感影像预处理：由于研究区范围较大，GoogleEarthEngine平台生成的遥感影像自动分块。因此将获取的遥感影像导入ArcGIS，进行拼接预处理。

一、数据获取及预处理在影像获取过程中发现，gee平台基于1年时间跨度的影像镶嵌生成的单独的B8、B11波段影像存在色差。分析原因：在遥感数据处理中，研究区不同范围拍摄影像可能因为云量、无数据等原因无法选择时间特别靠近的影像。而不同时间下由于气候等原因，地表辐射是不同的，这就造成了多景影像镶嵌后存在色差。解决：为改善色差和平衡云量，选择满足云量小于5的影像，时间跨度由原来的1年缩小至4个月（5-9月）。

二、北京市城市边界提取1.在arcmap中加载哨兵第8通道（近红外）和第11通道（中红外）的数据；2.利用arcmap中toolboxs中的栅格计算器输入公式（R11-R8）/（R11+R8）计算NDBI值；3.双击ndbi数据，点击符号系统，选择已分类，分为两类，以0值为断点；

二、北京市城市边界提取4.打开栅格计算器输入SetNull(“NDBI”0,1)，将NDBI0的像元设为NoData；5.创建渔网，利用toolboxs中的创建渔网进行，将输出要素类设置名为fish.shp，范围设置为和研究区一样，像素高度和宽度设置为0.01、0.01，几何类型更改为polygon；

二、北京市城市边界提取6.使用空间分析工具中的以表格显示分区统计工具，输入栅格为fish.shp，字段选择FID，输出表格；7.打开计算后的表格，添加新字段perc，在perc字段使用字段计算器输入AREA/(0.01*0.01)*100；

二、北京市城市边界提取8.利用创建的渔网进行连接表格，将表格属性赋予渔网fish图层；9.打开渔网的符号系统对perc字段进行分层设色，自然间断法分为两类，断点为21；

二、北京市城市边界提取10.打开渔网属性表，点击按属性选择，输入perc字段=21，右击渔网，选择导出数据；11.选择toolboxs中的融合（dissolve）工具，输入要素为上一步导出的数据，融合字段选择id，将创建多要素的勾选去掉；

二、北京市城市边界提取12.在融合之后的文件的属性表中添加字段area，使用计算几何计算面积，按降序排列，选中面积大于10km2的斑块导出数据.

三、提取结果分析计算北京城市边界的面积：2019年建设用地面积为5966.546平方千米，约占北京市面积的36.4%；2023年建设用地面积为10944.78平方千米，约占北京市面积的66.7%；4年间北京市建设用地面积扩张2倍，扩张区域主要体现在北京市的西部和北部城区。

三、提取结果分析计算城市扩张强度指数（UEI）：2019-2023年北京市城市扩张强度指数为20.859，体现北京市城市发展较快。计算紧凑度指数：2019年紧凑度为0.161165；2023年紧凑度为0.223032；表明北京市城区有集中紧凑发展的趋势。

三、提取结果分析生成城市轮廓的最小凸包：利用凸包属性MBG_Width（最短轴长）和MBG_Length（最长轴长）计算形状比2019年形状比为2052.458，延伸率为1.5670962023年形状比为2860.682，延伸率为1.538946说明19年城市的带状特征比23年明显，北京市以环状发展。2019年城市轮廓最小凸包2023年城市轮廓最小凸包