基于气象卫星资料的地类面积计算方法.pdfVIP

第 2 期 ,总第 40 期 国 土 资 源 遥 感 No . 2 ,1999 1999 年 6 月 15 日 　REMOTE SENSIN G FOR LAND RESOURCES Jun . ,1999 基于气象卫星资料的 地类面积计算方法 纪 　平 　　易浩若 ( 中国林业科学研究院资源信息研究所 ,北京 　100091) 摘 　要 　提出了二种用 NOAA 卫星数据不通过几何校正 ,就能快速 、准确计算各地类面积的 思想 , 同时给出了计算方法并通过实例进行了验证 。 关键词 　面积计算方法 　气象卫星资料 分类号 　TP79 : S29 0 　引言 气象卫星以其运行周期短 、覆盖范围广 ,资料被广泛应用于遥感的各个应用领域 。如主要 农作物估产 ,旱情 、洪涝灾害 、森林火灾 、雪灾等重大自然灾害的监测与评估 ,海冰监测 ,沙漠化 动态监测等 。如何快速 、准确地测算出各地类面积 ,是上述各项应用研究所必须解决的基本问 题之一 ,也是国内外同行所致力追求的目标 。以往的普遍做法是 ,利用NOAA 卫星星下点像元 ( 2) 面积 1. 1 ×1. 1 km 与受灾面所包含的像元个数相乘 ,得到受灾面积 。但是 ,越来越多的实践 表明 ,这种方法所计算出来的面积与实际受灾面积相比,有较大的出入 : 一般而言 ,在星下点 附近 ,该方法计算出来的受灾面积比实际受灾面积大 ; 而在远离星下点的地方 ,则比实际受灾 面积小 。为此 ,本文提出两种利用气象卫星已知信息 ,在不进行几何校正的情况下 ,快速 、准确 地计算各地类面积的方法 ,即: 行列号法和经纬度法 。 1 　行列号法 目前所使用的NOAA 气象卫星数据 ,其扫描带覆盖地表宽度约为 2 800 km ,星下点像元分 辨率为 1. 1 km ,扫描角度为 ±55. 4°,扫描采样的瞬时视场角为 0 . 0 013 rad ,而扫描时的步进为 0 . 000 944 rad 。这也就是说 ,在扫描方向上和步进方向上像元之间存在着重叠 ,这种重叠使得 每个像元所代表的面积比实际的面积要大 。同时 , 由于地球曲率的影响 ,像元大小随着扫描角 的增大而变大 ,在图像边缘部分的像元面积可以达到星下点像元面积的4～5 倍 。根据资料显 示 ,沿卫星步进方向的像元变形比扫描方向的变形要小的多 〔1〕。由于本文计算的是像元覆盖 的地面实际面积 ,因此 ,不考虑步进方向的覆盖 ,只以扫描时的步进角为分割来计算 。总的计 算思想是 ,计算每个被识别像元的面积 ,再将其求和 。而对于每个像元先计算其扫描线方向上 收稿 日期 : 1998 - 10 - 29 ,1998 - 12 - 30 改回。 ·56 · 　　　　国　土 　资 　源 　遥 　感 1999 年 在地球表面扫过的弧长 d , 再与步进弧长相乘 。下面推导弧长 d 的计算公式 。 如图 1 所示 , O 为地球中心, O1 为卫星扫描中心, A 点为第N 个像元的扫描起始点与地球 表面的交点, B 点为第N 个像元的扫描结束点与地球表面的交点, 也是第 N + 1 个像元的扫描 起始点与地球表面的交点, θ 为扫描辐射计在扫到第 N 个像元扫描起始点时已转过的扫描 N 角, Ψ 为 OO 与 OA 之间的夹角 。设某地类的像元共有 k 个, 每个像元的列号分别为 m , …, N 1 1 mk 。则 : 55. 4 　　　　　 θ = N