生成X，Y，Z与meshgrid坐标并保存到.csv

Question

我正在尝试创建一个均匀间隔的3D PointCloud/网格，并将其导出到CSV。

我设法找到了解决此问题的解决方案（请参阅链接和代码bellow），但可悲的是没有解释，我正在挠头，试图理解循环3的原因。

资料来源：我如何导出x，y，y，y，y，y，y，y数组？

import numpy as np
import csv

sp=(30)
x=np.arange(313000, 313120, sp)
y=np.arange(6220000,6220120, sp)
z=np.repeat(15,4)
x_mesh, y_mesh, z_mesh=np.meshgrid(x,y,z)

coords = []
for a, b, c in  zip(x_mesh, y_mesh, z_mesh):
    for a1, b1, c1 in zip(a, b, c):
        for a2, b2, c2 in zip(a1, b1, c1):
            coords.append((a2, b2, c2,))
print(coords)

with open('coords.csv','w') as f:
    out = csv.writer(f, delimiter=',')
    out.writerows(zip(*coords))
    f.close()

Answer 1

您有3个循环，因为您使用的是3个轴。

Nunpy Meshgrid创建ND坐标数组。例如，简单的2D数组在每个方向上有2个点：

x = ['x0', 'x1']
y = ['y0', 'y1']

x_mesh, y_mesh = np.meshgrid(x, y)

print('x')
print(x_mesh)
print('y')
print(y_mesh)
print('---')

将打印出：

x
[['x0' 'x1']
 ['x0' 'x1']]
y
[['y0' 'y0']
 ['y1' 'y1']]
---

因此，要获得单个坐标（例如'x0，y0'），您需要将它们加入其中。由于是2D，您只需要2个循环。

coords = []

for a, b in  zip(x_mesh, y_mesh):
    for a1, b1 in zip(a, b):
        coords.append((a1, b1))

print('coords')
print(coords)
print('---')

会打印出来：

coords
[('x0', 'y0'), ('x1', 'y0'), ('x0', 'y1'), ('x1', 'y1')]
---

您有3个维度，因此Meshgrid具有额外的嵌套数组。因此3循环。