I tried to reproduce the given code, and I have a few suggestions that can improve the code and create 3D cube porous material. Here are a few of them:
1. Relative Pore Radii:
You have set fixed values for the pore radii. This approach works well if the pore radii are always smaller than the matrix dimension, but it could lead to problems if the pore radii exceed the matrix dimension. To make the code more flexible, set the pore radii relative to the matrix size. The code snippet below illustrates this:
dimension = 20; % Size of the 3D matrix
porosity = 0.02; % Porosity (e.g., 2%)
minPoreRadius = 0.05 * dimension; % Relative to the matrix size
maxPoreRadius = 0.1 * dimension; % Relative to the matrix size
2. Avoiding Complete Matrix Dilation:
When you use “imdilate” function to create pores, the entire “materialMatrix” can become filled with ones by the end of the loop, which is not ideal for visualization. Instead, you can store the positions of the pores separately and then visualize those. This approach avoids cluttering the matrix and provides a clearer visualization. You can achieve this with the following code snippet:
p = [];
for i = 1:numPores
pore = strel('sphere', poreRadii(i));
materialMatrix(positionsX(i), positionsY(i), positionsZ(i)) = 1;
materialMatrix = imdilate(materialMatrix, pore);
[cpx, cpy, cpz] = ind2sub(size(materialMatrix), find(materialMatrix));
cp = [cpx, cpy, cpz];
p = [p; cp]; % Store the pore locations
p = unique(p, 'rows'); % Remove duplicate points
materialMatrix = zeros(dimension, dimension, dimension); % Reset the matrix
end
The above modifications to the code will help you create a 3D cube porous material.
For more information on “imdilate” function, you may refer to the MathWorks documentation mentioned below:
I hope this is helpful!
