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

3 views (last 30 days)
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:amjadkhan_cs@yahoo.com "*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=omread('A001.jpg'); 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)
if true
% code
end

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!