“资源一号 ”02C卫星影像土地利用变化检测能力分析

周 亦 ¹ ，吕 从 ² ，王慧敏 ¹

(1 .北京中色测绘院有限公司 ,北京 100012 ；2 .河北省基础地理信息中心 ，河北 石家庄 050031 )

摘要 ：“资源一号”02C 卫星自成功发射以来，被广泛应用于土地资源、矿产资源、地质环境调查，以及国土资源、地质灾害应急监测等应用领域。“资源一号”02C 卫星是一颗填补我国高分辨率遥感数据空白的卫星，是根据国土资源主体业务需求定制的第一颗国产高分辨率业务卫星。本文以吉林省和浙江省某试验区“资源一号”02C 卫星数据为例，对“资源一号”02C 卫星数据融合影像进行了土地利用遥感监测变化信息检测能力的方法试验。根据试验区数据变化信息提取试验结果，对“资源一号”02C 卫星影像在土地利用动态遥感监测中的变化检测能力和适用性开展了试验、分析和评价工作。

关键词 ：“资源一号”02C ；土地利用；变化检测能力

2011年12月22日11：26，我国在太原卫星发射中心用“长征四号乙”运载火箭，成功将“资源一号”02C卫星送入太空。“资源一号”02C卫星是一颗填补国内高分辨率遥感数据空白的卫星，是根据国土资源主体业务需求定制的第一颗国产高分辨率业务卫星^[1-2]。“资源一号”02C卫星具有高空间分辨率、宽覆盖、快重访等特点，搭载的两台2.36 m高分辨率相机和1台5 m/10 m全色多光谱相机幅宽54 km，能够实现3～5 d内对全球任一地点进行重复观测^[3]。

全国土地利用变更调查监测与核查遥感监测工作是原国土资源部重要的日常性业务之一。通过卫星遥感技术，监测特定时间段内以新增建设用地变化为主的土地利用变化情况。变化信息提取是土地利用动态遥感监测中的关键环节之一^[4]。因此，变化信息的检测能力是评价“资源一号”02C卫星数据在土地利用动态遥感监测中的适用性的基本因素之一。本文将对“资源一号”02C卫星影像进行变化检测能力的试验和评价工作。

1 “资源一号“02C卫星影像变化检测能力试验

以吉林省某试验区和浙江省某试验区为示范区，制作试验区02C卫星的HR融合遥感影像图，叠加2010年遥感影像提取新增建设用地图斑。按照历年《全国土地利用变更调查监测与核查遥感监测技术方案》要求，提取的疑似新增建设用地图斑主要有4种类型：①前时相影像解译为植被覆盖或明显非建设用地特征，后时相影像解译为明显建设痕迹；②上一年度影像解译为植被覆盖或明显非建设用地特征，本年度影像解译为建设用地的推土或填土痕迹特征；③上一年度影像解译为推图或填土特征，本年度影像解译为明显建设用地特征；④上一年度影像解译为植被覆盖或明显非建设用地特征，本年度影像解译为建设用地围墙或栅栏等圈占地特征^[5]。不同类型的图斑在02C影像上的表现如图1所示。

图1 不同类型的图斑示例

以本期QuickBird叠加上一年度遥感影像提取出的新增建设用地图斑为真值，将02C卫星的HR融合影像提取的变化图斑与真值对比，分析02C影像所提取图斑的正确率、遗漏率及面积准确率，最终分析02C数据在土地利用动态遥感监测中的应用潜力。

1.1 吉林试验区

吉林试验区选用的02C卫星影像对比Quick Bird与World View影像提取新增建设用地信息。吉林试验区提取新增建设用地图斑205块，2 747.3亩(1亩=666.67 m²)。提取的最小变化图斑面积为0.2亩。图斑分布情况如图2所示。

图2 吉林试验区遥感监测图斑分布

对随机抽取的15块面积介于1～5亩之间的图斑进行精度评价，单个图斑中误差为0.4亩，平均差异程度为11.8%。其中差异度最大的为37.5%，差异度最小的为0。

1.2 浙江试验区

