How to calculate the volume of a 3D plot defined by a matrix

Question

Wissem-Eddine KHATLA 2021-12-16

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1612585-how-to-calculate-the-volume-of-a-3d-plot-defined-by-a-matrix

评论： Wissem-Eddine KHATLA 2021-12-16

Hi everyone,

So basically, I am trying to calculate a volume. More precisely, I have a height field that is defined in every point of a 2D plane and nammed "h_real" in my code. This quantity is a 3D matrix because each height field is computed at every timestamp so for each time frame corresponds a given height field.

I would like to estimate the evolution of the volume associated to the height's progression during the considered time period (which is the length of the 3D matrix). I have tried this script but the result obtained is not accurate...

Could someone help me to find a new strategy to do so please ?

Thank you in advance

 for i = 1:length(h_real)
    
    volume(i) = sum(h_real(:,:,i),'all')*(scale^2)*(1e6);
    plot(volume)
    ylabel('Volume (ml)')
    xlabel('Number of images')
    title('Evolution du volume au cours du temps');
    
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Matt J 2021-12-16

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1612585-how-to-calculate-the-volume-of-a-3d-plot-defined-by-a-matrix#answer_856680

编辑：Matt J 2021-12-16

在 MATLAB Online 中打开

It looks fine, though I would do it this way:

    volume=reshape( sum(h_real,[1,2]), 1,[])*(scale^2)*(1e6);
  
    
    plot(volume)
    ylabel('Volume (ml)')
    xlabel('Number of images')
    title('Evolution du volume au cours du temps');

3 个评论
显示 1更早的评论隐藏 1更早的评论

Matt J 2021-12-16

编辑：Matt J 2021-12-16

在 MATLAB Online 中打开

So you think that if this method doesn't give me a relevant result it's because there is something wrong with my initial data "h_real" ?

Yes.

Don't you think that using the "sum" function is a little bit wrong in a certain way ?

Not if your data is smooth. You might be able to do a bit better with trapz().

volume=reshape( trapz(trapz(h_real,1),2), 1,[])*(scale^2)*(1e6);

Wissem-Eddine KHATLA 2021-12-16

Thank you for your answers : it helps me to reconsider again my initial experimental results.

Best regards.

请先登录，再进行评论。

How to calculate the volume of a 3D plot defined by a matrix

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

3 个评论
显示 1更早的评论隐藏 1更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

How to calculate the volume of a 3D plot defined by a matrix

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

3 个评论 显示 1更早的评论隐藏 1更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论