pore size distribution of a 1x100 surface

2 次查看(过去 30 天)
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!
  4 个评论
Matt J
Matt J 2023-7-4
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).

请先登录,再进行评论。

采纳的回答

Image Analyst
Image Analyst 2023-7-5
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.

更多回答(1 个)

Khushi Yadav
Khushi Yadav 2023-7-4
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.
Hope this addresses your question!
  3 个评论
Image Analyst
Image Analyst 2023-7-4
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
Adriana 2023-7-4
编辑:Adriana 2023-7-4
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!

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 2-D and 3-D Plots 的更多信息

产品


版本

R2019a

Community Treasure Hunt

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

Start Hunting!

Translated by