选择椭圆作为感兴趣区域（ROI）

Question

我使用imellipse选择椭圆作为我感兴趣的区域（ROI）。问题是我要选择的椭圆大约是45度，当我使用imellipse时，它似乎是水平或垂直的90度。

如何更改椭圆的方向？

感谢。

Answer 1

您需要旋转椭圆的坐标。像这样：

npts = 1e4;
t = linspace(0,2*pi,npts);
theta = pi/4;
aspect = [5 1]; % [x y]
x = aspect(1)*sin(t+theta);
y = aspect(2)*cos(t);
plot(x, y);

---

如果要使用imellipse在图像上绘制椭圆，可以提取顶点并对其进行转换：

figure, imshow('pout.tif');
h = imellipse;
exy = h.getVertices
theta = pi/12;
M = [cos(theta), sin(theta); -sin(theta), cos(theta)]
exy_centered = bsxfun(@minus,exy,mean(exy))
exyRot = bsxfun(@plus,exy_centered*M,mean(exy));
hold on
plot(exyRot(:,1),exyRot(:,2),'r') % orig: plot(exy(:,1),exy(:,2),'r')

要填写椭圆，创建蒙版，请使用roifill或roipoly：

w=getfield(imfinfo('pout.tif'),'Width');
h=getfield(imfinfo('pout.tif'),'Height');
bw = roipoly(zeros(h,w),exyRot(:,1),exyRot(:,2));