Paraview将地理坐标系转化成直角坐标系

由于CitcomCU的输出文件网格点坐标是（余纬度，经度，径向深度），属于地理坐标系，如果不做处理直接进行可视化的话，点位图是这样的：

要进行坐标转换首先要将Rectilinear Grid类型的数据转化为Structured Grid，可以添加一个Rectilinear Data to Point Set的Filter。

这样就会多出来一个RectilinearDatatoPointSet1。

然后继续添加一个Calculator的Filter，将地理坐标通过运算转化为直角坐标，具体运算公式如下：

iHat*coordsX*cos(coordsY-3.141592/2)*cos(coordsZ)+jHat*coordsX*cos(coordsY-3.141592/2)*sin(coordsZ)+kHat*coordsX*sin(coordsY-3.141592/2)

将这一段添加到Calculator的属性栏中，并勾选Coordinate Results选项，表示用计算结果更新网格点的坐标，最后点Apply应用：

可以看到现在的网格点已经变成了球盖形状：

到这一步转换的工作完成。

这篇文章只是拿CitcomCU的输出举例，实际使用可能需要自己改一下相关的运算公式。