Hi, I have curve fitted an histogram using the gaussian function. Now I would like to know how to extract a series of cell diameters from the gaussian fit.

Question

Akshay Kumar Pakala 2021-10-28

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1574093-hi-i-have-curve-fitted-an-histogram-using-the-gaussian-function-now-i-would-like-to-know-how-to-ex

评论： yanqi liu 2021-11-3

采纳的回答： yanqi liu

Histogram with curve fit (x = cell diamters y = frequency)

3 个评论
显示 1更早的评论隐藏 1更早的评论

Akshay Kumar Pakala 2021-10-28

In my histogram the x values are cell diamters

Akshay Kumar Pakala 2021-10-29

please let me know if you want me to add more comments to the code.

Also, I would like to know, how to create random 3d points based on my gaussian fit?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

yanqi liu 2021-11-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1574093-hi-i-have-curve-fitted-an-histogram-using-the-gaussian-function-now-i-would-like-to-know-how-to-ex#answer_821083

编辑：yanqi liu 2021-11-2

在 MATLAB Online 中打开

sir，please check the follow code to get some information

for more information,you can read the book 《计算机视觉与深度学习实战》，Thank you！

clc; clear all; close all;

BW = imread('https://www.mathworks.com/matlabcentral/answers/uploaded_files/783383/image.png');

[nr,nc,layers] = size(BW);

%% converting the image to 2d so the image can be labledd %%

if layers>1

BW = BW(:,:,1);

end

% end of lablling %%%%%%%%%%%%%%%%%

binaryimage = BW < 128;

bw = logical(binaryimage);

imwrite(bw, 'imagetoapp.jpg');

%% lablelling the image %%%%%%

[labeledimage, numberofblobs] = bwlabel(binaryimage,8);

props = regionprops(labeledimage, 'Equivdiameter');

x = numel(props);

y = zeros(x,1);

for i = 1:1:length(props)

y(i,1) = struct2array(props(i,1));

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%% LAbled image to rgb

colouredimage = label2rgb(labeledimage,'jet', [0,0,0.5]);

coloredLabelsImage = label2rgb (labeledimage, 'jet', 'k', 'shuffle');

imshowpair(labeledimage,colouredimage,'montage')

hold on

%% generating histogram and curve fitting it using gaussian function%%%%%%%%%

figure()

y1 = histogram(y,45);

values = y1.Values;

a=y1.BinEdges;

for i = 1 : length(a)-1

b(i) = mean(a(i:i+1));

end

[barheights, position] = hist(y,45);

positiontranspose = position';

valuestranspose = values';

gauss = 'd + a*exp(-(1/2)*((x-b)/c)^2)'

startpoints = [40 20 10 0.5];

myfit = fit(positiontranspose,valuestranspose,gauss,'start',startpoints);

hold on

plot(myfit)

z(:,1) = myfit(x);

zt=myfit(b);

%% plot properties

xlabel('cell diamters')

ylabel('frequency')

hold on; plot(b, zt, 'r*')

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%

5 个评论
显示 3更早的评论隐藏 3更早的评论

Akshay Kumar Pakala 2021-11-2

Hey Yanqui Liu

I did not understand what you meant by that. Can you please explain?

Thank you

yanqi liu 2021-11-3

such as two，make x----y

three，can as [x(1) x(2)]----y

请先登录，再进行评论。

Answer 2

yanqi liu 2021-10-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1574093-hi-i-have-curve-fitted-an-histogram-using-the-gaussian-function-now-i-would-like-to-know-how-to-ex#answer_819263

sir, may be upload the data mat file or the figure file to analysis

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Akshay Kumar Pakala 2021-10-29

please let me know if you want me to add more comments to the code.

Also, I would like to know, how to create random 3d points based on my gaussian fit?

请先登录，再进行评论。

Answer 3

Akshay Kumar Pakala 2021-10-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1574093-hi-i-have-curve-fitted-an-histogram-using-the-gaussian-function-now-i-would-like-to-know-how-to-ex#answer_819733

编辑：Walter Roberson 2021-11-1

在 MATLAB Online 中打开

clc
clear all
BW = imread('binary.PNG');
[nr,nc,layers] = size(BW);
%% converting the image to 2d so the image can be labledd %%
if layers>1
    BW = BW(:,:,1);
end
% end of lablling %%%%%%%%%%%%%%%%%
binaryimage = BW < 128;
bw = logical(binaryimage);
imwrite(bw, 'imagetoapp.jpg');
%% lablelling the image %%%%%%
[labeledimage, numberofblobs] = bwlabel(binaryimage,8);
props = regionprops(labeledimage, 'Equivdiameter');
x = numel(props);
y = zeros(x,1);
for i = 1:1:length(props)
    y(i,1) = struct2array(props(i,1));
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% LAbled image to rgb
colouredimage = label2rgb(labeledimage,'jet', [0,0,0.5]);
coloredLabelsImage = label2rgb (labeledimage, 'jet', 'k', 'shuffle');
imshowpair(labeledimage,colouredimage,'montage')
hold on
%% generating histogram and curve fitting it using gaussian function%%%%%%%%%
figure()
y1 = histogram(y,45);
values = y1.Values;
[barheights, position] = hist(y,45);
positiontranspose = position';
valuestranspose = values';
gauss = 'd +  a*exp(-(1/2)*((x-b)/c)^2)'
startpoints = [40 20 10 0.5];
myfit = fit(positiontranspose,valuestranspose,gauss,'start',startpoints);
hold on
plot(myfit)
z(:,1) = myfit(x);
%%  plot properties
xlabel('cell diamters')
ylabel('frequency')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Hi, I have curve fitted an histogram using the gaussian function. Now I would like to know how to extract a series of cell diameters from the gaussian fit.

3 个评论
显示 1更早的评论隐藏 1更早的评论

采纳的回答

5 个评论
显示 3更早的评论隐藏 3更早的评论

更多回答（2 个）

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

Hi, I have curve fitted an histogram using the gaussian function. Now I would like to know how to extract a series of cell diameters from the gaussian fit.

3 个评论 显示 1更早的评论隐藏 1更早的评论

采纳的回答

5 个评论 显示 3更早的评论隐藏 3更早的评论

更多回答（2 个）

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

3 个评论
显示 1更早的评论隐藏 1更早的评论

5 个评论
显示 3更早的评论隐藏 3更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

0 个评论
显示 -2更早的评论隐藏 -2更早的评论