Acreage estimation of mango orchards using object-oriented classification and remote sensing
RENChuanshuai1,3,, YEHuichun1,2, CUIBei1,2, HUANGWenjiang1,2 1. Key laboratory of Digital Earth Science,Institute of Remote Sensing and Digital Earth,Chinese Academy of Sciences,Beijing 100094,China2. Hainan Key Laboratory of Earth Observation,Sanya 572029,China3. University of Chinese Academy of Science,Beijing 100049,China 收稿日期:2017-02-10 修回日期:2017-07-3 网络出版日期:2017-08-20 版权声明:2017《资源科学》编辑部《资源科学》编辑部 基金资助:基金项目:海南省自然科学基金项目(2016CXTD015)海南省应用技术研发与示范推广专项(ZDXM2015102)海南省重大科技计划项目(ZDKJ2016021-02) 作者简介: -->作者简介:任传帅,男,山东济宁人,硕士,主要从事热带亚热带作物遥感监测研究。E-mail:rencs@radi.ac.cn
关键词:芒果林地;面向对象分类;植被覆盖度;数字高程模型;遥感提取 Abstract China has the second largest mango acreage in the word,and the planting area of mangos has followed a trend of sustained growth. Remote sensing has proved very useful in collecting information about the area of cereals and forest,but has rarely been applied to mango orchards. To utilize remote sensing in mango data collection,we attempted to estimate the acreage of mango orchards in China. Data from SPOT-6 and ASTER GDEM covering the study area was used. The boundary mask and sample points were used for acreage estimation. Vegetation coverage (FVC)and SLOPE factors were used to extract information about mango orchards using the object-oriented method. We found that mango orchard producer’s accuracy obtained by object-oriented classification and FVC and SLOPE data segmentation is as high as 92.81%;user accuracy is as high as 97.19%. Compared with maximum-likelihood classification and object-oriented classifi-cation without FVC and SLOPE data segmentation,object-oriented classification with FVC and SLOPE data segmentation has a higher accuracy for mango orchard extraction:producer accuracy increased by 19.28% and 3.62%,and user accuracy increased by 8.62% and 3.86%,respectively. Therefore,object-oriented classification with auxiliary data improves the accuracy of mango orchard extraction and this method provides an effective technical reference for other artificial orchards.
本研究选取海南省三亚市西部一块面积约45km2的区域作为研究区,区域内土地利用类型多样,主要包括芒果林地、耕地、建设用地、水体以及其它植被等。研究区地形复杂,主要包括山地、丘陵和平原等。研究区属于热带海洋性季风气候,全年平均气温26.6ºC,年极端高气温35.0ºC,年极端低气温14.6ºC,年平均降雨量1930.7mm[15]。近年来,三亚市芒果种植面积急速增长,据统计,20世纪90年代初三亚市芒果种植面积约800hm2,年产量约500t;2001年芒果种植面积增加到5407hm2,产量也增加到1.7万t;2013年,三亚市全境芒果种植面积已经达到10 065hm2,芒果年产量达12.4万t[15]。随着芒果产业的发展,三亚市芒果种植已经形成一定的区域布局,全市有梅山、三公里村、南山(抱古)、北岭、雅量等多个大型芒果种植基地。研究区具体位置见图1。 显示原图|下载原图ZIP|生成PPT 图1研究区地理位置及地面调查样点分布 -->Figure 1The study area and investigation sites -->
2.3.1 分类方法 (1)最大似然法分类 。最大似然法是最常用的遥感影像监督分类方法之一[26,27]。该方法的基本原理是通过对遥感数据像元统计分析,假定各类地物均为正态分布,通过训练样本,根据判决准则构建非线性判别函数,通过判定每个像元对各类别的归属度,把各像元分配到相应地物中去[28]。 (2)面向对象分类。多尺度分割是面向对象分类的关键,分割的合适与否直接影响到信息提取的精度[29,30]。图像分割即首先选定一个像元作为基准,根据相邻像元间的对比度和光谱异质性、光谱同质性以及异质性阈值对图像像元进行分割和合并,再利用形状异质性对产生的多边形进行修正,最终得到同质的、相互连接的不同斑块[31,32]。分割过程中为了充分利用研究数据的纹理信息,在对SPOT-6影像进行分割时,同时对FVC和SLOPE数据进行分割。经过多次对比分析,在分割尺度为50,合并尺度为80的情况下,芒果林地大多被分割到同一多边形,且零碎地块较少。并应用基于知识规则定义地物类别的面向对象分类方法,通过对实地考察样本的统计分析,建立分类规则对研究区进行分类。 2.3.2 研究方案 分别利用面向对象方法和最大似然法对研究区的芒果林地进行识别提取,通过对分类结果的精度验证,得到芒果林最优提取方法。在面向对象分类中,为了验证FVC和SLOPE因子参与的多尺度分割是否能进一步提高芒果林提取精度,将面向对象分类按FVC和SLOPE是否参与分割分为两种模式。整个研究方案分为以下三种模式: (1)基于最大似然法分类的研究区芒果林遥感提取,记为M1。 (2)基于面向对象分类,在影像分割时仅仅对SPOT影像进行多尺度分割,FVC和SLOPE数据不做分割处理的研究区芒果林遥感提取,记为M2。 (3)基于面向对象分类,在影像分割时对SPOT影像、FVC和SLOPE数据均进行多尺度分割的研究区芒果林遥感提取,记为M3。 具体技术路线如图2所示。 显示原图|下载原图ZIP|生成PPT 图2本研究的技术路线 -->Figure 2The technical route of this study -->
3 结果及分析
3.1 地物特征的选择和提取
通过对研究区的实地调查,发现研究区主要植被覆盖类型有芒果林地、其它林地及农田,其它土地利用类型包括水体、建筑物和裸地。由于本研究主要讨论芒果林地提取方法,而且建筑用地、裸地和水体的光谱信息明显异于植被,利用FVC可直接掩膜去除,故将建筑用地、水体和裸地合并为非植被地类,最终将分类类别设置为:非植被地类、耕地、其它林地和芒果林地四大类。通过对不同地类样点数据统计,生成芒果林地与其它林地FVC分布概率密度柱状图(图3)以及芒果林地和耕地坡度分布概率密度柱状图(图4)。 显示原图|下载原图ZIP|生成PPT 图32013年12月研究区芒果林地与其它林地植被覆盖度(FVC)的概率密度对比 -->Figure 3The probability density histogram of fractional vegetation cover (FVC)of mango fields and other forests in December 2013 -->
从图3的统计结果可以看出,大部分芒果林地的FVC值集中在0.3~0.8区间内。芒果林地概率密度分布形式符合正态分布,其中在0.5~0.6区域达到最大值,之后随FVC值的增加而逐渐减小。相比而言,其它林地FVC值多集中在0.7~1.0之间。芒果林地与其它林地FVC值在0.7~0.8区间内虽有部分交集,但所占比例较小,故利用FVC设定阈值,可较好地区分芒果林地与其它林地。从图4统计结果中可以看出,大部分耕地的坡度值集中在0~2°区间内,芒果林地坡度范围较大,但多集中在4°~25°区间内,两者之间交叉范围较小,故利用坡度因子设定阈值,可较好的区分芒果林地和耕地。由此可生产以下分类规则。 显示原图|下载原图ZIP|生成PPT 图42013年12月研究区芒果林地与耕地在不同坡度的概率密度对比 -->Figure 4The probability density histogram of slope of mango fields and cultivated lands in December 2013 -->
(1)水体、建设用地和裸地等非植被地类利用FVC进行掩膜处理直接去除。 (2)耕地多种植在地势平缓的平原地带,耕地类型主要分为有农作物覆盖的耕地及收割后的耕地两大类。对于收割后耕地,由于其光谱特性更接近于裸露地表,故将其归为裸地。由于耕地和芒果林地坡度相差大,故可利用坡度信息来区分有作物覆盖耕地与林地,根据样点统计数据将坡度小于2.125°的植被分为耕地。 (3)其它林地主要包括橡胶园地及热带雨林。通过对研究区调查发现,其它林地植被覆盖度较高,在融合后的SPOT-6(1.5m分辨率)影像中较难分辨出林间裸露地表,林地郁闭度较高。芒果林属于人工种植的乔木林,作为经济林的一种,芒果林地有明显区别于其它林地的特征。芒果林种植密度不易过大,一般株距3.5m,行距4.5m[33],故在遥感影像上有明显的线性纹理特征。芒果林地由于郁闭度较低,其FVC低于其它林地,故可利用FVC来区分芒果林地与其它林地,根据统计数据将0.773设定为芒果林地与其它林地的分割阈值。 根据上述可分性研究以及分类规则,利用统计数据得到的合理阈值对研究区进行分类。最终得到M3模式下芒果林地信息提取结果。同时应用M1模式和M2模式对研究区进行分类,以作为对照。3种分类方法的芒果林地提取结果如图5所示。 显示原图|下载原图ZIP|生成PPT 图5基于不同模式的研究区土地利用类型分类结果 -->Figure 5The land use classification results of this study area based on different classification modes -->
3.2 不同方法的分类结果对比分析
通过对比三种不同分类模式,可以明显发现面向对象分类方法所获得的结果精度明显高于最大似然监督分类法的结果。为进一步量化分类精度,根据实地考察数据从研究区内选取验证样本对分类结果进行验证。三种分类方法对应芒果林地提取精度见表2。 通过表2可以看出,M1模式提取芒果林地生产者精度为73.53%,用户精度为88.57%。根据实际调查,研究区内耕地集中在平原地带,而应用最大似然法分类,部分山区丘陵地带的芒果林地被错分为耕地,且芒果林地与其它林地存在严重的混分现象。这说明仅仅依靠光谱信息进行计算机自动分类并不能很好地提取芒果林地。M2模式提取芒果林地的精度有显著提高,相较于M1模式生产者精度提高了15.66%,用户精度提高了4.76%,有效地避免了芒果林地被错分做其它植被的现象。而将FVC和SLOPE因子参与到影像分割的M3模式能进一步提高芒果林地的提取精度,相较于M2模式生产者精度提高3.62%,用户精度提高了3.86%,该方法能很好地解决芒果林地的错分和漏分问题。 Table 2 表2 表2基于不同分类方法的芒果林地分类精度 Table 2Classification accuracy of mango based on different classification methods
利用遥感影像监督分类从植被物种复杂区域提取出单一植被物种仍比较困难。在人工园林提取方面,综合考虑人工园林的各种特征(如种植密度、坡度等)对于提高人工园林的提取精度有很大帮助。本文基于高分辨率卫星SPOT-6影像数据和面向对象分类方法,通过引入能反映芒果林地种植密度的植被覆盖度因子和能反映芒果林地地形特征的坡度因子作为辅助数据,提取芒果林地的生产者精度达92.81%,用户精度达97.19%,较最大似然法分类以及传统面向对象分类结果在精度上均有了较大提高。说明在人工园林遥感信息提取方面,引入目标地物自身特征及环境因素作为辅助分类信息,一定程度上可以提高遥感信息的提取精度。当然,本研究仍存在很多不足,仅利用了光谱、植被覆盖度、地形等因子来提取芒果林地信息,对于芒果林地在遥感影像上所表现出的特有的纹理、形状、面积等特征仍考虑不足。因此,如何结合目标地物的光谱、纹理、形状、面积等特征以及环境信息来进行分类参数的优选将会是今后研究的重点。 The authors have declared that no competing interests exist.
[ Huang GD,Su MH,Wang CT.Present situation and development countermeasures of mango standardization in China [J]. ,2010,(2):19-22.] [本文引用: 1]
[2]
JulienY,Sobrino JA,Jiménez-Muñoz J C. Land use classi-fication from multitemporal Landsat imagery using the Yearly Land Cover Dynamics(YLCD)method [J]. ,2011,13(5):711-720. [本文引用: 1]
[Ma YY,Zhang CX,Zhang JC,et al. Research on object-oriented classification method assisted with NDVI/SLOPE in extracting cassava:taking Wuming county for example [J]. ,2015,31(1):49-53.] [本文引用: 1]
[Shen JL,Yu WY,Qi XP,et al. The ecosystem research on Momoge conservation zone based on remote sensing data and object-oriented classification [J]. ,2008,(6):71-75.] [本文引用: 1]
[7]
GaoY,Mas JF,Maathuis B H P,et al. Comparison of pixel-based and object-oriented image classification approaches-a case study in a coal fire area,Wuda,Inner Mongolia,China [J]. ,2006,27(18):4039-4055.
[Wang HJ,Le CF.Extraction of paddy field using objects-oriented fuzzy classification approach based on texture characteristic [J]. ,2008,24(5):40-43.] [本文引用: 1]
[12]
San YJ,Tong QX,Qin QM.The object-oriented method for wetland information extraction [J]. ,2008,28(1):79-82. [本文引用: 1]
[13]
WrightC,GallantA.Improved wetland remote sensing in Yellowstone National Park using classification trees to combine TM imagery and ancillary environmental data [J]. ,2007,107(4):582-605.
[14]
Frohn RC,ReifM,LaneC,et al. Satellite remote sensing of isolated wetlands using object-oriented classification of Landsat-7 data [J]. ,2009,29(3):931-941. [本文引用: 1]
[Wang XP.Geometric rectification and effect analysis of remote sensing image [J].,2008,25(9):102-105.] [本文引用: 1]
[18]
Chen SH,SuH,ZhangR,et al. Fusing remote sensing images using à trous wavelet transform and empirical mode decom-position [J]. ,2008,29(3):330-342. [本文引用: 1]
[19]
AiazziB,BarontiS,SelvaM.Enhanced gram-schmidt spectral sharpening based on multivariate regression of MS and pan data [J]. ,2006,6:3806-3809. [本文引用: 1]
[WangL,He XB,Bao YH,et al. Error analysis on gradient information of sloping land derived from the low-resolution DEM [J].,2010,38(7):3644-3646.] [本文引用: 1]
[Feng LL,He ZM,Liu XF,et al. Dynamic monitoring of the fractional vegetation cover in Jilin Province based on MODIS-NDVI data [J]. ,2014,31(4):492-499.] [本文引用: 1]
[QiaoF,Zhang KB,Zhang SY,et al. RS monitoring on the dynamic change of vegetation coverage in a farming-pastur-ing ecotone-a case study in Yanchi county,Ningxia Hui Auto-nomous Region [J]. ,2006,23(2):283-288.] [本文引用: 1]
[26]
Otukei JR,BlaschkeT,WoldaiT,et al. Land cover change assessment using decision trees,support vector machines and maximum likelihood classification algorithms [J]. ,2010,12(1):S27-S31. [本文引用: 1]
[27]
Soheili MajdM,SimonettoE,PolidoriL.Maximum likelihood classification of high-resolution SAR images in Urban Area [J]. ,2012,(4):395-407. [本文引用: 1]
[Xu WB,LinN,Lu WH,et al. Maximum likelihood method for the extraction of monitoring and surveying information on Diaoyu Dao by aerial remote sensing [J]. ,2013,32(5):548-552.] [本文引用: 1]
[FanL,Cheng YZ,Wang LG,et al. Estimation of winter wheat planting area using object-oriented method based on multi-scale segmentation [J]. ,2010,31(6):44-51.] [本文引用: 1]
[QiaoT,Zhang HQ,Chen YF,et al. Extraction of vegetation information based on NDVI Segmentation and object-oriented method [J]. ,2013,28(4):170-175.] [本文引用: 1]
[ Bao HY,LiY,ZhaoP.The research of object-oriented classification method in remote sensing image with texture analysis [J]. ,2009,(4):33-37.] [本文引用: 1]
[Wang HJ,Le CF.Extraction of paddy field using objects-oriented fuzzy classification approach based on texture characteristic [J]. ,2008,24(5):40-43.] [本文引用: 1]