landsat5大气校正

LandsatTM辐射定标和大气校正步骤2011-01-04来源：小木虫作者：zxf200315309LandsatTM辐射定标和大气校正不一、数据准备从usgs网站或者马里兰大学下载TM原始数据，usgs网站下载的数据是原始数据，在envi软件File–OpenExternalFile–Landsat–Geotiffwithmeta中只需打开***********_MTL.txt即可打开所有波段数据（除band6）u马里兰大学网站下载的数据有可能不是原始数据，在envi软件File–OpenExternalFile–Landsat–Geotiffwithmeta中只需打开***********.met即可打开所有波段数据（除band6）二、辐射定标1.由于ENVI4.7中有专门进行辐射定标的模块。将原始TM的影像打开以后，选择BasicTools–Preprocessing–CalibrationUtilities–LandsatCalibration2.进入下一步参数选择：根据传感器类型选择Landsat4,5或者7。从遥感影像的头文件中获取DataAcquisition的时间，Sunelevation。如果你是用File–OpenExternalFile–Landsat–Geotiffwithmeta(Fast)的方法打开的话，sunelevation就已经填好了。这里CalibrationType注意选择为Radiance。输出文件，定标就完成了。三、大气校正简单一点的大气校正可以采用ENVI的FLAASH模块，以下就是FLAASH操作的步骤：1.FLAASH模块的进入方法是Spectral–FLAASH，或者是BasicTools–Preprocessing–CalibrationUtilities–FLAASH。2.FLAASH模块的操作界面分为三块：最上部设定输入输出文件；中间设定传感器的参数；下部设定大气参数。3.首先设定输入输出文件。FLAASH模块要求输入辐亮度图像，输出反射率图像。之前我们进行了辐射定标，得到辐亮度图像，在这里要把BSQ格式的图像转换为BIL或者BIP格式的图像，然后再InputRadianceImage中选择转换格式后的图像。(BasicTools–ConvertData(BSQ,BIL,BIP))。这里注意，当输入图像后，程序会让你选择ScaleFactor，即原始辐亮度单位与ENVI默认辐亮度单位之间的比例。ENVI默认的辐亮度单位是μW/cm2?sr?nm，而之前我们做辐射定标时单位是W/m2?sr?μm，二者之间转换的比例是10，因此在下图中选择Singlescalefactor，填写10.000。4.此外，如果TM影像的头文件中没有波段的信息，在这里也要求你提供一个.txt文件以包含此信息。那么，准备好一个.txt文件，其中含有一列TM每个波段中心波长的信息。如果，我们打开的是Geotiffwithmeta，就不用填写波段信息了。5.在OutputReflectanceFile和OutputDirectoryforFLAASHfiles里面设定输出文件的文件名和位置。6.设定传感器参数。首先是SceneCenterLocation，即遥感图像中心的坐标，以及FlightDate,FlightTimeGMT，这三者都可以在TM的头文件中找到，填入即可。7.在SensorType菜单中选择LandsatTM5。此时Sensoraltitude自动填上为705km。而PixelSize填为30m。8.根据遥感影像研究区实际情况，填写GroundElevation，比如呼和浩特市为1.05km。9.最关键的为大气参数部分:a)AtmosphericModel(大气模式):共有Sub-ArcticWinter(SAW)，Mid-LatitudeWinter(MLW)，U.S.Standard(US)，Sub-ArcticSummer(SAS),Mid-LatitudeSummer(MLS)和Tropical(T)。根据经纬度和时间可以选定研究区的大气模式，见ENVIHelp。b)AerosolModel（气溶胶模式）：有Rural,Urban,Maritime和Tropospheric四种选择。根据实际情况选择即可。关于此四种模式的解释见ENVIHelp。c)当我们选择TM时，可选的参数还有AerosolRetrieval和InitialVisibility。这两个参数对最后的结果又相当重要的影像，因此最好能调查到当地的InitialVisibility。此外，AERONET在全世界各地有测定AOD（AtmosphericOpticalDepth）的站点，可以查询AOD以后转换为消光系数，通过消光系数估算能见度，此步骤比较繁琐，在此不予详述。如果采用AerosolRetrieval中的K-T算法计算Visibility，且能够计算出结果的话，则采用K-T算法的能见度，否则采用InitialVisibility所指定的能见度。d)关于AerosolRetrieval。如果选择了下拉菜单中的K-Tmethod，那么需要在MultispectralSettings中设定参数,在AssignDefaultValuesBasedonRetrievalConditions中选择Over-landRetrievalStandard(660:2100nm)即可。根据不同的研究区可以设定不同的模式。其他设定可以不改变。Apply即可注：根据网上资料整理。