cosi-corr是一个可以加载到ENVI的软件包，由加利福尼亚理工学院的François Ayoub, Sébastien Leprince, and Lionel Keene开发，并且提供源代码，该软件受到了NSF的资助。软件的主要功能是可以满足各种影像的正射校正和影像配准。
针对下载下来的软件包参考《cosi-corr_guide.pdf》，进行配置使用具体详细过程如下：
（本机环境介绍：win7 64位，ENVI4.7，软件包cosi-corr_pak16dec09）
一、材料
1.软件包：http://www.tectonics.caltech.edu/slip_history/spot_coseis/download_software.html（尊重原创，请注册下载）
2.数据：http://pan.baidu.com/s/1i3xIsXf
二、配置
1.论坛：http://www.tectonics.caltech.edu/forum/
2.作者邮箱：leprincs@caltech.edu
3.小编邮箱：surveymofan@163.com
4.配置
（1）将下载文件夹中的cosi_corr.sav文件放到D:\Program Files(x86)\ITT\IDL71\products\envi47\save_add 文件夹中。
（2）将NLMeansFilterDll.def、NLMeansFilterDll.dll、NLMeansFilterDll.dlm文件放到D:\Program Files (x86)\ITT\IDL71\bin\bin.x86文件夹中。
（3）http://www.microsoft.com/en-us/download/details.aspx?id=29下载Visual C++ runtime，便于编译成功
（4）调整ENVI缓存大小ENVI->File->Preferences–>Miscellaneous对话框的右下角有cache size和image title size分别设置为200和40.
（5）重启ENVI会得到如下界面
三、使用（卫星影像）（针对sample data）
1.进行变化检测的操作流程：
Ancillaryfile -->Topographic modeling（Topographic菜单）-->Select Tie Points:Image to Image -->Tie Points to GCPs -->Optimization -->Orthorectification/Resampling -->Correlation
2.整体流程图如下：
3.生成*.anc文件
输入文件：
SPOT Leader/ Dimap File:lead_01.dat
输出文件:
Ancillary Data File :lead_4_1998
记录sun elevation：63.9、sun azimuth：137.4
同理对另一幅图像如此操作，得到一个lead2_2000.anc文件（太阳高度角和太阳方位角就不用记录了，这个在生成shaded DEM时用，而只用生成一个shaded DEM就可以）。
4.生成shaded DEM文件,点击Topographic菜单中的Topographic modeling
输入文件：
选择SRTM文件DEM_UTM_SRTM_HectorMine_interp_sinc_10m
输出文件：
Output：shaded DEM
输入相应的太阳高度角和方位角，输出文件的名字为shaded DEM，点击OK，可得到shaded DEM文件。
5.打开SPOT4影像（SPOT4定义为1影像，SPOT2定义为2影像）和shaded DEM
6.以shaded DEM为基准影像，1影像为辅助影像进行选取同名点。点击Select Tie Points: Image to Image
最后保存文件为：ICP_shadedDEM_SPOT4.pts
7. 生成GCPs，点击The Points to GCPs，生成GCP_shadedDEM_SPOT4.pts文件
8.优化GCPs，点击GCPs optimization
输入文件：
raw image：imag_01.dat
reference image：shaded DEM
DEM file：DEM_UTM_SRTM_HectorMine_interp_sinc_10m
输出文件：
Optimized GCPS：GCP_optimization_shadedDEM_SPOT4.txt
9.正射纠正，点击Orthorectification/Resampling
输入文件：
ancillary file：spot4_1998.anc
GCPs file：GCP_optimization_shadedDEM_SPOT4.txt
DEM：DEM_UTM_SRTM_HectorMine_interp_sinc_10m
输出文件：
Mapping Matrices：Mapping_spot4_Matrices.mat
点击from raw image
点击蓝色按钮，点击OK
10.重采样
输入数据：
image：imag_01.dat
mapping matrices：Mapping_spot4_Matrices.mat
输出数据：
resampled image：resampled_spot4_image
11.以1影像的正射影像为基准影像纠正2影像，选取纠正后的1影像和影像2同名点
12.将ICP转化为GCP，点击Tie Points to GCPs
输入文件：
Tied Points File：ICP_SPOT4_SPOT2.pts
Reference image：resampled_spot4_image
DEM File：DEM_UTM_SRTM_HectorMine_interp_sinc_10m
输出文件：
GCPS File：GCP_SPOT4_SPOT2.pts
13.优化控制点，点击GCPS Optimization
输入文件：
Raw image：imag_02.dat
Reference image：resampled_spot4_image
DEM File：DEM_UTM_SRTM_HectorMine_interp_sinc_10m
Ancillary Data File：spot2_2000.anc
GCPS/Tied Points/ICP File：GCP_SPOT4_SPOT2.pts
输出文件：
Optimized GCPS：GCP_optimization_SPOT4_SPOT2.txt
14.正射纠正2影像，点击orthorectification/resampling
（备注：只勾选Orthorectification，即先进行旋转矩阵的生成）
输入文件：
Ancillary Data File：spot2_2000.anc
GCPS File：GCP_optimization_SPOT4_SPOT2.pts
DEM：DEM_UTM_SRTM_HectorMine_interp_sinc_10m
From the raw：imag_02.dat
（点击蓝色按钮）
输出文件：
Mapping Matrices：Mapping_spot2_Matrices.mat
15.重采样，点击点击orthorectification/resampling
（备注：只勾选Resampling，即进行重采样）
输入文件：
Image：imag_02.dat
Mapping Matrices：Mapping_spot2_Matrices.mat
输出文件：
Resampling image：resampled_spot2_image
16.利用影像相关进行变化监测，点击correlation
输入数据：
Pre-Event image：resampled_spot4_image
Post-Event image：resampled_spot2_image
输出文件：
Correlation File：Correlation_Image
17.最后得到如下结果
第一幅是东西方向的变化，第二幅是南北方向的变化，第三幅是信噪比。
（本帖供大家参考交流，不涉及任何商业用途，请读者遵守软件开发作者的规则与要求）
