Login  Register

Re: Convert macro from NIH Image to ImageJ

Posted by dinurf on Jun 20, 2016; 10:00am
URL: http://imagej.273.s1.nabble.com/Convert-macro-from-NIH-Image-to-ImageJ-tp5015304p5016692.html

Dear Gabriel,

Thank you for your explanation. I am a bit confused since you explained
about dilation dimension. Are you referring to my first&previous question
about dilation macro? Because my recent&new question is about Euclidean
Distance Map macro.

But from your explanation about dilation dimension, I tried it as well. I
calculated the perimeter (length) first by dividing Area(epsilon)/epsilon
and also plotted the log radius vs log perimeter as you mentioned. I found
the fractal dimension for Koch Snowflake using that method was 1,2283,
since the slope is - 0,2283. That value is quite far from the theoretical
value, which is 1,2685.

I don't used an image with end points, I used the dilation macro/dilation
method for closed contours.

Thank you and thank you for the paper as well!

Regards,
Dini


On Mon, Jun 20, 2016 at 4:31 PM, Gabriel Landini <[hidden email]>
wrote:

> > I downloaded a macro from ImageJ website a long time ago, the macro is
> > called FractalEDM (attaced) to compute an object's fractal dimension.
> > Because I am new using ImageJ, I don't really understand what has written
> > within the macro. So if you don't mind, Could you help me please to let
> me
> > know what step by step within the macro or general step? I cannot contact
> > the main author, her email was no longer used.
> >
> > For example, for the first line, (requires "1.30p") and I see that they
> > used "run(Invert)" which is I don't understand why.
> >
> > Besides that, I am also trying to used menu Distance Map. The steps I
> used:
> >
> > 1. Make an image into binary
> > 2. Menu --> Process --> Binary --> Distance map
> > 3. Menu --> Analyze --> Histogram , to get the list of data of values and
> > counts.
>
> It is not the area of the Minkowski 'sausage', but its length, what you
> use to
> compute the 'dilation dimension'.
> The length is related to the diameter of the disc used for the dilation:
> Length(epsilon)=Area(epsilon)/epsilon, where epsilon is the diameter of the
> dilation disc,
>
> You can do that with the Maximum filer as you know exactly the radius, so
> you
> divide the area of the sausage at a given radius by
> (2*radius)+1 to get the length.
> Then plot log((2*radius)+1) vs. the log of the length of the 'sausage' for
> that radius. Then 1-slope should give you the fractal dimension.
>
> Be aware that this procedure only works for closed contours. If you are
> computing this for a tree or network with end points, you overestimate the
> length of the 'sausage' at each free end. So you need to compensate for
> that
> (which is not straightforward to do).
>
> See fir example Eins S. An improved dilation method for the measurement of
> fractal dimension. Acta Stereologica 1995;14(2):169-178.
> http://popups.ulg.ac.be/0351-580X/index.php?id=831&file=1&pid=825
>
> Maybe computing the box dimension would be easier as it does not suffer
> from
> that problem.
>
> Cheers
>
> Gabriel
>
> --
> ImageJ mailing list: http://imagej.nih.gov/ij/list.html
>

--
ImageJ mailing list: http://imagej.nih.gov/ij/list.html