基于GF-7卫星多角度特征进行作物识别

基于GF-7卫星多角度特征进行作物识别 Crop recognition by multi-angle features of GF-7 satellite 2 first-author 孙智虎，研究方向为农业遥感。E-mail： 1045376795@qq.com 孙智虎，研究方向为农业遥感。E-mail： 1045376795@qq.com 1045376795@qq.com 1 3 4 corresp 介绍：张锦水，主要从事深度学习与遥感大数据挖掘研究。E-mail： zhangjs@bnu.edu.cn 介绍：张锦水，主要从事深度学习与遥感大数据挖掘研究。E-mail： zhangjs@bnu.edu.cn zhangjs@bnu.edu.cn 2 1 3 4 5 遥感科学国家重点实验室 北京师范大学, 北京100875 遥感科学国家重点实验室 北京师范大学, 北京100875 State Key Laboratory of Remote Sensing Science , Beijing Normal University, Beijing 100875, China State Key Laboratory of Remote Sensing Science , Beijing Normal University, Beijing 100875, China 中国地质大学(北京)土地科学技术学院, 北京 100083 中国地质大学(北京)土地科学技术学院, 北京 100083 School of Land Science and Technology, China University of Geosciences, Beijing 100083, China School of Land Science and Technology, China University of Geosciences, Beijing 100083, China 北京市陆表遥感数据产品工程技术研究中心, 北京100875 北京市陆表遥感数据产品工程技术研究中心, 北京100875 Beijing Engineering Research Center for Global Land Remote Sensing Products, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China Beijing Engineering Research Center for Global Land Remote Sensing Products, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China 北京师范大学地理科学学部遥感科学与工程研究院, 北京 100875 北京师范大学地理科学学部遥感科学与工程研究院, 北京 100875 Institute of Remote Sensing Science and Engineering, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China Institute of Remote Sensing Science and Engineering, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China 北京工业职业技术学院, 北京 100042 北京工业职业技术学院, 北京 100042 Beijing Polytechnic College, Beijing 100042 Beijing Polytechnic College, Beijing 100042 多角度遥感对地观测能够提供更加丰富、多方向的遥感特征，提高地类之间的可区分性，为地物覆盖的精确识别打下坚实的数据基础。GF-7是我国继ZY-3号卫星后的首颗亚米级测绘卫星，这为利用多角度特性解决“异物同谱”的问题，提高作物的识别精度带来了机遇。本文利用GF-7号前视、后视全色及后视多光谱数据，各种特征组合输入到支撑向量机分类器进行分类，相对于光谱、纹理等特征，分析多角度特征对作物识别精度的作用。结果表明，较仅应用光谱特征，光谱与角差特征组合使用大蒜和冬小麦的制图精度分别提高了4.07%和3.15%，用户精度分别提高了6.73%和2.12%；较应用光谱与纹理特征，光谱、纹理与角差特征组合使用大蒜和冬小麦的制图精度分别提高了3.14%和1.01%，用户精度分别提高了5.11%和0.67%。通过McNemar检验分析，这种分类精度的提高是稳定的，角差特征使用能有效提高作物的识别精度。究其原因，多角度特征对不同作物类型在多角度观测时的光谱响应具备特有的差异性，这种差异提高了作物之间的可分性，从而保证作物遥感识别的精度。 Multi-angle remote sensing can provide richer and multi-directional features for ground objects observation, improve the distinguishability between land types, and lay a solid data foundation for the accurate identification of ground cover. GF-7 is the first domestic sub meter surveying and mapping satellite after ZY-3 satellite, which brings an opportunity to solve the problem of "foreign matter homospectrum" using multi-angle characteristics and improve the identification accuracy of crops. This paper use GF-7 forward-looking, backward looking panchromatic and backward looking multi-spectral data, and various features combinations are input to the SVM (support vector machine) classifier for classification, used to analyze the influence of multi-angle features on crop recognition accuracy relative to spectral and texture features. The results show that, compared with only spectral features, add the angle difference feature, the production accuracy of garlic and winter wheat is increased by 4.07% and 3.15%, and the user accuracy is increased by 6.73% and 2.12%, respectively. Compared with the combination of spectral and texture features, add the angle difference feature, the production accuracy of garlic and winter wheat is increased by 3.14% and 1.01%, and the user accuracy increased by 5.11% and 0.67%, respectively. Through the analysis of McNemar test, the improvement of classification accuracy is stable, angle difference feature can effectively improve the identification accuracy of crops. Tracing it to its cause, the multi-angle characteristics of GF-7 satellite have unique differences in the spectral response of different crop types during multi-angle observation. The difference improves the separability between crops so as to ensure the accuracy of crop remote sensing mapping. GF-7 支撑向量机 角差 遥感 冬小麦 大蒜 GF-7 SVM angle difference remote sensing winter wheat garlic 1　引 言 1 2019年11月7日我国成功发射的GF-7是继ZY-3卫星后又一多角度立体测绘卫星，空间分辨率达到亚米级，具有1：10000比例尺立体测绘能力，能观测到更为丰富的地物立体特征并可用于高精度的地形数字产品生产和土地覆盖制图，在农业监测、林业管理、城市规划等方面具有广阔的应用前景。 基于高分辨率遥感影像光谱和高清特性进行高精度农作物制图一直是遥感应用领域的重点研究方向。一般而言，作物分类都是需要利用遥感图像的光谱、空间等特征进行作物的精准识别。如，宋德娟等（2018）基于GF-2影像构建WWRSE深度学习模型，准确提取出冬小麦的种植分布； 吴锦玉等（2019） 吴锦玉等（2019） 利用GF-2号的光谱特征、植被指数与多种纹理特征能够对塑料大棚进行识别； Zhang等（2020） Zhang等（2020） 利用GF-1等米级影像构建MPSPNet网络模型进行大尺度耕地分布提取。这些研究均表明米级/亚米级影像空间特征是能够准确表达作物的独有特征，为作物的精准识别提供基础。 高分辨率立体测绘卫星多角度对地观测能够反演地形数据作为辅助信息支撑地物识别，或构建多角度遥感特征参与分类。如， Tian等（2013） Tian等（2013） 利用Cartosat-1影像生成的DSM（Digital Surface Model）进行工业区、林区的变化检测提取。 Li等（2016） Li等（2016） 利用ZY-3卫星影像生成地形数据，进行矿区和农业制图。 Huang等（2014） Huang等（2014） 利用ZY-3提取的DSM及形态学等特征进行城市的变化检测识别。 Matasci等（2015） Matasci等（2015） 利用WorldView-2多角度序列影像，分析影像拍摄角度对城市地区分类影响。这些研究是通过DSM表达出地表要素的高度差异作为独立特征变量支撑信息的提取。另外， Huang等（2018） Huang等（2018） 利用ZY-3影像，构建多级角差特征支撑城市场景分类，用于提高城市地区的分类精度。通过这些研究证实，通过立体像对提取DSM或利用多角度观测结果生成的空间结构特征能够提高地物的可辨识能力，是光谱、纹理等特征的有力补充。 Khatami等（2016） Khatami等（2016） 总结了在光谱分类的基础上纹理特征能将总体分类精度提高12.1%，多角度特征平均提高8%，表明多角度特征对地物识别的贡献。因此，本文针对GF-7卫星，研究多角度遥感特征对农作物精确识别的贡献，为促进GF-7影像多角度特征应用到作物遥感制图中积累研究基础。 本研究利用GF-7影像多光谱、前后视影像进行冬小麦、大蒜的识别，重点探讨三个方面的内容： 1. 分析GF-7卫星影像光谱、纹理和角差特征对地物的表达能力和可分性，重点分析角差特征与其他特征的区别。 2. 利用支撑向量机（Support Vector Machine， SVM）分类器进行作物识别，分析不同遥感特征对作物识别的贡献，尤其是角差特征对作物识别精度的作用。 3. 分析作物分类结果的不确定性，评估不同特征作物空间误差的表达，探讨作物专题分类结果不确定性与分类误差的内在关系，探讨角差特征与分类结果不确定性之间的关系。 2　研究区与数据 1 2.1　研究区概况 2 研究区位于华北平原中部，山东与河南的交界处，如 图 1 图 1 所示。该地区属于温带大陆性季风气候，夏热冬冷，降水集中在夏季，越冬作物主要是冬小麦和大蒜，夏季作物以玉米为主。对于冬小麦和大蒜，属同期作物，单用光谱等特征进行识别存在相混，这种现象在我国华北地区常见，本文探讨利用GF-7多角度特征解决这一难题。 2.2　GF-7影像介绍及预处理 2 GF-7卫星搭载1台双线阵相机和1台激光测高仪，其中双线阵相机具有优于0.8米空间分辨率的前视、后视全色和优于3.2米空间分辨率的后视多光谱成像能力，激光测高仪配置有2波束的1064nm激光用于获取地面30m光斑范围内的高程。GF-7单景遥感影像幅宽20km，重访周期为60天。 卫星拍摄时太阳、前视影像、后视影像的方位角和天顶角关系如 图 2 图 2 所示。 本研究影像（ 图 1 图 1 ）获取时间为2019年11月13日。以Google Earth影像为地理参考影像，目视选取控制点，将后视全色影像进行几何配准，随后将后视多光谱影像和前视全色影像均配准到后视全色影像上，配准误差小于0.5个像素。 李霖等（2014） 李霖等（2014） 、 刘会芬等（2014） 刘会芬等（2014） 指出Gram-Schmidt融合方法（ Clayton等，1971 Clayton等，1971 ）能够有效地继承多光谱影像的光谱信息和全色影像的空间信息，有利于地物的分类，因此本文选用Gram-Schmidt方法进行后视全色与后视多光谱影像数据的融合。本文后面作物分类所采用的光谱和纹理特征均来自于融合影像。 2.3　野外调查 2 为掌握研究区的作物种植情况，2019年11月30日-12月1日进行了野外调查，采用艾尔思公司开发的“AW拍”手机APP（iOS、安卓版本可从应用市场下载）采集了682个样本点和对应的现场照片，其中冬小麦296个，大蒜226个，其他作物34个，林地26个，草地20个，裸地58个，其他22个。基于采样点位置信息及照片信息，共数字化出36个野外调查样方（见 图3 图3 ），其中351个冬小麦、216个大蒜、842个其他地块（包括建筑、裸地、其他作物、林地、草地和水域），地块分布见 图 1 图 1 。以矢量地块范围和属性为基础，构建作物分类的检验样本，样本数量和类型见 表1 表1 。 分类类别 描述 训练样本像元个数 检验样本像元个数 冬小麦 越冬农作物，一般在10月份种植，种植跨度时间较长 9000 5443404 大蒜 与冬小麦种植时间相近，地表有塑料薄膜覆盖 7000 982025 其他 包括建筑、水域、林地、草地、裸地、蔬菜等 25800 2738733 通过野外调查了解到，该研究区内的冬小麦播种时间跨度时间较长，绝大部分冬小麦处于分蘖期（如 图 3 图 3 （b）、（c）），仍有少部分冬小麦处于出苗期。大蒜是种植面积仅次于冬小麦的农作物，在调查期内其表面覆有白色地膜，如 图3 图3 （a）所示。从调查图片可以看出，大蒜在这一时期的高度与冬小麦相差不多，不同之处在于大蒜垄间间隔较冬小麦宽，背景裸土信息暴露更多，且有地膜附着地面。 3　研究方法 1 3.1　技术流程 2 本文技术流程如 图 4 图 4 所示，主要包括：影像预处理、特征和训练集测试集构建、不同特征组合的SVM分类、精度评价以及分类结果不确定性的定量评价。 3.2　样本集构建 2 该研究区主要地物类型有冬小麦、大蒜、蔬菜（大葱、花椰菜、甘蓝等）、建筑（居住地、道路、蔬菜大棚、厂房等）、裸地（裸露的地表和休耕地）、水域（河流和池塘等）和草地，各类型在影像上呈现的特征如 图 5 图 5 所示。 已有研究表明，遥感制图精度高低与分类器和训练样本相关（ Khatami等，2016 Khatami等，2016 ）。为保证样本的代表性，在外业实测样本的支持下，我们以10*10像素框在影像上均匀地选择出解译单元进行地物类型的识别，构建训练样本，其像元数量见 表1 表1 。分类体系定义为冬小麦、大蒜和其他（包含其他作物、林地、草地、裸地、水域和建筑）。 3.3　特征数据集 2 本研究中综合考虑光谱、纹理和多角度特征进行作物分类，分析不同特征对作物分类精度的影响。 3.3.1　角差特征 3 各种土地覆盖在太阳光照射下会呈现出一定程度的光谱反射各向异性（ Su等， 2007 Su等， 2007 ）。GF-7卫星从前、后视角度观测地表，依据双向反射分布函数（bidirectional reflectance distribution function，BRDF）（ Jean-Louis等，1992 Jean-Louis等，1992 ； Schaepman-Strub等，2006 Schaepman-Strub等，2006 ）的原理可知，多角度影像能够反映出地物的双向反射差异特征差异，以及地物的结构信息（ Chopping等，2008 Chopping等，2008 ），这一特征信息称为角差特征。角差特征是基于同种地物具有稳定的BRDF，而不同空间结构的地物BRDF有所不同，与光谱特征一样，是地物区分其他地物特有的特性（ Huang等，2018 Huang等，2018 ）。因此，本文将GF-7号卫星影像提取的角差作为一个衡量地物的空间维度特征，见 公式（1） 公式（1） ： <math display="block" id="M1"><mi>A</mi><mi>n</mi><mi>g</mi><mo>_</mo><mi>D</mi><mi>i</mi><mi>f</mi><mo>=</mo><mi>F</mi><mi>W</mi><msub><mrow><mi>D</mi></mrow><mrow><mi>p</mi><mi>a</mi><mi>n</mi></mrow></msub><mo>-</mo><mi>B</mi><mi>W</mi><msub><mrow><mi>D</mi></mrow><mrow><mi>p</mi><mi>a</mi><mi>n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28331941&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28331940&type= 42.41799927 4.57200003 其中， FWD pan 为前视全色影像， BWD pan 为后视全色影像。 3.3.2　纹理特征 3 灰度共生矩阵是由观测窗口、像元移动步长和方向、分布概率的灰度等级定义的一种统计分析方法（ Lan和Liu，2018 Lan和Liu，2018 ），是纹理特征的量化指标（ Haralick，1979 Haralick，1979 ）。 本文为消除特征维度冗余带来的数据维度灾难（ 杨闫君等，2018 杨闫君等，2018 ），对融合影像进行主成分分析，第一主成分波段分量占影像信息的78.05%，因此这一分量信息能够代表整个图像，进一步提取第一主成分波段的8种独立纹理特征，包括均值（Mean）、方差（Variance）、均一性（Homogeneity）、对比度（Contrast）、相异性（Dissimilarity）、信息熵（Entropy）、二阶矩（Second Moment）和相关度（Correlation）。 朱秀芳等（2019） 朱秀芳等（2019） 指出纹理窗口对地物纹理特征表达影响较大。本文设置了3*3、5*5、7*7、9*9、11*11、13*13、15*15的不同窗口提取纹理特征，进行SVM构建和作物识别，并使用检验样本评估冬小麦和大蒜的精度（ Congalton等，1991 Congalton等，1991 ）。采用两种作物制图精度与用户精度之和的均值作为评价指标来确定最优窗口。经测试，3*3窗口的纹理特征是用于作物分类最优的（见 图 6 图 6 ）。 3.3.3　特征组合 3 光谱、角差和纹理是从不同维度表达作物特征， 表2 表2 列出了本文所采用的所有特征。采用不同的特征组合进行作物分类以验证不同特征对作物分类结果的影响是常用的策略（ 于婉婉等，2021 于婉婉等，2021 ）。本文以光谱为基础特征，分别分析纹理特征和角差特征对分类精度的影响，特征组合包括：A、A+B、A+C和A+B+C。 代码 分类 特征 A 光谱特征 蓝光、绿光、红光、近红外 B 角差特征 角差 C 纹理特征 均值、方差、均一性、对比度、相异性、信息熵、二阶矩、相关性 3.4　SVM参数优化 2 SVM常应用于图像识别领域，通过线性（linear）、多项式（polynomial）、sigmoid、rbf（Gaussian）等核函数（Chang和Lin，2011）将数据映射到高维空间，对数据集进行线性\非线性划分，构建几何间隔最大的超平面，找到各类别之间的最优分界线以期最大程度划分数据集，在农作物遥感分类中得到广泛的应用（ Khatami等，2016 Khatami等，2016 ）。 本文基于Python scikit-learn库进行模型训练和预测推断，采用网格搜索和10折交叉验证方法（ Hsu等，2003 Hsu等，2003 ）来寻找RBF核函数的最优参数cost和gamma（ Li等，2016 Li等，2016 ； Duro等，2012 Duro等，2012 ），结果见 表3 表3 。 特征集 cost gamma 总体精度 A 32768 8 90.13% A+B 32768 8 93.49% A+C 32768 2 94.09% A+B+C 32768 2 95.1% 3.5　精度评价 2 以数字化的野外调查样方结果作为检验样本（ 表1 表1 ），与分类结果进行对比，构建混淆矩阵（ Congalton等，1991 Congalton等，1991 ），评估不同特征组合的SVM分类结果的总体精度、Kappa系数，冬小麦和大蒜的制图精度、用户精度。 3.6　McNemar检验 2 McNemar检验是基于自由度为一的卡方检验（ Waldner等，2015 Waldner等，2015 ），比Z-test具有更高的精度和敏感性，适用于相同数据（ Alshehhi 等，2017 Alshehhi 等，2017 ）检验来检验实验结果准确性差异是否具有统计学意义（ Waldner等，2015 Waldner等，2015 ），也就是要回答“方法A识别精度比方法B精度高，是否通过显著性检验”，而不是简单利用精度指标的高低评价方法之间的优劣。 McNemar检验是基于不同特征组合分类器的分类结果对各像素的一致性假设，即判断各分类器是否将同一单个像素划分为同一类别，计算公式如下： <math display="block" id="M2"><msubsup><mrow><mi>χ</mi></mrow><mrow><mn mathvariant="normal">0.95,1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo>=</mo><mfrac><mrow><msup><mrow><mo stretchy="false">(</mo><mfenced open="|" close="|" separators="|"><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">12</mn></mrow></msub><mo>-</mo><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">21</mn></mrow></msub></mrow></mfenced><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">12</mn></mrow></msub><mo>+</mo><msub><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">21</mn></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28331810&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28331949&type= 37.50733566 10.15999985 针对对应的每个像元，其中， n 12 代表组合1分类器正确分类但组合2分类器错误分类的像元个数， n 21 正好相反。 本文将McNemar检验用于在95%置信度下检验不同特征组合的农作物识别结果是否有显著性差异，从而证实方法间分类精度差异的可行度。 3.7　分类不确定性表达 2 上面的精度评价方法能够表达出专题图的整体分类精度（ Foody，2002 Foody，2002 ），但无法体现出空间误差的分布。信息熵可以综合不同来源误差，用一个数字表达不确定性的分布和程度（ Dehghan 和 Ghassemian，2006 Dehghan 和 Ghassemian，2006 ），信息熵被用来量化表达遥感分类结果的不确定性（ Foody G M和Atkinson P M，2002 Foody G M和Atkinson P M，2002 ； Dehghan 和 Ghassemian， 2006 Dehghan 和 Ghassemian， 2006 ）。本文针对不同分类特征组合的不确定性结果，从空间异质性上分析各种特征对分类结果精度的影响，从中分析作物精度影响的内在机理。 分类器对于最终类别的归属概率可以描述分类的不确定，但这忽略了其他类别的归属概率（ Waldner等，2015 Waldner等，2015 ）。 Dehghan和Ghassemian（2006） Dehghan和Ghassemian（2006） 提出利用熵来综合反映分类结果的不确定性。按照 公式（3） 公式（3） 计算每一个分类像元x的熵： <math display="block" id="M3"><mi>E</mi><mi>n</mi><mi>t</mi><mi>r</mi><mi>o</mi><mi>p</mi><mi>y</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>n</mi></mrow></msubsup><mrow><mo>-</mo><mi>p</mi><mi>i</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mi>*</mi></mrow><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><msub><mrow><mi mathvariant="normal">g</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mi>p</mi><mi>i</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28331812&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28331811&type= 55.20266724 4.82600021 其中， n 代表分类的类别数， pi（x） 代表像元为第 i 种地物类型的概率。 本文进一步采用 Waldner等（2015） Waldner等（2015） 提出的归一化不确定性（Normalized Uncertainty Criteria， NUC）计算方法来表达像元的不确定性，如 公式（4） 公式（4） 所示： <math display="block" id="M4"><mi>N</mi><mi>U</mi><mi>C</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo><mo>=</mo><mn mathvariant="normal">1</mn><mo>-</mo><mfrac><mrow><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><msub><mrow><mi mathvariant="normal">g</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">(</mo><mi>M</mi><mo stretchy="false">)</mo><mo>-</mo><mi>E</mi><mi>n</mi><mi>t</mi><mi>r</mi><mi>o</mi><mi>p</mi><mi>y</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><msub><mrow><mi mathvariant="normal">g</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">(</mo><mi>M</mi><mo stretchy="false">)</mo></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28331814&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=28331813&type= 56.72666931 9.14400005 公式（4） 公式（4） 中 M 代表分类规则所划分类别的总数。 NUC （ x ）值越高，分类结果越不确定。 4　结果分析 1 4.1　特征差异分析 2 4.1.1　角度差异分析 3 依据各地物类型角差值的概率分布图（ 图 7 图 7 ），分析角差特征对冬小麦、大蒜、其他作物、林地、草地、裸地、建筑和水域等地物的表达能力。大蒜和水体角差值几乎全部大于0，这是卫星前视观测接收到太阳光的反射要高于后视观测，大蒜和水体表面有一定的镜面反射特性。 刘小珊（2015） 刘小珊（2015） 指出，对于土壤，当探测器方向与光源方向相同，探测器接收的能量是光源直接照射部分，土壤颗粒之间的阴影很少，后向反射能量高于前向反射能量，因此其角差值为负，这对于冬小麦、其他作物、林地、草地的角差值大部分为负是同理的。这样的多角度反射差异与太阳入射方向、卫星接收高度角以及地物类型所特有的表面结构特性有关。 大蒜的角差值远高于冬小麦、其他作物、林地、草地和裸地，大蒜与这些地物类型的角差值分布几乎没有重叠，这是由于大蒜种植区域覆有白色地膜具有较强的镜面反射特性，从而使得大蒜种植区域与其他地物具有显著的角差反射差异，可见通过角差特征可有效地区分大蒜与这几种地物。大部分情况下，冬小麦的角差值大于裸地、其他作物和草地，小于大蒜和水域，在分布上没有或仅有少部分重叠，这是通过角差特征将冬小麦与其他几种地物有效区分的基础。 大蒜角差值分布与水域和建筑有较多的重叠，冬小麦角差值与林地、建筑有较多的重叠部分，若仅通过角差特征识别农作物，则这些地物之间会产生一定的混淆。 冬小麦和大蒜角差值分布集中，说明研究区冬小麦和大蒜具有相对稳定的BRDF，类内稳定的角差特征有利于作物识别。林地、建筑和水域角差值分布离散，是由于其类内形态有较大差异，特征信息表达存在波动性。 裸地的角差值最小，与其余地物存在明显的不同。其他作物与草地角差值分布大部分都有重叠，角差特征难以区分这两种地物。建筑区域的角差值集中在0值附近分布，且多为正值，这可能是建筑物三维特征更为明显导致的。 4.1.2　空间特征差异分析 3 依据min-max标准化方法计算特征归一化值，统计得到地物类型的特征均值、标准差图（ 图 8 图 8 ），以分析不同地物空间特征的可分性。 图 8 图 8 中箱型图横线值为均值，上下端与中间横线的差距表示标准差，图形上下幅度差异越大异质性越强。其他作物、林地、水体和建筑在各特征上标准差较大，是由于它们的类内形态具有较大的差异。 冬小麦与大蒜在蓝光、绿光和红光波段上有显著的差异，说明冬小麦与大蒜光谱上具有一定的可分性；冬小麦的均值特征大部分小于大蒜和其他作物，卫星接收到冬小麦种植区域的地物反射能量较小；近红外特征可反映农作物的生长旺盛程度，冬小麦的近红外特征小于其他作物；冬小麦与林地在光谱特征上混淆较大，而在均一性、相异性、二阶矩和信息熵特征有显著差异。 图 8 图 8 显示出大蒜在绿光特征与林地、草地仅有少部分重叠，绿光特征可将大蒜与林地、草地区分；大蒜与各类土地覆盖物在8种纹理特征上大都有较多的重叠，这种混淆是和当地作物种植环境及作物生长阶段相关，因此加入纹理特征会造成信息冗余将会影响大蒜的识别精度。 4.2　分类结果及精度评价 2 图 9 图 9 是不同特征组合分类结果及其不确定性评估结果，用以整体显示不同特征尤其是角差特征对分类结果和不确定性的影响。 图 9 图 9 中不确定性评估结果呈现出明显的边界效应，这可能是由于本文研究的GF-7影像处于在轨测试时期，其数据存在部分问题。 表 4 表 4 为四种特征组合（见2.3.3章节）分类的混淆矩阵，并以此计算出分类的总体精度、Kappa系数以及冬小麦和大蒜的制图精度、用户精度，得到 图 10 图 10 分类精度对比图。 13 特征集SVM分类结果 3 A 3 A+B 3 A+C 4 A+B+C 冬小麦 大蒜 其他 冬小麦 大蒜 其他 冬小麦 大蒜 其他 冬小麦 大蒜 其他 3 真实地类属性 冬小麦 4398238 133642 911693 4569971 86554 787048 4641864 154040 647669 4696752 98100 648721 大蒜 83734 823347 74963 43388 863350 75306 57718 818289 106037 33002 849151 99891 其他 323289 67218 2343841 266230 41054 2427064 191909 53666 2488773 184827 53312 2496209 制图精度 80.80% 83.84% 85.72% 83.95% 87.91% 88.76% 85.27% 83.33% 91.02% 86.28% 86.47% 91.29% 用户精度 91.53% 80.39% 70.38% 93.65% 87.12% 73.78% 94.90% 79.76% 76.76% 95.57% 84.87% 76.93% 对冬小麦识别结果，A+C特征比A特征具有更高的识别精度，制图精度提升了4.47%，用户精度提升了3.37%。从混淆矩阵看，精度提升主要是冬小麦与其他类别之间错入错出的减少，这说明纹理特征能显著增强冬小麦与其他地物类型之间的差异。在所有组合中，A+B+C特征集能够取得最高的作物分类制图精度和用户精度，这说明角差特征是冬小麦识别的一个重要判别因素，角差特征与纹理特征的结合更加有利于改进冬小麦识别精度。分析发现冬小麦的制图精度略低，是由于研究区种植的冬小麦大部分刚播种，部分冬小麦光谱特征与裸地的光谱特征较为相近。 对大蒜识别结果，A+B特征比A+C特征取得了更高的识别精度，制图、用户精度分别提升了4.58%、7.36%，是由于大蒜的双向反射差异与其余地物显著不同（ 图 7 图 7 ），角差特征对于大蒜这种覆有地膜情形能够有效地进行表达。在A+B特征集基础上增加纹理特征会导致大蒜的用户精度降低2.25%，冗余的纹理特征加入不利于大蒜的识别。 4.3　McNemar检验分析 2 为进一步分析不同组合分类精度的显著性，针对A、A+B、A+C、A+B+C特征集构建的SVM分类器，以验证样本的SVM分类结果为统计对象，分别两两之间进行McNemar检验。 对不同特征组合的分类器分类结果进行对比得到检验结果（见 表 5 表 5 ），其中f12代表精度评价样本中分类器1错误分类而分类器2正确分类的像元数量，f21则相反。各组合特征的McNemar检验结果X 2 值均大于10.83（p < 0.001），均通过显著性检验，说明不同特征组合的分类结果彼此间具有显著性差异，新特征的加入对作物分类精度的提升是稳定且可信的。 分类器组合vs f12 f21 X 2 p-value A vs A+B 577142 282183 101242.51 <0.001 A vs A+C 828854 445354 115421.88 <0.001 A vs A+B+C 937273 460587 162554.61 <0.001 A+B vs A+C 621576 533035 6789.59 <0.001 A+B vs A+B+C 613105 431378 31617.88 <0.001 A+C vs A+B+C 288678 195492 17934.70 <0.001 4.4　分类结果不确定分析 2 Waldner等（2015） Waldner等（2015） 指出了正确分类的像素倾向于显示相对较低的不确定性值，表达出空间误差的分布情况。在研究区东北和西南地区随机选取1个子区（子区位置见 图 1 图 1 ），分析不同特征组合分类结果及其不确定性，见 图 11 图 11 。 在Ⅰ、Ⅱ子区内，特征集A、A+B分类结果存在把冬小麦地块内的部分区域错分为其他类别的情况，这些区域不确定性值基本大于0.5，而A+C和A+B+C特征组合改善这种错分情况，降低了这些区域的不确定性。 在冬小麦地块区域，基于特征集A+C识别的不确定性略低于特征集A+B，发现特征集A+C能够取得更好的分类结果。特征集A+B+C的分类不确定性值绝大部分小于0.5，其对冬小麦地块的划分更加完整，角差特征与纹理特征结合可以有效提高冬小麦的识别精度，这也与定量化评价相吻合。 对大蒜的分类结果分析，特征集A+B和A+B+C在大蒜区域的不确定性值基本小于0.6，相对于特征集A和A+C具有更低的分类不确定性，大蒜识别结果也更加准确；A+B和A+B+C特征组合在大蒜地块的少部分区域产生错分现象，将大蒜错分为其他地物，是由于这部分大蒜种植较晚导致其光谱信息与裸地相似，这些错分区域不确定性值大于0.6。 在子区Ⅰ、Ⅱ，特征集A分类结果显示林地地块区域分类结果表现出较高的不确定性，同时林地过多的被分为冬小麦，在特征集A+C、A+B+C分类结果中，林地的不确定性显著降低，林地区域仅有极少部分划分为了冬小麦。说明纹理特征所提取出的地物水平结构的差异是区分林地与冬小麦的关键，而角差特征刻画出的地物垂直结构差异可作为纹理特征的补充，为作物准确识别提供可靠的信息。 5　结论与讨论 1 本文通过组合GF-7号卫星影像的光谱特征、纹理特征和角差特征对冬小麦、大蒜识别的影响，证明了GF-7多角度特征对于农作物识别的作用。主要结论如下： （1） GF-7卫星影像的角差特征能显著增加地物之间的可区分度。较仅用光谱特征，光谱与角差特征组合使用大蒜和冬小麦的制图精度分别提高了4.07%和3.15%，用户精度分别提高了6.73%和2.12%；较光谱+纹理特征组合，光谱+纹理+角差特征组合使用大蒜和冬小麦的制图精度分别提高了3.14%和1.01%，用户精度分别提高了5.11%和0.67%，证明了融入角差特征在一定程度上能提高作物的识别精度。 （2）通过McNemar检验不同特征组合的分类结果，证实各特征组合的分类结果差异具有统计学意义，角差特征和纹理特征的加入对农作物识别结果的影响在95%置信度内是可靠的，为利用多角度遥感影像选择不同特征组合进行作物识别提供了依据。 （3）基于分类器熵值的不确定性结果能直观反映出冬小麦和大蒜识别结果的空间误差分布情况，光谱特征+纹理特征+角差特征组合的分类结果整体不确定性值最低，分类结果可靠。 本文研究结果为推进利用GF-7卫星多角度特征提高作物识别能力打下了基础。今后研究将围绕以下方面开展工作：当农作物生长中期阶段，探讨GF-7号多角度影像提取作物高度对作物识别的影响；纹理特征为大蒜的分类带来了混淆，究其原因以及解决办法仍需要深入地探讨；进一步从多时相、不同作物生长农业环境的角度出发，分析多角度影像在整个物候期表现出的空间结构特征变化，这将进一步丰富应用多角度特征进行作物识别的理论基础。 参考文献（ References ） Alshehhi R , Marpu P R , Woon W L and Mura M D . 2017 . S-imultaneous extraction of roads and buildings in remote se-nsing imagery with convolutional neural networks . ISPRS -Journal of Photogrammetry and Remote Sensing , 130 : 139 - 49 [ doi: 10.1016/j.isprsjprs.2017.05.002 http://dx.doi.org/10.1016/j.isprsjprs.2017.05.002 ] S-imultaneous extraction of roads and buildings in remote se-nsing imagery with convolutional neural networks Chang C C . and Lin C J . 2011 . LIBSVM: A Library for support v-ector machines(Article). ACM Transactions on Intelligent Syste-ms and Technology , 2(3): 27: 1 - 27 : 27 [ doi: 10.1145/1961189.1961199 http://dx.doi.org/10.1145/1961189.1961199 ] Lin C J Chopping M , Moisen G G , Su L H , Laliberte A , Albert R , John V M and Debra P . C P . 2008 . Large area mapping of southwestern forest crown cover, canopy height, and biomass using the NAS-A Multiangle Imaging Spectro-Radiometer. Remote Sensing of -Environment , 112 ( 5 ): 2051 - 63 [ doi: 10.1016/j.rse.2007.07.024 http://dx.doi.org/10.1016/j.rse.2007.07.024 ] C P Clayton D G . 1971 . Gram Schmidt Orthogonallization [J]. AppliedStatistics , 20 ( 3 ): 335 - 338 [ doi: 10.2307/2346771 http://dx.doi.org/10.2307/2346771 ] Gram Schmidt Orthogonallization Congalton R G . A review of assessing the accuracy of classifications of remotely sensed data . Remote Sensing of Environment , 1991 , 37 ( 1 ): 270 - 9 . [ doi: 10.1016/0034-4257(91)90048-B http://dx.doi.org/10.1016/0034-4257(91)90048-B ] A review of assessing the accuracy of classifications of remotely sensed data Dehghan H and Ghassemian H . 2006 . Measurement of uncertainty b-y the entropy: application to the classification of MSS data . Int-ernational Journal of Remote Sensing , 27 ( 18 ): 4005 - 14 [ doi: 10.1080/01431160600647225 http://dx.doi.org/10.1080/01431160600647225 ] Measurement of uncertainty b-y the entropy: application to the classification of MSS data Duro D C , Franklin S E , Dubé and Monique G . 2012 . A comparis-on of pixel-based and object-based image analysis with selected machine learning algorithms for the classification of agricultural landscapes using SPOT-5 HRG imagery . Remote Sensing of En-vironment , 118 : 259 - 72 [ doi: 10.1016/j.rse.2011.11.020 http://dx.doi.org/10.1016/j.rse.2011.11.020 ] A comparis-on of pixel-based and object-based image analysis with selected machine learning algorithms for the classification of agricultural landscapes using SPOT-5 HRG imagery Foody G M . 2002 . Status of land cover classification accuracy asse-ssment . Remote Sensing of Environment , 80 ( 1 ): 185 - 201 .[ doi: 10.1016/S0034-4257(01)00295-4 http://dx.doi.org/10.1016/S0034-4257(01)00295-4 ] Status of land cover classification accuracy asse-ssment Foody G M and Atkinson P M . 2002 . Uncertainty in remote sensing and GIS . John Wiley & Sons , [ doi: 10.1080/13658810310001620861 http://dx.doi.org/10.1080/13658810310001620861 ] Uncertainty in remote sensing and GIS Haralick R . M . 1979 . Statistical and structural approaches to texture. Proceedings of the IEEE , 67 ( 5 ): 786 - 804 [ doi: 10.1109/proc.1979.11328 http://dx.doi.org/10.1109/proc.1979.11328 ] M Hsu C W , Chang C C and Lin C J . 2003 . A practical guide to sup-port vector classification . 臺北市 : 國立臺灣大學資訊工程學 . [ doi: 10.1111/j.1365-3016.1995.tb00168.x http://dx.doi.org/10.1111/j.1365-3016.1995.tb00168.x ] A practical guide to sup-port vector classification Huang X , Chen H J and Gong J Y . 2018 . Angular difference featur-e extraction for urban scene classification using ZY-3 multi-ang-le high-resolution satellite imagery . ISPRS Journal of Photogram-metry and Remote Sensing , 135 : 127 - 41 [ doi: 10.1016/j.isprsjprs.2017.11.017 http://dx.doi.org/10.1016/j.isprsjprs.2017.11.017 ] Angular difference featur-e extraction for urban scene classification using ZY-3 multi-ang-le high-resolution satellite imagery Huang X , Lu Q K and Zhang L P . 2014 . A multi-index learning a-pproach for classification of high-resolution remotely sensed im-ages over urban areas . ISPRS Journal of Photogrammetry and -Remote Sensing , 90 : 36 - 48 [ doi: 10.1016/j.isprsjprs.2014.01.008 http://dx.doi.org/10.1016/j.isprsjprs.2014.01.008 ] A multi-index learning a-pproach for classification of high-resolution remotely sensed im-ages over urban areas Khatami R , Mountrakis G and Stehman S V . 2016 . A meta-an-alysis of remote sensing research on supervised pixel-based land-cover image classification processes: General guideli-nes for practitioner and future research . Remote Sensing of Environment , 177 : 89 - 100 [ doi: 10.1016/j.rse.2016.02.028 http://dx.doi.org/10.1016/j.rse.2016.02.028 ] A meta-an-alysis of remote sensing research on supervised pixel-based land-cover image classification processes: General guideli-nes for practitioner and future research Lan Z Y and Liu Y . 2018 . Study on multi-scale window deter-mination for GLCM texture description in high-resolution -remote sensing image geo-analysis supported by GIS and -domain knowledge . ISPRS International Journal of Geo-Inf-ormation , 7 ( 5 ): 175 [ doi: 10.3390/ijgi7050175 http://dx.doi.org/10.3390/ijgi7050175 ] Study on multi-scale window deter-mination for GLCM texture description in high-resolution -remote sensing image geo-analysis supported by GIS and -domain knowledge Li L , She M Y and Luo H . Comparison on fusion algorithms of ZY-3 panchromatic and multi-spectral images . 2014 . Tr-ansactions of the Chinese Society of Agricultural Engineer-ing , 30 ( 16 ): 157 - 165 Comparison on fusion algorithms of ZY-3 panchromatic and multi-spectral images 李霖 ， 佘梦媛 ， 罗恒 . 2014 . ZY-3卫星全色与多光谱影像融合方法比较 . 农业工程学报 , 30 ( 16 ): 157 - 65 [ doi: 10.3969/j.issn.1002-6819.2014.16.021 http://dx.doi.org/10.3969/j.issn.1002-6819.2014.16.021 ] ZY-3卫星全色与多光谱影像融合方法比较 Li X J , Chen W T , Cheng X W and Wang L Z . 2016 . A Co-mparison of Machine Learning Algorithms for Mapping of Complex Surface Mined and Agricultural Landscapes Using ZiYuan-3 Stereo Satellite Imagery . Remote Sensing , 8 ( 6 ): 514 [ doi: 10.3390/rs8060514 http://dx.doi.org/10.3390/rs8060514 ] A Co-mparison of Machine Learning Algorithms for Mapping of Complex Surface Mined and Agricultural Landscapes Using ZiYuan-3 Stereo Satellite Imagery Liu H F , Yang Y B , Yu S , Kong L T and Zhang Y . 2014 . A-daptability evaluation of different fusion methods on ZY-3 and Landsat8 images . Remote Sensing For Land and Reso-urces , 26 ( 4 ): 63 - 70 A-daptability evaluation of different fusion methods on ZY-3 and Landsat8 images 刘会芬 ， 杨英宝 ， 于双 ， 孔令婷 ， 章勇 . 2014 . 遥感图像不同融合方法的适应性评价——以ZY-3和Landsat8图像为例 . 国土资源遥感 , 26 ( 4 ): 63 - 70 [ doi: 10.6046.gtzyyg.2014.04.11 http://dx.doi.org/10.6046.gtzyyg.2014.04.11 ] 遥感图像不同融合方法的适应性评价——以ZY-3和Landsat8图像为例 Liu X S . 2015 . Study on soil bidirectional reflectance characte-ristics and application of soil standard spectral library . Hu-azhong Agricultural University . Study on soil bidirectional reflectance characte-ristics and application of soil standard spectral library 刘小珊 . 2015 . 土壤二向反射特性研究与标准光谱库的应用 . 华中农业大学 , [ doi: 10.7666/d.Y2803877 http://dx.doi.org/10.7666/d.Y2803877 ] 土壤二向反射特性研究与标准光谱库的应用 Matasci G , Longbotham N , Pacifici F , Kanevski M and Tuia D . 2015 . Understanding angular effects in VHR imagery and their significance for urban land-cover model portability: A study of two multi-angle in-track image sequences . ISP-RS Journal of Photogrammetry and Remote Sensing , 107 : 99 - 111 [ doi: 10.1016/j.isprsjprs.2015.05.004 http://dx.doi.org/10.1016/j.isprsjprs.2015.05.004 ] Understanding angular effects in VHR imagery and their significance for urban land-cover model portability: A study of two multi-angle in-track image sequences Jean-Louis R , Marc L and Pierre-Yves D . 1992 . A bidirectional reflectance model of the Earth's surface for the correction of remote sensing data. Journal of Geophysical Research -Atmospheres, 97 (D 18 ): 20455 - 68 . A bidirectional reflectance model of the Earth's surface for the correction of remote sensing data. Journal of Geophysical Research -Atmospheres, 97 Schaepman-Strub G , Schaepman M E , Painter T H , Dangel S and Martonchik J V . 2006 . Reflectance quantities in optical remote -sensing—definitions and case studies . Remote Sensing of Envir-onment , 103 ( 1 ): 27 - 42 [ doi: 10.1016/j.rse.2006.03.002 http://dx.doi.org/10.1016/j.rse.2006.03.002 ] Reflectance quantities in optical remote -sensing—definitions and case studies Song D J , Zhang C M , Yang X X , Li F , Han Y J , Gao S and Do-ng H Y . 2020 . Extracting winter wheat spatial distribution infor-mation from GF-2 image . Journal of Remote Sensing , 24 ( 5 ): 596 - 608 Extracting winter wheat spatial distribution infor-mation from GF-2 image 宋德娟 ， 张承明 ， 杨晓霞 ， 李峰 ， 韩颖娟 ， 高帅 ， 董海燕 . 2020 . 高分二号遥感影像提取冬小麦空间分布 . 遥感学报 , 24 ( 5 ): 596 - 608 [ doi: 10.11834/jrs.2020828 http://dx.doi.org/10.11834/jrs.2020828 ] 高分二号遥感影像提取冬小麦空间分布 Su L H , Mark J. C , Albert R , John V . M and Debra P . C . P. 2007 . Support vector machines for recognition of semi-arid vegetation types using MISR multi-angle imagery. Remote Sensing of Env-ironment, 107 ( 1-2 ): 299 - 311 [ doi: 10.1016/j.rse.2006.05.023 http://dx.doi.org/10.1016/j.rse.2006.05.023 ] M Tian J , Reinartz P . d’Angelo P and Ehlers M . 2013 . Region based -automatic building and forest change detection on Cartosat-1 st-ereo imagery. ISPRS Journal of Photogrammetry and Remote S-ensing , 79 : 226 - 39 [ doi: 10.1016/j.isprsjprs.2013.02.017 http://dx.doi.org/10.1016/j.isprsjprs.2013.02.017 ] d’Angelo P and Ehlers M Waldner F , Canto G S and Defourny P . 2015 . Automated annual cr-opland mapping using knowledge-based temporal features . ISPRS Journal of Photogrammetry and Remote Sensing , 110 : 1 - 13 [ doi: 10.1016/j.isprsjprs.2015.09.013 http://dx.doi.org/10.1016/j.isprsjprs.2015.09.013 ] Automated annual cr-opland mapping using knowledge-based temporal features Wu J Y , Liu X L , Shi Z T and Fu Z . 2019 . Plastic greenhouse rec-ognition based on GF-2 data and multi-texture features . Transac-tions of the Chinese Society of Agricultural Engineering , 35 ( 4 ): 106 - 13 Plastic greenhouse rec-ognition based on GF-2 data and multi-texture features 吴锦玉 ， 刘晓龙 ， 柏延臣 ， 史正涛 ， 付卓 . 2019 . 基于GF-2数据结合多纹理特征的塑料大棚识别 . 农业工程学报 , 35 ( 12 ): 173 - 183 [ doi: 10.11975/j.issn.1002-6819.2019.12.021 http://dx.doi.org/10.11975/j.issn.1002-6819.2019.12.021 ] 基于GF-2数据结合多纹理特征的塑料大棚识别 Yang Y J , Tian Q J , Zhan Y L , Tao B and Xu K J . 2018 . Effects of Spatial Resolution and Texture Features on Multi-spectral R-emote Sensing Classification , Journal of Geo-information Scienc-e , 20 ( 1 ): 99 - 107 Effects of Spatial Resolution and Texture Features on Multi-spectral R-emote Sensing Classification 杨闫君 ， 田庆久 ， 占玉林 ， 陶波 ， 徐凯健 . 2018 . 空间分辨率与纹理特征对多光谱遥感分类的影响 . 地球信息科学学报 , 20 ( 1 ): 99 - 107 [ doi: 10.12082/dqxxkx.2018.170177 http://dx.doi.org/10.12082/dqxxkx.2018.170177 ] 空间分辨率与纹理特征对多光谱遥感分类的影响 Yu W W , Xu K J , Zhao P , Shen P J and Zhao Y J . 2021 . Influen-ce of red-edge spectrum of sentinel-2 image on identification of dominant tree species in different growing periods . Geography and Geo-Information Science , 37 ( 3 ): 9 Influen-ce of red-edge spectrum of sentinel-2 image on identification of dominant tree species in different growing periods 于婉婉 ， 徐凯健 ， 赵萍 ， 申鹏举 ， 赵月娇 . 2021 . Sentinel-2影像红边谱段对不同生长期区域优势树种识别的影响 . 地理与地理信息科学 , 37 ( 3 ): 9 [ doi: 10.3969/j.issn.1672-0504.2021.03.007 http://dx.doi.org/10.3969/j.issn.1672-0504.2021.03.007 ] Sentinel-2影像红边谱段对不同生长期区域优势树种识别的影响 Zhang D J , Pan Y Z , Zhang J S , Hu T G , Zhao J H , Li N and Chen Q . 2020 . A generalized approach based on convolutional n-eural networks for large area cropland mapping at very high re-solution . Remote Sensing of Environment , 247 : 111912 [ doi: 10.1016/j.rse.2020.111912 http://dx.doi.org/10.1016/j.rse.2020.111912 ] A generalized approach based on convolutional n-eural networks for large area cropland mapping at very high re-solution Zhu X F , Li S B and Xiao G F . 2019 . Method on extraction of area and distribution of plastic-mulched farmland based on UAV images . Transactions of the Chinese Society of Agricultural Engineering , 35 ( 4 ): 106 - 113 Method on extraction of area and distribution of plastic-mulched farmland based on UAV images 朱秀芳 ， 李石波 ， 肖国峰 . 2019 . 基于无人机遥感影像的覆膜农田面积及分布提取方法 . 农业工程学报 , 35 ( 4 ): 106 - 13 [ doi: 10.11975/j.issn.1002-6819.2019.04.013 http://dx.doi.org/10.11975/j.issn.1002-6819.2019.04.013 ] 基于无人机遥感影像的覆膜农田面积及分布提取方法 10.11834/jrs.20221644 TP701 A 孙智虎，张锦水，洪友堂，杨珺雯，朱爽.XXXX.基于GF-7卫星多角度特征进行作物识别.遥感学报，XX（XX）： 1-12 孙智虎，张锦水，洪友堂，杨珺雯，朱爽.XXXX.基于GF-7卫星多角度特征进行作物识别.遥感学报，XX（XX）： 1-12 SUN Zhihu，ZHANG Jinshui，HONG Youtang，YANG Junwen，ZHU Shuang. XXXX. Crop recognition by multi-angle features of GF-7 satellite. National Remote Sensing Bulletin， XX（XX）：1-12 SUN Zhihu，ZHANG Jinshui，HONG Youtang，YANG Junwen，ZHU Shuang. XXXX. Crop recognition by multi-angle features of GF-7 satellite. National Remote Sensing Bulletin， XX（XX）：1-12 国家重点研发计划“粮食丰产增效科技创新专项”子课题（编号：2017YFD0300402-6）；国家高分辨率对地观测系统重大专项（民用部分）（编号：20-Y30F10-9001-20/22） Supported by Supported by National key research and development program (2017YFD0300402); Major Project of High Resolution Earth Observation System (Civil Part) (20-Y30F10-9001-20/22) 2021-10-08 2022-05-25 遥感学报