《遥感影像相对辐射校正实验》由会员分享，可在线阅读，更多相关《遥感影像相对辐射校正实验（5页珍藏版）》请在金锄头文库上搜索。
1、相对辐射校正打开待进行校正的同一地区的两个时相的影像.本实验以突泉县为例1. 打开两幅影像，观察是否数据无误（包括波段顺序是否正确,影像有无较大偏 移）.2. 观察两幅影像的目视效果，或者通过ENVI菜单Basic Tools中的statistics的 compute statistics来获取两幅影像的像元DN值分布.将两幅中DN值差异大, 或者更接近真实值的作为主影像;将DN值差异小,或者偏离真实值的影像作为 从影像，来进行校正.注意:一定记清楚哪幅为主影像,哪幅为从影像.（突泉县主影像为2007年影像, 从影像为2010年影像）3. 获取两时相影像的pif点.Pif点即伪不变特征点，通常
2、选取影像中的深水面，无植被覆盖的沙地，面积较大 的建设用地等.PIF点选取应当广泛一些.首先将两幅影像Link起来.利用ENVI的ROI在主影像上选取pif点.选取的时候 每个PIF点应当面积小，准确.在主影像选取时，应当link着从影像，以便观察选取的PIF点是否发生变化.如果发生变化，应当弃选.祯I TjuIFile R0I_TjTe Opti ons HelpWiradow： Image ( Scroil * Zoom1- OffROI目孙归C olcrPiselsGoto Stats Grow PixelDalet eFM皿 #1 R2 笛6/：36 土#1 Scroll (0.Sei
3、虹t AIL I | Hid: ROIw | Show BO选好之后，将ROI存储为ASCII.存储的时候提示以哪幅影像为基准，那么通过两次操作选取不同基准的影像，就得到了两幅影像的PIF点数据.（请将生成的数据命名时添加县名和年份等必要信息）#1 ROI ToolI i=i | 回4, 获得相对辐射校正系数.这一步利用最小二乘法处理数据.将获得的数据.创建一个EXCEL空表格.为了方便查找,最好现在就把它重命名一下（我将我的 表格命名为突泉数据）.打开表格.将主影像生成的txt也打开，将里面的从；ID 开始的文本拷贝到表格中.注意粘贴后，下方的粘贴选项.选择使用文本导入 向导，持续下一步，最
4、后完成.这时，每个数据元分别放在了单独的表格空格-i0字沌对齐方式中.ABcDEF23272531112OC5072355. 0045. 70CC.D?121. C232S27122113D2025. 0050393B5. 0045. 499636121. 时29*221 ,R日矩.Mw冲姑一 nn45. i9932121 . f30蓉的，UL50393B5. 0045. 493&27ML E32312和2212392025. 005039355. 0045. 493366121. 333227222123X：G5. OC5039355.0045. 499357121. C3133237322
5、12舞 1995, 005039355. 0045. 499362121. E*、34A 4-椅3B925. 00祁.m。 nn45. 27SSG1121. E命上.yiZ55SUL5025975. 004b. ytbyi11. t379F0C5025945. 0045. 373321121. E敏字占 1，ID Y hap 基33堕33妨4142将不必要的数据而膝，只留下必要的ID和6个波段即可，下图所示.（这其中有必要自己稍微调整一下表头部分，但并不麻烦，此处不赘述.）霁响事-于理甜TTCTFik| iSX子_jABcDEFGHLIDZOQTElZ007E22007E32Q07E42007
6、E52D07E6215023ia1286S/23IF1 s9nII Jocrs-Q1 1 oac我们将这个sheetl重命名为突泉2007”.同样的方法，将从影像的pif数据写入sheet2.并命名为”突泉2010”.进行数据处理时，我们需要一个波段一个波段的来.拿第一波段BI为例.一. 我们将两个时相的B1数据拷贝到sheet3中（当然最好将sheet3改名为”B1” 之类的）.二. 做散点图.选中B1页的数据（包括表头），点击”插入”中的”图表”之”散点图”, 第一类图即可.三. 线性回归.打开”数据”右端的”数据分析”.如下图所示.四.分析工具中选择”回归”.确定后弹出回归表.方若件杵；
7、T舌日右Fl束井析 相羌泰数枯若亮睡光订指数平滑?-宜拉率林办方至 冲EI汁廿析1折打一且-需要选中或填写的有：X、Y值的输入区域（），标志，置信度（95%），新 工作表组，残差，线性拟合图。确定取消帮助（H）本实验中,X值为从影像的值;Y的值为主影像的值.回归输入Y值输入区域（X）:x值输入区域QD:v标志田）M置信度（1）输出区域（jo）： 佥，新工作表组史）:新工作薄地）残差0残差残差图但）标准残差（T）丁正态分布正态概率图但）注意：选中数据“标志”和不选“标志”，X、Y值的输入区域是不一样的:前者包括数据标志:E007BLE010B1后者不包括数据标志.(前面实验的步骤显然是包括了数据
8、标志，即表头).五.确定后得到分析结果.如下图所示：|也、W 4TIp4 - V突泉农魅世kx - MiCfSKoftExcT_flX-7TR3、 瓦纣 *目阅ifejiH XUS (? 一 M CL 61747S L7111T02AECDEFGKIJXLNI.副丽ARE OUTPUT2010Q1 Line Fit Plot1回怛统计4 iltlple (L 95BB35 K SQUiTt (L 9212T96 Adjusted (L918B&3 T折谁误差3.1617658溉爵值LC方差左折11 df归差计一目我6O3 4ss ns fmL fl ed:己 f9畦 7+y 9977.749
9、397?2.曲-网明9.E盼8 9M粕琵9 431 巳 &39L517 Interc&pl 15-7&C513沱DM日15&旎胜标推误左 1占T&t fvalue Lwer 9瑞11耸就勺牌!7眼9#傅上限沮QW 盅组眺的?b nem s. iBE-oa 11.2952 20- sees il 翌5网 脸 知龊 队茹：!葛JA9H5g 42E-扪 蛟浏。熠1构覆也5演华，0,饥电7辿EO13DJ.5D2*O1QC130212223jtEsmukL cnrrrr啾值252627232949- 97K.4凯 4151243uSDET349. 07K552. EZHX残差0. 02M A 5B0E3
10、4 Lu 1412&91. 0ZO3 L 77E332L 77623LVI 57614232E： S3.34161 -L 2151司5| 52.7Blg| 0. 21B573闷1H 奖杲河7T 突呆哀11。 凤源口，in 3的Ml ，扯*3 石El UM6b ：-： 不过我们需要的只是两个数据，即第三个表中的Intercept(记作b)和2010B1(b).相对辐射校正公式为y=ax+b;同样的，我们可以获取其他波段校正需要的a和b.5, 波段计算并生成影像在bandmath中输入公式,需要注意的是公式外需要加上uint().如uint(0.56386b1+15.6)，这样生成波段的数值为整型.将6个波段全部计算完成后,生成一幅新的影像.最后再利用矢量边界数据生成掩膜，将界外区域的DN值变成0.实验完成.
