Radial averaging of 2-d tif image

Question

Kevin Wang 2016-2-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266546-radial-averaging-of-2-d-tif-image

回答： Sergey Loginov 2021-11-5

I am looking to compute a radial average of a 2-d tif microscopy image.

I want to start from a given pixel as the center, go progressively radially and record the intensities of all the pixels lying on the locus of the circle of that radius. Out of all these, I want to compute a graph of radius vs. Intensity.

Any ideas? I am looking to create a .m file to do this. Thank you!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Image Analyst 2016-2-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266546-radial-averaging-of-2-d-tif-image#answer_208527

编辑：Image Analyst 2019-9-12

average_radial_profile.m

See my attached demo. Feel free to adapt, like to change the center or whatever. It gives the average radial profile within the blobs but you can use the same concept to do it over the whole image. The whole-image version is shown in the Answer below this one. Of course you could always just use ginput() to identify the center, compute the radiii, and use accumaray. Or just use a double for loop. So, there are several ways.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Image Analyst 2016-2-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266546-radial-averaging-of-2-d-tif-image#answer_208541

编辑：Image Analyst 2019-9-12

average_radial_profile_2.m

OK, so maybe that demo was too hard for you to adapt. Here is an easier, more straightforward demo using two for loops. It's easy to understand (as compared to accumarray which would probably be faster). See attached script below the image.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Sergey Loginov 2021-11-3

Dear ImageAnalyst,

Thanks for mentioning accumarray function here! At image sizes greater than 1000x1000 it is much much faster than this circle drawing stuff!!

Regards

请先登录，再进行评论。

Answer 3

Hugo Trentesaux 2019-2-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266546-radial-averaging-of-2-d-tif-image#answer_359977

编辑：Hugo Trentesaux 2019-9-11

在 MATLAB Online 中打开

My version of the function :

function profile = radialAverage(IMG, cx, cy, w)
    % computes the radial average of the image IMG around the cx,cy point
    % w is the vector of radii starting from zero
    [a,b] = size(IMG);
    [Y, X] = meshgrid( (1:a)-cx, (1:b)-cy);
    R = sqrt(X.^2 + Y.^2);
    profile = [];
    for i = w % radius of the circle
        mask = (i-1<R & R<i+1); % smooth 1 px around the radius
        values = (1-abs(R(mask)-i)) .* double(IMG(mask)); % smooth based on distance to ring
        % values = IMG(mask); % without smooth
        profile(end+1) = mean( values(:) );
    end
end

3 个评论
显示 1更早的评论隐藏 1更早的评论

Hugo Trentesaux 2019-9-11

You're right.

Sergey Loginov 2021-11-3

Dear Hugo,

I have seen your are have answered quite a lot o these questions about the radial profiles with the same code. I have myself implemented it in a similar manner recently. However it is really much much faster to use accumarray function from MatLab!!

Regards

SL

请先登录，再进行评论。

Answer 4

Sergey Loginov 2021-11-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/266546-radial-averaging-of-2-d-tif-image#answer_824839

Using accumarray is indeed so much faster!

function [Tics,Average]=radial_profile(data,radial_step)

%main axii cpecified:

x=(1:size(data,2))-size(data,2)/2;

y=(1:size(data,1))-size(data,1)/2;

% coordinate grid:

[X,Y]=meshgrid(x,y);

% creating circular layers

Z_integer=round(abs(X+1i*Y)/radial_step)+1;

% % illustrating the principle:

% % figure;imagesc(Z_integer.*data)

% very fast MatLab calculations:

Tics=accumarray(Z_integer(:),abs(X(:)+1i*Y(:)),[],@mean);

Average=accumarray(Z_integer(:),data(:),[],@mean);

end

Sergey Loginov (2021). Very Fast Radial Profile (https://www.mathworks.com/matlabcentral/fileexchange/101480-very-fast-radial-profile), MATLAB Central File Exchange. Retrieved November 5, 2021.