摘要: 一、Landsat7的读入步骤： 1. 把Erdas Imagine8.5安装目录下的D:Program FilesIMAGINE 8.5etcsphoid.tab的投影Krasovsky改为Krassovsky。 2. 修改TM影像的头文件，按照删除一个字符...
一、Landsat7的读入步骤：
1. 把Erdas Imagine8.5安装目录下的D:Program FilesIMAGINE 8.5etcsphoid.tab的投影Krasovsky改为Krassovsky。
2. 修改TM影像的头文件，按照删除一个字符，增加一个空格和增加一个字符，删除一个空格的原则，进行修改，以保证该头文件的总bytes保持不变。例如：bands=123457改为bands=2345后面加两个空格。并且把1和7波段的相应的文件名删去，例如：FILENAME=L71125027_02720000701_B30.FST改为FILENAME = 中间有29个空格，最好把1和7波段相对应的GAINS AND BIASES IN ASCENDING BAND NUMBER ORDER也按照上述原则删去。
3. 修改中央经线：修改其中央经线的经纬度数。在同一个带的5和7的投影参数进行对调，例如把117000000修改为2.042035000000000。各带的参数不同。例如：下列是修改后的一个例子。
BANDS PRESENT =2345
FILENAME = FILENAME =L71122026_02620000509_B20.FST
FILENAME =L71122026_02620000509_B30.FST FILENAME =L71122026_02620000509_B40.FST
FILENAME =L71122026_02620000509_B50.FST FILENAME =
REV L7A
GAINS AND BIASES IN ASCENDING BAND NUMBER ORDER
-6.000000000000000 0.817254878025429
-4.500000000000000 0.639607867072610
-4.500000000000000 0.939215686274510
-1.000000000000000 0.128470589132870
GEOMETRIC DATA MAP PROJECTION =TM ELLIPSOID =KRASSOVSKY DATUM =KRASSO
USGS PROJECTION PARAMETERS = 6378245.0000000000000 6356863.0187999997000
1.0000000000000 0.0000000000000 2.146755000000000000000
0.0000000000000 500000.0000000000000 0.0000000000000
0.0000000000000 0.0000000000000 0.0000000000000
0.0000000000000 0.0000000000000 0.0000000000000
0.0000000000000
USGS MAP ZONE = 0
UL = 1200802.1991E 495056.4839N 293884.179 5528093.264
UR = 1225723.5268E 492612.6448N 496847.458 5478308.251
LR = 1222317.1513E 475221.4484N 454223.597 5304539.517
LL = 1193856.1625E 481623.2975N 251260.318 5354324.530
CENTER = 1211660.0000E 485160.0000N 374053.888 5416316.390 3482 2981
OFFSET = 0 ORIENTATION ANGLE =-13.78
SUN ELEVATION ANGLE =55.1 SUN AZIMUTH ANGLE =147.9
4. 导入Erdas Imagine8.5，采用TM Landsat7 Fast-L7A EROS类型即可，保存使用*.img格式的图像。
5. 进行投影转换。目标文件的投影类型为Gauss_krug，椭球体为Krassovsky，Datum为undefined，中央经线的垂直比例尺为1.00000 ，向东偏移量为215000000（其它带不同）米。中央经线各带之间不同。其它参数默认即可。
6. .img到.tif的转换，用img2tif转换程序：img2tif –I f:12028.img –o f:12028.tif –rgb 3 2 1即可。
二、Landsat5的读入步骤：
1. 不要修改Erdas安装目录下的投影文件，保持Krasovsky不变。
2. 修改头文件中的Krassovsky为Krasovsky，即删除一个s在后面添加一个空格
3. 修改TM影像的头文件，按照删除一个字符，增加一个空格和增加一个字符，删除一个空格的原则，进行修改，以保证该头文件的总bytes保持不变。例如：bands=1234567改为bands=234后面加四个空格。
4. 修改中央经线参数：把原来的2.042035000000000修改为1170000.000000000，注意中间不要改变其子节数。
5. 修改中央经线129.000000为1290000.0000即在129后面添加4个0，然后删除4个0 例如：下面是修改好的一个头文件：
PROJECTION PARAMETERS = 6378245.000000000000000 6356863.018773047300000
1.000000000000000 0.000000000000000 1230000.000000000
0.000000000000000 500000.000000000000000 0.000000000000000
0.000000000000000 0.000000000000000 0.000000000000000
0.000000000000000 0.000000000000000 0.000000000000000
0.000000000000000 EARTH ELLIPSOID =Krasovsky
SEMI-MAJOR AXIS =6378245.000 SEMI-MINOR AXIS =6356863.019
PIXEL SIZE =30.00 PIXELS PER LINE= 6523 LINES PER IMAGE= 5728
UL 1231757.4392E 453223.3317N 21523375.630 5045070.272
UR 1254436.9378E 451051.3063N 21715644.928 5008802.481
LR 1251644.3237E 434018.5874N 21683798.004 4839969.873
LL 1225339.6486E 440115.2910N 21491528.706 4876237.664
BANDS PRESENT =2345 BLOCKING FACTOR = 2 RECORD LENGTH =13046
SUN ELEVATION =60 SUN AZIMUTH =133
CENTER 1241746.9434E 443614.8336N 21602918.094 4941855.922 3244 2890
OFFSET= 43 REVB
6.导入Erdas Imagine8.5，采用TM Landsat EOSAT Fast Format类型即可，保存使用*.img格式的图像。
7.进行投影转换。目标文件的投影类型为Gauss_krug，椭球体为Krasovsky，Datum为undefined，中央经线的垂直比例尺为1.00000，向东偏移量为215000000（其它带不同）米。中央经线各带之间不同。其它参数默认即可。
8 _ . img到 _ . tif的转换，用img2tif转换程序：
img2tif –I f:12028.img –o f:12028.tif –rgb 3 2 1即可。
三、LandSat5和LandSat7中央经线对调的相互关系：
50-2.042035000000000 -1170000.
51-2.146755000000000 -1230000.
52-2.251475000000000 -1290000.
53-2.356194000000000 -1350000.