一、Landsat_5 TM 参数介绍

注：以下信息参考百度百科
Landsat-5卫星主要轨道特性参数如下：
近极近环形太阳同步轨道
倾角：98.2。
轨道周期：99分钟
穿越赤道时间：上午9点45分+/-15分钟
扫描带宽度：185km
重访周期：16天
景覆盖范围：184×185.2km
Landsat一5的专题制图仪(Thematic Mapper，TM)分为7个波段，各波段的参数如下所示：
在这里插入图片描述

二、Landsat_5 TM下载文件概况

在这里插入图片描述

三、 Landsat_5 TM 头文件信息预览

group 1 元数据文件信息

FILE_DATE = 2017-02-21T17:42:17Z >>>>>>数据生产日期
DATA_CATEGORY = "NOMINAL"       >>>>>>数据类型

group 2 产品元数据

ELEVATION_SOURCE = "GLS2000"    >>>>>>>>高程基准
OUTPUT_FORMAT = "GEOTIFF"        >>>>>>>>输出格式
EPHEMERIS_TYPE = "PREDICTIVE" >>>>>>>>星历表类型
SPACECRAFT_ID = "LANDSAT_5"     >>>>>>>>航空平台名称
SENSOR_ID = "TM"                             >>>>>>>>传感器名称
SENSOR_MODE = "SAM"                   >>>>>>>>传感器扫描模式
WRS_PATH = 127                               >>>>>>>>全球影像标记符号系统的影像编号
WRS_ROW = 039                                >>>>>>>>全球影像标记符号系统的影像w位置    
CORNER_UL_LAT_PRODUCT = 31.26054      >>>>>>>影像的角点纬度(一般为四个)
CORNER_UL_LON_PRODUCT = 106.38116   >>>>>>>影像的角点经度(一般为四个)
CORNER_UL_PROJECTION_X_PRODUCT = 631500.000   >>>>>>>>影像的角点X坐标
CORNER_UL_PROJECTION_Y_PRODUCT = 3459300.000  >>>>>>>影像的角点Y坐标
FILE_NAME_BAND_1 = "LT05_L1TP_127039_19860802_20170221_01_T1_B1.TIF"    >>>>>>>波段全称(一般有八个)
GROUND_CONTROL_POINT_FILE_NAME = "LT05_L1TP_127039_19860802_20170221_01_T1_GCP.txt"    >>>>>>>> 地面控制点文件
METADATA_FILE_NAME = "LT05_L1TP_127039_19860802_20170221_01_T1_MTL.txt"    >>>>>元数据文件

group 3 图像属性

CLOUD_COVER = 7.00    >>>>>>云量
CLOUD_COVER_LAND = 7.00   >>>>>>>云覆盖范围
SUN_AZIMUTH = 104.71729242   >>>>>> 太阳入射方位角
SUN_ELEVATION = 57.37237100  >>>>>>太阳入射高度角
EARTH_SUN_DISTANCE = 1.0148616    >>>>>>>>日地距离
SATURATION_BAND_1 = "Y"     >>>>>>各波段的色彩饱和度(一般有7个值)
GROUND_CONTROL_POINTS_VERSION = 4   >>>>>>地面控制点版本
GROUND_CONTROL_POINTS_MODEL = 206   >>>>>>地面控制点模型
GEOMETRIC_RMSE_MODEL = 3.719       >>>>>>RMSE回归模型
GROUND_CONTROL_POINTS_VERIFY = 2883   >>>>>>地面控制点验证
GEOMETRIC_RMSE_VERIFY = 0.118      >>>>>>>>RMSE回归模型验证

group 4 最大最小辐照度

RADIANCE_MAXIMUM_BAND_1 = 169.000   >>>>>>该波段最大辐照度
RADIANCE_MINIMUM_BAND_1 = -1.520       >>>>>该波段最小辐照度

group 5 最大最小反射率

REFLECTANCE_MAXIMUM_BAND_1 = 0.281290  >>>>>>最大反射率
REFLECTANCE_MINIMUM_BAND_1 = -0.002530  >>>>>>最小反射率

group 6 最大最小像素DN值

QUANTIZE_CAL_MAX_BAND_1 = 255   >>>>>>>>最大像素DN值
QUANTIZE_CAL_MIN_BAND_1 = 1        >>>>>>>> 最小像素DN值

group 7 产品参数
group 8 辐射测量重新缩放值

RADIANCE_MULT_BAND_1 = 6.7134E-01  >>>>>>该波段增值（辐照度）7个值
RADIANCE_ADD_BAND_1 = -2.19134  >>>>>>>>该波段偏差值 （辐照度）7个值
REFLECTANCE_MULT_BAND_1 = 1.1174E-03  >>>>>>>该波段增值（反射率）7个值
REFLECTANCE_ADD_BAND_1 = -0.003647   >>>>>>>>>该波段偏差值 （反射率）7个值

group 9 热常数

K1_CONSTANT_BAND_6 = 607.76    >>>>>>>波段6的K1常数值
K2_CONSTANT_BAND_6 = 1260.56  >>>>>>>波段6的K2常数值

group 10 投影参数

MAP_PROJECTION = "UTM"    >>>>>>>>影像投影
DATUM = "WGS84"    >>>>>>>>  参考基准
ELLIPSOID = "WGS84"   >>>>>>  参考椭球
UTM_ZONE = 48      >>>>>>>UTM投影所在带号
GRID_CELL_SIZE_REFLECTIVE = 30.00   >>>>>>>网格大小
ORIENTATION = "NORTH_UP"     >>>>>>>方位传感器 
RESAMPLING_OPTION = "CUBIC_CONVOLUTION"    >>>>>>>重采样

四、利用MATLAB读取显示Landsat_5 TM影像

注：
1）由于读取未进行拉伸，所以显示略显暗淡，若需要显示完美，可适当进行拉伸；
2）下图为各波段裁剪一部分后在matlab中显示的结果。

clc
clear all
band_data_file='C:\Users\稳魂\Desktop\学习\landsat4\';
tiff_data=dir(strcat(band_data_file,'*.tif'));
len=length(tiff_data);
image={};
for i=1:len-1image_name = tiff_data(i).name;image{i} =  imread(strcat(band_data_file,image_name));subplot(4,2,i);imshow(image{i});title("land-sat4" + "-Band" + num2str(i)+'.tiff');subplot(4,2,8);image_name = tiff_data(8).name;image{8} =  imread(strcat(band_data_file,image_name));imshow(image{8})title('landsat-4-BAQ.tiff');
end