A lot depends on the characteristics of your binary images, and what you want the result to look like. Here's something simple that you may be able to build on:
A = rand(50, 50, 50) < 0.01; % synthetic data
[x y z] = ind2sub(size(A), find(A));
plot3(x, y, z, 'k.');
