matlab计算栅格数据逐像元趋势,基于Matlab的栅格数据一元线性回归及显著性检验(slope趋势分析)...
%by yinlichang3064@163.com
在进行长时间序列的栅格数据分析时,如NDVI,fvc,LAI,NPP,GPP,需要知道每个格点的长期趋势。
如果再arcgis中进行一元回归计算,需要将整个公式展开,并且容易出错,而matlab提供了强大的回归
计算功能,通过逐像元的迭代,可以得到每个格点的趋势及显著性栅格图像,话不多说,见以下代码
[a,R]=geotiffread('D:\日降水插值\年pre2000.tif');%先导入某个图像的投影信息,为后续图像输出做准确
info=geotiffinfo('D:\日降水插值\年pre2000.tif');
[m,n]=size(a);
years=16; %表示有多少年份需要做回归
data=zeros(m*n,years);
k=1;
for year=2000:2015 %起始年份
file=['D:\日降水插值\年prec',int2str(year),'.tif'];%注意自己的名字形式,这里使用的名字是年prec2000.tif,根据这个可修改
bz=importdata(file);
bz=reshape(bz,m*n,1);
data(:,k)=bz;
k=k+1;
end
xielv=zeros(m,n);p=zeros(m,n);
for i=1:length(data)
bz=data(i,:);
if max(bz)>0 %注意这是进行判断有效值范围,如果有效范围是-1到1,则改成max(bz)>-1即可
bz=bz';
X=[ones(size(bz)) bz];
X(:,2)=[1:years]';
[b,bint,r,rint,stats] = regress(bz,X);
pz=stats(3);
p(i)=pz;
xielv(i)=b(2);
end
end
name1='D:\simulation\黄河流域降水0.1度\year\趋势\年prec一元线性回归00-15趋势值.tif';
name2='D:\simulation\黄河流域降水0.1度\year\趋势\年prec一元线性回归00-15_P值.tif';
geotiffwrite(name1,xielv,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
geotiffwrite(name2,p,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
%一般来说,只有通过显著性检验的趋势值才是可靠的
xielv(p>0.05)=NaN;
name1='D:\simulation\黄河流域降水0.1度\year\趋势\通过显著性0.05检验的年prec一元线性回归00-15趋势值.tif';
geotiffwrite(name1,xielv,R,'GeoKeyDirectoryTag',info.GeoTIFFTags.GeoKeyDirectoryTag);
在用上述代码的时候请引用以下文献:
Yin et al. 2020. Irrigation water consumption of irrigated cropland and its dominant factor in China from 1982 to 2015. Advance in water resources. 143. doi: 10.1016/j.advwatres.2020.103661
Yin, L.; Wang, X.; Feng, X.; Fu, B.; Chen, Y. A Comparison of SSEBop-Model-Based Evapotranspiration with Eight Evapotranspiration Products in the Yellow River Basin, China. Remote Sens. 2020, 12, 2528.
更多需求,请查看个人介绍
matlab计算栅格数据逐像元趋势,基于Matlab的栅格数据一元线性回归及显著性检验(slope趋势分析)...相关推荐
- MATLAB计算英文符号信源熵,基于matlab的英文文本信息量的统计-正确.doc
基于matlab的英文文本信息量的统计-正确.doc 实验一英文文本信息量的计算一. 实验目的及要求a) 实验目的1. 通过本实验熟悉Matlab软件编程环境2. 编写M文件实现对英文文本信息量的统计 ...
- matlab计算结构自振频率,基于Matlab的多自由度系统固有频率及振型计算
基于Matlab的多自由度系统固有频率及振型计算 阅读:25212010-04-13 21:38 标签:杂谈 可参考文涛,基于Matlab语言的多自由度振动系统的固有频率及主振型计算分析,2007 对 ...
- 基于matlab 的燃油喷雾图像处理方法,基于MATLAB的燃油喷雾图像处理方法.doc
基于MATLAB的燃油喷雾图像处理方法 基于MATLAB的燃油喷雾图像处理方法 摘要:提出了基于MATLAB的燃油喷雾图像处理方法.通过对喷雾图像的采集, 分割和滤波处理, 将多幅图像转换融合为一幅 ...
- matlab x(n)16点DFT,[基于MATLAB的数字信号处理实例分析.doc
[基于MATLAB的数字信号处理实例分析 湖北文理学院理工学院 学生结业论文 课程名称:MATLAB教程 结业论文名称:基于MATLAB的数字信号处理实例分析 专业名称:通信工程 班级:1011 学号 ...
- matlab 线性回归 参数显著性,基于Matlab的栅格数据一元线性回归及显著性检验(slope趋势分析)...
%by yinlichang3064@163.com 在进行长时间序列的栅格数据分析时,如NDVI,fvc,LAI,NPP,GPP,需要知道每个格点的长期趋势. 如果再arcgis中进行一元回归计算, ...
- 基matlab的水果识别的应用,基于MATLAB的水果识别的数字图像处理
基于MATLAB的水果识别的数字图像处理 图像处理 ( 报告 ) 题目 基于 MATLAB 的 水果识别的数字图像处理 指导教师 职称 教授 学生姓名 学号 专 业 院(系) 完成时间 2016 年 ...
- matlab speex的语音处理模块_基于MATLAB的条形码识别系统[GUI,可识别几十个图片]...
1.1 应用MATLAB识别条形码总体设计 本设计为基于MATLAB的条形码识别系统,带有一个GUI界面. GUI设计界面 1.1.1 程序总体设计思路 在上节中我们简单的介绍了MATLAB的发展以及 ...
- matlab两轮自平衡小车,基于MATLAB的两轮自平衡小车系统模型辨识.pdf
基于MATLAB的两轮自平衡小车系统模型辨识 学兔兔 第1期 (总第170期) 机 械 工程 与 自动 化 NO.1 2012年 2月 MECHANICAL ENGINEERING & AUT ...
- matlab正交gold码的相关性,基于Matlab的CDMA通信系统仿真
基于Matlab的CDMA通信系统仿真 1 绪 论 1.1课题背景及目的 20世纪60年代以来,随着民用通信事业的发展,频带拥挤问题日益突出.CDMA(Code Diveision Multiple ...
最新文章
- 迷宫出路代码_如何在软件开发的迷宫中找到自己的出路
- 【HAOI2010】订货
- 2017-9-15Opencv 杂
- HDU 4763 Theme Section(KMP+枚举公共前后缀)
- SpringBoot快速开发利器:Spring Boot CLI
- 操作系统角度谈测试管理和自动化测试
- ubuntu+php环境下的Memcached 安装方法
- JBoss模块示例–模块化Web应用程序
- js定时器和linux命令locate
- LeetCode 1291. 顺次数(模拟)
- 50 个 Bootstrap 插件
- 家卫士扫地机器人好吗_扫地机器人哪个牌子好?市场最全智能扫地机器人品牌解析_扫地机器人...
- MySQL5.7的多源复制
- 20145120 《Java程序设计》实验一实验报告
- 浅谈AlphaGo背后所涉及的深度学习技术
- 如何更换您的Gmail地址
- Python爬虫知网
- Android音频可视化
- 【数据结构】栈的实现与简单应用
- 9、图解门禁系统电锁分类及安装方法