ENVI扩展工具：偏相关系数计算器

1 偏相关系数

本节内容来自百度百科。

在多要素所构成的系统中，当研究某一个要素对另一个要素的影响或相关程度时，把其他要素的影响视作常数（保持不变），即暂时不考虑其他要素影响，单独研究两个要素之间的相互关系的密切程度，所得数值结果为偏相关系数。

在多元回归分析中，在消除其他变量影响的条件下，所计算的某两变量之间的相关系数。在多元相关分析中，简单相关系数可能不能够真实的反映出变量X和Y之间的相关性，因为变量之间的关系很复杂，它们可能受到不止一个变量的影响。这个时候偏相关系数是一个更好的选择。

2 工具安装

因为调用了新版接口，所以本工具只能运行在 ENVI5.5 及以上版本。

可访问 envi.geoscene.cn/envi_license 申请最新版ENVI试用许可。

方法1：使用 ENVI App Store 进行 ENVI 扩展工具的安装与管理。

App Store 下载地址：envi.geoscene.cn/appstore

方法2：手动下载 https://envi.geoscene.cn/appstore/parleieent

下载 zip 压缩包并解压，将得到的 extensions 和 custom_code 文件夹拷贝到如下 ENVI 安装路径，覆盖同名文件夹即可：

• ENVI 5.5 及以上 - C:\Program Files\Harris\ENVI5x\

重启 ENVI 即可使用。

注：本工具支持 ENVI 5.5 及更高版本，使用 ENVI 5.6 测试。

3 工具使用

工具采用分块技术，对输入图像大小没有限制。

工具位于Toolbox/Extensions/Partial Correlation Coefficient

双击启动，弹出如下主界面：

输入参数：

· Independent Variable Raster：指定自变量数据。

· Dependent Variable Raster：指定因变量数据。

· Controlled Variable Rasters：指定要移除其影响的控制变量数据。可指定一个或多个栅格图像。

工具要求和特性：

· 输入的所有栅格必须具备相同的波段数；

· 当有栅格行列数不同时，将自动进行波段组合使其行列数一致后再计算偏相关系数；

· 结果背景值为NaN。

特殊说明：

核心功能是调用的IDL函数P_CORRELATE，在计算过程中个别像元有可能会报如下错误：

Inversionfailed due to singular array or small pivot element.

为了完成计算，修改源代码，使其在报此错误时返回0，即计算结果中此像元值为0。

4 接口调用

因为采用ENVITask开发方式，可用在ENVI Modeler、ENVI Server、IDL中。

下面是IDL中调用语法，其中INPUT_RASTER3可输入1个raster或raster数组。

task = ENVITask('PartialCorrelationCoefficient')
task.INPUT_RASTER1= raster1
task.INPUT_RASTER2= raster2
task.INPUT_RASTER3= raster3
task.OUTPUT_RASTER_URI= output_uri
task.Execute
outraster= task.OUTPUT_RASTER