generation of 2D array of circular ring
조회 수: 7 (최근 30일)
이전 댓글 표시
채택된 답변
DGM
2023년 7월 8일
편집: DGM
2023년 7월 8일
In this case, I'm going to do an antialiased image instead of a binary image.
sz = [300 400]; % image size [y x]
c = [150 200]; % circle center [y x]
dmin = 175; % diameters
dmaj = 250;
maskinner = drawcircle(sz,c,dmin/2);
maskouter = drawcircle(sz,c,dmaj/2);
annmask = maskouter.*(1-maskinner);
imshow(annmask,'border','tight')
% draw smooth circle
function circ = drawcircle(sz,c,r)
xx = 1:sz(2);
yy = (1:sz(1)).';
circ = sqrt((xx-c(2)).^2 + (yy-c(1)).^2); % no quick tricks this time
circ = min(max((1+r-circ)/2,0),1);
end
댓글 수: 0
추가 답변 (1개)
Jayant
2023년 7월 8일
Here are the steps how you can write a function which takes D, d, x, y, width and height as its input and gives a 2D array as output.
- Create a ringArray of (width, height) dimesion.
- Traverse the ringArray. For each element (i, j), calculate the euclidean distance from the centre(x,y).
- If the distance>=d or <=D, the assign 1, else assign 0.
- The function returns the ringArray which is the required 2D array.
댓글 수: 1
DGM
2023년 7월 8일
편집: DGM
2023년 7월 8일
No loops necessary.
sz = [300 400]; % image size [y x]
c = [150 200]; % circle center [y x]
dmin = 175; % diameters
dmaj = 250;
xx = 1:sz(2);
yy = (1:sz(1)).';
rr = (xx-c(2)).^2 + (yy-c(1)).^2; % squared euclidean distance
annmask = (rr <= (dmaj/2)^2) & (rr >= (dmin/2)^2); % compare against square radii
imshow(annmask,'border','tight')
The bit with using squared distance simply avoids doing an expensive sqrt() on the entire array, so it's an easy way to make simple comparisons like this faster.
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!