TM的大气纠正-代码

1、 TM 的大气校正：在没有地面实测数据或者其它气象数据的情况下，通常用波段回归分析，进行波段间相对校正，其原理很多遥感书中都有，一般是以第 7 波段为校正波段，将其它波段一一与之比较，确定各波段的偏移量。将以下部分拷贝，另存为*.gmd 文件，在 erdas 的 model 里运行该模型即可PAGESIZE 8.66667, 11.5444 INCHES;CELLSIZE MINIMUM;PRINTERPAGESIZE 8.5, 11;MARGINS 0.5, 0.5, 0.5, 0.5;ORIENTATION PORTRAIT;PRINTSCALE 100;WINDOW UNION;PROJ

2、ECTION DEFAULT;AOI NONE;OPTIMIZE YES;SCALAR ID 4;TITLE n4_Float;POSITION 0.9, 5.94445;DATATYPE DOUBLE;VALUE 0;SHOW;CHILD 13, 15;SCALAR ID 6;TITLE n6_Float;POSITION 2.14444, 5.95556;DATATYPE DOUBLE;VALUE 0;CHILD 13, 15;SCALAR ID 8;TITLE n8_Float;POSITION 3.84444, 5.95556;DATATYPE DOUBLE;VALUE 0;CHILD

3、 13;SCALAR ID 10; TITLE n10_Float;POSITION 5.54444, 5.95556;DATATYPE DOUBLE;VALUE 0;CHILD 13;SCALAR ID 12;TITLE n12_Float;POSITION 7.04444, 5.95556;DATATYPE DOUBLE;VALUE 0;CHILD 13, 15;FUNCTION ID 13;TITLE $n10_Float - ;POSITION 2.13333, 7.45556;VALUE ($n10_Float - ($n6_Float * $n12_Float) / $n4_Flo

4、at) / ($n8_Float - ($n6_Float * $n6_Float) / $n4_Float);AREA UNION;CHILD 14;SCALAR ID 14;TITLE n14_Float;POSITION 2.12222, 8.7;DATATYPE DOUBLE;VALUE 0;SHOW;CHILD 15;FUNCTION ID 15;TITLE $n12_Float / $n4_Float - $n14_Float * ;POSITION 6.02222, 8.71111;VALUE $n12_Float / $n4_Float - $n14_Float * ($n6_

5、Float / $n4_Float);AREA UNION;CHILD 16;SCALAR ID 16;TITLE n16_Float;POSITION 6.02222, 10.0444; DATATYPE DOUBLE;VALUE 0;SHOW;TEXT ID 17;TITLE Linear Regression;POSITION 4.21111, 0.155556;FONT new century schoolbook;SIZE 24;STYLE BOLD;STYLE ITALICS;STYLE UNDERLINED;TEXT ID 18;TITLE X = Band1;POSITION

6、3.54445, 0.766671;FONT new century schoolbook;SIZE 24;TEXT ID 20;TITLE n;POSITION 0.844444, 5.95556;FONT new century schoolbook;SIZE 24;TEXT ID 21;TITLE Sum X;POSITION 2.14444, 5.95556;FONT new century schoolbook;SIZE 24;TEXT ID 22;TITLE Sum X*2;POSITION 3.84444, 5.95556;FONT new century schoolbook;

7、SIZE 24;TEXT ID 23; TITLE Sum X*Y;POSITION 5.56666, 5.94445;FONT new century schoolbook;SIZE 24;TEXT ID 24;TITLE Sum Y;POSITION 7.04444, 5.95556;FONT new century schoolbook;SIZE 24;TEXT ID 25;TITLE a;POSITION 2.12222, 8.7;FONT new century schoolbook;SIZE 24;TEXT ID 26;TITLE b;POSITION 6.03333, 10.02

