基于动态时间弯曲的时序sar影像耕地提取方法
【专利摘要】本发明涉及一种基于动态时间弯曲的时序SAR影像耕地提取方法，步骤包括：构建时序SAR影像；提取耕地参考时间序列；计算待分类像元时间序列与耕地参考时间序列之间的动态时间弯曲距离；计算结果阈值分割，待分类像元归类为耕地与非耕地；分割结果空域滤波，滤除孤立的耕地像元，填补连片耕地之间的缝隙，得到耕地的最终提取结果。本发明考虑到耕地时间序列特有的“时间轴弯曲”现象，使用动态时间弯曲距离（DTW）作为相似性度量标准，从而实现耕地像元与非耕地像元的划分，解决了传统方法无法适应时间轴畸变的时间序列相似性度量这一问题，提高了耕地的提取精度。本发明方法适应性强，提取精度可达82%以上，能够满足实际生产的需要。
【专利说明】基于动态时间弯曲的时序SAR影像耕地提取方法
【技术领域】
[0001]本发明涉及一种遥感影像信息提取方法，特别涉及一种基于动态时间弯曲的时序SAR影像耕地提取方法。
【背景技术】
[0002]耕地是粮食生产的基本条件，是确保粮食安全的物质基础，掌握实时的耕地信息对国家农业政策的制定具有重要意义。遥感技术经过几十年的发展，已经形成了 “星——空一地” 一体的对地观测网络，并产生了海量堆积数据，因此使用遥感技术进行大范围耕地信息的快速提取领域具有独到的优势。
[0003]耕地表面附着的作物拥有随季节变化十分明显的物候特征，不同生长周期(如播种期、营养生长期、生殖生长期和成熟期等)在遥感影像上有着不同的表现，因此使用单一时间节点的遥感影像往往不能有效提取出耕地信息。时间序列遥感影像(以下简称时序遥感影像)是卫星传感器对同一地区不同时间进行观测所获得的数值序列，通过构建时序遥感影像，可以识别观测数值的时间变化特征，并利用该特征进行耕地信息的有效提取。然而受观测条件多变、农作物物候差异较大等因素的影响，现阶段利用时序遥感影像反演耕地的工作仍存在不少问题，大致可归纳为两个方面:
[0004]I)所使用遥感数据的不足。当前用于耕地反演的时序遥感数据大多数为SPOT/VEGETATION、NOAA/AVHRR、MODIS等影像数据构建的NDV1、EVI等植被指数产品(Sakamoto, 2005 ；ffardlow, 2008 ；Dheeravath, 2010)，它们拥有较高的时间分辨率，然而空间分辨率很低(最高仅250m)，复合像元问题十分严重；同时，卫星传感器的光学观测过程中常受到云层遮挡，得到的时间序列存在大量噪声；此外，NDVI本身会受到土壤背景和含水量等因素较大的影响，且在高植被覆盖区域容易达到饱和(Baush，1993 ;唐怡，2006)。合成孔径雷达(SAR)为主动式传感器，发射的微波可以穿透云层，能够全天时不受天气影响完成对地观测任务。随着技术的发展，SAR传感器的空间分辨率以及时间分辨率都有了极大的提升，欧空局的ENVISAT ASAR数据空间分辨率最高可达10m，中国HJ-1C星SAR数据空间分辨率最高可达4m，德国TerraSAR-X、加拿大RADARSAT-2的SAR数据空间分辨率最高可达Im ;同时，SAR可侧视成像的特性保证了影像的高时间分辨率，许多卫星传感器能够提供时间间隔小于10天的SAR影像数据(TerraSAR-X对全球95%的地区可实现2天重访)。然而尽管SAR影像应用范围越来越广，却少有将其构建时间序列进行耕地信息提取的研究。
[0005]2)所采用方法的不足。利用时序遥感影像提取耕地的关键技术是时间序列相似性度量，即定量评估待分类像元时间序列与耕地参考时间序列的相似性程度。对大区域耕地而言，其种植的农作物种类多样，不同种类农作物的物候期开始时间不一，持续长短不同，致使耕地时间序列在时间轴上存在一定的平移和弯曲(Calendar Shift)。而当前进行遥感时间序列相似性度量的方法大体上可分为三类:距离法(欧氏距离、马氏距离等)、代数法(作差法、作商法等)和相关法(Pearson相关系数、Spearman相关系数等)，它们均无法适用于“时间轴移位”的时间序列。动态时间弯曲(Dynamic time warping, DTW)是一个在计算机科学中常用的时间序列相似度算法，它可以在一定程度上自动寻找发生平移和弯曲的时间节点，并对序列的时间轴进行动态弯曲以实现自适应匹配，从而计算动态弯曲距离，作为两条时间序列的相似度，从机理上看十分适于耕地时间序列的相似性分析，然而相应的研究或应用却并不多见。
【发明内容】
[0006]本发明要解决的技术问题是:克服现有技术缺点，提出一种基于动态时间弯曲的时序SAR影像耕地提取方法。
[0007]本发明基于动态时间弯曲的时序SAR影像耕地提取方法，包括以下步骤:
[0008]第一步、时序SAR影像构建一对研究区研究时段内的原始SAR影像数据进行预处理和配准，并按照时间先后顺序构建时序SAR影像；
[0009]第二步、耕地参考时间序列提取一选择研究区内的若干耕地像元作为训练样本，针对所述时序SAR影像中的每一幅影像，分别计算训练样本中耕地像元的后向散射系数的均值和标准差，若所述标准差之和小于设定的可靠性阈值，则利用所述均值依据对应影像的成像时间构建耕地参考时间序列；否则，重新选择训练样本并重复本步骤，直至标准差小于可靠性阈值；
[0010]第三步、时间序列相似度计算一针对每个待分类像元，构建其后向散射系数时间序列，并计算所述后向散射系数时间序列与耕地参考时间序列之间的动态时间弯曲距离，作为两时间序列之间的相似度；
[0011]第四步、计算结果阈值分割——针对相似度计算结果，使用预先设定的分割阈值进行二值分割，将动态时间弯曲距离大于分割阈值的对应待分类像元归类为非耕地像元，其余归为耕地像元，获得阈值分割结果；
[0012]第五步、分割结果空域滤波——使用空间操作算子对所述阈值分割结果进行空域滤波，滤除孤立的耕地像元，填补连片耕地之间的缝隙，得到耕地的最终提取结果。
[0013]本发明基于动态时间弯曲的时序SAR影像耕地提取方法，还具有如下改进:
[0014]1、所述第一步中，预处理至少包括轨道参数更新、辐射校正、地形校正、几何校正和去噪滤波。
[0015]2、所述第一步中，影像的配准误差不大于0.5像素。
[0016]3、所述第二步中，构成训练样本的耕地像元通过人工识别提取，训练样本中的耕地像元数量不低于9个。
[0017]4、第二步中的可靠性阈值Thl通过公式Thl =NXcXR计算获得，其中，N为时序SAR影像中的影像数量，c为比率系数，取值范围为[0.3,0.5]，R为所有影像的训练样本中耕地像元的后向散射系数均值的最大值与最小值之差。
[0018]5、所述第三步中，利用待分类像元在每一幅影像中的后向散射系数按照对应影像的成像时间构建待分类像元的后向散射系数时间序列。
[0019]6、所述第四步中，分割阈值Th2通过公式Th2 = C1XX计算获得，其中，C1为调节系数，其取值范围为0.3?I, Σ为第二步中所有标准差的求和。
[0020]7、所述第五步中，空间操作算子为众数滤波算子、数学形态学开闭算子、中值滤波算子中的一种，空间操作算子的尺寸size根据下式获得:
【权利要求】
1.一种基于动态时间弯曲的时序SAR影像耕地提取方法，包括以下步骤: 第一步、时序SAR影像构建——对研究区研究时段内的原始SAR影像数据进行预处理和配准，并按照时间先后顺序构建时序SAR影像； 第二步、耕地参考时间序列提取一选择研究区内的若干耕地像元作为训练样本，针对所述时序SAR影像中的每一幅影像，分别计算训练样本中耕地像元的后向散射系数的均值和标准差，若所述标准差之和小于设定的可靠性阈值，则利用所述均值依据对应影像的成像时间构建耕地参考时间序列；否则，重新选择训练样本并重复本步骤，直至标准差小于可靠性阈值； 第三步、时间序列相似度计算一针对每个待分类像元，构建其后向散射系数时间序列，并计算所述后向散射系数时间序列与耕地参考时间序列之间的动态时间弯曲距离，作为两时间序列之间的相似度； 第四步、计算结果阈值分割——针对相似度计算结果，使用预设的分割阈值进行二值分割，将动态时间弯曲距离大于分割阈值的对应待分类像元归类为非耕地像元，其余归为耕地像元，获得阈值分割结果； 第五步、分割结果空域滤波——使用空间操作算子对所述阈值分割结果进行空域滤波，滤除孤立的耕地像元，填补连片耕地之间的缝隙，得到耕地的最终提取结果。
2.根据权利要求1所述的基于动态时间弯曲的时序SAR影像耕地提取方法，其特征在于:所述第一步中，预处理至少包括轨道参数更新、辐射校正、地形校正、几何校正和去噪滤波。
3.根据权利要求1所述的基于动态时间弯曲的时序SAR影像耕地提取方法，其特征在于:所述第一步中，影像的配准误差不大于0.5像素。
4.根据权利要求1所述的基于动态时间弯曲的时序SAR影像耕地提取方法，其特征在于:所述第二步中，构成训练样本的耕地像元通过人工识别提取，训练样本中的耕地像元数量不低于9个。
5.根据权利要求1所述的基于动态时间弯曲的时序SAR影像耕地提取方法，其特征在于:第二步中的可靠性阈值Thl通过公式Thl = NXcXR计算获得，其中，N为时序SAR影像中的影像数量，c为比率系数，取值范围为[0.3,0.5]，R为所有影像的训练样本中耕地像元的后向散射系数均值的最大值与最小值之差。
6.根据权利要求1所述的基于动态时间弯曲的时序SAR影像耕地提取方法，其特征在于:所述第三步中，利用待分类像元在每一幅影像中的后向散射系数按照对应影像的成像时间构建待分类像元的后向散射系数时间序列。
7.根据权利要求1所述的基于动态时间弯曲的时序SAR影像耕地提取方法，其特征在于:所述第四步中，分割阈值Th2通过公式Th2 = C1XX计算获得，其中，C1为调节系数，其取值范围为0.3~I, Σ为第二步中所有标准差的求和。
8.根据权利要求1所述的基于动态时间弯曲的时序SAR影像耕地提取方法，其特征在于:所述第五步中，空间操作算子为众数滤波算子、数学形态学开闭算子、中值滤波算子中的一种，空间操作算子的尺寸size根据下式获得:
【文档编号】G06Q50/02GK104008552SQ201410266813
【公开日】2014年8月27日 申请日期:2014年6月16日 优先权日:2014年6月16日
【发明者】钟礼山, 朱明媛, 李满春, 程亮, 王结臣, 王亚飞, 姜鹏辉, 孙越凡, 马磊, 潘航 申请人:南京大学, 南京市国土资源信息中心