自然资源遥感, 2022, 34(1): 115-126 doi: 10.6046/zrzyyg.2021064

技术方法

基于集成学习和多时相遥感影像的枸杞种植区分类

史飞飞,1,2,3,4,5, 高小红,1,3,4,6, 肖建设2,5, 李宏达1,3,4, 李润祥1,3,4, 张昊1,3,4

1.青海师范大学地理科学学院,西宁 810008

2.青海省气象科学研究所,西宁 810008

3.青海省自然地理与环境过程重点实验室,西宁 810008

4.青藏高原地表过程与生态保育教育部重点实验室,西宁 810008

5.青海省防灾减灾重点实验室,西宁 810008

6.高原科学与可持续发展研究院, 西宁 810008

Classification of wolfberry planting areas based on ensemble learning and multi-temporal remote sensing images

SHI Feifei,1,2,3,4,5, GAO Xiaohong,1,3,4,6, XIAO Jianshe2,5, LI Hongda1,3,4, LI Runxiang1,3,4, ZHANG Hao1,3,4

1. School of Geographical Sciences, Qinghai Normal University, Xining 810008, China

2. Institute of Qinghai Meteorological Science Research, Xining 810008, China

3. Key Laboratory of Physical Geography and Environmental Process of Qinghai Province, Xining 810008, China

4. Ministry of Education Key Laboratory of Tibetan Plateau Land Surface Processes and Ecological Conservation, Xining 810008, China

5. Key Laboratory of Disaster Prevention and Mitigation of Qinghai Province, Xining 810008, China

6. Academy of Plateau Science and Sustainability, Xining 810008, China

通讯作者: 高小红(1963-),女,博士,博士生导师,主要研究方向为遥感应用与地理数据空间分析。Email:xiaohonggao226@126.com

责任编辑: 李瑜

收稿日期: 2021-03-10   修回日期: 2021-07-20  

基金资助: 青海省自然科学基金项目“基于GEE云平台与Landsat卫星长时间序列数据的湟水流域30多年土地利用/土地覆被时空变化研究”资助编号(2021-ZJ-913)

Received: 2021-03-10   Revised: 2021-07-20  

作者简介 About authors

史飞飞(1991-),男,博士研究生,工程师,主要研究方向为遥感应用与地理数据空间分析。Email: shifeifei1203@126.com

摘要

利用遥感技术对柴达木盆地枸杞种植区进行精准提取对当地政府开展市场管理与调控具有重要意义。以典型枸杞种植区诺木洪农场为例,选取Landsat8 OLI和GF-1 WFV影像构建作物生长期内时序NDVI/EVI数据,并采用4种新颖的集成学习分类器(LightGBM,GBDT,XGBoost,RF)和2种应用广泛的机器学习分类器(SVM,MLPC)对枸杞种植区进行分类。研究结果表明: ①LightGBM(90.4%),GBDT(90.4%),XGBoost(89.31%)和RF(86.96%)分类器能获得较高的分类精度,并以LightGBM+EVI的总体分类精度最高,达到了91.67%,Kappa系数为0.90; ②EVI指数在枸杞生长中后期表现更为灵敏,并在同一分类器下使用EVI时序数据能获得更好的枸杞作物制图效果; ③利用GBDT,XGBoost和RF分类器的特征重要性评分方法进行枸杞种植区分类时相特征优选,能够在获取高分类精度的同时进一步降低数据冗余。

关键词: 作物分类; 枸杞; NDVI/EVI时序; 集成学习

Abstract

It is significant for the market management and regulation of local government to accurately extract wolfberry planting areas in the Qaidam Basin using remote sensing technology. Taking the Nuomuhong Farm, a typical wolfberry planting area, as an example, this study selected Landsat8 OLI and GF-1 WFV images to construct the time-series NDVI/EVI data of the crop growth period. Then, this study employed four novel ensemble learning classifiers (i.e., LightGBM, GBDT, XGBoost, and RF) and two widely used machine learning classifiers (SVM and MLPC) to classify wolfberry planting areas. The results show that: ① Relatively high classification accuracy were obtained using LightGBM (90.4%), GBDT (90.4%), XGBoost (89.31%), and RF (86.96%). Most especially, LightGBM-EVI yielded the highest overall classification accuracy (91.67%), with a Kappa coefficient of 0.90; ② Enhanced vegetation index (EVI) is more sensitive in the middle-late stage of the wolfberry growth period. For the same classifier, better mapping effects of wolfberry planting areas can be obtained when time-series EVI data were used; ③ Data redundancy can be further reduced while obtaining high classification accuracy by determining the optimal temporal features of NDVI/EVI classification using the feature importance scores of the GBDT, XGBoost, and RF classifiers.

Keywords: crop classification; wolfberry; NDVI/EVI time series; ensemble learning

