%%作者?鄭州大學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（1n）;
N=zeros（1n）;%預定義格網對應長度
parfor?i=1:1:n
????len（i）=i;%設置網格大小
end
%%?主函數
tStart?=?tic;
load（‘kernel.mat‘）;
edge=extraedge（imagekernel）;
tCls1????=?toc（tStart）;
disp（[‘計算所用時間‘sprintf（‘time=?%.1f%‘tCls1）]）
for?i=1:n
????N（i）=gridlength（edgelen（i））;%得到每種網格邊界長度
end
tCls2???=?toc（tStart）;
disp（[‘計算所用時間‘sprintf（‘time=?%.1f%‘tCls2）]）
%方程擬合
lenS=zeros（1n）;NS=zeros（1n）;
parfor?i=1:n
????lenS（i）=log10（len（i））;
????NS（i）=log10（N（i））;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1015??2019-11-01?19:45??海岸線維數計算\Algorithms\extraedge.m

?????文件????????999??2019-11-01?20:30??海岸線維數計算\Algorithms\gridlength.m

?????文件????????185??2019-11-01?17:31??海岸線維數計算\Algorithms\kernel.mat

?????文件??????21492??2019-10-25?20:52??海岸線維數計算\Images\allgrd.tif

?????文件?????228372??2019-10-26?17:13??海岸線維數計算\Images\allgrd2.tif

?????文件????1558984??2019-11-01?18:52??海岸線維數計算\Images\中國海岸線.tif

?????文件??375573770??2019-11-01?18:46??海岸線維數計算\Images\中國海岸線da.tif

?????文件???????1233??2020-01-12?16:27??海岸線維數計算\main.m

?????文件??????68071??2019-11-01?21:37??海岸線維數計算\中國海岸線結果.png

?????文件??????64507??2019-11-01?21:19??海岸線維數計算\隨便一張圖片的結果.png

?????目錄??????????0??2019-11-01?19:49??海岸線維數計算\Algorithms

?????目錄??????????0??2019-11-01?21:27??海岸線維數計算\Images

?????目錄??????????0??2019-11-01?21:45??海岸線維數計算

-----------?---------??----------?-----??----

????????????377518628????????????????????13