CN104504279B - 检测蝗灾的方法 - Google Patents
检测蝗灾的方法 Download PDFInfo
- Publication number
- CN104504279B CN104504279B CN201410852062.5A CN201410852062A CN104504279B CN 104504279 B CN104504279 B CN 104504279B CN 201410852062 A CN201410852062 A CN 201410852062A CN 104504279 B CN104504279 B CN 104504279B
- Authority
- CN
- China
- Prior art keywords
- mrow
- normalized
- index
- locust
- msub
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 206010035148 Plague Diseases 0.000 title claims abstract description 64
- 241000607479 Yersinia pestis Species 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 24
- 241000238814 Orthoptera Species 0.000 title claims abstract description 11
- 230000010287 polarization Effects 0.000 claims abstract description 81
- 238000003384 imaging method Methods 0.000 claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims description 40
- 230000001427 coherent effect Effects 0.000 claims description 8
- 238000001514 detection method Methods 0.000 description 9
- 230000003287 optical effect Effects 0.000 description 4
- 238000010606 normalization Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 241000238631 Hexapoda Species 0.000 description 2
- 241000254023 Locusta Species 0.000 description 2
- 241000254032 Acrididae Species 0.000 description 1
- 241000254022 Locusta migratoria Species 0.000 description 1
- 241000346164 Locusta migratoria manilensis Species 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 241000894007 species Species 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供一种检测蝗灾的方法。所述方法包括：(A)获取被检测区域的多源遥感数据，其中，所述多源遥感数据包括：全极化合成孔径雷达数据、中分辨率成像光谱仪温度和专题制图仪植被指数；(B)通过全极化合成孔径雷达数据获得平均极化散射角和极化散射熵；(C)通过平均极化散射角、极化散射熵、中分辨率成像光谱仪温度和专题制图仪植被指数获得差值蝗灾指数；(D)根据差值蝗灾指数来确定被检测区域是否发生蝗灾。
Description
技术领域
本发明涉及农业领域，更具体地讲，涉及一种检测蝗灾的方法。
背景技术
东亚飞蝗是诱发中国北方蝗灾最常见的昆虫，以禾本植物为食，是我国北方的优势种，也是严重破坏我国北方生态系统的昆虫之一。由于蝗灾发生后大量的蝗虫会吞食禾田，使农产品遭到完全破坏，引发严重的经济损失以致因粮食短缺而发生饥荒。因此，准确地检测蝗灾对于蝗灾危害的及时遏制至关重要。
目前，主要利用下面两种方法来检测蝗灾：利用气象数据、中分辨率成像光谱仪(MODIS)数据计算归一化植被指数来检测蝗灾；或者，通过专题制图仪(TM)影像获取各种植被指数数据进行东亚飞蝗的遥感反演来检测蝗灾。然而，这两种方法都是利用单一的影像数据来检测蝗灾，检测的精度不高，难以获得准确的检测结果。
因此，现有的检测蝗灾的方法不能准确地检测蝗灾。
发明内容
本发明的示例性实施例在于提供一种检测蝗灾的方法。所述方法能够提高对蝗灾的检测精度，从而能够准确地检测蝗灾。
根据本发明示例性实施例，提供一种检测蝗灾的方法，包括：(A)获取被检测区域的多源遥感数据，其中，所述多源遥感数据包括：全极化合成孔径雷达数据、中分辨率成像光谱仪温度和专题制图仪植被指数；(B)通过全极化合成孔径雷达数据获得平均极化散射角和极化散射熵；(C)通过平均极化散射角、极化散射熵、中分辨率成像光谱仪温度和专题制图仪植被指数获得差值蝗灾指数；(D)根据差值蝗灾指数来确定被检测区域是否发生蝗灾。
可选地，获取的全极化合成孔径雷达数据为极化散射矩阵，其中，步骤(B)包括：(B1)将所述极化散射矩阵转化为极化相干矩阵；(B2)获得所述极化相干矩阵的特征向量和特征值；(B3)将所述极化相干矩阵转化为旋转矩阵、单位矩阵和旋转矩阵的共轭矩阵的乘积，通过所述旋转矩阵和所述特征向量获得极化散射角，并通过获得的极化散射角计算平均极化散射角；(B4)通过所述特征值获得极化散射熵。
可选地，步骤(C)包括：(C1)分别将平均极化散射角、极化散射熵、中分辨率成像光谱仪温度和专题制图仪植被指数归一化，获得归一化的平均极化散射角、归一化的极化散射熵、归一化的中分辨率成像光谱仪温度和归一化的专题制图仪植被指数；(C2)通过归一化的平均极化散射角、归一化的极化散射熵、归一化的中分辨率成像光谱仪温度和归一化的专题制图仪植被指数获得差值蝗灾指数。
可选地，在步骤(A)中，获取的多源遥感数据包括：第一时刻获取的多源遥感数据和第二时刻获取的多源遥感数据，其中，第一时刻早于第二时刻，其中，在步骤(C)中，通过分别与第一时刻和第二时刻相对应的平均极化散射角、极化散射熵、中分辨率成像光谱仪温度和专题制图仪植被指数获得差值蝗灾指数，在步骤(D)中，根据差值蝗灾指数来确定被检测区域在第二时刻是否发生蝗灾。
可选地，在步骤(D)中，当差值蝗灾指数大于第一预定阈值时，确定被检测区域发生蝗灾，当差值蝗灾指数不大于所述第一预定阈值时，确定被检测区域未发生蝗灾。
可选地，步骤(D)包括：根据差值蝗灾指数和归一化的专题制图仪植被指数来确定被检测区域是否发生蝗灾。
可选地，在步骤(D)中，当差值蝗灾指数大于第一预定阈值且分别与第一时刻和第二时刻相应的归一化的专题制图仪植被指数均大于第二预定阈值时，确定被检测区域在第二时刻发生蝗灾，否则，确定被检测区域在第二时刻未发生蝗灾。
可选地，所述差值蝗灾指数DPLI用下面的公式来表示：
其中，NLIafter为与第二时刻相应的归一化的中分辨率成像光谱仪温度，NDVIafter为与第二时刻相应的归一化的专题制图仪植被指数，NHafter为与第二时刻相应的归一化的极化散射熵，Nαafter为与第二时刻相应的归一化的平均极化散射角，NLIpre为与第一时刻相应的归一化的中分辨率成像光谱仪温度，NDVIpre为与第一时刻相应的归一化的专题制图仪植被指数，NHpre为与第一时刻相应的归一化的极化散射熵，Nαpre为与第一时刻相应的归一化的平均极化散射角。
在根据本发明示例性实施例的检测蝗灾的方法中，通过包含两个光学数据和一个微波数据的多源遥感数据来检测蝗灾，提高了蝗灾检测的精度，从而实现了对蝗灾的准确检测。
附图说明
通过下面结合示例性地示出实施例的附图进行的描述，本发明示例性实施例的上述和其他目的和特点将会变得更加清楚，其中：
图1示出根据本发明的示例性实施例的检测蝗灾的方法的流程图；
图2示出根据本发明的示例性实施例的检测蝗灾的方法中获得平均极化散射角和极化散射熵的步骤的流程图。
具体实施方式
以下，将参照附图更充分地描述本发明的示例性实施例，示例性实施例在附图中示出。然而，可以以许多不同的形式实施示例性实施例，并且不应被解释为局限于在此阐述的示例性实施例。相反，提供这些实施例从而本公开将会彻底和完整，并将完全地将示例性实施例的范围传达给本领域的技术人员。
图1示出根据本发明的示例性实施例的检测蝗灾的方法的流程图。
在步骤S100，获取被检测区域的多源遥感数据。其中，所述多源遥感数据包括：全极化合成孔径雷达(SAR)数据、中分辨率成像光谱仪(MODIS)温度和专题制图仪(TM)植被指数。这里，获取的全极化SAR数据为微波数据，获取的MODIS温度和TM植被指数均为光学数据。
作为示例，为了比较被检测区域在不同时刻的差异，获取的多源遥感数据可包括：第一时刻获取的多源遥感数据和第二时刻获取的多源遥感数据，其中，第一时刻早于第二时刻。
另外，由于全极化SAR数据、MODIS温度和TM植被指数都有其特定的数据格式和投影方式，可对全极化SAR数据、MODIS温度和TM植被指数进行预处理(例如，进行数据的拼接、裁剪、提取、投影以及数值转换)，使全极化SAR数据、MODIS温度和TM植被指数具有一致的数据格式和投影方式。例如，可借助POLSARPRO(SAR数据处理工具)对全极化SAR数据进行预处理、借助MRT(MODIS数据处理工具)对MODIS温度进行预处理、借助ENVI(遥感图像处理平台)对TM植被指数进行预处理，使全极化SAR数据、MODIS温度和TM植被指数具有一致的数据格式和投影方式。
在步骤S200，通过全极化SAR数据获得平均极化散射角和极化散射熵。
获取的全极化SAR数据可为极化散射矩阵，可通过极化散射矩阵及其特征向量和特征值来获得平均极化散射角和极化散射熵。
图2示出根据本发明的示例性实施例的检测蝗灾的方法中获得平均极化散射角和极化散射熵的步骤S200的流程图。
在步骤S210，将所述极化散射矩阵转化为极化相干矩阵。
在步骤S220，获得所述极化相干矩阵的特征向量和特征值。
作为示例，可由下面的公式获得所述极化相干矩阵的特征向量和特征值：
其中，T为所述极化相干矩阵，λi为所述极化相干矩阵T的特征值，ei为所述极化相干矩阵的特征向量，为ei的伴随矩阵，i的值域是{1,2,3}。
在步骤S230，将所述极化相干矩阵转化为旋转矩阵、单位矩阵和旋转矩阵的共轭矩阵的乘积，通过所述旋转矩阵和所述特征向量获得极化散射角，并通过获得的极化散射角计算平均极化散射角。
作为示例，可用下面的公式来获得所述极化散射角αi：
其中，ei为极化相干矩阵的特征向量，sinαi、cosαi、sinβi和cosβi均为所述旋转矩阵中的元素，和均为所述极化相干矩阵中的元素。
并且，可通过下面的公式来计算平均极化散射角
其中，i的值域是{1,2,3}。
在步骤S240，通过所述特征值获得极化散射熵。
作为示例，可用下面的公式来获得所述极化散射熵H：
其中，λi、λj均为所述极化相干矩阵的特征值，i和j的值域均是{1,2,3}。
再次参照图1，在步骤S300，通过平均极化散射角、极化散射熵、MODIS温度和TM植被指数获得差值蝗灾指数。
这里，通过由微波数据获得的平均极化散射角和极化散射熵，以及相关的光学数据：MODIS温度和TM植被指数可获得用于检测蝗灾的差值蝗灾指数。
作为示例，可先分别将平均极化散射角、极化散射熵、MODIS温度和TM植被指数归一化，获得归一化的平均极化散射角、归一化的极化散射熵、归一化的MODIS温度和归一化的TM植被指数。
具体地说，由于平均极化散射角、极化散射熵、MODIS温度和TM植被指数的量纲存在差异，不同数据的值域之间相差巨大。为避免较大值对较小值的遮盖，可将平均极化散射角、极化散射熵、MODIS温度和TM植被指数进行归一化，使所有数据值域均为[0,1]，以消除量纲的影响。
例如，可用下面的公式来归一化：
其中，Xi'表示归一化后的数据，Xi为原始数据(需要归一化的数据)，Xmax为需要归一化的所有数据中的最大值数据，Xmin为需要归一化的所有数据中的最小值数据。
接下来，可通过归一化的平均极化散射角、归一化的极化散射熵、归一化的MODIS温度和归一化的TM植被指数获得差值蝗灾指数。
作为另一示例，当获取的多源遥感数据包括第一时刻获取的多源遥感数据和第二时刻获取的多源遥感数据时，可通过分别与第一时刻和第二时刻相对应的平均极化散射角、极化散射熵、MODIS温度和TM植被指数获得差值蝗灾指数。在这种情况下，可先将分别与第一时刻和第二时刻相应的平均极化散射角、极化散射熵、MODIS温度和TM植被指数归一化，获得分别与第一时刻和第二时刻相应的归一化的平均极化散射角、归一化的极化散射熵、归一化的MODIS温度和归一化的TM植被指数。然后，通过分别与第一时刻和第二时刻相应的归一化的平均极化散射角、归一化的极化散射熵、归一化的MODIS温度和归一化的TM植被指数获得差值蝗灾指数。
所获得的差值蝗灾指数DPLI可用下面的公式来表示：
其中，NLIafter为与第二时刻相应的归一化的MODIS温度，NDVIafter为与第二时刻相应的归一化的TM植被指数，NHafter为与第二时刻相应的归一化的极化散射熵，Nαafter为与第二时刻相应的归一化的平均极化散射角，NLIpre为与第一时刻相应的归一化的MODIS温度，NDVIpre为与第一时刻相应的归一化的TM植被指数，NHpre为与第一时刻相应的归一化的极化散射熵，Nαpre为与第一时刻相应的归一化的平均极化散射角。
由于发生蝗灾时，归一化的MODIS温度会增大，归一化的TM植被指数、归一化的极化散射熵和归一化的平均极化散射角均会减小。因此，当蝗灾发生时，差值蝗灾指数DPLI会增大，所以，差值蝗灾指数DPLI可反映蝗灾情况。
在步骤S400，根据差值蝗灾指数来确定被检测区域是否发生蝗灾。
可将差值蝗灾指数与预定的阈值进行比较，根据比较结果来确定被检测区域是否发生蝗灾。
在一示例中，当差值蝗灾指数大于第一预定阈值时，确定被检测区域发生蝗灾，当差值蝗灾指数不大于所述第一预定阈值时，确定被检测区域未发生蝗灾。
在另一示例中，当获取的多源遥感数据包括第一时刻获取的多源遥感数据和第二时刻获取的多源遥感数据时，可根据差值蝗灾指数来确定被检测区域在第二时刻是否发生蝗灾(即，确定较晚的时刻是否发生蝗灾)。
作为示例，当差值蝗灾指数大于第一预定阈值时，确定被检测区域在第二时刻发生蝗灾，当差值蝗灾指数不大于所述第一预定阈值时，确定被检测区域在第二时刻未发生蝗灾。
作为优选示例，根据差值蝗灾指数和归一化的TM植被指数来确定被检测区域是否发生蝗灾。
这里，当只需要检测农田区域的蝗灾情况时，可首先排除荒漠、草地及湿地的影响。由于荒漠、草地及湿地的TM植被指数通常比较低，因此，可通过设置TM植被指数的阈值来排除荒漠、草地及湿地。
具体地说，当差值蝗灾指数大于第一预定阈值且分别与第一时刻和第二时刻相应的归一化的TM植被指数均大于第二预定阈值时，确定被检测区域在第二时刻发生蝗灾，否则，确定被检测区域在第二时刻未发生蝗灾。
这里，第一阈值和第二阈值可根据试验的方式来确定。例如，可根据不同的环境特点，通过试验的方式来确定。
在根据本发明示例性实施例的检测蝗灾的方法中，通过包含两个光学数据和一个微波数据的多源遥感数据来检测蝗灾，提高了蝗灾检测的精度，从而实现了对蝗灾的准确检测。
尽管已经参照其示例性实施例具体显示和描述了本发明，但是本领域的技术人员应该理解，在不脱离权利要求所限定的本发明的精神和范围的情况下，可以对其进行形式和细节上的各种改变。
Claims (6)
1.一种检测蝗灾的方法，包括：
(A)获取被检测区域的多源遥感数据，其中，所述多源遥感数据包括：全极化合成孔径雷达数据、中分辨率成像光谱仪温度和专题制图仪植被指数；获取的多源遥感数据包括：第一时刻获取的多源遥感数据和第二时刻获取的多源遥感数据，其中，第一时刻早于第二时刻，
(B)通过全极化合成孔径雷达数据获得平均极化散射角和极化散射熵；
(C1)分别将与第一时刻和第二时刻相对应的平均极化散射角、极化散射熵、中分辨率成像光谱仪温度和专题制图仪植被指数归一化，获得分别与第一时刻和第二时刻相对应的归一化的平均极化散射角、归一化的极化散射熵、归一化的中分辨率成像光谱仪温度和归一化的专题制图仪植被指数；
(C2)通过分别与第一时刻和第二时刻相应的归一化的平均极化散射角、归一化的极化散射熵、归一化的中分辨率成像光谱仪温度和归一化的专题制图仪植被指数获得差值蝗灾指数；
(D)根据差值蝗灾指数来确定被检测区域在第二时刻是否发生蝗灾。
2.如权利要求1所述的方法，其中，获取的全极化合成孔径雷达数据为极化散射矩阵，
其中，步骤(B)包括：
(B1)将所述极化散射矩阵转化为极化相干矩阵；
(B2)获得所述极化相干矩阵的特征向量和特征值；
(B3)将所述极化相干矩阵转化为旋转矩阵、单位矩阵和旋转矩阵的共轭矩阵的乘积，通过所述旋转矩阵和所述特征向量获得极化散射角，并通过获得的极化散射角计算平均极化散射角；
(B4)通过所述特征值获得极化散射熵。
3.如权利要求1所述的方法，其中，在步骤(D)中，当差值蝗灾指数大于第一预定阈值时，确定被检测区域在第二时刻发生蝗灾，当差值蝗灾指数不大于所述第一预定阈值时，确定被检测区域未发生蝗灾。
4.如权利要求1所述的方法，其中，步骤(D)包括：根据差值蝗灾指数和归一化的专题制图仪植被指数来确定被检测区域在第二时刻是否发生蝗灾。
5.如权利要求4所述的方法，其中，在步骤(D)中，当差值蝗灾指数大于第一预定阈值且分别与第一时刻和第二时刻相应的归一化的专题制图仪植被指数均大于第二预定阈值时，确定被检测区域在第二时刻发生蝗灾，否则，确定被检测区域在第二时刻未发生蝗灾。
6.如权利要求1所述的方法，其中，所述差值蝗灾指数DPLI用下面的公式来表示：
<mrow>
<mi>D</mi>
<mi>P</mi>
<mi>L</mi>
<mi>I</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>NLI</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>NDVI</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>NH</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>N&alpha;</mi>
<mrow>
<mi>a</mi>
<mi>f</mi>
<mi>t</mi>
<mi>e</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>NLI</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>NDVI</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>NH</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>N&alpha;</mi>
<mrow>
<mi>p</mi>
<mi>r</mi>
<mi>e</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中，NLIafter为与第二时刻相应的归一化的中分辨率成像光谱仪温度，NDVIafter为与第二时刻相应的归一化的专题制图仪植被指数，NHafter为与第二时刻相应的归一化的极化散射熵，Nαafter为与第二时刻相应的归一化的平均极化散射角，NLIpre为与第一时刻相应的归一化的中分辨率成像光谱仪温度，NDVIpre为与第一时刻相应的归一化的专题制图仪植被指数，NHpre为与第一时刻相应的归一化的极化散射熵，Nαpre为与第一时刻相应的归一化的平均极化散射角。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410852062.5A CN104504279B (zh) | 2014-12-31 | 2014-12-31 | 检测蝗灾的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410852062.5A CN104504279B (zh) | 2014-12-31 | 2014-12-31 | 检测蝗灾的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104504279A CN104504279A (zh) | 2015-04-08 |
| CN104504279B true CN104504279B (zh) | 2017-09-15 |
Family
ID=52945676
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410852062.5A Active CN104504279B (zh) | 2014-12-31 | 2014-12-31 | 检测蝗灾的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104504279B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111368716B (zh) * | 2020-02-28 | 2023-04-28 | 重庆市规划和自然资源调查监测院 | 一种基于多源时空数据的地质灾害灾毁耕地提取方法 |
| CN114419431B (zh) * | 2021-12-23 | 2024-07-19 | 深圳先进技术研究院 | 一种蝗灾潜在高风险区识别方法、装置、设备及存储介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1580764A (zh) * | 2003-08-08 | 2005-02-16 | 中国科学院遥感应用研究所 | 基于大规模爆发性农害昆虫的生育周期的虫灾监测方法 |
| CN102937574A (zh) * | 2012-07-20 | 2013-02-20 | 北京农业信息技术研究中心 | 一种基于卫星影像的病虫害信息提取方法 |
| CN103034910A (zh) * | 2012-12-03 | 2013-04-10 | 北京农业信息技术研究中心 | 基于多源信息的区域尺度病虫害预测方法 |
| CN103955606A (zh) * | 2014-04-23 | 2014-07-30 | 北京大学 | 一种基于遥感技术的草原蝗灾渐进式预测方法 |
-
2014
- 2014-12-31 CN CN201410852062.5A patent/CN104504279B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1580764A (zh) * | 2003-08-08 | 2005-02-16 | 中国科学院遥感应用研究所 | 基于大规模爆发性农害昆虫的生育周期的虫灾监测方法 |
| CN102937574A (zh) * | 2012-07-20 | 2013-02-20 | 北京农业信息技术研究中心 | 一种基于卫星影像的病虫害信息提取方法 |
| CN103034910A (zh) * | 2012-12-03 | 2013-04-10 | 北京农业信息技术研究中心 | 基于多源信息的区域尺度病虫害预测方法 |
| CN103955606A (zh) * | 2014-04-23 | 2014-07-30 | 北京大学 | 一种基于遥感技术的草原蝗灾渐进式预测方法 |
Non-Patent Citations (4)
| Title |
|---|
| Temporal filtering of successive MODIS data in monitoring a locust outbreak;Y. ZHA et al.;《International Journal of Remote Sensing》;20051220;第26卷(第24期);5665-5674 * |
| 基于Landsat8卫星数据的蝗虫遥感监测方法研究;黄健熙等;《农业机械学报》;20141225;1-13 * |
| 极化与干涉信息融合的SAR影像林地信息提取;谢雯君;《中国优秀硕士学位论文全文数据库 农业科技辑》;20140315;D049-141 * |
| 集成GIS和RS技术的作物病虫害监测预报研究进展;宫彦萍等;《农业机械学报》;20080930;第24卷(第增刊2期);331-336 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104504279A (zh) | 2015-04-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pölönen et al. | Hyperspectral imaging based biomass and nitrogen content estimations from light-weight UAV | |
| Srivastava et al. | Tests for covariance matrices in high dimension with less sample size | |
| US9373052B2 (en) | Shape and photometric invariants recovery from polarisation images | |
| CN105352895B (zh) | 高光谱遥感数据植被信息提取方法 | |
| Koukal et al. | Evaluation of semi-empirical BRDF models inverted against multi-angle data from a digital airborne frame camera for enhancing forest type classification | |
| CN110412573A (zh) | 极化sar图像分解方法及存储介质 | |
| CN104504279B (zh) | 检测蝗灾的方法 | |
| CN104316491A (zh) | 基于同步-异步二维近红外相关谱检测牛奶掺尿素的方法 | |
| CN109298402B (zh) | 基于通道融合的极化特征提取方法 | |
| CN114637001B (zh) | 雷达的目标检测方法、装置、设备及存储介质 | |
| CN106772371A (zh) | 基于极化干涉sar分类应用的极化定标参数需求分析方法 | |
| CN102914774A (zh) | 一种综合孔径辐射计图像反演方法 | |
| Xu et al. | Comparison of NDVI data between GF-1 PMS1 and ZY-3 MUX sensors | |
| Zhao et al. | A channel imbalance calibration scheme with distributed targets for circular quad-polarization SAR with reciprocal crosstalk | |
| Van Niel et al. | A simple method to improve field-level rice identification: toward operational monitoring with satellite remote sensing | |
| Yang et al. | An improved EM algorithm for remote sensing classification | |
| CN108710816B (zh) | 一种基于概率融合的极化sar图像建筑物检测方法 | |
| CN108008378A (zh) | 一种sar虚假目标能量比值确定方法及装置 | |
| CN109374652B (zh) | 基于毫米波信号的湿度感知与检测方法及系统 | |
| CN111076817A (zh) | 基于光学多维信息一体化感知系统的异常目标检测方法 | |
| CN102879770B (zh) | 基于sal回波数据的目标振动检测方法 | |
| Ferro-Famil et al. | Forest mapping and classification at L-band using Pol-inSAR optimal coherence set statistics | |
| Cloude et al. | Forest structure estimation using polinsar | |
| CN105184753B (zh) | 一维信号处理引导下的遥感影像条带噪声快速滤除方法 | |
| CN102708544B (zh) | 用于高光谱遥感图像逐波段修正负相关的光谱角填图方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |