Rayleigh and Gaussian gaussian and rayleigh distribution MRI Late Enhancement for cardiac images

Hello

i found this code on Mathworks.com, when i run this code for my cardiac LV image then i face withe this type of problem. plz i need help here

my id:[e-mail]

Error using  .* 
Integers can only be combined with integers of the same class, or scalar doubles.

Error in linspace (line 30)
    y = d1 + (0:n1).*(d2 - d1)/n1;

Error in hist (line 94)
        edges = linspace(miny,maxy,x+1);

Error in dd (line 13)
[c1,c2]=hist(I(:)',50);"

clear;clear all;close all;
I=dicomread('IM_0001.dcm');
I=uint8(round(double (I)));
I=I(:)',[muhat, sigmahat]=normfit(I); 
mean_value=round(muhat); 
sigma_value=round(sigmahat);

gauss(i)=100*exp( - (x(i)-mean_value)^2/(2*sigma_value^2)/(sigma_value*sqrt(2*pi)));

%5end
%%x=0:1:29952; for i=1:size(x')

rayl(i)=500*exp( - (x(i)^2/(2*sigma_value^2)*(x(i)/sigma_value^2)));

  %%end
  %%figure(3),plot(x,gauss, x,rayl),hold on imhist(I)
  [c1,c2]=hist(I(:)',50);
  f=fit(c2', c1', 'gauss1'); 

gauss=f.a*exp( - (c2-f.b).^2/ (2*f.c^2) );

hist(I(:),50),hold on, plot(c2,gauss,c2,rayl)

Answer 1

[Walter Roberson]

Change the

hist(I(:),50)

to

hist(double(I(:)),50)

[A-J KHAN]

First i say thanks "Walter Roberson" after change the hist(double(I(:)),50) i face with this problem can you plz help me here.

Error using fit>iFit (line 353) NaN computed by model function, fitting cannot continue. Try using or tightening upper and lower bounds on coefficients.

Error in fit (line 108) [fitobj, goodness, output, convmsg] = iFit( xdatain, ydatain, fittypeobj, ...

Error in dd (line 14) f=fit(c2', c1', 'gauss1');

[Walter Roberson]

After you do the

I=dicomread('IM_0001.dcm');

please show

max(I(:))
min(I(:))

I suspect that your original data is not in the range 0 to 255, so your uint8() is probably badly distorting the shape of the data, resulting in data that cannot be numerically be fit by a single Gaussian.

In particular CT data is in the range -1600 to +30000

[A-J KHAN]

this images are MRI images, and i am new in matlab so i don't know to how to create fit.

[Walter Roberson]

After you do the

I=dicomread('IM_0001.dcm');

please show

max(I(:))
min(I(:))

[A-J KHAN]

dear sir i upload the code and image here you can check plz

clear;clear all;close all;
I=imread('A001.jpg');
max(I(:))
min(I(:))
hist(double(I(:)),50)

I=uint8(round(double (I)));  
I=I(:)',[muhat, sigmahat]=normfit(I); mean_value=round(muhat); sigma_value=round(sigmahat);
x=0:1:29952; for i=1:size(x')
gauss(i)=100*exp( - (x(i)-mean_value)^2/(2*sigma_value^2)/(sigma_value*sqrt(2*pi)));
end
x=0:1:29952; for i=1:size(x')
  rayl(i)=500*exp( - (x(i)^2/(2*sigma_value^2)*(x(i)/sigma_value^2)));
end
figure(3),plot(x,gauss,x,rayl),hold on imhist(I)

[Walter Roberson]

I cannot proceed on this until you post max(I(:)) and min(I(:)) right after reading the .dcm file. .jpg files are not the same as .dcm files.

By the way: you need a comma before the imhist(I) .

I also firmly recommend that you use semi-colon instead of comma after I=I(:)'

[A-J KHAN]

Dear sir Walter Roberson, i say thanks to u and plz send me u r email id i say personal thanks to you, its running already....i posted the code for other student if some one need.

clear;clear all;close all;
I=dicomread('IM_0001.dcm');
max(I(:),0);
min(I(:),1);
hist(double(I(:)),50);

I=uint8(round(double (I)));  
I=I(:)',[muhat, sigmahat]=normfit(I); mean_value=round(muhat); sigma_value=round(sigmahat);
x=0:1:29952; for i=1:size(x')
gauss(i)=100*exp( - (x(i)-mean_value)^2/(2*sigma_value^2)/(sigma_value*sqrt(2*pi)));
end
x=0:1:29952; for i=1:size(x')
  rayl(i)=500*exp( - (x(i)^2/(2*sigma_value^2)*(x(i)/sigma_value^2)));
end
figure(3),plot(x,gauss,x,rayl),hold on,imhist(I)

Answer 2

[A-J KHAN]

Dear sir i still in trouble because i am new in matlab