如何从三维灌溉数据中分割二维网格？

Question

我有一些三维域点形式的土壤含水量（Theta）数据（列的CSV文件x, y, z, Theta）。我想从3D域的某些特定位置获取横截面（图中的截面ABCD）。我想计算横截面上5*5网格中θ的值，但是网格中每个节点周围的点与未知点不共面。我以前在python中对2D域做过这样的工作，但是3D域对我来说似乎更复杂。我发现了plotly can make something like that in its virtual environment，但我希望它输出一个numpy数组或pandas DataFrame，以在jupyter笔记本中绘制出一个轮廓。

我知道寻找网格需要通过插值或网格化从相邻点中找出图形中每个点的值，比如P0，然后使用matplotlib绘制横截面，但是我不知道怎么做。在

相关的问题是，matplotlib或类似的库中是否提供从三维网格切片二维网格的功能？在

谢谢你的帮助。在

Answer 1

你试过和侏儒玩吗？它是一个专门用于处理非结构化网格的库，所以听起来可能对您有一些用处。看看这些example notebooks。

Answer 2

根本问题是三维插值。有很多包可以完成这种类型的工作，或者您可以自己编写（使用，例如KDE，这基本上只是一种平滑/分块的类型）。关于这个话题有很多材料，比如

• 这个答案https://stackoverflow.com/a/15753011/230468

• scipy docs

• 这个扩展的set of option on scicomp.stack

• 和this blog post (with some good examples)