一种利用雷达遥感数据监测小麦倒伏的方法及装置制造方法
【专利摘要】本发明涉及一种利用雷达遥感数据监测小麦倒伏的方法及装置，包括：获取监测区域的全极化合成孔径雷达遥感影像；处理所述合成孔径雷达遥感影像，提取不同极化通道的后向散射系数；提取监测区域内小麦种植区域的地块单元；计算每个地块单元所有像素在不同极化通道上的平均后向散射系数，并利用所述平均后向散射系数构建极化指数；将所述极化指数作为判断指标，对小麦的倒伏状况进行判断。采用该方法及装置，实现了大面积、快速准确监测小麦倒伏状况，掌握小麦倒伏的发生范围和严重程度，是灾情诊断、及时防控、损失评估的关键；雷达遥感数据不受天气影响，且雷达遥感观测对结构变化十分敏感，利用全极化合成孔径雷达数据监测倒伏更具优势。
【专利说明】一种利用雷达遥感数据监测小麦倒伏的方法及装置
【技术领域】
[0001]本发明涉及农业灾害及雷达遥感应用【技术领域】，尤其涉及一种利用雷达遥感数据监测小麦倒伏的方法及装置。
【背景技术】
[0002]小麦倒伏是农业生产的常见灾害之一。小麦一旦发生倒伏，植株水分、养分的运转以及光合作用都会降低，还会诱发各种病虫害，严重影响籽粒灌浆过程，最终影响到小麦产量与籽粒品质的形成，严重倒伏时产量损失可达27%。此外由于倒伏不利于机械收获，人力收割成本的增加也会加重农田收益的损失。大面积、快速监测小麦倒伏状况是掌握灾情、及时防控、评估损失的关键，对于农业部门及时获取农田小麦生长信息具有重要价值。
[0003]遥感技术的发展为倒伏的宏观快速监测提供了一种有效手段。目前，不少研究者已经尝试用遥感技术来监测倒伏，刘良云等(遥感学报，2005)发现冬小麦冠层光谱反射率随倒伏角度的增加而增加，并利用Landsat ETM影像倒伏前后归一化植被指数(NDVI)的变化，成功监测了小麦倒伏的发生程度；刘占宇等(红外与毫米波学报，2009)利用可见光/近红外地面光谱仪数据，通过支持向量机方法对倒伏水稻进行了识别；张杰、胡宗杰等(安徽农业科学，2011)研究了灌浆期小麦倒伏冠层光谱变化的原因和特征；张竞成等(SensorLetters, 2012)利用地面高光谱数据和连续小波技术评估了倒伏对玉米品质的影响；吴尚蓉等(中国农业资源与区划，2013)基于图像和光谱技术评估了倒伏对冬小麦产量的影响
坐寸ο
[0004]可以发现，目前这些方法主要利用了光学遥感的光谱特性:倒伏后光谱反射率发生变化，进而基于变化的光谱反射率判断倒伏发生的程度。但是光学遥感的光谱技术有其局限性，倒伏引起的光谱变化比较微弱，且往往淹没在复杂多变的混合光谱之中，实际上导致光谱变化的因素较多，如农田环境和其他胁迫如病虫害、水肥胁迫等，因此难以将倒伏的弱信息从众多影响因素中提取出来；而且发生倒伏时，往往伴随着灾害天气，光学遥感数据受阴雨天气的影响无法保证数据的及时获取。
【发明内容】
[0005](一)要解决的技术问题
[0006]本发明所要解决的技术问题是:基于光学遥感的光谱技术监测小麦倒伏状况存在一定的局限性，目前大面积快速监测小麦倒伏状况手段匮乏的问题。
[0007](二)技术方案
[0008]为此目的，本发明提出了一种利用雷达遥感数据监测小麦倒伏的方法，包括以下步骤:
[0009]S1:获取监测区域的全极化合成孔径雷达遥感影像；
[0010]S2:处理所述合成孔径雷达遥感影像，提取不同极化通道的后向散射系数；
[0011]S3:提取监测区域内小麦种植区域的地块单元；[0012]S4:计算每个地块单元所有像素在不同极化通道上的平均后向散射系数，并利用所述平均后向散射系数构建极化指数；
[0013]S5:将所述极化指数作为判断指标，对小麦的倒伏状况进行判断。
[0014]优选地，所述方法还包括:
[0015]实时采集监测区域的全极化合成孔径雷达遥感影像。
[0016]优选地,所述步骤S2具体包括:
[0017]对所述合成孔径雷达遥感影像的四个极化通道进行辐射定标；
[0018]将辐射定标后每个极化通道的复数据生成强度数据；
[0019]对每个极化通道的所述强度数据进行地理编码；
[0020]对地理编码后的图像数据进行几何精校正和斑点噪声的去除处理。
[0021]优选地，所述将辐射定标后每个极化通道的复数据生成强度数据具体为:
[0022]获取每个极化通道上每一个像素的复数据；
[0023]计算每个极化通道上每一个像素的强度A，公式为:
[0024]A=a2+b2
[0025]其中,A为每个极化通道上任一个像素的强度,a为该像素对应的实部数据，b为该像素对应的虚部数据。
[0026]优选地,所述步骤S3具体包括:
[0027]利用已有的土地利用数据库和/或根据面向对象的分类方法得到农作物种植区划地图；
[0028]对所述农作物种植区划地图进行矢量化，得到小麦种植区域的地块单元。
[0029]优选地,所述步骤S4具体包括:
[0030]以地块单元为分析单元计算每个地块单元内所有像素在不同极化通道上的平均后向散射系数；所述平均后向散射系数包括:地块单元内所有像素的HH极化平均后向散射
系数<，地块单元内所有像素的VV极化平均后向散射系数c<v，地块单元内所有像素的HV极化平均后向散射系数.7
[0031]利用所述平均后向散射系数构建极化指数，公式为:
[0032]piI = σΗΗ /σνν ;
[0033]ΡΙ2 = CThh /crHV ；
[0034]其中，PI1为第一极化指数，PI2为第二极化指数。
[0035]优选地,所述步骤S5具体包括:
[0036]对比所述极化指数与正常小麦地块的极化指数；
[0037]获取对比的差异程度，并根据所述差异程度利用预设阈值对小麦的倒伏状况进行判断。
[0038]此外，本发明还提出了一种利用雷达遥感数据监测小麦倒伏的装置，包括:
[0039]获取模块，用于获取监测区域的全极化合成孔径雷达遥感影像；
[0040]处理模块，用于处理所述合成孔径雷达遥感影像，提取不同极化通道的后向散射系数；[0041]提取模块，用于提取监测区域内小麦种植区域的地块单元；
[0042]计算模块，用于计算每个地块单元所有像素在不同极化通道上的平均后向散射系数，并利用所述平均后向散射系数构建极化指数；
[0043]判断模块，用于将所述极化指数作为判断指标，对小麦的倒伏状况进行判断。
[0044]优选地，所述装置还包括:
[0045]采集模块，用于实时采集监测区域的全极化合成孔径雷达遥感影像。
[0046]优选地，所述处理模块包括:
[0047]辐射定标单元，用于对所述合成孔径雷达遥感影像的四个极化通道进行辐射定标；
[0048]强度单元，用于将辐射定标后每个极化通道的复数据生成强度数据；
[0049]地理编码单元，用于对每个极化通道的所述强度数据进行地理编码；
[0050]优化单元，用于对地理编码后的图像数据进行几何精校正和斑点噪声的去除处理。
[0051](三)有益效果
[0052]通过采用本发明公开的一种利用雷达遥感数据监测小麦倒伏的方法及装置，实现了大面积、快速准确监测小麦倒伏状况，掌握小麦倒伏的发生范围和严重程度，是灾情诊断、及时防控、损失评估的关键；此外，与传统的光学遥感方法相比，雷达遥感数据不受天气影响，且雷达遥感观测对结构变化十分敏感，利用全极化合成孔径雷达数据监测倒伏将更具优势。
【专利附图】
【附图说明】
[0053]通过参考附图会更加清楚的理解本发明的特征和优点，附图是示意性的而不应理解为对本发明进行任何限制，在附图中:
[0054]图1是本发明一种利用雷达遥感数据监测小麦倒伏的方法流程图；
[0055]图2是本发明实施例中第一极化指数/ Wv随生育进程的变化示意图；
[0056]图3是本发明实施例中第二极化指数Wh / Wv随生育进程的变化示意图；
[0057]图4是利用第一极化指数获得的小麦倒伏监测结果示意图；
[0058]图5是利用第二极化指数获得的小麦倒伏监测结果示意图；
[0059]图6是利用HH单极化后向散射系数获得的小麦倒伏监测结果示意图；
[0060]图7是利用VV单极化后向散射系数获得的小麦倒伏监测结果示意图；
[0061]图8是利用HV单极化后向散射系数获得的小麦倒伏监测结果示意图；
[0062]图9是本发明一种利用雷达遥感数据监测小麦倒伏的装置模块图。
【具体实施方式】
[0063]下面将结合附图对本发明的实施例进行详细描述。
[0064]本发明实施例提出了一种利用雷达遥感数据监测小麦倒伏的方法，如图1所示，包括以下步骤:
[0065]S1:获取监测区域的全极化合成孔径雷达遥感影像；[0066]S2:处理所述合成孔径雷达遥感影像，提取不同极化通道的后向散射系数；具体包括:对所述合成孔径雷达遥感影像的四个极化通道进行辐射定标；将辐射定标后每个极化通道的复数据生成强度数据；对每个极化通道的所述强度数据进行地理编码；对地理编码后的图像数据进行几何精校正和斑点噪声的去除处理。
[0067]其中，所述将辐射定标后每个极化通道的复数据生成强度数据具体为:
[0068]获取每个极化通道上每一个像素的复数据；
[0069]计算每个极化通道上每一个像素的强度A，公式为:
[0070]A=a2+b2
[0071]其中,A为每个极化通道上任一个像素的强度,a为该像素对应的实部数据，b为该像素对应的虚部数据。
[0072]S3:提取监测区域内小麦种植区域的地块单元；具体包括:利用已有的土地利用数据库和/或根据面向对象的分类方法得到农作物种植区划地图；对所述农作物种植区划地图进行矢量化，得到小麦种植区域的地块单元。
[0073]S4:计算每 个地块单元所有像素在不同极化通道上的平均后向散射系数，并利用所述平均后向散射系数构建极化指数；具体包括:以地块单元为分析单元计算每个地块单元内所有像素在不同极化通道上的平均后向散射系数；所述平均后向散射系数包括:地块
单元内所有像素的HH极化平均后向散射系数，地块单元内所有像素的VV极化平均后
向散射系数，地块单元内所有像素的HV极化平均后向散射系数;利用所述平均后向散射系数构建极化指数，公式为:
[0074]P】i = crHH j σνν.[0075]ΡΙ2 = σ*π /σ?°ιν.[0076]其中，PI1为第一极化指数，PI2为第二极化指数。
[0077]S5:将所述极化指数作为判断指标，对小麦的倒伏状况进行判断；具体包括:对比所述极化指数与正常小麦地块的极化指数；获取对比的差异程度，并根据所述差异程度利用预设阈值对小麦的倒伏状况进行判断。
[0078]其中，本发明实施例公开的利用雷达遥感数据监测小麦倒伏的方法还包括:
[0079]实时采集监测区域的全极化合成孔径雷达遥感影像。
[0080]本发明实施例以内蒙古自治区额尔古纳市上库力农场为研究区，其地处中国最东北部，大兴安岭西北麓，呼伦贝尔草原北端，50.28°~50.39° N, 120.76°~120.89° E，属于典型的东北高寒地区黑土地春小麦种植带，大陆性气候显著。该农场为农垦系统国有农场，农场面积约2，800hm2，地形起伏较小，种植结构相对简单，机械化程度高，地块整齐均匀，主要种植小麦、油菜和大麦3种作物，均为大面积种植，平均每地块17hm2。由于冬季时间长，温度低，农作物一年一季，作物生育期从5月中旬至9月初。2013年7月底，由于大雨、大风等天气因素，该农田小麦发生较大面积倒伏。实地调查了研究区31个地块的倒伏状况，发现13地块在8月3日时处于不同程度的倒伏状态，而另18个地块一直正常生长。
[0081]本发明实施例中具体的实施步骤包括:
[0082]获取研究区域2013年8月3日的全极化Radarsat_2遥感影像,此时研究区小麦正处于抽穗灌浆期。Radarsat-2为加拿大新一代合成孔径雷达卫星，性能较为稳定，辐射精度较高，幅宽25kmX 25km, C波段,使用其SLC (single look complex)格式产品。
[0083]需要说明的是，对本发明实施例还获取了其它生育期的4景Radarsat-2影像，如表1所示，以辅助证明本发明方法的有效性，在其它时间没有倒伏，极化指数并无太大变化。该5景Radarsat-2影像成像模式、成像入射角等参数完全一致，如表2所示:
[0084]表1 5景雷达影像获取日期和对应生育期
【权利要求】
1.一种利用雷达遥感数据监测小麦倒伏的方法，其特征在于，该方法包括以下步骤: S1:获取监测区域的全极化合成孔径雷达遥感影像； S2:处理所述合成孔径雷达遥感影像，提取不同极化通道的后向散射系数； S3:提取监测区域内小麦种植区域的地块单元； S4:计算每个地块单元所有像素在不同极化通道上的平均后向散射系数，并利用所述平均后向散射系数构建极化指数； S5:将所述极化指数作为判断指标，对小麦的倒伏状况进行判断。
2.根据权利要求1所述的方法，其特征在于，所述方法还包括: 实时采集监测区域的全极化合成孔径雷达遥感影像。
3.根据权利要求1所述的方法，其特征在于，所述步骤S2具体包括: 对所述合成孔径雷达遥感影像的四个极化通道进行辐射定标； 将辐射定标后每个极化通道的复数据生成强度数据； 对每个极化通道的所述强度数据进行地理编码； 对地理编码后的图像数据进行几何精校正和斑点噪声的去除处理。
4.根据权利要求3所述的方法，其特征在于，所述将辐射定标后每个极化通道的复数据生成强度数据具体为: 获取每个极化通道上每一个像素的复数据； 计算每个极化通道上每一个像素的强度A，公式为:
A=a2+b2 其中，A为每个极化通道上任一个像素的强度，a为该像素对应的实部数据，b为该像素对应的虚部数据。
5.根据权利要求1所述的方法，其特征在于，所述步骤S3具体包括: 利用已有的土地利用数据库和/或根据面向对象的分类方法得到农作物种植区划地图； 对所述农作物种植区划地图进行矢量化，得到小麦种植区域的地块单元。
6.根据权利要求1所述的方法，其特征在于，所述步骤S4具体包括:以地块单元为分析单元计算每个地块单元内所有像素在不同极化通道上的平均后向散射系数；所述平均后向散射系数包括:地块单元内所有像素的HH极化平均后向散射系数 ,地块单元内所有像素的VV极化平均后向散射系数,地块单元内所有像素的HV极化平均后向散射系数； 利用所述平均后向散射系数构建极化指数，公式为:
OO
PJ = ijHH /crVV *
OO
PJ2 — crHH /aHV * 其中，PI1为第一极化指数，PI2为第二极化指数。
7.根据权利要求1所述的方法，其特征在于，所述步骤S5具体包括: 对比所述极化指数与正常小麦地块的极化指数；获取对比的差异程度，并根据所述差异程度利用预设阈值对小麦的倒伏状况进行判断。
8.一种利用雷达遥感数据监测小麦倒伏的装置，其特征在于，包括: 获取模块，用于获取监测区域的全极化合成孔径雷达遥感影像； 处理模块，用于处理所述合成孔径雷达遥感影像，提取不同极化通道的后向散射系数； 提取模块，用于提取监测区域内小麦种植区域的地块单元； 计算模块，用于计算每个地块单元所有像素在不同极化通道上的平均后向散射系数，并利用所述平均后向散射系数构建极化指数； 判断模块，用于将所述极化指数作为判断指标，对小麦的倒伏状况进行判断。
9.根据权利要求8所述的装置，其特征在于，所述装置还包括: 采集模块，用于实时采集监测区域的全极化合成孔径雷达遥感影像。
10.根据权利要求8 所述的装置，其特征在于，所述处理模块包括: 辐射定标单元，用于对所述合成孔径雷达遥感影像的四个极化通道进行辐射定标； 强度单元，用于将辐射定标后每个极化通道的复数据生成强度数据； 地理编码单元，用于对每个极化通道的所述强度数据进行地理编码； 优化单元，用于对地理编码后的图像数据进行几何精校正和斑点噪声的去除处理。
【文档编号】G01S7/41GK103969632SQ201410116033
【公开日】2014年8月6日 申请日期:2014年3月26日 优先权日:2014年3月26日
【发明者】赵春江, 杨浩, 杨贵军, 杨小冬, 徐新刚, 顾晓鹤, 宋晓宇, 张竞成, 董彦生, 冯海宽 申请人:北京农业信息技术研究中心