％百分比在度之间递增量％griddeginc = 0.5;％程度％用于纬度扫描的度数以学位的％％%%％。％geoc_lon = -180：griddeginc：180;％程度％扫描的％％大地纬度值。％geod_lat = -90：griddeginc：90;％程度％％百分比转换为地理尺度以获得大量地高度。％重新= 6378137;％f = 1/298.256415099;％geoc_lat = geod2geoc（geod_lat，0，f，Re）;％％％循环通过每个纬度阵列的经度值-90：90。LONIDX =大小（GEOC_LON，2）： -  1：1％％％的经度必须与纬度阵列相同％lon = geoc_lon（lonidx）* on（1，numlatte）;％程度％GeoidResults（1：结束，LonIdx）= GeotioHeight（Geoc_lat，Lon，'none'）;％％ 结尾

Geoidfilename =.'geoidresults_05deg_180.mat';加载（Geoidfilename）;海岸=负载（'海岸'）;

％使用| MESHM |创建2-D图H2D =数字;设置（H2D，'位置'，[20 75 700 600]，'工具栏'那'数字'）;用于将大地带高度映射到GLAT / LON的大量参考矩阵。rrr = makerefmat（'rastersize'，大小（地别的滤），......'latlim'，[-90 90]，'lonlim'，[-180 180]）;Ast2dgeoidplot（RRR，Geoidresults，Coast，Griddeginc）％使用VR帆布查看大猩猩高度www2d = vrworld（'Astgeoidheights.wrl'）;打开（www2d）％实际大地大麻高度供参考Geoidgrid = Vrnode（www2d，'egm96_grid'）;stapthheights = getfield（地理群岛，'高度'）;％＃OK 对于滑块控制，％初始化高度为0Geoidgrid.height = 0 *实际亮度;用于绘制和设置参数的％大小画布Geoidcanvas2d = Vr.Canvas（www2d，'父母'，H2D，......'抗锯齿'那'上'那'navspeed'那'非常慢'那......'navmode'那'检查'那'单位'那'标准化'那......'观点'那'看法'那'位置'，[。15 .04 .7 .42]）;％创建滑块Slid = Astgeoidslider（Geoidcanvas2d）;

H3D =数字;设置（H3D，'位置'，[20 75 700 600]）;％设置轴hmapaxis = axesm（'地球'那'网格'那'上'）;设置（hmapaxis，'位置'，[。5 .8 .4]）视图（85,0）轴离开3-D全球的％绘图数据Meshm（GeoidResults，RRR）％情节陆地大规模轮廓Plotm（Coast.lat，Coast.long，'颜色'那'K'）Colormap（'喷射'）;％绘图标题标题（{'egm96大猩猩高度'; ['网格增量：'，num2str（griddeginc），'度;高度单位：米']}）彩色杆;％3-D全球：使用VR帆布的大猩猩高度www3d = vrworld（'Astgeoidsphere.wrl'）;打开（www3d）％位置帆布Geoirocanvas3d = VR.Canvas（www3d，'父母'，H3D，......'抗锯齿'那'上'那'navspeed'那'非常慢'那......'navmode'那'检查'那'单位'那'标准化'那......'位置'，[。15 .04 .7 .4]）;Vrdrawnow;

关闭（H2D，H3D）关闭（www2d）;关闭（www3d）;删除（www2d）;删除（www3d）;