Radial Distribution of particle centers

Hi everyone, I have the following aggregate structure of different colored particles. I want to compute the radial distribution of particle centers for the aggregate (agglomerate). This structured is formed with the help of voxel information from binary images obtained through X-ray tomography. Any suggestion would be highly appreciated, thanks in advance.

Answer 1

[Image Analyst]

Just create an array. Like if you have an array of xCenters and yCenters and zCenters

distances = sqrt(xCenters.^2 + yCenters .^ 2 + zCenters .^ 2);
radialDistribution = zeros(1, length(distances));
for k = 1 : length(xCenters)
    thisDistance = round(distances(k)); % Convert distance into an index.
    radialDistribution(thisDistance) = radialDistribution(thisDistance) + 1;
end
bar(radialDistribution, 'BarWidth', 1);
grid on;

[mejdi rogbist]

Thank you for this reply.

I had a position file.xyz so how can i do to load this data from file ?

also i try to run this code but i fond an error in "radialDistribution(thisDistance) = radialDistribution(thisDistance) + 1;" <-- Index exceeds matrix dimensions.

[Image Analyst]

Try

thisDistance = floor(distances(k)); % Convert distance into an index.

[mejdi rogbist]

Thank You. But the same error still again

[Image Analyst]

Attach an mat file with the centers vectors in it.

[mejdi rogbist]

I just try with five values of my data to test this code, after i will introduce my file (xyz extension).

[Image Analyst]

I will not be typing that all in, from that picture. Would you if you were me, when you can simply attach the m-file with the paper clip icon?

[mejdi rogbist]

Thank's a lot Image Analyst, you find attached the m file. I couldn't attach my data position (file.xyz). I will be happy to find a solution of my problem.