ENVI/IDL编程：批量使用伪不变目标法进行相对辐射校正

最新推荐文章于 2025-06-19 21:57:49 发布

银木荷

最新推荐文章于 2025-06-19 21:57:49 发布

阅读量1.7k

点赞数 4

CC 4.0 BY-SA版权

分类专栏： ENVI/IDL 文章标签： ENVI/IDL

本文链接：https://blog.youkuaiyun.com/qq_41300185/article/details/87374427

这篇博客介绍了如何使用ENVI/IDL编程批量进行伪不变目标法的相对辐射校正。通过导出ROI，计算一元线性回归方程的系数，然后对多个待校正影像进行校正。文中提供了四种不同的校正思路，并验证了代码结果与已知方法的一致性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

原理：参考博客http://blog.sina.com.cn/s/blog_764b1e9d0100pwby.html

可直接利用的工具：http://blog.sina.com.cn/s/blog_764b1e9d0102y3b6.html

IDL矩阵运算参考：http://blog.sciencenet.cn/home.php?mod=space&uid=346157&do=blog&quickforward=1&id=886512

矩阵运算：http://homepages.see.leeds.ac.uk/~lecjm/Teaching/IDL_course/Notes/notes/node22.html

既然都已经有了工具，为什么还要写这篇博文呢?当然是想要批量批量批量！假设已有伪不变目标（.roi格式），使用该roi，以时刻T1为参考影像，其余多个时序为待校正时序，难到要一个一个的操作吗？哈哈哈哈哈哈哈哈哈！当然也是可以的，如果有代码，岂不是更好呢。本着一贯折磨自己的原则，踏上了代码的不归路。。。

主要的一元线性回归方程 y=ax+b的系数计算如下：

源代码如下：

pro zy3relativenormpif
;2019/2/15 采用PIF相对辐射校正方法，批量操作
;数据：1.伪不变目标（地物）： ROI格式 2.参考影像 3.待校正影像（与前者的地理位置应该相同，否则ROI找不到位置
;方法：参考博文 http://blog.sina.com.cn/s/blog_764b1e9d0100pwby.html
;主要步骤：1.分别导出各影像的roi为坐标点 2.计算回归系数 3.波段运算校正系数
COMPILE_OPT idl2
e=envi()
DIRPATH='待校正的影像路径'
ROIPATH='不变特征roi路径'
RESPATH='输出路径'
filelist=FILE_SEARCH(DIRPATH,'*.hdr',count=countall)
refername='参考时序的影像名称'
referpath=FILEPATH(refername,ROOT_DIR=DIRPATH)
rroispath=FILEPATH('referroi.c