Do not plot the shape in the middle and polar transformation

1 次查看（过去 30 天）
Alex Perrakis
Alex Perrakis 2021-11-18
评论： Alex Perrakis 2021-11-18
Hello Guys, i have the above seen plot and i would like not to plot the shape seen within the circle and the do a polar tranformation. So that that circle appears to be a line.
  6 个评论
显示 4更早的评论
Alex Perrakis
Alex Perrakis 2021-11-18
I tried to attach em as tiff but it does not accept em. so i attached em as png

请先登录，再进行评论。

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

采纳的回答

Adam Danz
Adam Danz 2021-11-18
编辑：Adam Danz 2021-11-18
Ran in:
This demo uses a cleaned up version of your code and then does the following,
  1. Computes center of noisy circle
  2. Shifts noisy circle to (0,0)
  3. Converts shifted image to polar coordinates
  4. finds the spike in rho values (radii) based on thresholding
  5. Removes the spike.
If you want a smooth circle, you could fit the final (x,y) values to a cirlce (see example) or just use the median rho values in polar coordinates.
% Your code cleaned up
pic01=imread("1.png");
pic050=im2gray( imread("50.png") );
newpicture=imsubtract(pic01,pic050);
%imshow(newpicture)
% level = graythresh(newpicture);
level = 0.1405;
newbinpic = imbinarize(newpicture,level);
%imshowpair(newpicture,newbinpic,'montage')
newbinpic2=bwpropfilt(newbinpic,'perimeter',1);
% imshowpair(newbinpic2,newpicture,'montage');
boundaries = bwboundaries(newbinpic2,'noholes');
% Plot raw (x,y) values
x = boundaries{1}(:,2);
y = boundaries{1}(:,1);
figure()
subplot(2,2,1)
plot(x,y, 'r-')
title('raw values')
axis equal tight
hold on
% Assuming (x,y) is a noisy circle, find center
cnt = [min(x)+range(x)/2, min(y)+range(y)/2];
xline(cnt(1))
yline(cnt(2))
% Shift to (0,0) and compute polar coordinates
[theta, rho] = cart2pol(x-cnt(1), y-cnt(2));
% Plot in polar coord
subplot(2,2,2)
polarplot(theta, rho)
title('Polar coord.')
% Define a reasonable threshold to isolate the spike in RHO values
% 90% of the median looks OK
subplot(2,2,3)
plot(rho)
threshold = median(rho)*.9;
yline(threshold, '-k', 'Threshold')
ylabel('rho'); xlabel('index')
title('rho values')
% Find start and end of peak
peakIdx = find(rho < threshold);
peakIdx(2:end-1) = [];
xline(peakIdx(1)) % plot start of peak
xline(peakIdx(2)) % plot end of peak
% Remove the peak
theta(peakIdx(1):peakIdx(2)) = []; % you can do the same with x
rho(peakIdx(1):peakIdx(2)) = []; % you can do the same with y
% Plot the circle without the interior object
subplot(2,2,4)
polarplot(theta,rho)
title('Inner removed')
  4 个评论
显示 2更早的评论
Adam Danz
Adam Danz 2021-11-18
编辑：Adam Danz 2021-11-18
It looks like range is in the stats toolbox.
Here's an anonymous function you can include at the top of the code that will replace range()
range = @(x)max(x(:))-min(x(:));

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Weather and Atmospheric Science 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by