[转帖]tm的大气校正

tm的大气校正（回归分析法）

---------------------------------------------------------------------

发信人: ahdou0 (ahdou一到，众口不再难调), 信区: GIS
标 题: tm的大气校正（回归分析法）
发信站: 南京大学小百合站 (Thu Jun 3 23:02:46 2004)

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;
PROJECTION 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 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_Float) / ($n8_Float - ($n
6_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_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 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";
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.0222;
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 UNSIGNED8;
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.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 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_Integer) )) / (GLOBAL MEAN ( $n28_lani
er($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;
}
FUNCTION {
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_Output ) ";
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 century schoolbook";
SIZE 24;
}
SCALAR {
ID 46;
TITLE "n46_Integer";
POSITION 3.47778, 0.877778;
DATATYPE SIGNED32;
VALUE 1;
CHILD 35, 32, 29, 42;
}
SCALAR {
ID 47;
TITLE "n47_Integer";
POSITION 5.36667, 0.833333;
DATATYPE SIGNED32;
VALUE 2;
CHILD 42, 39;
}

<img src="images/post/smile/dvbbs/em02.gif" />

<P>有什么用呢？？</P><P>大侠介绍一下呗</P>

<a href="http://www.gisempire.com/bbs/showimg.asp?BoardID=58;filename=2004-9/20049217374720.jpg" target="_blank" ><img src="http://www.gisempire.com/bbs/showimg.asp?BoardID=58;filename=2004-9/20049217374720.jpg"></A>

没有用到大气参数,能做好校正,还是比较难的.