摘要
水稻是我国南方地区种植面积最大的粮食作物,准确及时地获取水稻的面积及空间分布信息对于加强农业生产管理、保障国家粮食安全具有重要意义。以往利用卫星遥感技术进行水稻分类研究中影像的红边波段类型没有确定,影响着水稻面积监测的精度和效率。
文章选取中国南方典型的水稻种植区——浙江省德清县为研究区,利用覆盖水稻主要生长期的5个时相(2019年5月24日、8月2日、8月17日、10月31日、11月15日)Sentinel-2卫星影像,采用随机森林算法进行水稻分类研究。分析研究区水稻及其余6种典型地物在不同波段下的反射光谱特征,采用J-M距离定量评价各时相下有无红边波段参与的各种地物光谱可分离性,比较不同红边波段下的水稻分类精度,优选水稻分类的关键时相和波段组合。
相较于红、绿、蓝波段,7种典型地物在红边、近红外波段光谱差异更为明显。相比于无红边波段,有红边波段参与时水稻与其余6种地物的J-M距离更大,可分离性更强。5个时相中,在8月17日(水稻拔节期)时,水稻与其他地物可分离性最大。当3个红边波段全部参与分类时,水稻与其余6种典型地物的总体分类精度最高为91.58%,Kappa系数为0.89。Sentinel-2影像参与分类的7个波段,按重要性排序是:近红外波段>红边波段2>红边波段3>红波段>红边波段1>绿波段>蓝波段。
关键词
我国水稻种植面积居世界第一,水稻产量占全国粮食总产量的34%,在我国粮食生产中占有至关重要的地
Sentinel-2影像相对MODIS、Landsat、GF-1 WFV影像来说空间分辨率更高,对地块破碎、种植复杂地区作物识别能力更强,另外还包含红边波段,对监测植被健康信息非常有效,红边波段已在农作物监测、农作物面积提取中广泛应
针对上述问题,文章选取浙江省德清县为研究区,首先分析各个时相上地物的光谱特征,然后比较有无红边波段条件下水稻与其他地物的可分性以及分类精度,分析引入红边波段对水稻识别能力的影响,最后比较不同类型红边波段参与分类情况下水稻提取精度,通过随机森林重要性度量验证结果。优选出适合水稻提取的红边波段类型及时相,为改善南方地区水稻提取精度提供依据。
为评价红边波段对水稻识别能力影响,该研究利用Sentinel-2遥感影像进行水稻提取研究。首先分析各个时相上地物的光谱特征,然后采用随机森林监督分类方法,获取有无红边波段、不同红边波段组合条件下分类精度,并计算J-M距离作为评价水稻与其他地物间可分离性指标。比较红边波段条件下水稻识别能力的提升作用,借此分析引入红边波段对水稻提取精度的影响,以及Sentinel-2影像红边波段的重要性。
该研究选取中国南方典型地区德清县作为研究区。德清县位于浙江省北部,地处30°26′E~30°42′E,119°45′N~120°21′N(

图1 研究区地理位置及野外样点分布
研究区水稻物候自移栽后,依次经历返青期、分蘖期、拔节孕穗期、抽穗期、成熟期5个生长阶段,水稻各个物候期对应的具体日期
数据源选择Sentinel-2卫星遥感影像,该卫星携带一枚多光谱成像仪(MSI),飞行高度为786km,可覆盖13个光谱波段(包括3个红边波段),幅宽达290 km。从可见光、近红外到短波红外,具有不同的空间分辨率,空间分辨率分别为10m、20m和60m,一颗卫星的重访周期为10d,两颗卫星互补,重访周期为5d。Sentinel-2卫星各波段详细信息如
Sentinel-2多光谱影像预处理步骤包括辐射定标、大气校正、图像拼接与裁剪、图像融合处理。采用ENVI5.3软件对Sentinel-2影像进行辐射定标,将遥感影像记录的原始DN值转换为辐射亮度值。大气校正是通过FLASSH大气校正模块将辐射定标后影像辐亮度转换为地表反射率,目的是为了消除大气和光照等因素对地物反射的影响。为了提高红边波段的空间分辨率,将20m的红边波段与10m的近红外波段进行融合,采用的是ENVI5.3模块Gram-Schmidt算法。对融合后的Sentinel-2影像进行拼接处理覆盖整个研究区,然后利用德清县矢量边界进行裁剪,得到研究区影像数据。
在德清县境内开展地面调查,为水稻提取和精度验证提供参考依据,选取的样点在研究区内尽量均匀分布。在德清县境内开展为期3d的野外调查,调查德清县主要土地利用类型,并结合高分辨率的Google影像对研究区内的建筑、水稻、水体等地物进行GPS定位。地面采集地物类型样本数如表3所示,其中2/3作为训练样本,1/3作为验证样本,同时利用相机对野外环境进行拍照记录,此次共采集野外样点分布如
地物可分性判定依据有J-M距离、欧式距离、光谱距离、离散度
(1) |
(2) |
随机森林模型整合了多棵决策树,是一种较为实用的集成学习方法,随机森林模型有两个重要参数,分别是决策树棵数以及分裂结点个数。随机森林是利用多棵决策树对数据进行训练、分类和预测的方
随机森林是由多个决策树构成,其中每一棵决策树之间是没有关联的。采用基尼指数对分类过程中每棵决策树的每一个节点进行纯度判断并进行最优属性划分,最终使得每个节点样本尽可能属于同一类别,随着划分过程的不断进行,结点的类别纯度越
为了分析研究区各个时相不同波段下水稻、林地、建筑、水体、其他植被、茭白、大棚7种典型地物的光谱差异,

图2 2019年5个时相不同波段下典型地物反射光谱变化
Sentinel-2影像包括3个红边波段,该研究分别基于全部7个波段及不包括红边的4个波段计算J-M距离以及分类精度。为了分析红边波段对水稻可分离性及识别能力的影响,
有红边波段分类精度均优于无红边波段,5月24日总体分类精度由87.38上升到88.36%,增加了0.02个百分点,水稻制图精度由69.80%上升到74.17%,上升了4.37个百分点;8月2日总体分类精度由81.83%上升到87.94%,增加了6.11个百分点,水稻制图精度由68.89%上升到80.11%,上升了11.22个百分点;8月17日总体分类精度由81.44%上升到91.58%,增加了10.14个百分点,水稻制图精度由83.18%上升到92.89%,上升了9.71个百分点;10月31日总体分类精度由87.73%上升到89.90%,增加了2.17个百分点,水稻制图精度由78.91%上升到81.00%,上升了2.09个百分点;11月15日总体分类精度由88.26%上升到89.11%,增加了0.85个百分点,水稻制图精度由82.66%下降到81.61%,降低了1.05个百分点。引入红边波段后,5个时相总体分类精度呈现上升的趋势,水稻制图精度除了11月15日其他时相均增加。从总体分类精度来看,8月17日最高为91.58%,8月2日最低为87.94%;从水稻制图精度来看8月17日最高为92.89%,5月24日最低为74.14%。最低的原因可能是水稻处于返青期,在遥感影像上水体与水稻光谱特征相似,水稻和水体出现严重错分的情况。
综上所述,8月17日影像水稻与其他地物间可分离性最好,总体分类精度、水稻制图精度最高,即在水稻生长拔节期能较好的识别水稻。原因是此时水稻位于水体低反射率与其他植被高反射率之间,与水体与其他植被均能区分开。
为了分析不同红边波段组合对水稻提取精度的影响,
由

图3 Sentinel-2卫星各波段在水稻分类中的重要性排序
8月17日单时相遥感影像德清县随机森林分类混淆矩阵如

图4 3个红边波段参与的7种典型地物分类
为了分析红边波段对水稻分类结果的影响,该文以Sentinel-2多光谱数据作为数据源,采用随机森林方法提取研究区内水稻分布,得到以下结论。
(1)7类典型地物的光谱特征在Sentinel-2影像上,红、绿、蓝波段各地物光谱差异小,难以区分,而红边波段、近红外波段各地物反射率变化则较为明显,为区分地物类型提供依据。
(2)与无红边波段相比,有红边波段的5个时相影像上水稻与其他地物之间的J-M距离均增加,其中8月2日林地—水稻、茭白—水稻J-M距离分别为1.976和2.000,可分离能力最优,8月17日其他植被—水稻、水稻—水体、建筑—水稻、大棚—水稻J-M距离均最大。从总体分类精度来看,8月17日最高为91.58%,8月2日最低为87.94%;从水稻制图精度来看8月17日最高为92.89%,5月24日最低为74.14%。由此得到8月17日影像水稻与其他地物间可分离性最好,总体分类精度、水稻制图精度也为最高,即在水稻生长拔节期时做适合进行水稻分类研究。
(3)Sentinel-2影像中3个红边波段不同组合参与分类,水稻识别能力Band6(红边波段2)>Band7(红边波段3)优于Band5(红边波段1);3个红边波段均参与分类时,总体分类精度最高,为91.58%,Kappa系数为0.89,水稻制图精度为92.89%;随机森林重要性排序结果与不同红边条件下分类精度比较的结果相吻合,即利用Sentinel-2影像提取水稻时,红边波段2、红边波段3发挥重要作用,需要着重关注。
(上接第143页)
目制定个性化的培养方案,并组织相关专家开展论证,在具体的课程操作中实现理论课与实践课的结合与衔接,将大学生培养成高技能的创新创业者。其次,借助特色小镇良好的创业氛围和优惠的配套政策和服务搭建创新创业实践平台,制定大学生创新创业训练计划,为大学生提供小微创业机会。各地特色小镇建设有政府相关优惠政策和创意经费、平台的支持,学校可对接特色小镇创客空间,为学生提供真正的实践机会和相应的创业支持,配备有实操经营管理经验的指导老师,帮助学生在创业中提升创新能力。同时各大院校可结合自身的办学理念和人才培养目标依托当地企业项目为不同专业学生创业寻找突破口,将学生自身的专业优势转化为小镇发展优势,实现学生、学校、企业及特色小镇建设的共赢。
推动创新创业教育,提升学生的创新创业能力是当前我国高等教育的现实需求,特色小镇的建设与大学生创新创业能力的培养相辅相成,相得益彰。《创新创业成功之道——基于梦想小镇发展实践的研究》一书介绍了我国经济转型升级背景下特色小镇建设的特点和价值,为各地特色小镇的建设指引了方向,也为大学生创新创业教育提供了新思路。
参考文献
胡琼, 吴文斌, 宋茜, 等. 农作物种植结构遥感提取研究进展. 中国农业科学, 2015, 48(10): 1900-1914. [百度学术]
陈仲新, 任建强, 唐华俊, 等. 农业遥感研究应用进展与展望. 遥感学报, 2016, 20(5): 748-767. [百度学术]
Salmon J M, Friedl M A, Frolking S, et al. Global rain-fed, irrigated, and paddy croplands: A new high resolution map derived from remote sensing, crop inventories and climate data. International Journal of Applied Earth Observation &Geoinformation, 2015, 38: 321-334. [百度学术]
王迪, 周清波, 陈仲新, 等. 基于合成孔径雷达的农作物识别研究进展. 农业工程学报, 2014, 30(16): 203-212. [百度学术]
Wardlow B, Egbert S, Kastens. Analysis of time-series MODIS 250 m vegetation index data for crop classification in the U.S. Central Great Plains. Remote Sensing Environment, 2007, 108: 290-310. [百度学术]
刘哲, 刘帝佑, 朱德海, 等. 作物遥感精细识别与自动制图研究进展与展望. 农业机械学报, 2018, 49(12): 1-12. [百度学术]
Shibayama M, Akiyama T.Seasonal visible, near-infrared and mid-infrared spectra of rice canopies in relation to LAI and above-ground dry phytomass. Remote Sensing of Environment, 1989, 27(2): 119-127. [百度学术]
Yang Y, Huang Y, Tian Q, et al. The extraction model of paddy rice information based on GF -1 satellite WFV images. Spectroscopy and Spectral Analysis, 2015, 35( 11) : 3255-3261. [百度学术]
Ozdogan M. The spatial distribution of crop types from MODIS data: Temporal unmixing using independent component analysis. Remote Sens. Environ . 2010, 114: 1190-1204. [百度学术]
Villa P, Stroppiana D, Fontanelli G, et al. In-season mapping of crop type with optica land x-band SAR data: A classification tree approach using synoptic seasonal features. Remote Sensing, 2015, 7: 12859-12886. [百度学术]
王建勋, 华丽, 邓世超, 等. 基于GF-1与MODIS时空融合的南方丘陵区水稻提取研究. 中国农业资源与区划, 2019, 40(5): 37-46. [百度学术]
Kontgis C, Schneide R A, Ozdogan M. Mapping rice paddy extent and intensification in the Vietnamese Mekong River Delta with dense time stacks of landsat data.Remote Sensing of Environment, 2015, 169: 255-269. [百度学术]
许青云, 杨贵军, 龙慧灵, 等.基于MODIS NDVI 多年时序数据的农作物种植识别. 农业工程学报, 2014, 30( 11) : 134-144. [百度学术]
张沁雨, 李哲, 夏朝宗, 等.高分六号遥感卫星新增波段下的树种分类精度分析.地球信息科学学报, 2019, 21(10): 1619-1628. [百度学术]
Christopher D E, Chen Z K. Comparison of broad-band and narrow-band red and near-infrared vegetation indexes. Remote Sensing of Environment, 1995, 54(1): 38-48. [百度学术]
Li X J, Chen G, Liu J Y, et al. Effects of rapid eye imagery's red-edge band and vegetation indices on land cover classification in an arid region. Chinese Geographical Science, 2017, 27(5): 827-835. [百度学术]
Kim H, Yeom J. Effect of red-edge and texture features for object-based paddy rice crop classification using RapidEye multi-spectral satellite image data. International Journal of Remote Sensing, 2014, 35(19): 7046-7068. [百度学术]
Yeom J M. Effect of red-edge and texture features for object-based paddy rice crop classification using RapidEye multi-spectral satellite image data. International Journal of Remote Sensing, 2014, 35(19): 7046-7068. [百度学术]
杨闫君, 占玉林, 田庆久, 等. 基于GF -1 /WFVNDVI 时间序列数据的作物分类.农业工程学报, 2015, 31( 11) : 155 -161. [百度学术]
刘佳, 王利民, 滕飞, 等.RapidEye卫星红边波段对作物面积提取精度的影响.农业工程学报, 2016, 32( 13) : 140-148. [百度学术]
谭玉敏, 夏玮.基于最佳波段组合的高光谱遥感影像分类.测绘与空间地理信息, 2014, 37(4): 19-22. [百度学术]
Pan Y, Li L, Zhang J, et al. Winter wheat area estimation from MODIS-EVI time series data using the Crop Proportion Phenology Index. Remote Sensing of Environment, 2012, 119: 232-242. [百度学术]
贺原惠子, 王长林, 贾慧聪, 等.基于随机森林算法的冬小麦提取研究.遥感技术与应用, 2018, 33(6): 1132-1140. [百度学术]
Radoux J, CHOMé G, Jacques D C, et al. Sentinel-2’s potential for sub-pixel landscape feature detection. Remote Sensing, 2016, 8( 6): 488. [百度学术]
Ghosh A, Sharma R, Joshi P K. Random forest classification of urban landscape using Landsat archive and ancillary data: Combining seasonal maps with decision level fusion. Applied Geography, 2014, 48( 2) : 31-41. [百度学术]
王娜, 李强子, 杜鑫, 等. 单变量特征选择的苏北地区主要农作物遥感识别. 遥感学报, 2017, 21(4) : 519-530. [百度学术]