Scipy/Numpy Reimann总高度

Question

我正在研究模拟物体在平面表面上的轨迹的可视化。目前，我已经提供的算法使用一个简单的轨迹函数（其中提供了速度和重力）以及Runge-Kutta积分来检查沿着曲线的n个点，其中速度变为0的点。我们正在打折任何大气相互作用。

我想要做的事情是从数字地形模型（光栅）中引入一个非平面曲面。我的想法是计算每个像素的黎曼和，并确定平面表面的偏移量是否等于或小于平面表面下方地形的偏移量。

使用numpy或scipy可以计算Reimann矩形的高度吗？相反，矩形的面积（中点很好）会工作，因为我知道宽度nd可以计算高度。

Answer 1

对于计算雷曼和可以看看numpy.cumsum（）。我不确定你是否可以用这种方法做一个曲面或只有一个数组。但是，您可以随时循环遍历地形的所有行，并将每行存储为二维数组。留下你所有的地形高度阵列。