海岸线维数计算.rar

%%作者 郑州大学2019级测绘科学与技术 %%核心公式来源：海岸线分维数计算方法及其比较研究 %%公式：lgN（ε）=-Dlgε+A(ε为len，N（ε）为N（len）) clc clear %% 将必要信息加载到环境路径中 add path in p = fileparts(mfilename('fullpath')); addpath(fullfile(p,'Algorithms')); addpath(fullfile(p,'Images')); %% 加载图像 image = double(imread('allgrd2.tif'));%中国海岸线.tif、allgrd2.tif %% 选用的网格大小 n=20; len=zeros(1,n); N=zeros(1,n);%预定义格网对应长度 parfor i=1:1:n len(i)=i;%设置网格大小 end %% 主函数 tStart = tic; load('kernel.mat'); edge=extraedge(image,kernel); tCls1 = toc(tStart); disp(['计算所用时间',sprintf('time= %.1f%',tCls1)]) for i=1:n N(i)=gridlength(edge,len(i));%得到每种网格边界长度 end tCls2 = toc(tStart); disp(['计算所用时间',sprintf('time= %.1f%',tCls2)]) %方程拟合 lenS=zeros(1,n);NS=zeros(1,n); parfor i=1:n lenS(i)=log10(len(i)); NS(i)=log10(N(i)); end myfunc=inline('beta(2)-(beta(1).*lenS)','beta','lenS'); beta=nlinfit(len,NS,myfunc,[0 0]); D=beta(1);A=beta(2); tCls3 = toc(tStart); disp(['计算所用时间',sprintf('time= %.1f%',tCls3)]) figure;%结果显示 subplot(1,1,1);imagesc(uint8(image));title(sprintf('此图的海岸线分维数为%.2f%', D))