如何命令梯度的重量
2视图(30天)
显示旧的评论
你好,我83帧图像dicom。我试着做梯度分割。
但是为什么只有一帧出现(如附)?
我假设函数imshow3D将显示所有83帧。
任何人都可以帮我吗?
%获得文件夹中所有文件的列表和所需的文件名称模式。
myFolder = (“C: \用户下载\阿克毛\ \智商NEMA GA68 51 (1) 2020 - 12 - __studies \ \智商NEMA GA68 5.1 _iq NEMA GA68 5.1 _pt_2020-12-02_114711_spine_pet.pelvis.dyn.q600_n83__00000);
filePattern = fullfile (myFolder,‘* .dcm‘);你需要什么模式变化百分比。
theFiles = dir (filePattern);
为L = 1:长度(theFiles)
baseFileName = theFiles (L) . name;
fullFileName = fullfile (theFiles (L)。文件夹,baseFileName);
流(1,“现在阅读% s \ n”,fullFileName);
%现在做任何你想做的这个文件的名字,
%,如阅读它作为图像数组与imread ()
RZ (:,: L) = dicomread (fullFileName);
结束
σ= 1.5;
W = gradientweight (RZ (:,: L),σ,“RolloffFactor”,0.01,“WeightCutoff”1);
R = 149;C = 138;
持有在;
情节(C R“r”。,“线宽”,1.5,“MarkerSize”15);
标题(与种子位置的原始图像)
打= 0.009;
(BW, D) = imsegfmm (W C R,打);
图中,imshow3D (BW)
标题(“分割图像”)
持有在;
情节(C R“r”。,“线宽”,1.5,“MarkerSize”15);
0评论
答案(0)
翻译的