浙江试验区选用02C卫星遥感正射影像与SPOT 5遥感影像叠加，提取新增建设用地信息。提取新增建设用地图斑1887块，19 751.9亩。提取的最小变化图斑面积为0.3亩。图斑分布情况如图3所示。

图3 浙江试验区遥感监测图斑分布

表1 吉林试验区遥感监测图斑汇总统计 (按类型 ) 亩

表2 吉林试验区遥感监测图斑汇总统计 (按面积 ) 亩

对随机抽取的11块面积介于10～20亩之间的图斑进行精度评价，单个图斑中误差为1.89亩，平均差异程度为8.1%。其中差异度最大的为26.7%，差异度最小的为0。

表3 浙江试验区遥感监测图斑汇总统计 (按类型 ) 亩

表4 浙江试验区遥感监测图斑汇总统计 (按面积 ) 亩

另外，500亩以上的图斑有1个，面积为1 095.3亩。

根据试验区所提取的变化图斑，从其面积精度、属性精度、最小可检测图斑及图斑边界精度4个方面对02C的HR融合数据在土地利用动态遥感监测中的应用潜力进行分析。

2 试验结果分析

所谓“任务驱动”，就是在学习专业技术的过程中，学生在教师的帮助下，紧紧围绕一个共同的任务活动中心，在强烈的问题动机的驱动下，通过对学习资源的积极主动应用，进行自主探索和互动协作的学习，并在完成既定任务的同时，引导学生产生一种学习实践活动。这种教学方法特别适用于学生学习操作类的知识和技能，尤其适于学习工程机械驾驶与维修方面的知识和技能。因此，“任务驱动”教学法在铁道工程机械系统故障排除课程教学中已得到运用，就笔者的教学实践来看，要有效地实施任务驱动教学，应注重以下几个方面。

2.1 属性精度分析

经分析，试验区提取的2092块图斑，实地为新增建设用地的图斑2058块，占总图斑数的93.4%，面积为22 282.3亩，占总图斑面积的96.6%；伪变化图斑34块，占总图斑数的1.5%，面积为216.9亩，占总图斑面积的0.9%；遗漏图斑112块，占总图斑数的5.1%，面积为573.8亩，占总图斑面积的2.5%。试验区新增建设用地图斑情况统计见表5。

W=[WC1×[We11, We12, We13, We14], Wc2×[We21, We22, We23, We24], Wc3×[We31, We32, We33]]=[0.036, 0.138, 0.095, 0.061, 0.214, 0.182, 0.049, 0.085, 0.056, 0.024, 0.060]

要说呢，哪个工地经得住一查呐？你爸的工地上也的确存在许多漏洞，小记者不依不饶，扬言要给你爸曝光。你爸哪见得了这个啊，直接叫几个工人把他们揍了一顿，撵出去了。那不是还有些拍你爸马屁的人吗？你爸那时候是纳税大户，是地方上的财神爷，有人跟报社领导扯上了关系，反而让那个记者挨了处分。

按不同的图斑类型分析，第1类中伪变化图斑21块，面积120.1亩；第2类中伪变化图斑15块，面积148.9亩；第3类中无伪变化图斑。伪变化图斑主要集中在类型前两类上，影像上表现为建设或推平未建部分，实地为非建设用地。

1.对标工作体系，强化问题整改。及时制定并印发《宝胜集团贯彻落实航空工业“1122”党建工作体系实施计划》，召开学习宣贯会，对体系重点内容进行详细解读，通过发放问题整改通知单，督促基层党组织及时自查，列出整改清单，明确具体责任人，限期整改到位。

表5 试验区新增建设用地图斑情况统计

本次变化信息提取共遗漏图斑112块，面积为573.8亩。按不同的图斑类型分析，第1类中遗漏图斑42块，面积90.57亩；第2类中遗漏图斑53块，面积431.9亩；第3类中遗漏图斑2块，面积为14.1亩。根据影像纹理，结合实际情况，对遗漏图斑进行了分析，按不同的图斑面积大小，将遗漏图斑按面积划分为6个等级，分别为1亩以下、1～5亩、5～10亩、10～20亩、20～50亩、50亩以上。由于受影像分辨率和原始数据质量的影响，大部分遗漏图斑为前两类的图斑，且数量集中在1～5亩的区域，第2类的图斑遗漏面积占大部分。

