pore size distribution of a 1x100 surface

Hello!

Quick question, I've created an aggregate of 1000 spheres with this code:

t = readmatrix('Sim 1.txt');
[xn,yn,zn]=sphere;
i=1;
while i<=1000
    for i=i
        x=t(i,1);
        y=t(i,2);
        z=t(i,3);
        r=t(i,4);
        
        xi=xn*r+x;
        yi=yn*r+y;
        zi=zn*r+z;
        
        p(i)=surf(xi,yi,zi);
        hold on
        i=i+1;
    end
end
hold off;

How can I now measure the pore size distribution of my aggregate (1x1000 surface)?

Thank you in advance!

[Matt J]

Your code does not run. See above.

[Adriana]

Dear Matt, that's because you need the attached file to run it

t = readmatrix('Sim 1.txt');

Error using readmatrix

Unable to find or open 'Sim 1.txt'. Check the path and filename or file permissions.

[xn,yn,zn]=sphere;

i=1;

while i<=1000

for i=i

x=t(i,1);

y=t(i,2);

z=t(i,3);

r=t(i,4);

xi=xn*r+x;

yi=yn*r+y;

zi=zn*r+z;

p(i)=surf(xi,yi,zi);

hold on

i=i+1;

end

end

hold off;

Now it should run.

[Matt J]

Yes, it runs now but you should run it for us so we can see what it produces (but I've done it for you now).

[Adriana]

Thank you Matt!

Answer 1

[Image Analyst]

So it seems that you are saying the spheres represent solid matter and the rest is air space or vacuum. And you want to know the pore size between the solid spheres. For that I'd start with bwdist but be aware that the air space outside the envelope of your sphere cluster will be one gigantic pore and it will be connected to all the air space (pores) inside the cluster so you really need to think about how you define pore.

Answer 2

[Khushi Yadav]

Hi @Adriana,

I understand that you want to read data from "Sim 1.txt" using "readmatrix" but getting error indicating that the file cannot be found or opened.

The code you provided ran successfully in my system. You can try providing the full path of the file to the "readmatrix" function, if your "Sim 1.txt" file is not located in the same directory as your MATLAB script.

For more information, you can refer to the following documentation section: If the file is not in the current folder or in a folder on the MATLAB path, then specify the full or relative path name in filename.

Hope this addresses your question!

[Adriana]

Dear Khushi,

Thank you for your answer. That's not really my question, I'm able to read the file and plot the surfaces (see attached figure). But I need to get the pore size distribution of the agreggate generated. Do you know how to do that?

Thank you!

[Image Analyst]

I don't see how it doesn't answer your question. I offered two possibilities

1. the "pores" are the spheres (void regions) in a solid matter, or
2. the "pores" are the empty space between the solid matter.

So it's got to be one or the other. Either the spheres are pores, or they're not.

[Adriana]

I'm sorry, I dont know why, but i dont seem to see your answer. I can only see te part about correcting my path to the file, not anything about the pores. My case is your possibility number 2. Can you help me?

Thank you in advance!