MATLAB的答案

谁知道我的T2 *的计算是错误的吗?

2视图(30天)
威廉
威廉 2023年3月23日
评论道: 威廉2023年4月4日
%变量定义:
面具= niftiread (“Gated_Image_TE1_1.nii.gz”);
% T2 *计算,让我们看看如果这工作
fitfunc = cftool;
%回波时间矩阵
x = (0.086883, 0.15, 0.3, 0.7);
% =图像矩阵;
图像= 0 (128128128 4);
图像(:,::1)= niftiread (“Original_Gated_Image_TE1_1.nii.gz”);
图像(:,::2)= niftiread (“Original_Gated_Image_TE2_1.nii.gz”);
图像(:,::3)= niftiread (“Original_Gated_Image_TE3_1.nii.gz”);
图像(:,::4)= niftiread (“Original_Gated_Image_TE4_1.nii.gz”);
图片/ max =图像(图像(:));
% %
抽搐
= 0(大小(面具));
B = 0(大小(面具));
C = 0(大小(面具));
D = 0(大小(面具));
某人= 0;
% f0 = fitoptions(“起源”,[0]);
为i = 1:128
为j = 1:128
parfork = 1:128
如果面具(i, j, k) = = 1
试一试
y =挤压(图像(i, j, k,:));
% myfit =适合(x, y,“exp2”、“曾经繁荣”,(0,0,0,0));
myfit =适合(x, y,“* exp (- x / b) + c”,曾经繁荣的(1、0.2、0.001),“UpperLimit”(10000 2 . 05),“LowerLimit”,(0,0,无穷));
(i, j, k) = myfit.a;
B (i, j, k) = myfit.b;
C (i, j, k) = myfit.c;
% D (i, j, k) = myfit.d;
抓
% B (i, j, k) = -10;
% disp(东西打破的);
某人= + 1;
结束
结束
结束
结束
结束
toc
所以当我尝试运行面具(Gated_ImageTE1_1.nii.gz) /美金汇扫描(Original_Gated_Image_TE1_1.nii.gz)没有回来。当我unkill捕捉功能整个最终拥有一个价值-10的面具。所以我不确定发生了什么和为什么。如果任何人都可以让我知道发生了什么,将不胜感激!提前谢谢。

在回答这个问题。

接受的答案

催化
催化 2023年3月23日
编辑:催化 2023年3月23日
试的问题……抓is that if you have 任何 在try块编码错误,它永远不会得到执行。在这里,您使用的参数名称“UpperLimit”和“LowerLimit”,不存在。
最好将请求一个exitflag从健康,而不是使用try / catch -
[myfit、~、输出]=健康(___)
如果output.exitflag < = 0
继续
结束
1评论
威廉
威廉 2023年4月4日
这帮助!对不起,我花了一段时间送还给你。忘了回复和马克帮助

登录置评。

更多的答案(0)

类别

找到更多的在图像算法在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!