PDF (10347KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

史飞飞, 高小红, 肖建设, 李宏达, 李润祥, 张昊. 基于集成学习和多时相遥感影像的枸杞种植区分类[J]. 自然资源遥感, 2022, 34(1): 115-126 doi:10.6046/zrzyyg.2021064

SHI Feifei, GAO Xiaohong, XIAO Jianshe, LI Hongda, LI Runxiang, ZHANG Hao. Classification of wolfberry planting areas based on ensemble learning and multi-temporal remote sensing images[J]. Remote Sensing for Natural Resources, 2022, 34(1): 115-126 doi:10.6046/zrzyyg.2021064

0 引言

青海柴达木盆地作为我国枸杞的新兴种植区,近年来因枸杞种植规模的无序扩大以及市场开拓不足致使供需失衡严重,已显现出枸杞价格持续低迷和田间消极管理引发减产等问题,而实现枸杞种植面积监测与有序调控是保障枸杞市场平稳、持续发展的关键手段[1,2]。与传统地面调查方法相比,遥感技术具备宏观、经济和高效等优势,已广泛应用于作物的分类、种植结构提取和种植面积监测等方面[3]

作物遥感识别与分类是获取农业土地利用信息的基础,近年来采用遥感观测资料进行作物分类的研究较多,其大致可分为2大类[4,5]: 一种为选用作物生长关键期内的单时相影像进行作物分类; 另一种为选用作物生长期内的多时相影像,经过构建时序的植被指数进行作物分类。受我国作物种植结构复杂性和作物光谱相似性影响,在采用单时相影像开展作物分类时,存在“最佳识别期”不易选取以及易出现“同谱异物,同物异谱”的问题,限制了该方法的分类精度[6]。而在同一区域内因不同作物的生长周期普遍存在差异,采用多时相影像能够获取光谱与物候特征信息,具备获取更佳作物分类精度的潜力[7]。在该方法应用初期以采用具有高时间分辨率的MODIS数据为主,并对全球广泛种植的大宗作物开展了遥感监测 [8,9,10],但因MODIS数据空间分辨率较低,且在种植结构复杂的中小规模种植区易出现大量混合像元,致使分类精度难以提升[11,12]。伴随遥感技术的发展,近年来涌现出利用多时相的Landsat7/8[13,14,15]、GF-1[16,17]和Sentinel-2[18,19,20,21]等中高空间分辨率的影像对香梨、葵花、苜蓿等特色经济作物进行提取的案例。通过对比上述研究发现,时序高分影像数据能获得丰富的作物冠层光谱、纹理和物候等特征信息,提高了作物识别的准确性,但也面临如单一数据源易存在关键物候期影像缺失,采用决策树和神经网络等分类器存在分类特征利用不充分等问题[22,23]。针对上述问题,多源数据应用能较好地解决关键时相数据的缺失,并且已有研究对Landsat8 OLI,Landsat7 ETM+和Sentinel-2A MSI等传感器,进行了波段特性对比以及数据间差异消除方法的构建[24,25],能为多源数据综合提供方法参考,但同时尝试引入具有更强学习与泛化能力的分类器可能是提升作物分类精度的关键[23,26]。任何单一的强分类器都有其自身的优势与不足,集成学习(ensemble learning)方法能综合运用各分类器的优势,以提高分类精度,这也是当前研究的热点方向[27]。集成学习依照个体学习器的形成方式可划分为串行或并行的序列化,前者典型代表为提升法(boosting),而后者的代表是自举汇聚法(bootstrap aggregating,bagging)和随机森林(random forest,RF),其中RF对多维特征和数据噪声具有更强的鲁棒性,已广泛应用于作物分类并取得较好效果[28,29],而采用boosting算法进行作物分类鲜有尝试。

本文以柴达木盆地典型枸杞种植基地诺木洪农场为研究区,选取Landsat8 OLI和GF-1 WFV影像构建作物生长期内时序植被指数数据,并采用轻量级梯度提升机(light gradient boosting machine,LightGBM)、梯度提升决策树(gradient boosting decision tree, GBDT)、极端梯度提升(extreme gradient boosting,XGBoost)、RF、支持向量机(support vector machine,SVM)和多层感知分类器(multilayer perceptron classifier,MLPC)共6种分类器,旨在探讨集成学习方法在开展枸杞种植区分类时的适应性。

1 研究区概况

诺木洪农场位于柴达木盆地都兰县境内,生态环境十分纯净,是全国规模较大的枸杞种植基地,被誉为“中国枸杞之乡”。农场地理坐标E96°15'~96°35',N36°20'~36°30',海拔为2 745~2 790 m,总面积约200 km2,诺木洪农场及周边枸杞种植区分布状况见图1。农场属高原大陆性气候,辐射强,昼夜温差大,气候干燥,年平均气温约5 ℃,降雨量约51 mm,而蒸发量在2 200 mm以上。农场地处冲洪积扇前缘的细土平原带,土壤肥沃并被大量开垦进行枸杞种植,南部为山前荒漠戈壁,而北部为沼泽地,诺木洪河为柴达木流域内第六大河流,由南至北穿过农场为其生产生活提供充足水源[30]

图1

图1   枸杞种植区分布

Fig.1   Distribution map of wolfberry planting area


2 数据源与研究方法

2.1 数据采集与预处理

2.1.1 影像数据与预处理

为构建具备高时间分辨率的时序数据集,选取研究区内作物生长阶段4—11月份内的6景Landsat8 OLI影像和3景GF-1 WFV影像,其中Landsat8 OLI影像在谷歌地球引擎(Google Earth Engine,GEE)平台下载(https://code.earthengine.google.com),空间分辨率30 m,重访周期16 d,GF-1 WFV数据在中国资源卫星应用中心下载(http://www.cresda.com/),空间分辨率为16 m,波段数为4,重访周期为2 d。数据详情如表1所示。

表1   影像数据列表

Tab.1  List of image data

影像获
取时间
数据类型传感
空间分
辨率/m
影像质量数据用途
2017-04-25Landsat8OLI30良(5%云量)构建时序多源遥感数据
2017-05-11Landsat8OLI30优(无云)
2017-05-27Landsat8OLI30良(8%云量)
2017-06-28Landsat8OLI30优(无云)
2017-07-14Landsat8OLI30优(无云)
2017-08-02GF-1WFV16优(无云)
2017-08-14GF-1WFV16优(无云)
2017-09-27GF-1WFV16优(无云)
2017-11-03Landsat8OLI30优(无云)

①OLI为陆地成像仪(多光谱),WFV为多光谱宽幅相机。

新窗口打开| 下载CSV


GF-1 WFV数据利用ENVI5.5软件进行预处理,以邻近时期的Landsat8 OLI数据和30 m DEM数据为参考,选用自动采集控制点的正射校正工具(RPC orthorectification using reference image, RPC)进行校正,并将GF-1 WFV进行辐射定标和FLAASH(fast line-of-sight atmospheric analysis of spectral hypercubes)大气校正。Landsat8 OLI数据在GEE平台中可直接获取地表反射率数据,但仍需在GEE平台中通过编写脚本使用多通道掩模算法(C function of mask, CFMask)对影像中的云、阴影和雪进行处理[31]。最后将预处理后的GF-1 WFV和Landsat8 OLI影像空间分辨率重采样为30 m,并将投影统一转换为Albers等面积投影。

2.1.2 实地野外调查数据

为获取准确且均匀分布的地面样点调查数据,用于分类器的训练和精度验证,按1 km×1 km格网将研究区划分为264个单元,于2017年9月在诺木洪农场开展实地调查,通过随机选取格网以及参考调查地点的可通达性,共对101个格网开展地面调查,期间利用GPS对主要地物类型进行定位并拍照,使用ASD FieldSpec 4型地物光谱仪获取枸杞作物的冠层光谱。因研究区内分布大范围的沼泽地和荒漠戈壁致使部分格网单元无法开展实地调查,因而利用GF-1影像通过目视判识获取该格网内的地物类型信息。最终通过调查将研究区待分类对象划分为幼龄期枸杞地、壮龄期枸杞地、行道树、城镇、沼泽地、戈壁滩和公路共7种类型,并获得各类地物样点共1 378个,样点分布见图2

图2

图2   研究区样点分布

Fig.2   Distribution of sample points in the study area


2.1.3 物候期观测数据

2016—2018年在诺木洪农场选取试验田开展作物物候观测,观测地枸杞树龄为5 a左右,图3展示了枸杞树的生长周期,一般枸杞树从5月上旬开始进入萌芽期,5月下旬进入春梢生长期,6月下旬进入开花盛期,7月上旬—9月下旬进入果熟期,期间枸杞树不断开花结果,当地果农一般开展2~3茬摘果,10月初进入落叶始期,至11月初果树开始进入休眠期。

图3

图3   枸杞作物物候期

Fig.3   Phenology of wolfberry crops


2.2 植被指数选取

植被指数能增强光学遥感影像中的植被信息,并被广泛用于描述生物物理特性的变化[32]。其中归一化植被指数(normalized difference vegetation index, NDVI)是监测植被覆盖度和生长状况的最佳指数,但受计算方法的影响容易在高植被覆盖区产生饱和,而增强植被指数(enhanced vegetation index,EVI)加入了土壤背景调节参数和大气修正参数,能弥补NDVI指数的不足[33]。使用植被指数进行作物物候期监测时,如对地面植被类型及覆盖度状况不了解,一般应综合选取NDVI和EVI指数互为补充[33],2种植被指数的计算公式分别为:

NDVI=ρNIR-ρRedρNIR+ρRed,
EVI=2.5(ρNIR-ρRed)L+ρNIR+C1ρRed-C2ρBlue,

式中: ρNIR, ρRedρBlue分别为近红外、红光波段和蓝光波段的反射率值; L,C1C2分别为土壤背景调节参数、大气校正红光参数和大气校正蓝光参数,其中L=1,C1=6,C2=7.5。

2.3 时序植被指数重构

受云雾和水汽因素影响,利用Landsat8 OLI和GF-1 WFV影像构建的NDVI/EVI时间序列数据会夹杂大的波动和噪声,需进行平滑和去噪[34],本文使用SG(Savitzky-Golay)滤波 [18],其计算公式为:

VIi*=j=-mmCjVIi+12m+1,

式中: VIi*为NDVI/EVI滤波后的模拟值; VI为原始的NDVI/EVI值; Cj为滤波窗口第j个NDVI/EVI值的卷积系数; 2m+1为滤波窗口的滑动数组宽度。本文通过反复测试将滑动窗口数设置为5,多项式次数设置为3时能获得最佳效果。

2.4 分类器构建与精度评价

在Visual Studio Code软件中使用sklearn (Scikit-learn),LightGBM和XGBoost等第三方库,实现了LightGBM,GBDT,XGBoost,RF,SVM和MLPC分类器搭建。当前集成学习是机器学习领域的研究热点,它具有准确性高、速度快、参数少等优点,已在叶绿素反演和森林地上生物量估算等方面开展了应用[35,36],其中LightGBM,GBDT和XGBoost均是在Boosting集成学习方式基础上的改进算法,但在基分类器、列采样方式、决策树生长策略、残差拟合方法和并行策略等方面存在差异[37]; RF是以决策树作为基分类器,并采用bagging方式进行集成的,能够处理高维数据,并且无需大训练样本、反复的模型训练和参数调整便能获取较好的分类精度[37]; SVM和MLPC是典型的非参数化机器学习方法,相比传统的参数化方法具有较强的学习及泛化能力,适用于土地利用、作物分类和参量反演等场景[38]

使用sklearn库中的样本分割方法,按照2∶8将样本集随机划分为测试集和训练集。同时分类器训练时超参数的设定对于模型分类效果影响十分关键,一般在参数集合较大时采用随机格网搜索法能够高效选定候选超参数,文中设定10折交叉验证和10次重复来评估超参数值的所有可能组合。通过实验LightGBM的boosting类型设定为GBDT,并对学习率、最大叶子节点数、树深度和子模型树设定值范围并进行寻优; GBDT和RF均对学习率、树节点分裂最小样本数、树终节点最少样本数、树深度和子模型树参数进行寻优; XGBoost中设定类别数为7,对树深度、学习率和损失函数进行寻优; SVM的核函数设定为径向基函数,并对惩罚系数进行寻优; MLPC对隐含层树及节点数、激活函数、权重优化器和学习率进行寻优。

精度评价主要从像素尺度出发通过构建混淆矩阵,计算总体分类精度(overall accuracy,OA)、制图精度(producer accuracy,PA)、用户精度(user accuracy,UA)和Kappa系数进行分类精度评价[39]

2.5 分类特征重要性评估

利用sklearn中的元转换器并与LightGBM,GBDT,XGBoost和RF分类器的特征重要性评分方法配合,用于评估不同时相NDVI/EVI特征在作物分类中的重要性。而SVM和MLPC分类器缺乏特征重要性评分能力,利用单变量特征选择通过对单变量的统计测试来抉择较好的分类特征。

3 结果与分析

3.1 时序植被指数曲线特征

采用SG方法逐像元重建了时序NDVI/EVI图像,并利用地面调查样点提取了研究区内7种地物类型的平均时序NDVI/EVI曲线(图4),经过分析重建后的NDVI/EVI曲线更加平滑,消除了“锯齿”和“毛刺”等噪声,表明SG方法具有较好的滤波和特征保持能力; 壮龄/幼龄期枸杞在不同生长阶段,其NDVI/EVI曲线表现的物候特征变化趋势较为接近,但增减幅度存在差异,如壮龄期与幼龄期枸杞在休眠期NDVI/EVI值均为最低,5月上旬—中旬枸杞由芽开放期转至展叶盛期,该阶段NDVI/EVI曲线上升趋势迟缓,5月下旬—6月下旬,枸杞由春梢开始生长至春梢开花盛期,该阶段NDVI/EVI曲线均呈快速增长,但壮龄期枸杞增幅较大,7月中旬—8月中旬为夏果形成与成熟期,该阶段枸杞NDVI/EVI曲线均达到峰值并保持稳定,其中壮龄期枸杞峰值高于幼龄期。9月中旬—下旬枸杞由秋果成熟盛期逐渐过渡至落叶始期,枸杞NDVI/EVI曲线逐渐开始降低。而到11月上旬,枸杞进入落叶盛期,枸杞NDVI/EVI曲线持续快速下降,伴随气温骤降以及降雪过程的来临,枸杞进入休眠期。研究区内行道树和沼泽地植被NDVI/EVI曲线也均表现出明显的物候特征,其中在田间道路两旁种植的行道树主要以挺直高大的新疆杨为主,生长期内NDVI/EVI值均高于枸杞等其他植被类型,而沼泽地植被的NDVI/EVI值也略高于幼龄期枸杞; 研究区内公路、城镇和戈壁滩非植被类型的NDVI/EVI值均低于0.2,并且无明显峰谷特征。

图4

图4   NDVI和EVI曲线

Fig.4   NDVI and EVI curve


3.2 时序植被指数曲线可分性

对比分析时序NDVI/EVI曲线对于枸杞等不同植被类型间的区分程度发现,因NDVI和EVI指数的计算公式导致同一类型的时序NDVI曲线值较EVI值整体偏高[19]。7种类型的时序NDVI/EVI曲线特征均存在差异,但是壮龄期枸杞与行道树在时序EVI曲线中的差异更加明显,并主要表现在枸杞生长期后期(8月下旬—11月上旬)。为定量分析NDVI/EVI指数在枸杞作物不同生长期内表现的差异性,选取壮龄期枸杞的NDVI/EVI时序曲线,并计算枸杞下一生长期与当前时期NDVI/EVI指数的增加或减少幅度,结果如图5所示,分析发现在枸杞芽开放期至开花盛期(5月11日—6月28日,生长初期),EVI与NDVI指数增幅呈快速上升趋势,而后在枸杞夏果形成至成熟盛期(7月14日—8月2日,生长中期),EVI与NDVI指数增幅降低其数值达到生长期内最大,在枸杞秋果形成至落叶盛期(8月14日—11月3日,生长中后期),EVI和NDVI指数降幅逐渐呈上升趋势。总体而言,EVI指数曲线相较于NDVI指数曲线在枸杞物候期的前期和中期增幅均较大,两者差异并不明显,而在中后期EVI指数曲线的降低幅度明显高于NDVI指数曲线,两者差异逐渐显现。

图5

图5   NDVI及EVI指数增加与减少幅度

Fig.5   NDVI and EVI indices increase and decrease


3.3 分类结果与精度评价

3.3.1 分类结果

利用NDVI/EVI时序数据采用6种分类器进行枸杞种植区分类,结果如图6图7所示,在12种分类结果中各地物类型的细节展现总体较好,其中幼龄期枸杞种植区连片分布在农场南部,多为新开垦的枸杞幼苗种植区,壮龄期枸杞分布在农场大部地区,并被行道树分割呈现规则的棋盘状分布,上述3种类型边界提取完整且准确。沼泽地主要分布在农场的北部,南部为广阔的戈壁滩,而城镇主要散落分布在枸杞种植区内,上述4种类型在部分分类结果中出现了相互混淆。对农场南部国道G109的提取存在间断并易出现碎斑,其原因为道路宽度较窄采用30 m空间分辨率影像对其识别能力存在不足。

图6

图6   基于NDVI时序数据的分类结果

Fig.6   Classification results based on NDVI time series data


图7

图7   基于EVI时序数据的分类结果

Fig.7   Classification results based on EVI time series data


3.3.2 精度评价

对12种分类结果进行精度验证与对比,结果如图8图9所示,分析表明以LightGBM+EVI的OA最高达到91.67%,Kappa系数为0.90,其对幼龄期枸杞和壮龄期枸杞的PA分别为91.0%和98.0%,而MLPC+EVI的OA最低为83.33%,Kappa系数为0.81,但对幼龄期枸杞和壮龄期枸杞的PA也分别达到82.0%和89.0%,表明上述6种分类器与2种分类数据组合均能实现枸杞种植区的准确提取。文中分类结果按分类精度由高至低依次为: LightGBM+EVI,GBDT+EVI,XGBoost+EVI/GBDT+NDVI,LightGBM+NDVI,XGBoost+NDVI,RF+EVI,RF+NDVI,SVM+EVI,MLPC+NDVI,SVM+NDVI和MLPC+EVI。基于6种分类器对采用NDVI/EVI时序数据下的平均分类精度进行对比分析, LightGBM和GBDT的平均OA均为90.40%,Kappa系数均为0.89,而分类效果最差的MLPC,其平均OA也达到84.06%,Kappa系数为0.82,表明6种分类器对枸杞种植区及其他地物类型均有较强区分能力,按照分类器精度由高至低依次为: LightGBM/GBDT,XGBoost,RF,SVM和MLPC。对研究区不同地物类型的精度进行分析,在12种分类结果中幼龄期枸杞的平均PA为91.08%,壮龄期枸杞的平均PA为94.17%,表明选用的分类器均能实现枸杞种植区的精准提取,通过混淆矩阵分析也存在壮龄期枸杞种植区易于与沼泽地和行道树类型发生混淆,而幼龄期枸杞易于与戈壁滩产生错分。研究区内行道树、公路、戈壁滩和沼泽地的平均PA也较高,分别达到了86.58%,88.25%,85.00%和85.42%,而6种分类器对城镇的提取效果普遍不佳,平均PA仅为79%,农场城镇规模较小且散落分布,并且其图像特征也与公路、戈壁滩类型较为接近,当影像空间分辨率不足以及分类器区分能力较差时易发生上述类型边界提取模糊和错分现象。比照NDVI和EVI这2种时序数据的分类效果,在同一分类器下,使用EVI时序数据的OA将高出1.09~2.54个百分点,同时在对壮龄期枸杞、行道树和沼泽地等类型进行提取时,采用EVI时序数据的PA较采用NDVI时序数据时的PA平均高1.70个百分点。

图8

图8   制图精度和用户精度的热图

Fig.8   Producer accuracy and user accuracy heatmap


图9

图9   总体分类精度和Kappa系数的雷达图

Fig.9   Overall classification accuracy and Kappa coefficient Radar map


3.4 分类时相优选与评价

通过模型特征重要性评分方法分析了不同时相NDVI/EVI数据对于分类的贡献程度,并对优选后的部分时相特征进行分类和精度验证,结果见图10,其中NDVI-OA与NDVI-Kappa分别为分类器特征选取前模型的OA和Kappa系数,NDVI-OA-S与NDVI-Kappa-S分别为分类器特征选取后模型的OA和Kappa系数,分析发现GBDT+EVI,XGBoost+NDVI,SVM+NDVI,SVM+EVI和RF+EVI共5种分类器与数据组合方式下,优选5~6个关键时相的NDVI/EVI数据能获得与采用全部时相NDVI/EVI数据相同的分类精度,并且GBDT+EVI和XGBoost+NDVI在选取6个特征后的OA均达到89.85%,Kappa系数为0.88,对幼龄期枸杞和壮龄期枸杞的PA分别为94.0%和96.0%,表明通过特征优选能在保证分类精度的同时进一步降低数据冗余。LightGBM通过特征重要性评估能够选择更少的时相特征,但对比特征优选前后的OA下降了3~8个百分点。

图10

图10   时相特征选取前后模型的总体分类精度和Kappa系数对比

Fig.10   Comparison of overall classification accuracy and Kappa coefficient of models before and after the selection of temporal features


对分类器选取的时相次数进行统计分析(图11),其中6月下旬和9月下旬时相被选取次数最多,该时段分别处于枸杞开花盛期和枸杞秋果成熟盛期,应为枸杞种植区分类的关键时相。同时通过对比NDVI和EVI这2种数据源的时相特征选取分布规律,针对EVI数据选取的特征时相集中分布在6—9月,而NDVI时序数据选取的时相特征集中分布在5—8月。

图11

图11   时序NDVI/EVI的时相选择次数统计

Fig.11   Statistics of time-phase selection times of NDVI and EVI


4 讨论

中高分辨率遥感数据为作物识别提供了丰富数据源,同时蓬勃发展的集成学习方法将推动作物识别精度的提升。首先,文中综合GF-1 WFV和Landsat8 OLI数据进行互补以获取时序较为完整的植被指数数据,前期宋军伟等[40]已对上述2种卫星的反射率数据进行对比,表明其对应波段均存在显著相关,并且对典型地物类型的可分离性也十分接近。本文在此研究基础上进一步对其计算的植被指数进行差异分析,具体是采用相同目标区域的同日过境数据,并对比分析了影像预处理前后2种数据提取的NDVI/EVI相关程度(图12),结果表明经影像预处理后提取的GF-1 NDVI/EVI与Landsat8 NDVI/EVI间的相关性被大幅提高,其决定系数由0.73~0.74上升至0.97~0.98,表明上述2种传感器的植被指数具有较高的一致性,并且在数据互补使用过程中开展辐射定标、大气校正、几何精纠正和异常值剔除等影像预处理能够提高数据间的一致性,而在后续研究中可通过建立GF-1和Landsat8卫星对应波段间的转换关系以更好地消除数据源间的差异; 同时,为获取准确且灵敏的植被指数,选用同类研究中普遍使用的NDVI和EVI指数,通过对比发现EVI指数在枸杞生育期的中后期更为灵敏,并与行道树等高覆盖植被类型间的时序曲线差异更明显,这与白燕英等[14]认为在生育期高峰期和高植被覆盖下EVI指数更加灵敏的结论较为一致。并且为进一步分析NDVI和EVI这2种时序数据对作物分类精度的影响,采用6种分类器进行评估,结果表明在同一分类器下使用EVI时序数据的分类精度普遍高出1.09~2.54个百分点,EVI指数用于作物分类具有更好的区分性; 最后,选用LightGBM,GBDT,XGBoost,RF,SVM和MLPC 6种分类器和NDVI/EVI这2种时序数据开展枸杞种植区分类,各分类器的平均OA分别为90.4%,90.4%,89.31%,86.96%,85.14%和84.06%。表明4种集成学习分类器相较于SVM和MLPC分类器,可获取更高的作物识别精度。目前,集成学习已在土地利用[41]、叶绿素反演[35]和地上生物量估算[36]方面进行应用并获得较好精度,而本研究则表明LightGBM,GBDT,XGBoost和RF分类器在枸杞作物分类中也具有较好适用性。

图12

图12   预处理前后GF-1和Landsat8数据的NDVI及EVI指数点密度对比

Fig.12   Comparison of NDVI and EVI indices point density of GF-1 and Landsat8 data before and after preprocessing


5 结论

综合应用Landsat8 OLI和GF-1 WFV影像,构建作物生育期内时序NDVI/EVI数据,并采用6种分类器对枸杞种植区进行分类,主要结论如下:

1)LightGBM,GBDT,XGBoost和RF集成学习分类器在枸杞种植区分类中能够获得更高精度,6种分类器精度由高至低依次为: LightGBM/GBDT,XGBoost,RF,SVM和MLPC,在分类结果中以LightGBM+EVI的OA值最高达到91.67%,Kappa值为0.90,其对幼龄期枸杞和壮龄期枸杞的PA值均达到91%以上。

2)时序NDVI/EVI数据均能较好区分枸杞作物以及种植区内其他地物类型,但EVI指数在枸杞作物生长中后期更为灵敏,并在相同分类器下使用EVI时序数据能够获得更好的枸杞制图精度。

