春季是农作物出苗的时期，地表异常干燥松散，抗风蚀能力很弱，在有大风刮过时，就会造成农作物大面积倒伏、植株折断、籽粒或果实脱落，减产严重。2009 年 4 月 16 日，阿克苏市遭遇大风沙尘暴天气，持续时间近 3 小时，瞬间风速达到 24． 6m/s。这次大风沙尘天气风力强、范围广、持续时间长、影响大。传统的风沙灾害监测与评估方法主要通过大田抽样调查的方法估算受灾面积和其分布情况。这种方法需要大量的田间调查工作，很难在区域尺度上实现灾害的及时准确监测。因此，快速获取风沙灾害发生区域的范围、受灾程度等信息，客观地反映风沙灾害发生的真实情况，并制定一套实时、快捷、方便、有效的风沙灾害监测方法已迫在眉睫。
由于遥感数据能够快速、准确、大范围评估农作物的灾害情况，使得卫星资料监测自然灾害受到广泛关注，并已成为研究热点。Tapia － Silva 等利用农作物的归一化植被指数标准生长曲线和遥感数据分类方法估计农作物受洪水灾害的损失; Baojuan Zheng 等利用多时相 Landsat 数据计算归一化耕作指数分析了作物残茬覆盖情况; 王堃等通过定量分析收获季节内多时相 NDVI ( Normalized Differential VegetationIndex，NDVI) 的变化大小和变化过程，实时监测研究区内玉米种植地块的收获期; 南颖等通过对比分析台风前后长白山归一化植被指数，提取了长白山自然保护区内遭受台风袭击后的灾区范围; 董燕生等根据冻害发生的时间和冬小麦的生长特征选择冬小麦越冬前分蘖期和年后起身期作为冻害分析的关键生长期，基于多时相 HJ － 1 CCD 遥感数据和统计方法分析了冬小麦冻害影响范围和受灾程度的遥感评估。
基于以上灾情监测方法，运用多时相遥感数据的灾情评估指标变化分析是监测农作物灾情的重要方法。考虑到风沙灾害监测的时效性，中分辨率成像光谱仪( Moderate － resolution Imaging Spectroradiometer，MODIS) 具有高时间分辨率和数据空间覆盖范围广的特点，在农作物遥感监测方面具有独特的优势。文中以阿克苏市为例，选择多时相 MODIS 数据产品，通过建立风沙灾害评估指标，获得灾情影响的农作物种植面积，为防灾减灾机构提供有力的信息支持。
1 材料与研究方法
1． 1 研究区概况
阿克苏市位于新疆天山南麓、塔克拉玛干沙漠西北边缘，属暖温带干旱气候地区，降雨量稀少，蒸发量大，气候干燥，适宜各类农作物生长，是国家重要的商品粮、商品棉基地。2000 ～2010 年，阿克苏市发生沙尘天气共 582 天，8 级以上大风 82 次，其中以 4 月次数最多。风沙灾害已成为阿克苏市农作物出苗期主要自然灾害之一，严重制约了农作物生产的发展。
1． 2 数据与预处理
( 1) 研究时间选择。根据阿克苏市地面气象观测站提供的总云量资料，选取无云的时间，即 4 月 15日( 农作物风沙灾害前) 、4 月 17 日( 农作物风沙灾害后) 和 4 月 29 日( 农作物风沙灾害恢复期) 作为农作物风沙灾害研究时间。为了对风沙灾害指标阈值的确定，需要选择季节相近的研究区农作物正常生长晴空图，时间分别为 2010 年 4 月 22 日、23 日和 24 日。
( 2) MODIS 数据及预处理。根据研究时间，此次共使用 2009 年 4 月 15 日、17 日和 29 日，以及 2010年 4 月 22 日、23 日和 24 日的 250 m MODIS 日地表反射率产品( MOD09GQ) 影像 6 期，2009 年逐年合成500m 土地覆盖类型产品 MCD12Q1 影像 1 期。对这些遥感数据均采用 MODIS 重投影工具进行重采样和拼接处理，转换成 tif 格式。根据土地覆盖类型产品中第一类国际地圈生物圈计划( IGBP) 全球植被分类方案，选择第 12 类( 农作物) 提取阿克苏市的农作物种植区。
1． 3 研究方法
1． 3． 1 风沙灾害植被指数变化情况
利用遥感监测风沙对农作物的灾害，通常是利用地表植被状况、特征等指标监测作物风沙灾害的情况。参考文献中指出，NDVI 探测低密度绿色植被的能力最强，是植被生长状态及植被覆盖度的最佳指示因子。因此，文中将农作物的 NDVI 值作为农作物风沙灾害的评估指标。
( 1) 归一化植被指数。根据每日地表反射率数据( MOD09GQ) 的近红外和红波段，通过公式( 1) 计算归一化植被指数。【1】
式中: NDVI 为归一化植被指数; ＲNIＲ为近红外波段反射率; ＲＲED为红波段反射率。
( 2) 植被指数变化分析。从光谱物征与作物冠层之间的关系看，绿色植物叶面叶绿素在光照条件下发生光合作用，在可见光红光波段有很强的吸收特性，而对近红外波段具有很强的反射率、高的透射率和极低的吸收，即红光波段反射率包含了植物冠顶层叶片的丰富信息，近红外光反射率包含了整个植物冠层内叶片的大部分信息，这是植被遥感监测的物理基础。农作物生长初期，叶绿素含量增多，红波段吸收和近红外反射率增强，植被指数增大。当农作物受到风沙灾害影响后，会出现农作物活性减弱、死苗等现象，红波段吸收和近红外反射率减弱，植被指数减小。因此，不同日期的数值相减是为了观测两个日期植被的变化，并用来探测由于风沙灾害造成的农作物变化。
文中将相邻时期 NDVI 两两相减后，用公式( 2) 逐像元判断其是否发生变化，即: 若△NDVI( i，j)＞ a，表示正向变化; 若△NDVI( i，j)＜ a，像元负向变化【2】
式中: NDVId1( i，j)和 NDVId2( i，j)分别是像元( i，j) 在日期 d1 和 d2 的 NDVI 值; △NDVI( i，j)是像元( i，j) 在日期 d1 和 d2 的 NDVI 差值。
由于 NDVI 值的变化不仅受植物生长和地物变化的影响，还会受到太阳光照强度、大气状况、卫星拍摄角度等因素影响，相邻日期的 NDVI 值也会存在较小的波动，所以 NDVI 变化阈值 a 的确定是一件较为困难而又关键的问题。根据文献，选择与受灾时间同时相的遥感数据，计算研究区农作物春季正常生长 NDVI 变化，定义风沙灾害的阈值为公式( 3) :【3】
式中: △NDVI、σ 分别为同时期正常生长年份研究区农作物未受灾害的植被指数变化的均值和标准差。
1． 3． 2 农作物风沙灾害的遥感监测方法
基于多时相遥感数据，农作物风沙灾害的遥感监测方法分为 4 个步骤:
( 1) 选择遥感数据及预处理。根据灾害发生的时间，选择灾前、灾后和恢复期 3 个时相的遥感数据，以及同时期正常生长年份的遥感数据，并进行预处理。
( 2) 计算△NDVI 的均值和标准差。根据受灾年份的灾前、灾后和恢复期的遥感数据，以及正常年份相同时相的遥感数据，分别计算研究区域内农作物△NDVI 的均值和标准差。
( 3) 根据公式( 3) 确定农作物受灾阈值 a。当农作物的△NDVIi，j＜ a 时，表示此像元受到风沙灾害的影响。
( 4) 确定受灾面积。农作物风沙灾害受灾面积通过不同类别风沙灾害影响范围内像元的统计得到。农作物受灾面积计算公式: Ai= A0× Pi( 4)式中: Ai为第 i 类风沙灾情的作物受灾面积; A0为农作物面积; Pi为不同风沙灾害类别所占比例。
2 结果与分析
2． 1 风沙灾害前后植被指数的变化情况
从图 1 中可以看出，阿克苏农作物在风沙灾害前后和恢复期 NDVI 分布逐渐趋于扁平状。NDVI最大值分别为 0． 36、0． 64、0． 73，80% 的农作物 ND-VI 分别在 0． 13、0． 19、0． 26 以下，其中: 占总作物总面积比例最大的 NDVI 值逐渐增大，由 0． 08( 17．51% ) 、0． 09 ( 14． 42% ) 到 0． 11 ( 9． 63% ) 。在风沙灾害前后和恢复期，NDVI 面积变化最为明显的区间是 NDVI ＜0． 1 和0． 15 ＜ NDVI ＜0． 4。NDVI ＜0． 1的面积变化明显，分别为 61． 21%、43． 24%、20．96% ; 0． 15 ＜ NDVI ＜ 0． 4 的面积分别为 11． 93% 、30． 21% 、39． 83% 。【图1】
2． 2 像元植被指数的变化情况
由于不同年份农作物物候的变化，文中基于2010 年4 月22 日、23 日和24 日的遥感数据，对多天相邻遥感数据风沙灾害阈值取平均作为风沙灾害的阈值。利用风沙灾害阈值的确定方法，分别计算研究区农作物 a23 － 22和 a24 － 23，并取它们的均值 a = ( a23 － 22+ a24 － 23) /2 作为风沙灾害的阈值。因此，文中阿克苏市农作物风沙灾害阈值为 －0． 005，即△NDVI( i，j)＜ － 0． 005 的区域被划定为农作物受灾区。
在风沙灾害前、后和恢复期，若农作物植被指数始终保持正向变化，无风沙危害症状，则农作物未受到风沙灾害的影响，保持正常生长; 否则，为农作物异常变化，受到风沙灾害的影响。因此，根据阿克苏市农作物在风沙灾害前、后和恢复期植被指数的变化，将农作物风沙灾害划分为三类: 负向 － 正向型、正向 － 负向型和负向 － 负向型。【表1】
2． 2． 1 负向 － 正向型
这种风沙灾害类型在风沙灾害后，植被指数减小，其中: NDVI ＞0． 15 的农作物 NDVI 都减小 0． 02 以上，NDVI 在 0． 3 ～0． 35 的农作物减小最大，为 0． 037，NDVI 在 0． 05 ～0． 1 的农作物减小最少，为 0． 013，平均减小 0． 022; 而在风沙灾害恢复期，植被指数逐渐增大，由 NDVI 在 0． 05 ～0． 1 的农作物增长 0． 022 到NDVI 在 0． 30 ～ 0． 35 的农作物增长 0． 233，平均增长 0． 11。由此可见，负向 － 正向型农作物风沙灾害的NDVI 越大，风沙灾害后和其恢复期变化就越大。处于生长子叶或真叶期的农作物，此时抗风沙能力较弱，重者主茎上部被风刮断，轻者叶片被风刮碎，NDVI 减小。但这些农作物具有较强的恢复能力，在风沙灾害后不久就开始恢复生长，NDVI 逐渐增大。
2． 2． 2 正向 － 负向型
农作物受到风沙灾害后好像正常生长，NDVI 增加的范围从 0． 03 到 0． 189; 但随着时间的推移，农作物受风沙灾害的影响也逐渐显现出来，NDVI 开始减小，NDVI 减小的范围从 0． 02 到 0． 065。由于风力较大，已播种的农作物地膜不同程度被风刮起，根系尚未健全的农作物幼苗裸露在外，增加了 NDVI 值。当风沙吹打农作物幼苗后，幼苗子叶和幼茎上的伤口受到大风长时间的吹干，叶绿素降低，导致 NDVI 逐渐减小。
2． 2． 3 负向 － 负向型
NDVI 在 0 ～ 0． 25 的农作物，特别是 NDVI 在 0 ～ 0． 05 的农作物容易受到这种类型的风沙灾害影响。
农作物在这种类型的风沙灾害后，NDVI 在 0 ～0． 05 的农作物减小 0． 067，恢复期 NDVI 减小 0． 045。沙尘覆盖在麦苗、棉花等大田农作物叶面上，影响叶面正常的光合作用和呼吸作用，使作物生长和发育受到抑制或损害。此外，受伤的农作物体内水分会从伤口散失，其子叶、生长点和靠近子叶的幼茎失水干枯，ND-VI 不断减小。
2． 3 农作物受灾面积
根据使用图像的分辨率大小，利用植被指数变化的类别，分别统计风沙对农作物不同影响的面积。
2009 年 4 月 16 日阿克苏市农作物风沙灾害面积 12697． 07hm2，这与实际调查数据基本相符。
阿克苏市三类风沙灾害的受灾情况是: 负向 － 正向型农作物受灾面积 4464． 90hm2，占总受灾面积的35． 16% ; 正向 － 负向型农作物受灾面积最大，为8006． 77hm2，占总受灾面积的 63． 06%; 负向 － 负向型农作物受灾面积最小，为 225． 39hm2，占总受灾面积的 1． 78%。其中: 在负向 － 正向型风沙灾害中，NDVI 在 0． 05 ～ 0． 1 之间的农作物受灾面积是1969． 49hm2，占这类风沙灾害面积的 44． 11%; 在正向 － 负向型风沙灾害中，NDVI 在 0． 05 ～0． 1 之间的农作物，受灾面积是 5082． 05hm2，占这类风沙灾害面积的 63． 47%; 在负向 － 负向型风沙灾害中，NDVI在 0． 1 ～0． 15 之间的农作物受灾面积是 96． 60hm2，占这类风沙灾害面积的 42． 86%。因此，在春季风沙灾害中应重点关注 NDVI 在 0． 05 ～ 0． 15 之间的农作物，并根据不同风沙灾害类型进行抗灾救灾工作。【图2】
3 讨论
目前，对农作物的风沙灾害监测没有明确的方法，而研究农作物长势监测的常用参数是 NDVI，因此，文中以 2009 年 4 月 16 日阿克苏市农作物为例，提出了基于 NDVI 的农作物风沙灾害阈值确定方法，并将受灾农作物分为正向 － 负向型、负向 － 正向型和负向 － 负向型三种类型。农作物风沙灾害遥感监测过程中，通过遥感数据处理和统计分析，得出农作物的受灾情况，方法简单，可操作性强，可以解决区域尺度上的及时监测问题。遥感监测结果显示，阿克苏市农作物风沙灾害面积与实际受灾面积的提取精度为 98． 24%，表明文中提出的农作物风沙灾害遥感监测方法是可行的。
文中主要考虑了农作物风沙灾害快速评估的要求，然而，当农作物遭受风沙灾害后，对农作物生育进程、生物量、质量以及产量的影响等还需要进一步地分析。
4 结论
( 1) 春季风沙灾害前后和恢复期阿克苏市农作物 NDVI 值变化明显，农作物 NDVI 最大值从 0． 36 增大到 0． 73，占农作物总面积比例最大的 NDVI 值从 0． 08 增大到 0． 11，其中: NDVI ＜0． 1 的农作物面积从61． 21% 大幅减小到 20． 96% ，0． 15 ＜ NDVI ＜ 0． 4 的农作物面积从 11． 93% 增加到 39． 83% 。
( 2) 基于风沙灾害前后和恢复期农作物 NDVI 值的变化，将农作物风沙灾害划分为三类: 正向 － 负向型、负向 － 正向型和负向 － 负向型。不同类型风沙灾害的农作物 NDVI 区间值分析显示，负向 － 正向型和正向 － 负向型农作物 NDVI 值越大，风沙灾害后和其恢复期变化就越大; 负向 － 负向型农作物 NDVI 值越小，风沙灾害后和其恢复期变化就明显。
( 3) 阿克苏市受灾农作物面积 12697． 07hm2，其中: 负向 － 正向型农作物受灾面积为 4464． 90hm2，正向 － 负向型农作物受灾面积8006． 77hm2，负向 － 负向型农作物受灾面积225． 39hm2，预防农作物风沙灾害的重点是 NDVI 值为 0． 05 ～0． 15 的农作物。
参考文献
［1］Tapia － Silva，F O，Itzerott S，Frster S，Kuhlmann B，et al． Estimation of flood losses to agricultural crops using remote sensing［J］． Physicsand Chemistry of the Earth，2011，36: 253 － 265．
［2］Zheng B J，Campbell J B，DeBeurs K M． Ｒemote sensing of crop residue cover using multi － temporal Landsat imagery［J］． Ｒemote Sensing ofEnvironment，2012，117: 177 － 183．
［3］王堃，顾晓鹤，程耀东，等． 基于变化向量分析的玉米收获期遥感监测［J］． 农业工程学报，2011，27( 2) : 180 －186．
［4］南颖，刘志锋，朱卫红，等． 长白山的台风灾区范围及植被恢复动态 － 基于遥感技术的研究［J］． 自然灾害学报，2011，20( 6) : 131 －139．
［5］董燕生，陈洪萍，王慧芳，等． 基于多时相环境减灾卫星数据的冬小麦冻害评估［J］． 农业工程学报，2012，28( 20) : 172 －179．
［6］邬明权，牛铮，王长耀． 利用遥感数据时空融合技术提取水稻种植面积［J］． 农业工程学报，2010，26( S2) : 48 －52．
