云地分离中太阳天顶角对反射率阈
值法阈值的影响
摘 要
作为云检测的一个重要组成部分，本文利用FY-2C云图资料，分析了云图反射率与太阳天顶角及地表性质、云性质之间的关系，得到不同性质的地表、云区对应的云图反射率值与太阳天顶角余弦之间的拟合关系式,并给出点聚图。参考这些关系式和点聚图，可确定不同太阳天顶角下的反射率阈值，从而为反射率阈值法云检测提供重要依据。
关键词： 云检测 FY-2C云图 太阳天顶角 反射率
影响数值天气预报准确性的因素很多，其中一个最重要的因素就是人们对云特性尚缺乏足够的了解。因此近年来云特性的观测研究日益受到人们的广泛关注。
当大气中有云层存在时，云在可见光波段具有较大的反射率。另一方面，云顶温度通常低于地表温度，有云存在时向外层空间的热辐射小于无云时的热辐射。基于上述事实，人们发展了多种利用卫星资料检测云的方法，其中具有代表性的有阈值方法和空间相干方法和卷云监测方法。
空间相干法主要利用发射场的空间结构来确定像元是否对应晴空或有云并推导探测器视场部分被云覆盖所对应的云量。利用空间相干法估计云量存在的问题在于其需要假定云层水平分布是均匀的。对于卷云检测，一种常用方法是利用11微米和12微米的亮温差以及8微米和11微米的亮温差来检测卷云。
阈值方法的基本思想是将一个区域上空的可见光或热红外大气窗口卫星图像按照某一辐射率阈值分为2个部分。辐射率小于该阈值的区域为晴空无云区域，辐射率大于该阈值的区域为被云覆盖的区域。如何选择辐射率阈值是利用阈值方法确定云量所要解决的关键问题。阈值方法一般分为亮度温度阈值法、反射率阈值法和双通道差值法。简单的亮度温度阈值法主要用于检测云顶温度较低的云，最适用于夜间海洋上的云检测。
与下垫面相比，云具有高的反射率和低的温度。因此，简单的可见光和红外窗区的阈值就可提供相当不错的云检测方法。当下垫面为冰雪、云为薄卷云、夜间出现低的层云或小的积云时，云和下垫面的辐射相似，难以用简单的可见和红外窗区阈值区分云和下垫面时，可采用多通道或吸收带光谱法
1
（比如用水汽吸收带1.38?m通道检测高空卷云）进行云检测。
目前存在的问题主要是利用可见光云图进行云检测时其反射率阈值的设定受太阳高度角较大影响，尤其是检测低云。而对于不同的地表情况反射率阈值也并不相同。
本文的主要工作是利用FY-2C云图资料分析不同性质的地表、云区对应的云图反射率与太阳天顶角余弦之间的关系，为反射率阈值法云检测提供依据。
1 可见光云图的遥感原理及特点
1.1可见光云图观测原理
卫星在可见光大气窗区观测，大气中的散射辐射很小时，忽略大气散射作用，卫星接收到的辐射为
??s??sL?(0;?,?)???(?,???sun,?sun)E?(?)?sunexp{?}exp{?}?sun?sat 是双向反射率。通常大气在水平方向近似为均匀的，地面为朗伯面，则由式，在有阳光时不考虑大气散射的辐射传输方程写为
sat L?(0;?,?)???s?L(?)?E(?)?sunexp{??s}exp{?} ??sun?~?L(?)E(?)?sunT?(??s,?sun)T?(??s,?) ?~ =
式中?L(?) 是朗伯面反照率。T?(??s,?sun)、T?(??s,?)分别是太阳光入射方向和卫星观测方向的大气透过率。由于在大气窗区
~~T?(??s,?sun)?T?(??s,?)?1所以上式又可以写为
Lsat?(0;?,?)?~?L(?)E?(?)?sun ? 从上式可以看到，由于E?(?)是入射大气顶的辐射，通常可以当作定值，因此卫星观测到的辐射可以近似地认为与地面的反照率?L(?)和太阳的天顶角余弦?sun成正比。对此下面作进一步的分析讨论。
1.2 基本特点
可见光云图是卫星扫描辐射仪在可见光谱段，如 AvHRR 仪器的 cH1
2
（0.68—0.725微米）通道或静止卫星的（ 0 . 52 — 0 . 75 微米 ) ，测量来自地面和云面反射的太阳辐射，如果将卫星接收到的地面目标物反射的太阳辐射转换为图像，卫星接收到的辐射越大就用越白的色调表示；而接收到的辐射越小，则用越暗的色调表示，这就得到可见光云图。在可见光云图上，物像的色调决定于反射太阳辐射的强度。而卫星接收到的反射太阳辐射决定于入射到目标物上的太阳辐射，及目标物的反照率。入射至目标物的太阳辐射又与太阳高度角有关。因此，在可见光云图上物像的色调与其本身的反照率和太阳高度角有关。下面对此进一步说明。
( l ）反照率对可见光云图上色调的影响从（ 5 . 2 ）式可以看到，在一定的太阳高度角下，卫星接收到的辐射仅决定于物体的双向反射率，如果将地面看成朗伯面，则卫星接收的辐射仅取决于物体的反照率，物体的反照率越大，它的色调越白；反照率越小，色调越暗。表 1 给出了各种云和地面目标物体的反照率，从表中可以看出：
① 水面的反照率最低，厚的积雨云最大；
② 积雪与云的反照率十分接近，所以仅从可见光云图上的色调难以区
别云和积雪；
③ 薄卷云与晴天积云、沙地的反照率也很接近，也不易区别它们。
表1 可见光云图上主要目标物的色调
色调 1.黑色 2.深灰色 3.灰色 4.灰白色 5.白色 6.浓白色 目标物 海洋、湖泊、大的河流 陆地上大面积森林覆盖区、牧场、草地、耕地 陆地上晴天积云、塔里木沙漠、陆地上单独出现的卷云 大陆上的中高云 积雪、冰冻的湖泊和海洋、中等厚度的云（中云、积云和层积云） 大块厚云、积雨云团 ( 2 ）太阳高度角对可见光云图上色调的影响 太阳高度角决定了卫星观测地面时的照明条件，太阳高度角越大，光照条件越好，卫星接收到的反射太阳辐射也越大，否则越小。这就是目标物的色调还与每天卫星观测的时刻和季节有关，如在北半球冬季中高纬度地区，太阳高度角很低，照明差，图片色调十分灰暗。又如卫星在早晨或傍晚观测，太阳高度角也很低，图片色调也很暗。对于同一图片上的各个点的太阳高度角也不同，如是上午的云图，图片右半侧（东面一侧）的太阳高度角较高，色调明亮，而左半侧，太阳高度角低，色调较暗。反之也可以根据这一特点判断云图的观测时刻，是否是可见光云图。对于静止卫星中午的云图，整个观测区的光照条件较好，物像间的反差明显，图片明亮。
3
表2 一些主要云和地面目标的反照率
云和地面目标物 1.积雨云 主要特征 大而厚 反照率 92 云和地面目标物 10.层云 11.卷云 2.积云 小，云顶在6km左右 86 12.卷层云 13.晴天积云 4.积云,层积云 陆地上，云量大于80% 69 14.中云（高层，高积云） 14.沙地 主要特征 薄，洋面上 薄，单独出在陆地上 3.卷层云 厚，下面由衷地云和降水 74 单独在陆地上 陆地上，云量大于80% 5.层积云 6.层云 7.沙漠 8.层积云 陆地上，云量大于80% 厚，出现在洋面上，云厚约0.5km 白砂 洋面成片 旧雪，已有3—7天 9.积雪 ，大部分在森林地区 新雪 68 64 60 60 59 80 中等厚度 谷地、平原、坡地 15.沙地和矮树林 16.植被 18.海洋，湖泊， 河流 9（7） 针叶林 68 77 17 12 29 36 32 反照率 42 2 FY-2C卫星及云图资料特点
2.1我国静止气象卫星风云二号（FY-2C）
风云二号（FY—2C）气象卫星是我国自行研制，于2003年10月19日上午9时18分在西昌卫星发射中心成功发射。
这次发射的风云二号（FY—2C）静止业务气象卫星位于3．6万公里高空的地球静止轨道，定点于105°E，可以24小时对我国天气系统及演变过程进行观测。FY-2C卫星包括下面子系统：遥感子系统，数据传输和广播子系统，数据收集系统，遥测和指令子系统，天线子系统，姿态和轨道控制子系统，能源子系统，热控和远地点发动机子系统。
FY-2C卫星可以获取白天可见光云图、昼夜红外云图和水汽分布图；向卫星覆盖区转发宽云图和低速区域图；收集覆盖区内地面数据收集平台采集的环境监测数据；监测太阳活动和卫星所处轨道的空间环境。我国天气系统大多从西边过来，此次发射的风云二号（FY—2C）气象卫星监测着从西边进
4
入我国的大多数天气系统，这大大增强了实时监测我国上游地区天气系统发生发展的能力，从而提高天气、数值预报和气候预测的准确性。同时还能转发分布在地面、高山、海洋、木、沙漠等各种数据收集平台观测的地球环境资料，为认识和开发利用可再生的风能、太阳能，防治沙漠化扩大发挥重要作用。
静止轨道气象卫星资料已经成为天气分析、预报服务必不可缺的重要监测手段，在灾害性天气的监测，特别是对台风、暴雨、洪水等重大自然灾害的监测中发挥了重要的作用。它的成功发射，标志着我国静止气象卫星由实验阶段正式向业务阶段成功迈进。
与FY-2A和FY-2B两颗卫星相比，风云二号02批的三颗卫星有下列改进：
1. 可见光和红外自旋扫描辐射计的通道由三个增加到五个，较大的增强了卫星的观测能力；卫星数据格式也有相应改变。
2. 增加星上电源能力，在卫星地影星蚀期间仪器不间断工作，加强了卫星安全保护措施；
3. 将FY-2A/B卫星的WEFAX云图模拟广播改为数字广播。
表3 FY-2C/D/E VISSR主要特征
波长 可见 0.5-0.75μm 红外 3.5-4.0μm 5km 140μrad 2500 HgCdTe 水汽 红外 红外 11.5-12.5μm 5km 140μrad 2500 HgCdTe 6.2-7.6μm 10.3-11.3μm 5km 140μrad 2500 HgCdTe 5km 140μrad 2500 HgCdTe 分辨率 1.25km FOV 35μrad 扫描线 2500×4 探测器 Si-phsto-diode S/N=0.5(反照率噪声 2.5%) S/N=43(反照率95%) 8 bits NEDT=0.5(300K) NEDT =1K(300k) NEDT=0.5(300K) NEDT =0.5(300K) 10 bits 8 bits 10 bits 10 bits 2.2 云图资料特点 本文主要应用可见光（0.5-0.75μm）和IR2（10.3-11.3μm）两个通道的云图资料及太阳高度角资料，云图资料显示的是灰度值，必须通过灰度与反照率对照表转换为反照率。每幅云图覆盖范围为60°N—60°S，45°E—165°
5