网格人为干扰度计算方法的matlab代码实现

matlab计算网格人为干扰度方法
以土地利用类型数据为底图，利用ARCGIS创建渔网工具创建渔网，以每一个网格为评价单元，计算人为干扰度。并把人为干扰度赋值给各网格。

具体步骤详见：刘富强,吴涛,蒋国俊, 等.海岸线与海岸景观格局对人为干扰度的动态响应——以营口市南部海岸为例[J].生态学报,2017,37(22):7427-7437. DOI:10.5846/stxb201609211902

具体计算公式如下：

MATLAB文件：
叠加渔网后的土地利用类型数据；
利用土地利用类型裁剪好的渔网面数据
MATLAB代码如下：

clc;
clear;
MAP = shaperead('E:\60m\60mwanggemian1.shp');            %土地利用类型与渔网的叠加文件
RENWEI = shaperead('E:\60m\60mrenwenganraomian1.shp');   %渔网面文件
date.ID=[MAP.FID_c60_Fe];              %渔网FID
ID_RENWEI=[RENWEI.OBJECTID];          % 渔网FID
date.LC=[MAP.ZHEJ1995_I];           %土地利用类型数据
date.area=[MAP.Shape_Area];                  %每一个小块的面积
len1= length(ID_RENWEI);
len2= length(date.LC);
for i=1:len2                           %给斑块赋值t=date.LC(i);switch(t)case {41,42,43}date.coef(i)=1;             % 水域                           case {45,46}date.coef(i)=2;             % 滩涂case {21,22,23,24}date.coef(i)=3;             % 林地case {31,32,33}date.coef(i)=4;             % 草地case {111,112,113.121,122,123}date.coef(i)=5;             % 耕地case {52 53}date.coef(i)=6;             % 其它建筑用地case {51}date.coef(i)=7;% 城镇用地end
endfor j=1:len1sum = 0;                             % sum每一个评价单元的总面积n = 0;X=find( date.ID == ID_RENWEI(j));len3=length(X);for k=1:len3sum = sum + date.area(X(k));n=n+1;temp.coef(n) = date.coef(X(k));       %记录第J个评价单元内的各斑块的干扰指数temp.area(n) = date.area(X(k));       %记录第J个评价单元内的各斑块的面积大小endu = 0;                                     %u为干扰指数for m = 1:n   temp_HDI = (temp.area(m)/sum)*temp.coef(m);u = u + temp_HDI;endHDI(j,1) = u;                          %记下每个中心点的干扰指数temp.coef=[];temp.area=[]; temp_HDI=[];
end
HDI1=num2cell(HDI);
l4=length(HDI)
for m=1:l4
RENWEI(m).renwei=HDI(m,1);
end
shapewrite(RENWEI,'E:\60m\fuzhiyuwang.s

网格人为干扰度计算方法的matlab代码实现相关推荐

python怎么画地图空间分异图_基于地形梯度的赣南地区生态系统服务价值对人为干扰的空间响应...
南方山地丘陵是我国主体功能区划中"两屏三带"生态安全战略格局的重要组成部分, 是我国华南和西南地区的生态屏障带, 是长江流域与珠江流域的分水岭及源头区, 对于区域乃至全球的气候调节 ...
matlab判断文件是否损坏,检查 MATLAB 代码文件是否有问题
示例以下示例使用了 lengthofline.m,它是带有可改进代码的 MATLAB 示例文件.您可以在 matlabroot/help/techdoc/matlab_env/examples 中可 ...
基于蒙特卡诺和拉格朗日乘子法的电动车调度【有序、无序充放电】【Matlab代码】
目录 1 概述 2 蒙特卡洛模拟方法介绍 3 拉格朗日乘子法 4 规模化电动汽车充电负荷预测计算方法 5 Matlab代码实现 1 概述电动汽车EV(Electric Vehicle)具有清洁环保. ...
【故障诊断】基于粒子群和多元宇宙算法优化最大相关峭度解卷积(MCKD)实现信号处理附matlab代码
1 内容介绍 1.1研究背景及意义制造业的发展支撑着国家经济的稳定快速发展,是人民生活幸福和国家局势稳定的保障.发展制造业离不开科技创新.加强制造业创新.推动产业转型升级是提升我国国际地位.实现 ...
VMD分解，matlab代码，包络线，包络谱，中心频率，峭度值，能量熵，近似熵，包络熵，频谱图，希尔伯特变换，包含所有程序MATLAB代码，-西储大学数据集为例
目录 1.选取数据 2.VMD函数-matlab代码 3.采用matlab脚本导入数据并做VMD分解 4. VMD分解图 5.计算中心频率 6.画包络线 7. 画包络谱 8. 计算峭度值 9.计算能量 ...
时域特征值提取的MATLAB代码实现（均方根、峰值因子、脉冲因子、裕度因子、峭度因子、波形因子和偏度等）
时域计算代码相对简单,下列代码主要来自求助matlab时域分析程序-峰值.均方根.裕度.峭度_百度知道.不过原文中求"峭度因子"的方法不对,即"Kr = sum(y.^4 ...
块坐标下降法算法matlab代码,《最优化：建模、算法与理论》、《最优化计算方法》代码展示...
<最优化:建模.算法与理论>.<最优化计算方法>代码展示此系列页面为配套<最优化:建模.算法与理论>.<最优化计算方法>的代码展示, 包含完整的算法和 ...
多尺度结构元素形态学边缘检测算法的研究-含Matlab代码
目录一.引言二.数学形态学理论概述三.实验验证四.参考文献五.Matlab代码获取一.引言使用数字图像处理技术来解决计算机视觉.人工智能.生物遥感器视觉等领域所涉及到的图像问题时,最重要 ...
【信号分解】基于LMD算法和ELMD算法实现管道泄漏信号处理附matlab代码
1 内容介绍在科技水平相当发达的今天,互联网+.大数据慢慢渗透进人们的生活当中,但科技的进步不仅仅要体现在生活质量水平的提高.经济的快速发展,更应该体现在对社会资源的合理利用.自建国以来,我国管道 ...

网格人为干扰度计算方法的matlab代码实现

网格人为干扰度计算方法的matlab代码实现相关推荐

最新文章

热门文章