cosi-corr操作详细步骤

       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.作者邮箱：[e-mail]

     3.小编邮箱：[e-mail]

     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.最后得到如下结果

      第一幅是东西方向的变化，第二幅是南北方向的变化，第三幅是信噪比。

 

（本帖供大家参考交流，不涉及任何商业用途，请读者遵守软件开发作者的规则与要求）