你可以认为这是一个完美的圆吗?
选择特定区域的DICOM文件
5视图(30天)
显示旧的评论
我想选择的区域在两个圆之间,内内圆和外外的一切都不感兴趣的。这就是我试图“激活”和“禁用”我想要的像素。通过这两次,我能选择的一部分,左边和右边的一部分,然后把它们加起来。我想完成与圆之间的整个区域像素“上”。甚至这不是我应该做的,但如果有人推荐我很乐意听到!(不能上传dicom文件我将在zip压缩)
同时,我将选择区域的绝对值所以不担心它改变颜色从第一个到最后一个。
清晰;clc;clf;关闭所有;
我(::1)= dicomread (“MRIm1.dcm”);
P =我(:,:1);
图
次要情节(2、2、1),imshow (P, [])
% %
Pb = imbinarize (P, graythresh (P));
se = strel(‘行’,4,4);
Pb = imerode (Pb, se);
道具= regionprops(铅、“地区”、“PixelIdxList”);
[m,指数]= max ([props.Area]);
Pbm = 0(大小(P, 1),大小(P, 2));
Pbm(道具(指数).PixelIdxList) = 1;
P ~ (Pbm) = 0;
%相同但绝对值,颠倒和右侧
=我(:,:1);
= abs ();
Ab = imbinarize (A, graythresh (A));
% Ab = bwareaopen (Ab, 250);
% Ab = imfill (Ab,“漏洞”);
se = strel(‘行’,4,4);
Ab = imerode (Ab、se);
props1 = regionprops (Ab,“区域”,“PixelIdxList”);
[n, index1] = max ([props1.Area]);
反弹道导弹= 0(大小(A, 1),大小(A, 2));
反弹道导弹(props1 (index1) .PixelIdxList) = 1;
(~ Abm) = 0;
C = P +;%增加他们
次要情节(2 2 2),imshow (P, [])
次要情节(2,2,3),imshow (A, [])
次要情节(2,2,4),imshow (C, [])
我们的目标是完成右下角。
接受的答案
prasanth年代
2022年12月14日
使用中值滤波去除噪声
我= medfilt2 (P);
找到梯度
[Gmag, Gdir] = imgradient (IM);
然后应用圆检测或任何其他methds分离循环
4评论
里克
2022年12月15日
逻辑索引:
%加载装入的例子形象
负载核磁共振成像
我=挤压(D(:,: 1,装天花板(结束/ 2)));
我= im2double (im);
次要情节(1、3、1)
imshow (im)
标题(原始图像的)
(X, Y) = ndgrid (linspace(1, 1,尺寸(im, 1)), linspace(1, 1,尺寸(im, 2)));
R =√x ^ 2 + y ^ 2);
CircleMask = R > 0.6 & R < 0.8;
次要情节(1、3、2)
imshow (CircleMask)
标题(“面具”)
NewImage = 0.5 * 1(大小(im));
%与逻辑索引复制数据
NewImage (CircleMask) = im (CircleMask);
次要情节(1、3、3)
imshow (NewImage)
标题(“图像与旧值在面具”)
