用Erdas 8.5读入Landsat 5、Landsat 7的投影参数

一、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.