Radial averaging of 2-d tif image

조회 수: 22 (최근 30일)
Kevin Wang
Kevin Wang 2016년 2월 4일
답변: 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!

채택된 답변

Image Analyst
Image Analyst 2016년 2월 4일
편집: Image Analyst 2019년 9월 12일
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.

추가 답변 (3개)

Image Analyst
Image Analyst 2016년 2월 4일
편집: Image Analyst 2019년 9월 12일
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
Sergey Loginov
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

댓글을 달려면 로그인하십시오.


Hugo Trentesaux
Hugo Trentesaux 2019년 2월 7일
편집: Hugo Trentesaux 2019년 9월 11일
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
Hugo Trentesaux
Hugo Trentesaux 2019년 9월 11일
You're right.
Sergey Loginov
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

댓글을 달려면 로그인하십시오.


Sergey Loginov
Sergey Loginov 2021년 11월 5일
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.

카테고리

Help CenterFile Exchange에서 Biomedical Imaging에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by