问题描述
我是否可以用 ndarray 替换 pyvista_ndarray 的坐标？
我的目标是重新投影 vtu 文件（非结构化网格）的点坐标。 mesh 对象的当前坐标位于坐标参考系统 epsg 27672 中，我希望它们位于 epsg 4326 (wgs84) 中。 为此，我使用 pyvista 模块打开 vtu 文件：
mesh = pv.read("arg_Tm2__t0002.pvtu")
type(mesh.points)
>pyvista.core.pyvista_ndarray.pyvista_ndarray
因此，mesh.points 给出了 3 个空间坐标。然后，我使用 pyproj 模块将 3 个坐标重新投影到 epsg 4326。通过组合 3 个结果 x,y,z numpy.ndarray，我现在得到一个 NumPy 数组类似于 mesh.points 的形状和大小。
# set the pyproj transformer
crs1 = CRS.from_epsg(27562)
crs2 = CRS.from_epsg(4326)
reproj = Transformer.from_crs(crs1,crs2)
# Reprojection
reproj_dataY,reproj_datax,altitude = reproj.transform(mesh.points[:,0],mesh.points[:,1],2])
reprojData = np.column_stack((reproj_datax,reproj_dataY,altitude))
#compare objects
print('original Mesh points -> ',mesh.points)
print('original Mesh type: ',type(mesh.points))
print('Reprojected points-> ',reprojData)
print('Reprojected type: ',type(reprojData))
Original Mesh points -> [[958427.33 119680.95 2396.288549 ]
[957754.39 120023.85 2430.1833881 ]
[957256.56 120241.02 2112.22953263]
...
[963366.748527 115096.364632 3054.75408138]
[963401.840285 113351.753238 3024.50286566]
[963497.913738 113339.696062 3048.83674197]]
Original Mesh type: <class 'pyvista.core.pyvista_ndarray.pyvista_ndarray'>
Reprojected points-> [[ 6.96487903 45.9823843 2396.288549 ]
[ 6.95646936 45.98581994 2430.1833881 ]
[ 6.95021969 45.98803333 2112.22953263]
...
[ 7.02498443 45.93857775 3054.75408138]
[ 7.02409542 45.92289079 3024.50286566]
[ 7.02532248 45.92273099 3048.83674197]]
Reprojected type: <class 'numpy.ndarray'>enter code here
现在，是时候替换 vtu 对象的坐标了：
mesh.points = reprojData
最后，我检查了修改后的mesh：X bounds 和 Y bounds 已经被修改并且范围是正确的。但是，该图显示了一条点线，而不是一个漂亮的 3d 对象。 :(.
你知道出了什么问题吗？您是否看到了另一种管理重新投影的方法？
解决方法
变换后XY和Z的取值范围明显不同：
>>>np.array(mesh.bounds).reshape((3,-1)).ptp(axis=1)
array([1.22515302e-01,7.78657599e-02,2.47978788e+03])
XY 确实以度为单位，Z 仍然以米为单位。这些数据的可视化表示是不切实际的，应该对数据进行调整。