网格人为干扰度计算方法的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)