How to use voxels of different sizes to measure the volume of a 3D object using Matlab?

3 次查看（过去 30 天）
M.S. Khan
M.S. Khan 2020-9-3
评论： M.S. Khan 2020-9-7
Hi Dear Community members,
i want to measure the volume of a cone using voxels. i am new in this field of voxels. Using the first link, a sphre of radius 8 has been voxelized. After voxelization, how can we count the total number of voxels. how to fix the size of voxels?
Using this method of author, how can we voxelize the cone to measure its volume. The attached file has xyz - coordiantes for the cone.
Thanks in advance for all cooperation and guidance.
https://www.mathworks.com/matlabcentral/fileexchange/50802-voxelplotter
https://www.mathworks.com/matlabcentral/answers/55633-how-to-drawing-3d-voxel-data?s_tid=answers_rc1-2_p2_MLT
% This scrypt illustrates the use of VoxelPlotter function to visualize
% voxel data stored in a 3d matrix
clear all
close all
clc
%Generating sinthetic input
gridesize=16;
R=8;
VoxelMat=zeros(gridesize,gridesize,gridesize);
for i=1:gridesize
for j=1:gridesize
for k=1:gridesize
if (i-gridesize/2)^2+(j-gridesize/2)^2+(k-gridesize/2)^2<R^2
VoxelMat(i,j,k)=1;
end
end
end
end
[vol_handle]=VoxelPlotter(VoxelMat,1);
view(3);
daspect([1,1,1]);
set(gca,'xlim',[0 gridesize], 'ylim',[0 gridesize], 'zlim',[0 gridesize]);

请先登录，再回答此问题。

采纳的回答

Bruno Luong
Bruno Luong 2020-9-3
编辑：Bruno Luong 2020-9-4
Vol estimate = 1.2692e+06
xyz = load('Cone_20000_points.txt');
[minxyz,maxxyz] = bounds(xyz,1);
N = 30*ones(1,3); % resolution, select with care
ijk = min(max(ceil((xyz-minxyz)./(maxxyz-minxyz).*N),1),N);
C = accumarray(ijk,1);
b = C > 0;
p = mean(C(b))
if p < 2 % Sanity check
warning('Your bin is too small, inaccurate result');
end
Vol = sum(b,'all') * prod((maxxyz-minxyz)./N)
  11 个评论
显示 9更早的评论
M.S. Khan
M.S. Khan 2020-9-5
Greetings Bruno, sorry for late reply. i was out of station.
Yes true, i got wrong results. Could you pleasue give me some tips how to count the number of voxels and how to measure the size of voxel.
M.S. Khan
M.S. Khan 2020-9-7
i could not still figure out, how to measure the total number of voxels used. Please anyone can help me?

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 3-D Volumetric Image Processing 的更多信息

标签

产品


版本

R2020a

Community Treasure Hunt

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

Start Hunting!

Translated by