3)选用GBDT,XGBoost和RF分类器的特征重要性评分方法进行NDVI/EVI数据分类时相特征优选后,能够保证在分类精度不损失的同时降低数据冗余提高分类效率。

4)文中仅选用单一年份的4—11月影像构建时序植被指数数据,而在后期研究中应考虑使用多年份的时序数据以进一步提高分类结果的可靠性。

(责任编辑: 李瑜)

参考文献

徐常青, 刘赛, 徐荣, .

我国枸杞主产区生产现状调研及建议

[J]. 中国中药杂志, 2014, 39(11):1979-1984.

[本文引用: 1]

Xu C Q, Liu S, Xu R, et al.

Investigation of production status in major wolfberry producing areas of China and some suggestions

[J]. China Journal of Chinese Materia Medica, 2014, 39(11):1979-1984.

[本文引用: 1]

朱时佳, 谢韶璟, 余秦胤, .

边远地区农特产品供需失配问题研究——以青海枸杞为例

[J]. 现代商贸工业, 2020, 41(2):25-27.

[本文引用: 1]

Zhu S J, Xie S J, Yu Q Y, et al.

Research on the problem of supply and demand mismatch of special agricultural products in remote areas:Taking Qinghai wolfberry as an example

[J]. Modern Business Trade Industry, 2020, 41(2):25-27.

