Droplet size calculation of a spray image

Question

Jason 2023-2-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1909550-droplet-size-calculation-of-a-spray-image

回答： Image Analyst 2023-2-12

particle.jpg

Hi predecessors,

I have this kind of spray image and I want to calculate the droplet size and distribution of the particles, and wondering if there is a way to do that on MATLAB so that I can learn from it. Thanks!

Here are what I do now, I really appreciate you can help to point out the errors or mistakes,and if you are me, what code will you use to realize this function.

clc
clear all;
image=imread('E:\picture_test\heitu.tiff');  
jqimage=image(2000:2500,2000:2500);  
imshow(jqimage,[0 4000]);title('jiequ');
[gradthresh,numofiter]=imdiffuseest(jqimage);
diffusedImage=double(imdiffusefilt(jqimage,'GradientThresh',gradthresh,'NumberOfIterations',numofiter));
imshow(diffusedImage,[0 4000]);
[testx,testy]=size(diffusedImage);
array_x=zeros(501,501);  
blank_x=zeros(501,501);
y=zeros(1,501);   
for i=1:testx     
    testx_temp=diffusedImage(i,:); 
       for i1=2:testy-1      
           y(i1)=abs((testx_temp(i1+1)-testx_temp(i1-1))/2);   
       end    
    array_x(i,:)=y;  
    [pks,locs]=findpeaks(array_x(i,:),"MinPeakHeight",40);
    blank_x_i=blank_x(i,:);  
    for i2=1:length(locs);
        blank_x_i(locs(i2))=1; 
    end
    blank_x(i,:)= blank_x_i; 
end
array_y=zeros(501,501);   
blank_y=zeros(501,501);
x=zeros(501,1);   
for i=1:testy     
    testx_temp=diffusedImage(:,i);  
       for i1=2:testx-1      
           x(i1)=abs((testx_temp(i1+1)-testx_temp(i1-1))/2);   
       end    
    array_y(:,i)=x;  
 
    [pks,locs]=findpeaks(array_y(:,i),"MinPeakHeight",40);
    blank_y_i=blank_y(:,i);   
    for i2=1:length(locs);
        blank_y_i(locs(i2))=1;   
    end
    blank_y(:,i)= blank_y_i;  
end
figure;
imagesc(blank_x);
hold on;
imagesc(blank_y);
add_ab=blank_x+blank_y;
imagesc(add_ab);
figure;
binaryimage=im2bw(add_ab);
imshow(binaryimage);
binaryimage=imfill(binaryimage,'holes');
imshow(binaryimage);title('binaryimage');
figure
labeledimage=bwlabel(binaryimage,8);
imshow(labeledimage);
coloredlabels=label2rgb(labeledimage,'hsv','k','shuffle');
hold on;
imshow(coloredlabels);
blobmeasurements=regionprops(labeledimage,add_ab,'all');
numberofblobs=size(blobmeasurements,1); 
boundaries=bwboundaries(binaryimage);   
numberofboundaries=size(boundaries,1);
for k1=1:numberofboundaries
    thisboundary=boundaries{k1};
    plot(thisboundary(:,2),thisboundary(:,1),'g','LineWidth',2);
end 
for k2=1:numberofblobs
    blobarea=blobmeasurements(k2).Area;
    blobecd(k2)=sqrt(4*blobarea/pi);
      fprintf('%d\n',blobecd(k2))
end

2 个评论
显示无隐藏无

Image Analyst 2023-2-10

在 MATLAB Online 中打开

This

jqimage=image(2000:2500,2000:2500);

does not work with the particle.jpg image you attached. It's not big enough.

Jason 2023-2-12

Thanks professor, I know this problem, now it is only for testing.Could you please help write a program that can implement this function?it may be better to learn.Thanks a lot

请先登录，再进行评论。

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

Answer 1

prasanth s 2023-2-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1909550-droplet-size-calculation-of-a-spray-image#answer_1169040

在 MATLAB Online 中打开

using 'kmeans' for segmentation and 'regionprops' to get properties of particles,

image=imread('particle.jpg');  
jqimage=rgb2gray(image);  
[idx,C]=kmeans(jqimage(:),5);
c_image=reshape(idx,size(jqimage));
droplets=c_image==1;
stats=regionprops(droplets);
average_area=mean([stats.Area]);