LANDSAT-5辐射标定.docVIP

中国遥感卫星地面站LANDSAT-5 TM的辐射定标 王文 李子扬 唐伶俐 2004-12-23 近年来LANDSAT-5卫星的内部定标灯逐渐老化，基于星上内部定标灯的辐射定标算法变得越来越不准确，由于LANDSAT-7卫星的SLC不能正常工作，提高LANDSAT-5卫星的辐射定标精度就显得更有意义。为此，NASA建议LANDSAT-5也使用类似于LANDSAT-7的传感器增益算法，使得这两颗卫星的辐射定标数据更具有可比性。 为了使广大用户更好地使用LANDSAT-5卫星数据，中国遥感卫星地面站于2004年5月对LANDSAT-5 TM数据处理系统进行了升级，更新后的处理系统在原有的CCRS 和NASA 基础上增加了NASA CPF和CCRS CPF两种算法，一般产品生产系统缺省选择CCRS CPF算法，即2004年5月以后生产的产品，如果没有特殊需求，都使用CCRS CPF算法。2004年5月以后选用CCRS CPF或NASA CPF处理的产品需要使用新的辐射定标算法。使用老的算法（CCRS和NASA）处理的产品不能使用新算法所给定的固定值来进行计算辐射亮度值。（固定值见下表1） 新的辐射定标算法只适用于1-5和7波段，6波段仍然沿用老的基于星上内部定标灯的辐射定标算法，对于该波段的新辐射定标算法数据正在研究中。下表给出了USGS提供的新算法所使用的LMAX和LMIN固定值。 Band LMIN LMAX 1 -1.52 193.00 2 -2.84 365.00 3 -1.17 264.00 4 -1.51 221.00 5 -0.37 30.20 6 1.2378 15.3032（没有改变） 7 -0.15 16.50 表1：USGS提供的新算法所使用的LMAX和LMIN固定值 需要注意的是，历史上LANDSAT-5 TM的定标信息以为辐射率单位。为了与LANDSAT-7 ETM+保持一致性，升级后产品的辐射率单位为。 使用升级后的处理系统生产的LANDSAT-5的FAST B格式产品中，头文件与老的处理系统生产的产品一样，仍然包含LMAX/LMIN的信息，用户可以直接使用这些信息计算GAIN和BIAS的值，但是依然应该注意的是FAST B头文件中的LMAX/LMIN值的单位不含有，所以在使用时还需考虑波段的波谱宽度，如下表所示： 波段号 频谱范围μ 频谱宽度 1 0.45 – 0.52 0.07 2 0.52 – 0.60 0.08 3 0.63 – 0.69 0.06 4 0.76 – 0.90 0.14 5 1.55 – 1.75 0.20 6 10.40 – 12.5 2.10 7 2.08 – 2.35 0.27 因为使用表1所提供的固定值将不用进行单位换算和波谱宽度相关计算，并且头文件中的LMAX/LMIN值由于处理过程中的系统误差而变得不精确，因此，我们建议用户直接使用USGS提供的LMAX/LMIN值来计算2004年5月以后生产的数据产品的GAIN和BIAS。 增益GAIN LMAX - LMIN /255 偏置BIAS LMIN 辐射亮度 L GAIN * DN + BIAS 我们选取了2003年7月18日的LANDSAT-5数据，PATH/ROW为141/30，对这4种算法进行了比较。下表为第4波段图像中所选取的几个点： 检查点 CCRS（老） NASA（老） CCRS CPF（新） NASA CPF（新） 点1灰度值（DN） 11 11 11 11 点1辐射亮度 L 7.967059 7.967059 8.088471 8.088471 点2灰度值（DN） 65 66 65 66 点1辐射亮度 L 48.02553 48.76735 55.20824 56.08082 点3灰度值（DN） 129 130 129 130 点1辐射亮度 L 95.50224 96.24406 111.0539 111.9265 可以看出，使用NASA/CCRS两种算法所获得的结果差别不大，同样NASA CPF/CCRS CPF两种算法之间差别也不大，但新算法与旧算法之间还是有不小的差别。 采用新的LANDSAT-5辐射定标算法，可以消除由于内部定标器老化所带来的误差，由于新算法与LANDSAT-7 ETM+所采用的辐射定标算法类似，使得两颗卫星的数据更具有一致性，给同时使用这两颗卫星数据的遥感应用带来更多的便利。