Landsat系列辐射定标参数整理.docxVIP

Landsat系列辐射定标参数整理 辐射定标参数整理 1.亮温计算 亮度温度是一个常用的温度概念，是在卫星高度上传感器探测波段范围内普朗克黑体辐射函数与传感器响应函数乘积积分得到的辐射值．亮度温度包含有大气和地表对热辐射传导的影响，不是真正意义上的地表温度。 式中，lλ为传感器探头孔处的光谱辐射强度，即车载辐射值，以实现像素DN值到绝对辐射值的转换。 1.1.星上辐射亮度（lλ） 遥感图像的亮度值（DN值）是以8bit编码的量化和校正数字图像。为了准确地检索地物的特征，有必要将DN值转换为机载辐射值。 1.1.1.landsat8 lλ=ml*qcal+al 通过查看影像的头文件，可以获取偏差参数：ml(radiance_mult_band_x)和al(radiance_add_band_x)为图像的增益和偏置。 1.1.2. 陆地卫星5/7 qcal为经过辐射校正的图像灰度值即dn值；lmax为探测器可检测到的最大辐射亮度，也是最大灰度值所相应的辐射亮；lmin为探测器可检测到的最小辐射亮度，也是最小灰度值所相应的辐射亮度。 表1 landsat5tm的lmin和Lmax值 表2landsat7etm+的lmin和lmax值 Qcalmax是传感器接收到的最大灰度值，qcalmin是传感器接收到的最小灰度值。（1） 作为 果没有元数据信息，qcalmin默认值1（tm和etm+1)或者0(mss)；qcalmax取默认值255(tm和etm+)或者127(mss)。（2）如果有元数据信息，qcalmin取值如下：对于lpgsproducts（thelevel1productgenerationsystem）取值为1，对于nlapsproducts（nationallandsatarchiveproductionsystem）在04april2021之前取值为0，在04april2021之后取值为1；qcalmax取值为127（mss）,255（tm、etm）。 注：LPG和NLAP分别是从两个数据处理系统获得的产品。从2022年12月开始，LPtM+和L5由LPG系统处理，L4TM和MSS由NLAPS系统处理。 表3landsat5/7的qcalmin和qcalmax的值 1.2. 预设常数k k1和k2是发射前预设的常量，具体值如下表所示。 2.大气顶层反射率（视发射率） 6.1.landsat5/7（tm/etm） π?lλ?d2 ρ=其中：ρ——地面相对反射率；d——日地天文单位距离；lλ——传感器光谱辐射值，即大气顶层的辐射能量；esun——大气顶层的太阳平均光谱辐射，即大气顶层太阳辐照度；1 注：Landsat7热红外波段（波段6）在格式1中始终设置为低增益（6L），在格式2中始终设置为高增益（6h） θ——太阳天顶角（单位为弧度）。 相关参数可以通过不同的方式获得。式中：日地天文单位距离d=1-0.01674cos（0.9856×（jd-4）×π/180）；jd为遥感成像的儒略日，可通过儒略日算法模块获取；太阳天顶角=90°-太阳高度角；太阳高度角可以从遥感数据的头文件中获取；大气顶部的太阳辐照度（eSUN）来自遥感部门定期测量和发布的信息。具体计算方法：① 地球天文单位距离D： 儒略日计算：jd=1721103.5+int(365.25*y)+int(30.6*m+0.5)+d 计算太阳和地球之间的距离：D=1-0.01674*cos（0.9856*（jd-2451545）*π/180/36525），其中YMD分别为月、年和日 ②太阳天顶角＝90°－太阳高度角（单位：弧度）太阳高度角在头文件中（sunelevation） 在这个头文件中，它是：太阳天顶角=90°-太阳高度 ③大气顶层太阳辐照度（esun）可从遥感权威单位定期测定并公布的信息中获取； 大气顶部的太阳辐照度（eSUN） 6.2.landsat8（oli） ρλ=（mp*qcal+ap）/sin（θse） 为增益参数（reflectance_mult_band_x）, （反射率加上波段），θSe是太阳高度