python - 切片并找到体积

Question

我有四列，即 x,y,z,zcosmo 。 zcosmo的范围是0.0<zcosmo<0.5 。

对于每个 x,y,z，都有一个 zcosmo。

当 x,y,z 被绘制出来，这就是它们的样子。

我想求这个图形的体积。如果我将其切成 50 个部分(按 zcosmo 升序)，那么每个部分都类似于一个圆柱体，我可以将它们相加以获得最终体积。

切片圆柱体的体积为pi*r^2*h ，就我而言 r = z/2 & h = x

例如，切片就像， <强> x,z for 0.0<zcosmo<0.01 找到此卷V1。然后x,z for 0.01<zcosmo<0.02 找到这个卷V2，依此类推，直到zcosmo=0.5

我知道要手动执行此操作(这当然很耗时)，方法是:

r1 = z[np.logical_and(zcosmo>0.0,zcosmo<0.01)] / 2 #gives me z within the range 0.0<zcosmo<0.01
h1 = x[np.logical_and(zcosmo>0.0,zcosmo<0.01)] #gives me x within the range 0.0<zcosmo<0.01

V1 = math.pi*(r1**2)*(h1)

这里 r1 和 h1 应该是 r1 = ( min(z) + max(z) ) / 2.0和h1 = max(x) - min(x) ，即最大值和最小值，以便我为每个切片获得一个体积

我应该如何创建一个代码来计算 zcosmo 切片范围内的 50 个体积切片？

Answer 1

使用for循环:

volumes = list()
for index in range(0, 50):
    r = z[np.logical_and(zcosmo>index * 0.01, zcosmo<index * 0.01 + 0.01)] / 2
    h = x[np.logical_and(zcosmo>index * 0.01, zcosmo<index * 0.01 + 0.01)]

    volumes.append(math.pi*(r**2)*(h))

最后，volumes 将是包含 50 个气缸的体积的列表。

您可以使用volume = sum(volumes)来获取形状的最终体积。