2.2 面积精度分析

根据实际情况，随机选取了一定数量的图斑，按面积划分为6个等级进行面积精度评价，分别为1亩以下、1～5亩、5～10亩、10～20亩、20～50亩、50亩以上，对随机抽取的15块面积≤1亩的图斑进行精度评价，单个图斑中误差为0.1亩，平均差异程度为13.2%。其中差异度最大的为33.3%，差异度最小的为0。

经分析，该区域205块图斑中，变化类型为第1类的图斑131块，面积为1 705.4亩；变化类型为第2类的图斑55块，面积为736.8亩；变化类型为第3类的19块，面积为305.1亩；变化类型为第4类的图斑个数为0。具体情况见表1和表2。

如果可以把所有生长过的海拉细胞堆起来的话，它们可能重达5000万吨。5000万吨是什么概念呢？大概就是 100 幢纽约帝国大厦的重量。如果将所有生长过的海拉细胞从头到尾排列起来，它们可以绕地球至少三圈，相当于1亿多米。但你要知道，拉克斯本人的身高只有 1.5 米。

对随机抽取的15块面积介于5～10亩之间的图斑进行精度评价，单个图斑中误差为0.73亩，平均差异程度为10.6%。其中差异度最大的为34.2%，差异度最小的为0。

经分析，该区域1887块图斑中，变化类型为第1类的图斑962块，面积为9781亩；变化类型为第2类的图斑609块，面积为6 825.4亩；变化类型为第3类的302块，面积为2 963.8亩；变化类型为第4类的14块，面积为181.7亩。具体情况见表3和表4。

对随机抽取的11块面积为50亩以上的图斑进行精度评价，单个图斑中误差为8.67亩，平均差异程度为7.1%。其中差异度最大的为28.3%，差异度最小的为0.5%。

对随机抽取的15块面积介于20～50亩之间的图斑进行精度评价，单个图斑中误差为3.499亩，平均差异程度为11.2%。其中差异度最大的为37.6%，差异度最小的为0。

沥青混凝料的最佳碾压遍数确定原则为在满足设计密度和孔隙率等前提下的最小碾压遍数来确定。由表4、表6和图2综合分析，沥青混合料的最佳碾压遍数为10遍。

综上可知，随着图斑面积的增大，单个图斑中误差逐渐增大，图斑之间的差异程度逐渐减小。

2.3 最小可检测图斑

变化信息最小上图图斑面积要求：土地利用变化信息提取新增建设用地规定的最小上图图斑面积为1亩,本次监测的新增建设用地最小上图图斑面积为0.2亩，因此本次试验结果表明：02C卫星HR融合影像能够满足土地利用遥感监测变化信息提取的精度要求^[6-8]。

2.4 图斑边界精度

随机选取了25块图斑，通过两种数据对比量测，图斑偏移量偏移最大值为8.18 m，最小值为0.20 m，均值为2.77 m,新增建设用地图斑界线与量测的中误差为3.37 m，满足规定限差5 m的要求。

根据表6（1）回归结果，上市公司收到政府补助后总体上存在信号效应，Sub与Own回归系数为正且在1%的水平上显著。6（2）分析了应计盈余管理的调节作用，Sub∗Da系数为正且在5%的水平上显著，说明应计盈余管理在政府补助的信号效应中发挥了正向的调节作用。6（3）验证了真实盈余管理的调节作用，Sub∗Rem的回归系数为正且在1%的水平上显著，说明企业的真实盈余管理不能为社会投资者所识别，收到政府补助的企业通过真实盈余管理行为获得了更多的社会投资者跟进。

3 结 论

通过试验结果分析，在属性精度方面，试验区内实地为新增建设用地的图斑占总图斑数的93.4%，面积占总图斑面积的96.6%；伪变化图斑占总图斑数的1.5%，面积占总图斑面积的0.9%；遗漏图斑占总图斑数的5.1%，面积占总图斑面积的2.5%。伪变化图斑和遗漏图斑集中分布在第1类和第2类上，遗漏数量集中分布在1～5亩的面积档。第2类的图斑遗漏面积占大部分。

