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

原理：参考博客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='待