Landsat5TM数据辐射定标
- 格式:pdf
- 大小:170.67 KB
- 文档页数:5
应用6S模型进行LANDSAT_TM影像大气校正应用6S 模型进行LANDSAT TM 影像大气校正一、辐射校正1、 1、用定标系数将原始DN 值转换为大气层顶太阳辐亮度L ;rescale cal rescale B Q G L +?=λL 为大气层顶太阳辐亮度,Q 为记录的电信号数值,rescaleG ,为通道增益,rescale B 为偏移量,定标系数可以在头文件中获得。
表1 LANDSAT5 TM 数据定标系数2、由大气层顶太阳辐亮度L 转换为反射率。
s p ESUN d L θπρλλcos 2=其中:pρ: 行星反射率λL : 传感器口径的光谱辐射值 d: 日地距离(以天文为单位)λESUN :Mean solar exoatmospheric irradiances 平均太阳外大气层辐射值s θ : 太阳天顶角表2 TM 太阳外大气层光谱辐射值表3,日地距离(以天文为单位)二、大气校正经过辐射校正后,象元灰度值转换为了反射率,我们使用6S模型对可见光和近红外波段进行大气校正。
1、辐射校正完的反射率是0-1之间的值,然后把它转换为0-100之间的数值;2、再在ENVI中把第一步中的反射率(0-100)存为RAW格式;3、然后在inputfiles中填写大气条件的输入文件、大气条件的输出文件名、待大气校正的输入文件(RAW格式)、待大气校正的图像大小;4、4、第3步中的大气条件的输入文件需要填写以下的几项:Landsat5 geometrical conditionsmonth,day,hh.ddd,long.,lat.tropical atmospheric modecontinental aerosols modelvisibility in km (aerosol model concentration)target at 600 m above sea levelsensor on board of satellitethird band of Landsat5the image has values of reflectance, DN is percent (actual values only 0-100, not 0-255)(-1)number of pixels of the image=number of bytes以上的这些参数可以根据实际情况进行填写。
辐射定标参数整理1.亮度温度计算亮度温度是一个常用的温度概念,是在卫星高度上传感器探测波段范围内普朗克黑体辐射函数与传感器响应函数乘积积分得到的辐射值.亮度温度包含有大气和地表对热辐射传导的影响,不是真正意义上的地表温度。
计算公式:其中,Lλ为传感器探孔处光谱辐射强度,即星上辐射亮度值,实现像素DN值转化为绝对辐射亮度值。
1.1.星上辐射亮度(Lλ)遥感影像的亮度值(DN值)都是经过量化和纠正过的以8bit编码的数字影像,为了精确反演地物特性,有必要将DN值转化为星上辐射亮度值。
ndsat8Lλ= M L*Q cal + A L通过查看影像的头文件,可以获取偏差参数:M L(RADIANCE_MULT_BAND_x)和A L(RADIANCE_ADD_BAND_x)为图像的增益和偏置。
1.1.ndsat5/7QCAL为经过辐射校正的图像灰度值即DN值;L max为探测器可检测到的最大辐射亮度,也是最大灰度值所相应的辐射亮;L min为探测器可检测到的最小辐射亮度,也是最小灰度值所相应的辐射亮度。
表 1 Landsat5 TM的Lmin和Lmax值表 2 Landsat7 ETM+的Lmin和Lmax值QCAL max为传感器接收到的最大灰度值,QCAL min为传感器接收到的最小灰度值。
(1)如果没有元数据信息,QCAL MIN默认值1(TM和ETM+1)或者0(MSS);QCAL MAX取默认值255(TM 和ETM+)或者127(MSS)。
(2)如果有元数据信息,QCAL MIN取值如下:对于LPGS Products(The level 1 product generation system)取值为1,对于NLAPS Products(National Landsat Archive Production System)在04 April 2004之前取值为0,在04 April 2004之后取值为1;QCAL MAX 取值为127(MSS), 255(TM、ETM)。
ENVI-专题五LandsatTM辐射定标与大气纠正专题五Landsat TM辐射定标与大气纠正图像处理流程分为以下几个步骤:图像的配准、重采样、定标以及大气纠正。
其中,图像配准是做所有工作的前提,是图像的几何纠正。
在进行定性或者定量遥感时都要进行图像配准来确定我们所要研究的目标。
在进行完图像配准(几何纠正)之后,为了使得输出图像的配置与输入图像向对应,因此要进行重采样。
定标以及大气纠正则是进行辐射量纠正的重要的过程,是进行定量遥感不可少的步骤。
分类是图像处理的最后一步,按某种使用意图分类之后的图像,对实际应用有很大帮助。
辐射定标和大气纠正都属于图像的预处理,辐射定标的目的是把图像上的DN值转为辐亮度或者是反射率,大气纠正的目的是消除或减少大气对图像的干扰。
以下是图像预处理的流程:一、辐射定标要对图象进行辐射定标,将图象的DN值转化为表观反射率,该过程的实现是通过应用以下两个步骤来实现的:(1)首先将图象的DN值转化为辐亮度:radiance=gain*DN+offset (1)式(2)然后将图象的辐亮度转化为表观反射率:(reflectance) ρ=π*L*d2/(ESUN*cos(θ))(2)式其中ρ为表观反射率,L为表观辐亮度,d为日地距离,ESUN为太阳平均辐射强度,θ为太阳天顶角。
(3)将以上两个步骤结合得:ρ=π*(gain*DN+offset)* d2/(ESUN*cos(θ))(3)式①日地天文单位距离D:D=1 - 0.01674 cos(0.9856× (JD-4)×π/180);JD为遥感成像的儒略日(Julian Day)D = 1 + 0.0167 * Sin(2 * PI * (days - 93.5) / 365);days是拍摄卫片的日期在那一年的天数,如2004年5月21号,则days=31+29+31+30+21=142。
计算得:D=1.01250756ENVI中的具体实现(以Landsat 7 ETM+为例):采用简单的波段运算例如,我们把2002-5-22的一幅ETM图像第3波段的DN值转化为表观反射率。
辐射定标、辐射校正、大气校正、正射校正等相关概念作为初学者,容易将这几个概念搞混。
为了较好地理解这几个概念,先介绍一下相关的术语 terminology。
DN值(Digital Number ):遥感影像像元亮度值,记录地物的灰度值。
无单位,是一个整数值,值大小与传感器的辐射分辨率、地物发射率、大气透过率和散射率等相关。
反映地物的辐射率radiance地表反射率:地面反射辐射量与入射辐射量之比,表征地面对太阳辐射的吸收和反射能力。
反射率越大,地面吸收太阳辐射越少;反射率越小,地面吸收太阳辐射越多,表示:surface albedo表观反射率:表观反射率就是指大气层顶的反射率,辐射定标的结果之一,大气层顶表观反射率,简称表观反射率,又称视反射率。
英文表示为:apparent reflectance4、行星反射率:从文献“一种实用大气校正方法及其在TM影像中的应用”中看到“卫星所观测的行星反射率(未经大气校正的反射率)”;在“基于地面耦合的TM影像的大气校正-以珠江口为例”一文有“该文应用1998年的LANDSAT5 TM影像,对原始数据进行定标、辐射校正,求得地物的行星反射率”。
因此行星反射率就是表观反射率。
英文表示:planetary albedo,辐射校正VS. 辐射定标辐射校正:Radiometric correction 一切与辐射相关的误差的校正。
目的:消除干扰,得到真实反射率的数据。
干扰主要有:传感器本身、大气、太阳高度角、地形等。
包括:辐射定标,大气纠正,地形对辐射的影响辐射定标:Radiometric calibration 将记录的原始DN值转换为大气外层表面反射率(或称为辐射亮度值)。
目的:消除传感器本身的误差,确定传感器入口处的准确辐射值方法:实验室定标、机上/星上定标、场地定标不同的传感器,其辐射定标公式不同。
L=gain*DN+Bias在ENVI4.8中,定标模块:Basic Tools>Preprocessing>Calibration Utilities>模块(只能对一个波段进行辐射定标)大气校正:Atmospheric correction 将辐射亮度或者表面反射率转换为地表实际反射率目的:消除大气散射、吸收、反射引起的误差。
遥感数据辐射定标和大气校正-Landsat 5 Thematic Mapper为例北京农业信息技术研究中心 遥感技术部yanggj@2008-12-261 辐射定标1.1 基本原理当TM 获取影像以后(Level 0),会将其转化为32 位浮点型的绝对辐亮度。
之后进一步处理,将绝对辐亮度变为8 位的DN 值(Digital Number),这也就是我们购买后拿到的数据(Level 1)。
如果要将L1 的DN 值转化为传感器处的辐亮度值(at-sensor spectral radiance),需要通过下面这个公式:其中,是传感器处的辐亮度值,注意单位是;为像元的DN 值;为传感器处最小辐亮度值所对应的DN 值(一般取0); 为传感器处最大辐亮度值所对应的DN 值(一般取255); 为光谱辐亮度的最小值,单位同 ; 为光谱辐亮度的最大值,单位同 。
上面的这个公式也可以改为:其中,各个波段的 , 以及 和 见表1。
需要注意的是,上述参数在2003 年5 月5 日前后是不一致的,所以在操作时,一定要搞清楚影像获取的时间。
参考文献:Chander, G. and Markham, B. (2003). Revised Landsat-5 TM radiometric calibration procedures and postcalibration dynamic ranges. IEEE Transactions on geosciences and remote sensing. Vol. 41, No. 11, November: 2674-2677.1.2 操作流程由于ENVI 4.4 中有专门进行辐射定标的模块,因此实际的操作十分简单。
将原始TM影像打开以后,选择Basic Tools¬Preprocessing¬Calibration Utilities¬Landsat TM进入下一步参数选择:根据传感器类型选择Landsat 4,5 或者7。
专题五Landsat TM辐射定标与大气纠正图像处理流程分为以下几个步骤:图像的配准、重采样、定标以及大气纠正。
其中,图像配准是做所有工作的前提,是图像的几何纠正。
在进行定性或者定量遥感时都要进行图像配准来确定我们所要研究的目标。
在进行完图像配准(几何纠正)之后,为了使得输出图像的配置与输入图像向对应,因此要进行重采样。
定标以及大气纠正则是进行辐射量纠正的重要的过程,是进行定量遥感不可少的步骤。
分类是图像处理的最后一步,按某种使用意图分类之后的图像,对实际应用有很大帮助。
辐射定标和大气纠正都属于图像的预处理,辐射定标的目的是把图像上的DN值转为辐亮度或者是反射率,大气纠正的目的是消除或减少大气对图像的干扰。
以下是图像预处理的流程:一、辐射定标要对图象进行辐射定标,将图象的DN值转化为表观反射率,该过程的实现是通过应用以下两个步骤来实现的:(1)首先将图象的DN值转化为辐亮度:radiance=gain*DN+offset (1)式(2)然后将图象的辐亮度转化为表观反射率:(reflectance) ρ=π*L*d2/(ESUN*cos(θ))(2)式其中ρ为表观反射率,L为表观辐亮度,d为日地距离,ESUN为太阳平均辐射强度,θ为太阳天顶角。
(3)将以上两个步骤结合得:ρ=π*(gain*DN+offset)* d2/(ESUN*cos(θ))(3)式①日地天文单位距离D:D=1 - 0.01674 cos(0.9856× (JD-4)×π/180);JD为遥感成像的儒略日(Julian Day)D = 1 + 0.0167 * Sin(2 * PI * (days - 93.5) / 365);days是拍摄卫片的日期在那一年的天数,如2004年5月21号,则days=31+29+31+30+21=142。
计算得:D=1.01250756ENVI中的具体实现(以Landsat 7 ETM+为例):采用简单的波段运算例如,我们把2002-5-22的一幅ETM图像第3波段的DN值转化为表观反射率。
landsat5波段介绍Landsat一5的成像传感器TM 获取的数据属于光学类遥感数据,目前中国科学院中国遥感卫星地面站(以下简称中国遥感卫星地面站)所生产的Landsat一5数据产品一共有四个级别,分别是0级、1级、系统级纠正(Systematic Geocorrection)与精纠正(Precision Geocorrection)。
0级产品是指像素值没有经过处理的图像数据,1级产品是指对0级产品进行辐射纠正后的产品,系统级纠正产品是在1级产品的基础上进行系统几何纠正后的产品,精纠正产品是引入了控制点信息进行几何精纠正后的产品。
相对于0级和1级产品,系统级纠正与精纠正的产品可以增加高程纠正(elevation correction)功能,高程纠正有三种选择级别:base、Coarse DEM与Fine DEM,base是全球尺度的高程纠正,CoarseDEM 是1:1O0万比例尺的DEM 数据,Fine DEM是1:25万以上比例尺的DEM 数据。
0级与1级产品的数据格式有EOsAT FAST 与CCRS LG—SOwG两种,系统级纠正与精纠正产品的数据格式有EOSAT FAST、CCRS LGSOWG与Geo Tiff三种。
0级产品没有经过辐射纠正和几何纠正处理,这类产品主要是面向具有一定遥感卫星数据处理经验的高级用户,所以一般情况下不建议普通用户使用。
目前中国遥感卫星地面站给用户提供最多的是经过辐射纠正和几何纠正的系统级纠正产品。
用户如果对图像定位精度要求比较高,平原地区可选用精纠正产品,高程较高的区域建议选用经过高程纠正的精纠正产品,其几何定位精度误差在1-2个像元之内。
因此从网上下载的数据应该是辐射纠正和几何纠正,还需要进行正射校正-精校正-辐射定标—layer stacking-大气校正-裁剪辐射定标会将分辨率不一致的波段丢失。
Landsat 5 Thematic Mapper 辐射定标和大气校正杨曦北京师范大学资源学院,旱灾与风险小组yangxi@2008-10-261 辐射定标1.1 基本原理当TM 获取影像以后(Level 0),会将其转化为32位浮点型的绝对辐亮度。
之后进一步处理,将绝对辐亮度变为8位的DN 值(Digital Number ),这也就是我们购买后拿到的数据(Level 1)。
如果要将L1的DN 值转化为传感器处的辐亮度值(at-sensor spectral radiance ),需要通过下面这个公式:其中,是传感器处的辐亮度值,注意单位是2/W m sr m μ••;为像元的DN 值;为传感器处最小辐亮度值所对应的DN 值(一般取0);为传感器处最大辐亮度值所对应的DN 值(一般取255);为光谱辐亮度的最小值,单位同;为光谱辐亮度的最大值,单位同。
上面的这个公式也可以改为:其中,各个波段的,以及和见表1。
需要注意的是,上述参数在2003年5月5日前后是不一致的,所以在操作时,一定要搞清楚影像获取的时间。
参考文献:Chander, G. and Markham, B. (2003). Revised Landsat-5 TM radiometric calibration procedures and postcalibration dynamic ranges. IEEE Transactions on geosciences and remote sensing. Vol. 41, No. 11, November: 2674-2677.1.2操作流程由于ENVI 4.4中有专门进行辐射定标的模块,因此实际的操作十分简单。
将原始TM影像打开以后,选择Basic ToolsÆPreprocessingÆCalibration UtilitiesÆLandsat TM进入下一步参数选择:根据传感器类型选择Landsat 4,5 或者7。
Envi 对Landsat进行辐射定标将DN值定标为辐亮度(Radiance)或表观反射率(Reflectance)1、定标为辐亮度(Radiance)1)由DN定标为辐亮度(Radiance)公式:•QCAL为原始量化的DN值•LMINλ为QCAL = 0时的辐射亮度值•LMAXλ为QCAL = QCALMAX时的辐射亮度值,2)各参数的取值:A、最小和最大量化定标参数QCALMIN和QCALMAX的取值:如果没有元数据信息:QCALMIN为最小量化定标像素值取默认值1(TM和ETM+))或者0 (MSS);QCALMAX为最大量化定标像素值取默认值255 (TM 和ETM+)或者127 (MSS)。
如果有元数据信息:QCALMIN是最小量化定标像素值取值如下:1:LPGS产品1:04 April 2004之后NLAPS产品0:04 April 2004之前NLAPS产品注:LPGS 和NLAPS分别是两种数据处理系统得到的产品:the Level 1 Product Generation System (LPGS) 和the National Land Archive Production System (NLAPS),从2008年12月份开始,L7 ETM+ 和L5都是以LPGS系统处理,L4 TM和MSS以NLAPS系统处理。
QCALMAX根据元数据信息取值为127(MSS), 254(ETM), 255(TM)。
B、最小和最大量化定标参数QCALMIN和QCALMAX的取值:LMINλ和LMAXλ的值取自Chander, Markham, and Helder (2009)的研究成果。
注意:按照如上公式计算出来的作为结果的辐亮度(Radiance)单位为:u W/(cm2*sr*um)。
2)定标系数LPGS处理产品的ETM+数据定标参数2、定标为表观大气反射率(ρ):3、Envy下实现图3 Fast格式的L7定标界面4、手动进行辐射定标当没有元数据时可以通过Band Math进行运算例如,我们把2002-5-22的一幅ETM图像第3波段的DN值转化为表观反射率。
Landsat 5 Themati c Mapper辐射定标和大气校正120102105龚鑫烨操作流程辐射定标:1、加载原始图像H1。
2、Basic Tools——Preprocessing——data-specifi c utiliti es——Landsat TM——landsat calibration进入下一步参数选择:根据传感器类型选择Landsat 4,5 或者7。
从遥感影像的头文件中获取Data Acquisi tion 的时间,Sun elevati on。
如果你是用File?Open ExternalFile?Landsat?Fast 的方法打开header.dat 的话,sun elevati on 就已经填好了。
这里Calibra tion Type 注意选择为Radiance。
输出文件,定标就完成了。
大气校正1、Spectra l——FLAASH/B asic Tools——Preprocessing——Calibra tion Utiliti es——FLAASH。
首先设定输入输出文件。
FLAASH模块要求输入辐亮度图像,输出反射率图像。
之前我们进行了辐射定标,得到辐亮度图像,在这里要把BSQ 格式的图像转换为BIL 或者BIP 格式的图像,然后再Input Radiance Image 中选择转换格式后的图像。
(Basic Tools——Convert Data(BSQ,BIL,BIP))。
这里注意,当输入图像后,程序会让你选择Scale Factor,即原始辐亮度单位与ENVI 默认辐亮度单位之间的比例。
E NVI 默认的辐亮度单位是2/ W cm sr nm μ••,而之前我们做辐射定标时单位是2/ W m sr m μ••,二者之间转换的比例是10,因此在下图中选择Singlescale factor,填写10.000。
Landsat 5和IRS P6数据辐射定标说明Landsat 5和IRS P6数据是两种应用很广泛的遥感数据,在利用两种数据来进行定量遥感研究过程中,不少用户会对它们的辐射定标问题产生疑问,现就中国科学院对地观测与数字地球科学中心(中国遥感卫星地面站)生产的Landsat 5 TM 和IRS P6数据的辐射定标问题做一说明(针对系统级纠正的FAST 格式产品)。
1. 辐射亮度计算辐射亮度计算采用如下的通用计算公式:max min min min max min()L L L QCAL QCAL L QCAL QCAL λ-=⋅-+- (1) 其中,L λ为光谱辐射亮度(w/(m 2μmsr)),QCAL 为像元的DN 值,QCAL max 和QCAL min 分别为像元可以取得的最大和最小DN 值(对于中国科学院对地观测与数字地球科学中心(中国遥感卫星地面站)生产的Landsat 5和IRS P6数据,QCAL min =0),L max 和L min 分别为QCAL= QCAL max 和QCAL= QCAL min 时的光谱辐射亮度值。
对于Landsat 5 TM 数据,QCAL max =255对于 IRS P6 LISS 3和LISS 4数据,QCAL max =255,对于AWIFS 数据,QCAL max =10231.1 TM 数据L max 和L min 获取TM 数据头文件中“RAD GAINS/BIASES ”之后的部分给出的是各波段的L max 和L min (注意这里给出的并非GAIN 和BIAS ),单位为mw/(cm 2sr),没有考虑波段宽度,需要除以波段宽度将单位变为mw/(cm 2μmsr),然后乘以10将辐射亮度单位变为标准单位(w/(m 2μmsr)),再参与有关的定量遥感计算。
从头文件计算TM 数据各波段的辐射亮度要考虑波段宽度并进行单位换算,比较复杂,另外中国科学院对地观测与数字地球科学中心(中国遥感卫星地面站)TM 数据处理系统采用的辐射定标算法与USGS TM 数据处理系统采用的辐射定标算法接近,而采用USGS 给出的定标参数来计算辐射亮度不用考虑波段宽度和单位换算,比较方便和精确,因此建议直接利用USGS 的定标参数来进行TM 数据的辐射定标[1],具体如下:L λ=gain*DN+bias (2)其中,L λ为光谱辐射亮度(w/(m 2μmsr),gain 和bias 为各波段的增益和偏置(w/(m 2μmsr)。
TM数据各波段中心波长值/ 波谱宽度即wavelength / FWHM Landsat 1-5Wavelength FWHMTM1:0.485μm/0.066TM2:0.569μm或0.56 / 0.082TM3:0.660 μm/0.067TM4:0.840μm或0.83 /0.128TM5:1.676μm或1.65 /0.217TM7:2.223μm或 2.22 /0.252TM6: 11.4μmLandsat 1-5表3 美国USGS发布的Lmax与Lmin值波段Lmax Lmin1 -1.52 193.002 -2.84 365.003 -1.17 264.004 -1.51 221.005 -0.37 30.206 1.2378 15.30327 -0.15 16.50表4 Landsat5TM 数据头文件中的Lmin与Lmax值波段Lmax Lmin1 1.26880 -0.01002 2.98126 -0.02323 1.76186 -0.00784 2.81771 -0.01935 0.65277 -0.00806 3.20107 0.259947 0.44375 -0.0040表5 TM 数据波段对应波谱宽度波段号频谱宽度1 0.0662 0.0823 0.0674 0.1285 0.2177 0.252Landsat 7OLI陆地成像仪包括9个波段,空间分辨率为30米,其中包括一个15米的全色波段,成像宽幅为185x185km。
OLI传感器高度为761km,轨道高度为705km。
OLI包括了ETM+传感器所有的波段,为了避免大气吸收特征,OLI对波段进行了重新调整,比较大的调整是OLI Band5 (NIR)(0.845–0.885 μm),排除了0.825μm处水汽吸收特征;OLI全色波段Band8波段范围较窄,这种方式可以在全色图像上更好区分植被和无植被特征;此外,还有两个新增的波段:蓝色波段(band 1; 0.433–0.453 μm) 主要应用海岸带观测,短波红外波段(band 9; 1.360–1.390 μm) (卷云波段)包括水汽强吸收特征可用于云检测;近红外band5和短波红外band9与MODIS对应的波段接近。
本次实验以烟台地区TM5数据为例,向大家介绍6S和FLAASH大气校正模型。
ENVI4.7进行Landsat TM 5辐射定标和大气校正一、辐射定标使用ENVI4.7对TM5进行辐射定标,输入数据格式为ENVI标准格式,一次只能定标一个波段。
(以20080227为例)1.打开原始数据(以Fast格式打开头文件),File---Open External File---Landsat---Fast,选择原始数据的头文件header.dat打开Gyanesh Chander等修正后的定标系数选择输出文件路径,输入输出文件名,点击OK进行辐射定标,得到定标结果c1(辐亮度)。
3.将定标后的7个辐亮度图像合并为一个图像得到s2*******(注意波段顺序要一一对应),合并后图像中没有wavelength信息,在头文件s2*******.hdr最后加上波长信息wavelength = {0.48500,0.56000,0.66000,0.83000,1.65000,11.40000,2.220000}wavelength units = micrometers4.格式转换,上述操作得到的图像是BSQ格式,转换为FLAASH要求的BIL或者BIP格式。
Basic Tools–Convert Data(BSQ,BIL,BIP),打开对话框:二、大气校正简单一点的大气校正可以采用ENVI的FLAASH模块,以下就是FLAASH操作的步骤:1.FLAASH 模块的进入方法是Spectral–FLAASH,或者是Basic Tools–Preprocessing–Calibration Utilities–FLAASH。
2.FLAASH 模块的操作界面分为三块:最上部设定输入输出文件;中间设定传感器的参数;下部设定大气参数。
3.首先设定输入输出文件。
FLAASH 模块要求输入辐亮度图像,输出反射率图像。
在InputRadiance Image 中选择转换格式后的图像。
在新版本ENVI中(4.5版本以及更新),没有单独设立ASTER和MODIS定标的工具。
对于ASTER L1A/L1B和MODIS 02级数据,在打开数据时会自动完成对数据的定标。
如图1所示打开ASTER L1B的结果,在波段列表中,自动读取各个波段的中心波长信息,并按照波段范围信息(VNIR、SWIR、TIR)分组波段。
其中VNIR、SWIR自动定标为辐射亮度,单位是:W/m2/sr/μm;TIR数据定标为大气表观温度值,单位:开尔文。
打开其中一个数据,浏览像元值,可以看到已经定标为浮点型的辐射亮度值。
图1 ASTER L1B数据如图2为打开MODIS 02级1km数据,其中250米和500米的波段经过重采样为1km 加入这个数据集中。
ENVI根据各个波段的中心波长信息定标为三个类型数据:反射率数据(Reflectance)、辐射亮度值数据(Radiance)和发射率数据(Emissive)。
其中反射率和发射率为0~1无单位值,辐射亮度值单位是:W/m2/μm/sr。
图2 MODIS 02级数据如果打开原始的ASTER和MODIS的DN值数据,可以在ENVI主菜单中选择File->preferences,切换到Miscellaneous面板,将Auto-Correct ASTER/MODIS项设置为NO。
3.2 Landsat数据定标ENVI4.7版本改进了Landsat数据定标的功能,对于Landsat4/5数据可以手动选择以下两种定标公式:式中:•QCAL为原始量化的DN值•LMINλ为QCAL = 0时的辐射亮度值•LMAXλ为QCAL = QCALMAX时的辐射亮度值注:LMINλ和LMAXλ的值取自Chander, Markham, and Helder (2009)的研究成果。
•QCALMIN是最小量化定标像素值(与LMINλ类似)。
取值如下:1:LPGS产品1:04 April 2004之后的NLAPS产品0:04 April 2004之前的NLAPS产品注:如果没有元数据信息,QCALMIN取默认值1(TM和ETM+))或者0 (MSS)。