3

在此处输入图像描述

我有一张图像,我想围绕一个点裁剪感兴趣的圆形区域。我在 MATLAB 中做了以下操作:

  vessel=imread('vessel.jpg');
  imshow( vessel)
  t = 0:pi/20:2*pi;
  xc=230; % point around which I want to extract/crop image
  yc=79;
  r=20;   %Radium of circular region of interest
  xcc = r*cos(t)+xc;
   ycc =  r*sin(t)+yc;
   roimaskcc = poly2mask(double(xcc),double(ycc), size(vessel,1),size(vessel,2));
   pr_gccc = find(roimaskcc);
   roimean_cc= mean(vessel(pr_gccc));
  figure, imshow(roimaskcc)

roimaskcc 是正确的,但是当我执行以下操作时,它会给出 nX1 矩阵,而不是掩码下的感兴趣区域:

  vessel_undermask=vessel(roimaskcc==1);

任何人都可以请。帮助提取兴趣点(xc,yc)周围的圆形兴趣区域。谢谢

4

2 回答 2

4

你正在做的是从vessel面具外面提取所有东西。由于没有圆形矩阵这种东西,Matlab 的解决方案是输出掩码内所有值的向量

虽然这在技术上是正确的,但该向量可能很难使用。另一种解决方案是让数据矩阵保持正方形,并将掩码外的所有内容设置为NaN

% make a copy      
vessel_undermask = vessel;

% NaN everything outside the mask (in R, G and B layers)  
vessel_undermask(repmat(~roimaskcc,[1,1,3])) = NaN;

imshow(vessel_undermask)

这应该为您提供一个更易于使用的矩阵。

注意

vessel_undermask(~isnan(vessel_undermask)) ==  vessel(roimaskcc)
于 2013-04-15T13:18:19.853 回答
1

是的,我知道了。我做了这样的事情:

           vesseltry=vessel;
          vesseltry(~roimaskcc)=0;

现在,vesseltry 是我的新图像,带有圆形感兴趣区域...

于 2013-04-15T13:17:43.383 回答