％百分比在度之间递增量％griddeginc = 0.5;％程度％用于使用纬度扫描的度数％的％％值。% lon = -180:gridDegInc:180;％程度％% %要扫描的大地纬度值。％geod_lat = -90：griddeginc：90;％程度％％％循环通过每个纬度阵列的经度值-90：90。LONIDX = SIZE（LON，2）： -  1：1％％％的经度必须与纬度阵列相同％lon = lon（lonidx）*那些（1，numlate）;％程度% geoidResults(1:end,lonIdx) = geoidheight(geod_lat,lon，'None');％％ 结尾

geoidFileName =“GeoidResults_05deg_180.mat”;加载（Geoidfilename）;海岸=加载（'海岸线.mat'）;

％使用| MESHM |创建2-D图H2D =数字;设置（H2D，“位置”，[20 75 700 600]，'工具栏'那“图”）;用于将大地带高度映射到GLAT / LON的大量参考矩阵。存款准备金率= 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]）;％创建滑块滑= astGeoidSlider (geoidcanvas2D);

h3D =图;集(h3D,“位置”，[20 75 700 600]）;％设置轴hmapaxis = axesm（'地球'那'网格'那'上'）;设置（hmapaxis，“位置”，[。5 .8 .4]）视图（85,0）轴离开3-D全球的％绘图数据meshm (geoidResults,存款准备金率)%地块地块轮廓Plotm（海岸.Coastlat，Coast.coastlon，'颜色'那'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）;