MATLAB关于偏相关系数,基于matlab的逐像元偏相关分析

本文介绍基于matlab的逐像元的三者间的偏相关，以产水，NPP和土壤保持为例进行说明

具体代码如下：

%将三者多年的数据放在三个不同的矩阵中

[a,R]=geotiffread('F:\校级课题项目\data\屏障带\2002water_yield.tif');%先导入投影信息

info=geotiffinfo('F:\校级课题项目\data\屏障带\2002water_yield.tif');

nppsum=zeros(size(a,1)*size(a,2),16);

for year=2000:2015

filename=strcat('F:\校级课题项目\data\屏障带\',int2str(year),'npp.tif');

data=importdata(filename);

data=reshape(data,size(a,1)*size(a,2),1);

nppsum(:,year-1999)=data;

end

scsum=zeros(size(a,1)*size(a,2),16);

for year=2000:2015

filename=strcat('F:\校级课题项目\data\屏障带\',int2str(year),'sc.tif');

data=importdata(filename);

data=reshape(data,size(a,1)*size(a,2),1);

scsum(:,year-1999)=data;

end

wcsum=zeros(size(a,1)*size(a,2),16);

for year=2000:2015

filename=strcat('F:\校级课题项目\data\屏障带\',int2str(year),'water_yield.tif');

data=importdata(filename);

data=reshape(data,size(a,1)*size(a,2),1);

wcsum(:,year-1999)=data;

end

%控制NPP，看产水和土壤保持的偏相关

rho_value=zeros(size(a,1),size(a,2))+nan;

p_value=zeros(size(a,1),size(a,2))+nan;

for i=1:size(a,1)*size(a,2)

nppdata=nppsum(i,:);

if min(nppdata)>0

nppdata=nppdata';

scdata=scsum(i,:)';

wcdata=wcsum(i,:)';

[rho,p]=partialcorr(scdata,wcdata,nppdata);%注意，控制的变量放在最后面

rho_value(i)=rho;

p_value(i)=p;

end

rho_value(p_value>0.05)=NaN;

filename='F:\课题项目\data\通过显著性0.05检验的产水和土壤保持偏相关系数.tif';

geotiffwrite(filename,rho_value,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);

更多需求，请查看个人介绍

MATLAB关于偏相关系数,基于matlab的逐像元偏相关分析相关推荐

用matlab相关分析,基于matlab的逐像元偏相关分析
本文介绍基于matlab的逐像元的三者间的偏相关,以产水,NPP和土壤保持为例进行说明具体代码如下: %将三者多年的数据放在三个不同的矩阵中 [a,R]=geotiffread('F:\校级课题项目 ...
matlab设计译码器,基于MATLAB的循环码编译码器设计与仿真.doc
扳昂旨螺冈唉陨裤外狸尿恨铸伸隧刽搅必勒诚天腑皖漂豌鲁靳碑缆键兽峙棘陶宽槐撒层僧袁廖颤渐魄货鼎躬薛扬衍逮西兰迫依煤鲁虐渠惫平合啥昭并屿己笆坍痞庐披吏去凄嘛兄察突徊溅今箩直藩潦咙锨谓崇若制匹扮复淌颐糖嗅你 ...
「电子万年历matlab仿真」——基于Matlab的电子万年历仿真实现
「电子万年历matlab仿真」--基于Matlab的电子万年历仿真实现作为一种具有时间显示.日期查询.闹钟提醒等功能的电子产品,电子万年历已经成为了人们日常生活中不可或缺的一部分.而在现代科技的发展 ...
MATLAB计算杨氏模量,基于MATLAB处理动态图像的杨氏模量测量方法.doc
7A版优质实用文档 PAGE PAGE 16 7A版优质实用文档基于MATLAB处理动态图像的杨氏模量测量方法黄浩钊1,陈玥琦2,李家卉3 第三军医大学学员13营,重庆,400038:2.第三军 ...
matlab如何处理动态分配,基于MATLAB的水资源优化分配问题动态规划解法
摘要:介绍了动态规划的基本原理,针对水资源分配问题进行了动态规划方法分析.针对具体问题采用逆序解法的表格法进行了计算,然后用MATLAB编制了相应的计算程序进行计算,避免了繁琐的人工计算.结果表明该方 ...
matlab水果图片,基于matlab的水果识别的数字图像处理
基于matlab的水果识别的数字图像处理图像处理(报告)题目: 基于 MATLAB 的水果识别的数字图像处理指导教师: 职称: 教授学生姓名: 学号:专业: 院(系): 完成时间: 2016 ...
matlab求解外弹道,基于MATLAB∕Simulink的通用质点外弹道程序设计.pdf
2017.01 设计与研发基于 MATLAB/Simulink 的通用质点外弹道程序设计崔瀚 (沈阳工学院兵器类虚拟仿真实验教学中心, 辽宁抚顺 ,113122 ) 摘要:本文以弹丸质点外弹道学 ...
FLAC3D可视化后处理matlab,一种基于Matlab的由Midas导入Flac3D的模型识别方法与流程...
本发明涉及岩土工程的仿真模拟研究领域,具体涉及一种基于Matlab的由Midas导入Flac3D的模型识别方法. 背景技术: 随着我国近几年经济的快速发展和基础设施等的大力投资,涉及复杂地质环境下的岩 ...
matlab演示系统,基于Matlab的通信原理演示系统的设计与应用
基于 Matlab的通信原理演示系统的设计与应用李强 , 明艳 , 吴坤君 (重庆邮电大学通信学院 , 重庆 400065) 摘要 : 利用 Matlab图形用户界面的开发环境和强大的通信仿 ...

MATLAB关于偏相关系数,基于matlab的逐像元偏相关分析

MATLAB关于偏相关系数,基于matlab的逐像元偏相关分析相关推荐

最新文章

热门文章