Cross-sectional area of a cone(ish)

Question

jlt199 2016-8-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/301518-cross-sectional-area-of-a-cone-ish

回答： jlt199 2016-9-9

采纳的回答： John D'Errico

Morning all,

I have a surface that looks something like a lopsided squished cone, hollow in the middle.

I would like to take the cross-sectional area at different heights, but I need some help. I'm sure it shouldn't be too difficult but I can't figure it out today. I can use the countour function to plot the contours at different heights, but not find the area of the cross-section.

Can anyone help me please?

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

请先登录，再进行评论。

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

Answer 1

John D'Errico 2016-8-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/301518-cross-sectional-area-of-a-cone-ish#answer_233289

And polyarea won't suffice to compute the area of a polygon? It should.

Just generate the contours as desired. Then use polyarea. Easy.

6 个评论
显示 4更早的评论隐藏 4更早的评论

John D'Errico 2016-8-31

编辑：John D'Errico 2016-8-31

In that case, the contour is a noisy mess. You get what you get then. What did you expect? You asked for the area enclosed. If the surface is a noisy, bumpy mess, then the area enclosed will be difficult to measure, by ANY means. Computer algorithms are not magically gifted. They cannot read your mind, and the human eye is very good at seeing the signal the analyst wants to see, while removing the spurious dreck. (It would have helped if you stated originally this would be a noisy surface.)

1. Smooth the surface as necessary, BEFORE you compute the contours. This can introduce problems of course. Smoothing will change the surface. Too much smoothing may change it unacceptably. Too little and you get what you showed.

A simple scheme for smoothing a regular surface in 2 dimensions is to convolve the surface with a radially symmetric Gaussian blur kernel. conv2 will do this nicely, but you will need to choose the smoothing kernel. What will matter here is the standard deviation of the Gaussian.

Alternatively, you could use the equivalent of a Savitsky-Golay filter, which will implicitly fit a local 2-d surface at every point, then return only the predicted value of the surface at that point. Here the local polynomial order will matter, and the local span used for the fit. Again, this can be achieved using conv2, although I do not know of any implementations of the above scheme for 2-d smoothing. Not hard to write. :-P

I suppose you could also use my gridfit tool on the FEX, passing in the array as a list of (x,y,z) points in 3-dimensions. Gridfit allows you to control the smoothing parameter. Gridfit will be slower than other tools that are based on conv of course, so I'd first recommend one of the other schemes I mentioned.

There may be some 2-d smoothing tools on the FEX. If so, I cannot validate them for goodness of code. (Except for gridfit.)

Again, any smoothing scheme will do things like flatten and widen the peaks.

2. Break the contour polygons into their separate pieces, in case there are disjoint polygons. You will need to know which one is important to you then.

jlt199 2016-9-9

Thanks for your help John, I have just got back to this problem after having to leave it to write reports :(

I was as surprised by the level of noise as you, it's only in about 5 of about several hundred surfaces, so I missed it originally.

I'm having a hard time working out the structure of the matrix C which is output from the contour3 function, I have tried looking at the code for contour3, but really didn't understand what was going on. Can anyone help me with this? I would like to isolate polygons as suggested, but I can't understand the data structure.

Another thought I had was to just send the portion of the surface that contains the "horn" to the contour3 function, but I can't think how to do that either.

If anyone can help with either of these two problems I would be very grateful.

This is the surface I am currently working with:

Many thanks

jlt199 2016-9-9

Ok, I've figured out the structure of C now. Now I just need to figure out what to do with it...

请先登录，再进行评论。

Answer 2

Chad Greene 2016-8-31

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/301518-cross-sectional-area-of-a-cone-ish#answer_233285

编辑：Chad Greene 2016-8-31

在 MATLAB Online 中打开

I'd use my C2xyz function which is on File Exchange to easily get the x,y values of a contour line. Below I'm using the built-in peaks data as an example dataset and getting the area of the polygon bounded by the z=6 line.

[X,Y,Z] = peaks(1000); 
pcolor(X,Y,Z); 
shading interp
colorbar
zval = 6;
hold on
C = contour(X,Y,Z,zval*[1 1],'k');

[x,y,z] = C2xyz(C); 
A = polyarea(x{1},y{1})
A =
  0.6661

2 个评论
显示无隐藏无

Sean de Wolski 2016-8-31

How is C2xyz different than contourc?

Chad Greene 2016-9-9

Hey Sean, the C2xyz function simply converts the difficult-to-interpret C matrix into more intuitive x and y values.

请先登录，再进行评论。

Answer 3

jlt199 2016-9-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/301518-cross-sectional-area-of-a-cone-ish#answer_234316

Thanks, I've managed to get it working in most cases. Except where the noise gets horrendous.

Many thanks for your time

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

请先登录，再进行评论。

Cross-sectional area of a cone(ish)

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

采纳的回答

6 个评论
显示 4更早的评论隐藏 4更早的评论

更多回答（2 个）

2 个评论
显示无隐藏无

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

另请参阅

类别

标签

产品

Community Treasure Hunt

Cross-sectional area of a cone(ish)

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

采纳的回答

6 个评论 显示 4更早的评论隐藏 4更早的评论

更多回答（2 个）

2 个评论 显示 无隐藏 无

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

另请参阅

类别

标签

产品

Community Treasure Hunt

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

6 个评论
显示 4更早的评论隐藏 4更早的评论

2 个评论
显示无隐藏无

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