大地测量观测数据可视化MATLAB工具箱：(2)球面数据的显示

接着上一专栏，我们这里介绍如何在球面显示数据。其实球面的数据显示具有非常形象的特点，比如之前我写过的GRACE等效水高的球面3D表达。

（1）绘制球面某点等距离的直线

首先我们介绍circ4sph.m函数，其计算在任意参考点的球面上具有恒定球面距离和恒定方位角的直线。

需要添加以下函数：R2.m, R3.m, jump2nan.m。

结果图：

（2）skyplot.m 在地心坐标系统中显示(天体的)位置，即观测位置的方位角和仰角。

skyplot函数为：

结果图：

（3）绘制3D的球体

其中的Ynm为一个MxN的二维矩阵。

函数sphrwarp.m

函数shpere2.m

结果图：

参考资料：

（1）https://www.gis.uni-stuttgart.de/en/research/downloads/visbundle/

（2）The toolbox VISBUNDLE (version #/date #) was developed by N. Sneeuw, M. Weigelt,M. Roth, M. Antoni, M. Tourian et. al. and its latest version is provided via download from the Institute of Geodesy (GIS), University of Stuttgart .