How to compare the shape of lipid layers

Hi,

I have these two grayscale images of lipids forming multiple layers. The images looks very different, however I would like to find a procedure to quantitatively compare these two images in terms of curvature and/or circularity of the different layers. Any suggestion is welcome.

Claudio

Answer 1

[Image Analyst]

Try thresholding at a decreasing set of values and getting the area fraction or perimeter length at each value. Something like

thresholds = 255 : -32 : 1;
for k = 1 : length(thresholds)
    t = thresholds(k)
    mask = grayImage > t;
    imshow(mask);
    drawnow;
    props = regionprops(mask, 'Area', 'Perimeter');
    allAreas = [props.Area];
    allPerims = [props.Perimeter];
    areaFraction = sum(allAreas) / numel(mask)
    totalPerimeter = sum(allPerims)
end

Write back if you need more help.