8、22;FONT new century schoolbook;SIZE 24;TEXT ID 27;TITLE y = aX + b;POSITION 3.66666, 9.91111;FONT new century schoolbook;SIZE 24;RASTER ID 28;TITLE n28_lanier;POSITION 0.933334, 0.7;FILENAME d:/program files/imagine 8.4/examples/lanier.img;INTERPOLATION BILINEAR;ATHEMATIC;这里就是我说的默认存储路径DATATYPE UNSIG

9、NED8;DECLARE Integer;COMPRESSION UNCOMPRESSED;COORDINATES MAP; AOI NONE;RECODE NO;CHILD 32, 29, 35, 39, 42;FUNCTION ID 29;TITLE GLOBAL SUM ;POSITION 3.95555, 2.15556;VALUE GLOBAL SUM ( ($n28_lanier($n46_Integer) * $n28_lanier($n46_Integer);AREA UNION;CHILD 30;TABLE ID 30;TITLE n30_Output;POSITION 3.

10、86667, 3.51111;SIZE 0;DATATYPE DOUBLE;OUTPUT;CHILD 31;FUNCTION ID 31;TITLE GLOBAL SUM ;POSITION 3.86666, 4.82223;VALUE GLOBAL SUM ( $n30_Output ) ;AREA UNION;CHILD 8;FUNCTION ID 32;TITLE GLOBAL SUM ;POSITION 2.32222, 2.17778;VALUE GLOBAL SUM ( $n28_lanier($n46_Integer) );AREA UNION;CHILD 33;TABLE ID

11、 33;TITLE n33_Output;POSITION 2.26666, 3.48889;SIZE 0;DATATYPE DOUBLE; OUTPUT;CHILD 34;FUNCTION ID 34;TITLE GLOBAL SUM ;POSITION 2.23333, 4.84445;VALUE GLOBAL SUM ( $n33_Output ) ;AREA UNION;CHILD 6;FUNCTION ID 35;TITLE GLOBAL SUM ;POSITION 0.911111, 2.15556;VALUE (GLOBAL SUM ( $n28_lanier($n46_Inte

12、ger) ) / (GLOBAL MEAN ( $n28_lanier($n46_Integer) , IGNORE 0 ) );AREA UNION;CHILD 36;TABLE ID 36;TITLE n36_Output;POSITION 0.922222, 3.48889;SIZE 0;DATATYPE DOUBLE;OUTPUT;CHILD 37;FUNCTION ID 37;TITLE GLOBAL SUM ;POSITION 0.922222, 4.82222;VALUE GLOBAL SUM ( $n36_Output ) ;AREA UNION;CHILD 4;FUNCTIO

13、N ID 39;TITLE GLOBAL SUM ;POSITION 7.18889, 2.18889;VALUE GLOBAL SUM ( $n28_lanier($n47_Integer);AREA UNION;CHILD 40; TABLE ID 40;TITLE n40_Output;POSITION 7.22222, 3.44445;SIZE 0;DATATYPE DOUBLE;OUTPUT;CHILD 41;FUNCTION ID 41;TITLE GLOBAL SUM ;POSITION 7.08889, 4.85556;VALUE GLOBAL SUM ( $n40_Outpu

14、t ) ;AREA UNION;CHILD 12;FUNCTION ID 42;TITLE GLOBAL SUM ;POSITION 5.55555, 2.21111;VALUE GLOBAL SUM ($n28_lanier($n46_Integer) * $n28_lanier($n47_Integer) );AREA UNION;CHILD 43;TABLE ID 43;TITLE n43_Output;POSITION 5.5, 3.52223;SIZE 0;DATATYPE DOUBLE;OUTPUT;CHILD 44;FUNCTION ID 44;TITLE GLOBAL SUM ;POSITION 5.55555, 4.78889;VALUE GLOBAL SUM ( $n43_Output ) ;AREA UNION;CHILD 10; TEXT ID 45;TITLE Y= Band2;POSITION 5.33333, 0.766671;FONT new c