近年来,气象卫星在火情监测方面的应用得到了广泛开展,已成为监测森林、草原火情、农村秸秆焚烧等日常业务重要组成部分,它为防火部门提供的火情监测信息主要是火点位置和范围。
根据EOS/MODIS和气象卫星AVHRR光谱特征和通道特性分析监测秸秆焚烧、森林火点的面积、位置,国内许多省级气象部门已开展此项工作[1-4]。邹春辉等[5]采用AVHRR对秸秆焚烧进行全程动态监测,对火点进行判识并计算火点面积,实现了地理信息和遥感火点信息、遥感火点监测图像的集成显示与查询,用1:25×104地理信息资料建立河南省遥感防火地理信息平台。然而省级产品火点行政定位只能精确到市县,对火点的判识指标也是针对全省范围的,对小范围的地区来说,难免造成漏判和误判。实际工作中往往需要对火点进行更精确的行政定位,并且在火点的判识上由于下垫面性质不同,判识指标的精确性还需进一步提高。鹤壁市气象局2004年安装了DVBS接收系统,为地区级火情监测技术的提高提供了充足的数据源。通过对火点判识指标进行细化和验证,提高了火点判识的准确度,减少了漏判和误判。利用3S技术叠加大比例尺地图,对火点位置进行更准确定位,在2004年和2005年为地方政府服务方面取得了良好效果。
1 AVHRR光谱特征和火点监测原理 1.1 AVHRR光谱特征根据斯蒂芬-波尔兹曼定律:
其中
黑体的全波长辐射能N与温度的四次方成正比。也就是说黑体温度只要有很小的变化,就会引起辐射的很大变化。而高温热源的温度将会引起辐射的更急剧增加,这种变化十分有利于判识高温热源。
极轨气象卫星AVHRR红外扫描辐射计的第3、4、5通道的波长范围分别为3.55~3.95μm、10.3~11.3μm、11.5~12. 5μm。
1.2 火点监测原理根据维恩位移定律:
黑体温度T和辐射峰值波长λmax成反比,即温度愈高,辐射峰值波长愈小,常温(约300K)地表辐射峰值波长在通道4波长范围左右,当温度增高时,尤其髙于500K后,辐射峰值波长迅速移向通道3波长范围。火点温度一般在550K以上,其热辐射峰值波长在3~6μm之间,靠近通道3波长范围。由普朗克辐射公式,对于相同的温度增量,辐射率的增量在通道3波长范围与通道4是不同的,通道3辐射率增长高于通道4, 而且温度愈高,通道3的增长率愈大,尤其当物体从300K常温增高到500K以上高温时,通道3辐射率增高数百倍,而通道4仅增高十几倍。利用卫星这一特性可以很容易把高温热源(火点)从背景中判识出来。
2 资料处理流程和判识指标 2.1 资料来源与处理本文所用卫星遥感资料来自国家卫星气象中心通过DVBS下发AVHRR卫星资料,实际火点资料来自鹤壁市秸秆禁烧办公室2004、2005年全市秸秆焚烧统计资料。极轨卫星遥感资料处理流程分为数据接收、数据预处理、图像处理等,其中预处理包括数据定标、几何校正、投影变换、数据融合等。火点资料挑选每天卫星过境时间前后着火面积在0.067hm2以上火点作为实例验证。
2.2 火点判识指标的确定及验证 2.2.1 云检测进行火点判识之前首先要进行云检测。云顶在可见光波段反照率在30%~70%,远大于陆地的反照率(20%以下),在白天可以利用这一特性把云剔除出来。通道4一般反映的是地表放射辐射,可以给出地表温度变化,云区和陆地在通道4的亮温值有显著差异,在夜间可利用通道4大于某一临界值来区别云区和非云区,这些值的大小分别随太阳高度角和季节而变化。
白天云检测判识指标:CH1m <CH1c;
夜间云检测判识指标:CH4m>CH4c。
此处CH1m为被判识像元通道1的反射率,CH4m为被判识像元通道4的亮温,CH4c,CH4c分别为通道1和通道4的云判识门槛值。
2.2.2 门檻值的精细化确定判识火点的关键在于确定合适的门槛值。由于不同子区下垫面性质不同,从而引起在同一幅图像中各子区亮温平均值不同及距平值不同。因此,对不同子区判识火点的门檻值也应该不同。即便对于同一子区,在不同季节不同时次各通道距平值也可能不同。因此,需要进行统计分析,建立每一子区在不同季节不同时次判识火点的门槛值。
根据下垫面性质,把鹤壁地区分成3个子区:西部山地丘陵区为Ⅰ区,占总面积的35%;中部农田区为Ⅱ区,占总面积的60%;东部沙土区为Ⅲ区,占总面积的5%。3个子区在通道4和通道5中的亮温值差异较小,但在通道3中差值较大。在实际工作中,通过对不同的子区根据季节时次不同确定不同的门檻值,可大大提高火点的判识准确度。通过对2004和2005年秋季秸秆焚烧期间不同子区不同时次在遥感通道的不同特性进行统计分析,结合两年来实际火点进行验证,确定了着火区、亚像元区、过火区判识指标(见表 1)。
利用精细化的人机交互判识方式,可以识别微弱的火点信息。在实际工作中最小可以识别出0.067hm2大小的火点,大大提高了火点的判识精度。
由于亚像元效应,以上判识条件还不能完全排除对云边缘处像元的误判。因为在云边缘处的像元中,既含有云信息,也含有陆地信息。由于有云的信息,CH3资料中太阳辐射反射将使得其亮温值高于陆地像元; 而CH4资料中陆地可能占有较大比例,亮温仍较高,达不到被云剔除的门槛值;CH1也是同样情况,这样就使得有些云边缘处像元容易被识别为火点。小范围内,可以通过精细选择CH3背景亮温与CH4背景亮温差异值,以及CH1、CH4云判识门檻值,从而避开云边缘处带来的干扰,提高火点判识的精度。
2.3 火点面积的计算AVHRR为10比特数据,量化等级为0~1023, 辐射分辨率较高,而空间分辨率为1.1km。卫星观测到的像元辐射率是像元范围内所有各部分地物辖射率累加的加权平均值,即:
式中Nt为卫星观测到的像元辐射率,t为辐射率Nt对应的亮温,ΔSi为像元中第i个小子区面积,NTi为该子区的辐射率,为该子区温度,S为像元总面积。
在日常火情监测中,经常监测到数个或数十个像元的火点。如果以像元分辨率表示明火点面积,则明显夸大了明火的实际面积。AVHRR红外通道的亮温可以反映在一定范围内不同面积不同温度的火点。因此,火点的像元数与明火的实际面积并不一定成正比。
当地面出现火点时,含有火点像元的混合像元辐射率由下式表述:
其中:i= 3,4,5,表示AVHRR红外通道序号,Nim为混合像元辐射率,P为亚像元火点(即明火区)面积占像元面积百分比,Nih火点辐射率,Nib为火点周围背景辐射率。背景辐射率可由混合像元周围的非火点像元(背景像元)获得近似值。
AVHRR有3个红外通道,尤其是有中红外和远红外两个波段的通道。根据维恩定律,对于相同温度增量(在相同面积上),在中红外和远红外引起的辐射率增量不同。温度增量越大,辐射率差异越大,由此计算的亮温值也将明显不同。根据这一特点,建立火点面积和温度的二元非线性方程组,利用牛顿迭代法,可以估算出着火点的实际面积和温度。
2.4 火点的精确定位实际工作中,由于上级产品只精确到市县,防火部门在组织扑救时往往还需进行人工地图查找,对准确进行火点定位和及时扑救带来了困难。利用3S技术精细化火情监测产品具有非常大的实用效果。
利用ArcView软件叠加前面的火情监测结果,结合1:5×104鹤壁市大比例尺地图及1:25×104鹤壁市土地利用图,可以很快确定火点所在乡镇、村名称及土地利用类型(图略)。
3 应用效果检验2004、2005年鹤壁市利用精细判识指标监测到火点61处(其中上级监测到火点52处),纠正上级产品漏判火点11处,错判2处(见表 2)。错判的2处火点都位于东部的沙土区,时间均为10月上旬13—15时。由于晴天午后沙土地表升温迅速,造成在通道3上沙土区亮温比其它农田区高6K左右,从而对火点判识产生干扰造成误判。此外,火点定位精度准确,利用1:5×104大比例尺地图可以具体判别到村庄,误差一般在200m以内,为防火部门及时准确发现火点提供了快捷直观的依据,取得了良好的服务效果。
(1) 根据不同下垫面确定的判识指标与季节及时次有关。经过精细化的判识指标不仅可以提高判识火点的准确度,减少误判和漏判,还可以提高火点判识精度,甚至能识别出0.067hm2左右的微弱火点信息。
(2) 对火点面积和温度的估算反映的是卫星观测的瞬时状态,其后有可能因当地风力、地形、可燃物密度等条件发生变化,因此应用时需考虑这些情况。
致谢:本文承蒙河南省气象局李平高工指导,特此致谢!
[1] |
张树誉, 李登科, 李星敏, 等. 卫星遥感在秸秆焚烧监测中的应用[J]. 气象, 2005, 30(9): 83-86. DOI:10.7519/j.issn.1000-0526.2005.09.018 |
[2] |
赵红旭. 云南省卫星监测火点三级定位系统[J]. 气象, 1999(7): 10-12. |
[3] |
张春桂. 基于RS与GIS技术的福建省森林火灾监测研究[J]. 福建林学院学报, 2004(1): 32-35. |
[4] |
李近文, 钟儒祥. 利用气象卫星监测森林火情的方法探讨[J]. 广州林业科技, 2005(3): 19-22. |
[5] |
邹春辉, 赵学斌, 刘忠阳, 等. 卫星遥感技术在秸杆焚烧监测业务中的应用[J]. 河南气象, 2005(3): 24-26. |
[6] |
刘健, 许健民, 方宗义. 利用NOAA卫星AVHRR资料分析云的性质[G]. 气象卫星遥感反演和应用论文集(上册): 397-403.
|