[本文引用: 1]

董金玮, 吴文斌, 黄健熙, .

农业土地利用遥感信息提取的研究进展与展望

[J]. 地球信息科学学报, 2020, 22(4):772-783.

DOI:10.12082/dqxxkx.2020.200192      [本文引用: 1]

农业用地占到全球土地面积近一半,农业土地利用(包括耕地及作物分布、种植制度、土地管理等)变化直接影响到粮食安全、水安全、生态安全和气候变化。遥感已经成为土地利用信息获取的重要手段,近年来中分辨率遥感卫星如Landsat、Sentinel以及中国高分卫星等的免费开放为国内外农业土地利用信息提取提供了前所未有的机遇,取得了一系列重要研究进展。本文从耕地分布、作物类型识别、农业种植制度以及农业土地管理4个角度分析了土地利用信息提取的最新研究进展。结果发现:① 耕地分布产品已经由过去的粗分辨率提升到10~30 m,耕地现状数据较为丰富,但挖掘遥感数据实现耕地变化历史回溯的能力有待加强;② 作物分类方面多采用地面调查数据和卫星遥感(Landsat和Sentinel-2为主)相结合的方式进行,在北美和欧洲得到了业务化运行,但对作物种植面积早期监测的能力有待加强;③ 基于遥感的农业种植制度信息获取(如撂荒)研究多集中在东欧等地区,在中国由于经济和政策因素导致的撂荒、轮作、休耕等现象也十分普遍,但具有针对性的遥感监测研究目前还相对缺乏;④ 农业土地管理措施信息提取方面,区域灌溉面积产品取得了重要进展,但数据的可靠性和准确性仍有待提高。在此基础上,我们结合遥感大数据、深度学习算法、云计算平台的发展对未来农业土地利用信息提取研究进行了展望:① 融合多源数据形成更高维度空间、光谱和时间信息的遥感大数据,提升特征提取和数据挖掘能力;② 机器学习和深度学习算法等智能化方法与基于地理学和物候信息的专家知识方法的耦合;③ 遥感云计算和大数据挖掘等前沿遥感和计算技术的应用。

