multifractal-compute 用MATLAB编写的图象多重分形程序

用MATLAB编写的图象多重分形程序,可直接计算得到多重分形谱

%%
% Author: Tegy J. Vadakkan
% Date: 09/08/2009 
% input a binary image
% the multifractal spectra is calculated based on the ideas in the paper by
% Posadas et al. Soil Sci. Soc. Am. J. 67:1361-1369 2003

clc;
clear all;
close all;
%%
indata=inputdlg（{‘input filename‘}）;
a = imread（indata{1}）;
[rows cols] = size（a）;
figure;imshow（a）;
npix = sum（sum（a））;
%% calculates niL which is the number of pixels in the ith box of size L 
% ideas from boxcount.m by F. Moisy have been borrowed here
width = rows;
p = log（width）/log（2）;
max_boxes = power（rows2）/power（22）;
nL = double（zeros（max_boxesp））;
for g=（p-1）:-1:0
            siz = 2^（p-g）;
            sizm1 = siz - 1;
            index = log2（siz）;
            count = 0;
            for i=1:siz:（width-siz+1）
    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3175  2009-09-08 19:50  multifractal.m

----------- ---------  ---------- -----  ----

                 3175                    1