在面积精度方面，试验区HR融合影像变化图斑平均差异程度在10%左右，且随着图斑面积的增大，单个图斑中误差逐渐增大，图斑之间的差异程度逐渐减小。

02C卫星HR融合影像最小可监检测图斑面积能够满足土地利用遥感监测变化信息提取的精度要求。图斑边界精度满足方案规定的限差5 m的要求。

囿于本文应用试验仅利用了工程化生产中变化信息提取的常规方法对两个试验区数据进行了应用试验，未能从大区域、规模化的角度开展深入的试验与评价，因此所得出的结论不可避免地具有局限性和片面性^[9-10]。要对“资源一号”02C卫星在土地利用动态遥感监测中变化检测能力进行更客观、更全面的评价，仍需进一步选用更多方法和大量数据进行更加深入与广泛的试验研究。

参考文献 ：

[1] 马利刚，张乐平，邓劲松，等. 资源一号“02C”遥感影像土地利用分类[J].浙江大学学报(工学版), 2013, 47(8):1508-1516.

[2] 周伟，黄炜，王彦佐，等. 资源一号02C卫星数据管理与服务系统研建[J]. 国土资源遥感, 2014, 26(1):179-185.

[3] 刘顺喜, 尤淑撑，张定祥，等.土地资源管理对我国后续资源卫星数据空间分辨率的需求分析[J]. 国土资源遥感, 2003，15(4):6-8.

[4] 柴渊，李万东.土地利用动态遥感监测技术与方法[M].北京：地质出版社，2011.

[5] 中国土地勘测规划院.全国土地利用变更调查监测与核查项目遥感监测任务技术方案[Z]. 2010—2017.

[6] 周亦，武娟，李琦，等.“资源一号”02C卫星影像在土地利用变化信息发现中的试验与分析[J]. 矿产勘查, 2012, 3(5):688-694.

[7] 周亦，李琦. 土地利用动态遥感监测中不同数据源正射影像制作关键技术浅析[J].矿产勘查，2011，2(6)：827-833.

[8] 刘鹰.土地利用动态遥感监测中变化信息提取方法的研究[J]. 遥感信息, 1999(4):21-24.

[9] 吴海平,黄世存.自动提取土地利用变化信息后处理方法的研究[J]. 国土资源遥感, 2005, 17(1):41-44.

[10] 肖立强，姚月，李小华. “资源一号”02C卫星影像几何校正方法实验[J].科技资讯，2014(8)：43-45.

Land use change detection capability analysis with “ZY -1 ”02C data

ZHOU Yi¹，LÜ Cong²，WANG Huimin¹

(1. Beijing CNNC Institute of Surveying and Mapping Co., Ltd., Beijing 100012, China; 2. Hebei Geomatics Center, Shijiazhuang 050031, China)

Abstract ： On December 22, 2011, China launched the “ ZY-1” 02C satellite successfully. “ ZY-1” 02C satellite is a high resolution remote sensing data to fill the domestic blank satellite. And it is the first domestic high resolution satellite services based on the land resources the main business demand. This paper took “ZY-1” 02C in-orbit test data in the test area of Jilin and Zhejiang province as an example, carried out experiments on the method of land use change detection capability with “ZY-1” 02C fusion images. According to the results of change detection experiments in test areas, we carried out experiment, analysis and evaluation of land use dynamic remote sensing monitoring change detection capability with “ZY-1” 02C data.

Key words ： “ZY-1” 02C; land use; change detection capability

中图分类号 ：P237

文献标识码： A

文章编号： 0494 -0911 (2019 )01 -0097 -04

引文格式 ：周亦，吕从，王慧敏.“资源一号”02C卫星影像土地利用变化检测能力分析[J].测绘通报，2019(1)：97-100.DOI:10.13474/j.cnki.11-2246.2019.0020.

收稿日期 ：2018-04-25

作者简介 ： 周 亦(1980—)，女，高级工程师，主要从事测绘地理信息与土地遥感方面的研究。E-mail：zhouyi801223@126.com

标签：“资源一号”02C论文;  土地利用论文;  变化检测能力论文;  北京中色测绘院有限公司论文;  河北省基础地理信息中心论文;