Dong J W, Wu W B, Huang J X, et al.

State of the art and perspective of agricultural land use remote sensing information extraction

[J]. Journal of Geo-Information Science, 2020, 22(4):772-783.

[本文引用: 1]

Yi Z W, Jia L, Chen Q T.

Crop classification using multi-temporal Sentinel-2 data in the Shiyang River basin of China

[J]. Remote Sensing, 2020, 12(24):4052-4073.

DOI:10.3390/rs12244052      URL     [本文引用: 1]

Zhong L H, Gong P, Biging G S.

Efficient corn and soybean mapping with temporal extendability:A multi-year experiment using Landsat imagery

[J]. Remote Sensing of Environment, 2014, 140:1-13.

DOI:10.1016/j.rse.2013.08.023      URL     [本文引用: 1]

Conese C, Maselli F.

Use of multitemporal information to improve classification performance of TM scenes in complex terrain

[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 1991, 46(4):187-197.

DOI:10.1016/0924-2716(91)90052-W      URL     [本文引用: 1]

Pena M A, Brenning A.

Assessing fruit-tree crop classification from Landsat8 time series for the Maipo Valley,Chile

[J]. Remote Sensing of Environment, 2015, 171:234-244.

DOI:10.1016/j.rse.2015.10.029      URL     [本文引用: 1]

Lobell D B, Asner G P.

Cropland distributions from temporal unmixing of MODIS data

[J]. Remote Sensing of Environment, 2004, 93:412-422.

DOI:10.1016/j.rse.2004.08.002      URL     [本文引用: 1]

Richard M, Sankey T T, Congalton R G, et al.

MODIS phenology-derived,multi-year distribution of conterminous US crop types

[J]. Remote Sensing of Environment, 2017, 198:490-503.

DOI:10.1016/j.rse.2017.06.033      URL     [本文引用: 1]

许青云, 杨贵军, 龙慧灵, .

基于MODIS NDVI多年时序数据的农作物种植识别

[J]. 农业工程学报, 2014, 30(11):134-144.

[本文引用: 1]

Xu Q Y, Yang G J, Long H L, et al.

Crop information identification based on MODIS NDVI time-series data

[J]. Transactions of the Chinese Society of Agricultural Engineering, 2014, 30(11):134-144.

[本文引用: 1]

平跃鹏, 臧淑英.

基于MODIS 时间序列及物候特征的农作物分类

[J]. 自然资源学报, 2016, 31(3):503-514.

[本文引用: 1]

Ping Y P, Zang S Y.

Crop identification based on MODIS NDVI time-series data and phenological characteristics

[J]. Journal of Natural Resources, 2016, 31(3):503-514.

[本文引用: 1]

郭昱杉, 刘庆生, 刘高焕, .

基于MODIS 时序NDVI主要农作物种植信息提取研究

[J]. 自然资源学报, 2017, 32(10):1808-1818.

[本文引用: 1]

Guo Y S, Liu Q S, Liu G H, et al.

Extraction of main crops in Yellow River delta based on MODIS NDVI time series

[J]. Journal of Natural Resources, 2017, 32(10):1808-1818.

[本文引用: 1]

汪小钦, 邱鹏勋, 李娅丽, .

基于时序Landsat遥感数据的新疆开孔河流域农作物类型识别

[J]. 农业工程学报, 2019, 35(16):180-188.

[本文引用: 1]

Wang X Q, Qiu P X, Li Y L, et al.

Crops identification in Kaikong River basin of Xinjiang based on time series Landsat remote sensing images

[J]. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(16):180-188.

[本文引用: 1]

白燕英, 高聚林, 张宝林.

基于Landsat8影像时间序列NDVI的作物种植结构提取

[J]. 干旱区地理, 2019, 42(4):893-901.

[本文引用: 2]

Bai Y Y, Gao J L, Zhang B L.

Extraction of crop planting structure based on time-series NDVI of Landsat8 images

[J]. Arid Land Geography, 2019, 42(4):893-901.

[本文引用: 2]

Pareeth S, Karimi P, Shafiei M, et al.

Mapping agricultural landuse patterns from time series of Landsat8 using random forest based hierarchial approach

[J]. Remote Sensing, 2019, 11(5):601-615.

DOI:10.3390/rs11050601      URL     [本文引用: 1]

魏鹏飞, 徐新刚, 杨贵军, .

基于多时相影像植被指数变化特征的作物遥感分类

[J]. 中国农业科技导报, 2019, 21(2):54-61.

[本文引用: 1]

Wei P F, Xu X G, Yang G J, et al.

Remote sensing classification of crops based on the change characteristics of multi-phase vegetation index

[J]. Journal of Agricultural Science and Technology, 2019, 21(2):54-61.

[本文引用: 1]

刘雅清, 王磊, 赵希妮, .

基于GF-1/WFV时间序列的绿洲作物类型提取

[J]. 干旱区研究, 2019, 36(3):781-789.

[本文引用: 1]

Liu Y Q, Wang L, Zhao X N, et al.

Extraction of crops in oasis based on GF-1/WFV time series

[J]. Arid Zone Research, 2019, 36(3):781-789.

[本文引用: 1]

杜保佳, 张晶, 王宗明, .

应用Sentinel-2A NDVI时间序列和面向对象决策树方法的农作物分类

[J]. 地球信息科学学报, 2019, 21(5):740-751.

DOI:10.12082/dqxxkx.2019.180412      [本文引用: 2]

农作物种植结构是农业生产活动对土地利用的表现形式。及时精确地获取农作物的空间分布信息对指导农业生产、合理分配资源以及解决粮食安全问题等具有重要意义。目前农作物信息提取研究大多局限于中低分辨率遥感影像的NDVI时间序列,影响了作物空间分布信息提取的准确性。随着Sentinel-2A卫星成功发射,为高分辨率NDVI时间序列的构建提供了可能。本文以黑龙江省北安市为研究区,基于覆盖完整生育期的Sentinel-2A多光谱数据,构建10 m分辨率的NDVI时间序列数据集,利用 Savitzky Golay (S-G) 滤波器对 Sentinel-2A NDVI时间序列数据进行平滑。基于典型时相的多光谱数据和NDVI时间序列构建面向对象决策树分类模型进行作物类型遥感识别。通过对样本的NDVI时间序列曲线分析,可以得出NDVI时间序列能够清晰地区分作物物候差异。此外,本文还利用面向对象分类和支持向量机(Support Vector Machine, SVM)分类两种方法,对典型时相的多光谱数据进行了作物分类对比实验,并对结果进行了对比分析。研究结果表明:① 典型时相多光谱数据引入平滑重构后的NDVI时间序列能够更好地描述作物的物候特性,能够准确刻画研究区作物发育情况,有效区分各类作物;② 通过对比分类实验发现,典型时相多光谱数据引入NDVI时间序列特征,增强了不同作物之间的光谱差异,提高了作物分类精度,总体精度和kappa系数较典型时相多光谱数据进行分类的结果分别提高了7.7% 和0.055;③ 基于面向对象的决策树分类模型在作物分类的结果中精度最高,总体精度为96.2%,kappa系数为0.892。本研究的方法为其他大区域农作物的分类提供了重要参考和借鉴价值。

Du B J, Zhang J, Wang Z M, et al.

Crop mapping based on Sentinel-2A NDVI time series using object-oriented classification and decision tree model

[J]. Journal of Geo-Information Science, 2019, 21(5):740-751.

[本文引用: 2]

Ali N, Nicolas B, Mario M, et al.

A novel approach for mapping wheat areas using high resolution Sentinel-2 images

[J]. Sensors, 2018, 18(7):2089-2110.

DOI:10.3390/s18072089      URL     [本文引用: 2]

Vrieling A, Meroni M, Darvishzadeh R, et al.

Vegetation phenology from Sentinel-2 and field cameras for a Dutch Barrier island

[J]. Remote Sensing of Environment, 2018, 215:517-529.

DOI:10.1016/j.rse.2018.03.014      URL     [本文引用: 1]

谷祥辉, 张英, 桑会勇, .

基于哨兵2时间序列组合植被指数的作物分类研究

[J]. 遥感技术与应用, 2020, 35(3):702-711.

[本文引用: 1]

Gu X H, Zhang Y, Sang H Y, et al.

Research on crop classification method based on Sentinel-2 time series combined vegetation index

[J]. Remote Sensing Technology and Application, 2020, 35(3):702-711.

[本文引用: 1]

李亚飞, 董红斌.

基于卷积神经网络的遥感图像分类研究

[J]. 智能系统学报, 2018, 13(4):550-556.

[本文引用: 1]

Li Y F, Dong H B.

Classification of remote sensing image based on convolutional neural network

[J]. CAAI Transactions on Intelligent Systems, 2018, 13(4):550-556.

[本文引用: 1]

Masoud M, Bahram S, Mohammad R, et al.

Very deep convolutional neural networks for complex land cover mapping using multispectral remote sensing imagery

[J]. Remote Sensing, 2018, 10(7):1119-1140.

DOI:10.3390/rs10071119      URL     [本文引用: 2]

徐光志, 徐涵秋.

Sentinel-2A MSI和Landsat8 OLI两种传感器多光谱信息的交互对比

[J]. 遥感技术与应用, 2021, 36(1):165-175.

[本文引用: 1]

Xu G Z, Xu H Q.

Cross-comparison of Sentinel-2A MSI and Landsat8 OLI multispectral information

[J]. Remote Sensing Technology and Application, 2021, 36(1):165-175.

[本文引用: 1]

Mancino G, Ferrara A, Padula A, et al.

Cross-comparison between Landsat8(OLI)and Landsat7(ETM+)derived vegetation indices in a mediterranean environment

[J]. Remote Sensing, 2020, 12(2):291-311.

DOI:10.3390/rs12020291      URL     [本文引用: 1]

史飞飞, 雷春苗, 肖建设, .

基于多源遥感数据的复杂地形区农作物分类

[J]. 地理与地理信息科学, 2018, 34(5):49-55,2.

[本文引用: 1]

Shi F F, Lei C M, Xiao J S, et al.

Classification of crops in complicated topography area based on multisource remote sensing data

[J]. Geography and Geo-Information Science, 2018, 34(5):49-55,2.

[本文引用: 1]

Zhang H, Li Y, Jiang Y, et al.

Hyperspectral classification based on Lightweight 3-D-CNN with transfer learning

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2019, 57(8):5813-5828.

DOI:10.1109/TGRS.36      URL     [本文引用: 1]

Markus I, Clement A, Tatjana K.

Tree species classification with random forest using very high spatial resolution 8-band WorldView-2 satellite data

[J]. Remote Sensing, 2012, 4(9):2661-2693.

DOI:10.3390/rs4092661      URL     [本文引用: 1]

蔡耀通, 刘书彤, 林辉, .

基于多源遥感数据的CNN水稻提取研究

[J]. 自然资源遥感, 2020, 32(4):97-104.doi: 10.6046/gtzyyg.2020.04.14.

DOI:10.6046/gtzyyg.2020.04.14      [本文引用: 1]

Cai Y T, Liu S T, Lin H, et al.

Extraction of paddy rice based on convolutional neural network using multi -source remote sensing data

[J]. Remote Sensing for Land and Resources, 2020, 32(4):97-104.doi: 10.6046 /gtzyyg.2020.04.14.

DOI:10.6046 /gtzyyg.2020.04.14      [本文引用: 1]

崔亚莉, 刘峰, 郝奇琛, .

诺木洪冲洪积扇地下水氢氧同位素特征及更新能力研究

[J]. 水文地质工程地质, 2015, 42(6):1-7.

[本文引用: 1]

Cui Y L, Liu F, Hao Q C, et al.

Characteristics of hydrogen and oxygen isotopes and renewability of groundwater in the Nuomuhong alluvial fan

[J]. Hydrogeology and Engineering Geology, 2015, 42(6):1-7.

[本文引用: 1]

Foga S, Scaramuzza P L, Guo S, et al.

Cloud detection algorithm comparison and validation for operational Landsat data products

[J]. Remote Sensing of Environment, 2017, 194:379-390.

DOI:10.1016/j.rse.2017.03.026      URL     [本文引用: 1]

Villa P, Bresciani M, Braga F, et al.

Comparative assessment of broadband vegetation indices over aquatic vegetation

[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 7(7):3117-3127.

DOI:10.1109/JSTARS.4609443      URL     [本文引用: 1]

白燕英, 高聚林, 张宝林.

基于NDVI与EVI的作物长势监测研究

[J]. 农业机械学报, 2019, 50(9):153-161.

[本文引用: 2]

Bai Y Y, Gao J L, Zhang B L.

Monitoring of crops growth based on NDVI and EVI

[J]. Transactions of the Chinese Society for Agricultural Machinery, 2019, 50(9):153-161.

[本文引用: 2]

解毅, 张永清, 荀兰, .

基于多源遥感数据融合和LSTM算法的作物分类研究

[J]. 农业工程学报, 2019, 35(15):129-137.

[本文引用: 1]

Xie Y, Zhang Y Q, Xun L, et al.

Crop classification based on multi-source remote sensing data fusion and LSTM algorithm

[J]. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(15):129-137.

[本文引用: 1]

Huang G M, Wu L F, Ma X, et al.

Evaluation of CatBoost method for prediction of reference evapotranspiration in humid regions

[J]. Journal of Hydrology, 2019, 574:1029-1041.

DOI:10.1016/j.jhydrol.2019.04.085      URL     [本文引用: 2]

Su H, Yang X, Lu W F, et al.

Estimating subsurface thermohaline structure of the global ocean using surface remote sensing observations

[J]. Remote Sensing, 2019, 11(13):1598-1620.

DOI:10.3390/rs11131598      URL     [本文引用: 2]

周志华. 机器学习[M]. 北京: 清华大学出版社, 2016:171-196.

[本文引用: 2]

Zhou Z H. Machine learning[M]. Beijing: Tsinghua University Press, 2016:171-196.

[本文引用: 2]

李宏达, 高小红, 汤敏.

基于CNN 的不同空间分辨率影像土地覆被分类研究

[J]. 遥感技术与应用, 2020, 35(4):749-758.

[本文引用: 1]

Li H D, Gao X H, Tang M.

Land cover classification for different spatial resolution images from CNN

[J]. Remote Sensing Technology and Application, 2020, 35(4):749-758.

[本文引用: 1]

Pleoianu A I, Stupariu M S, Sandric I, et al.

Individual tree-crown detection and species classification in very high-resolution remote sensing imagery using a deep learning ensemble model

[J]. Remote Sensing, 2020, 12(15):2426-2448.

DOI:10.3390/rs12152426      URL     [本文引用: 1]

宋军伟, 张友静, 李鑫川, .

基于GF-1与Landsat8影像的土地覆盖分类比较

[J]. 地理科学进展, 2016, 35(2):255-263.

DOI:10.18306/dlkxjz.2016.02.012      [本文引用: 1]

高分一号(GF-1)卫星具有多种分辨率与大幅宽结合、重访周期短等优势,而Landsat-8卫星具有多波段、高辐射分辨率等优势。针对不同传感器参数特点,利用支持向量机分类器(SVM)对同区域同期两种数据进行土地覆盖分类对比研究。结果表明:两种传感器对应波段决定系数均大于0.92;典型样本的光谱趋势一致性良好,但在农田与林地、不透水面与裸土的典型样本可分离性方面,Landsat-8优于GF-1;GF-1与Landsat-8的分类总精度分别为90.38%和90.07%,但不同地物类型的分类精度存在差异,波谱响应函数的差异可能是导致Landsat-8对林地的分类精度高于GF-1的原因;此外,GF-1对零碎分布地物类型的分类精度高于Landsat-8,主要原因是GF-1具有更高的空间分辨率。

Song J W, Zhang Y J, Li X C, et al.

Comparison between GF-1 and Landsat8 images in land cover classification

[J]. Progress in Geography, 2016, 35(2):255-263.

[本文引用: 1]

Vasilakos C, Kavroudakis D, Georganta A.

Machine learning classification ensemble of multitemporal Sentinel-2 images:The case of a mixed mediterranean ecosystem

[J]. Remote Sensing, 2020, 12(12):2005-2030.

DOI:10.3390/rs12122005      URL     [本文引用: 1]

/

京ICP备05055290号-2
版权所有 © 2015 《自然资源遥感》编辑部
地址:北京学院路31号中国国土资源航空物探遥感中心 邮编:100083
电话:010-62060291/62060292 E-mail:zrzyyg@163.com
本系统由北京玛格泰克科技发展有限公司设计开发