% gunzip南旧金山DEM文件到一个临时目录文件名= gunzip (“sl3d_sanfranciscos.dem.gz”，tempdir）；demFilename=文件名{1}；%读取1:24 000 DEM文件%将缺少的值标准化为NaNdeminfo = georasterinfo (demFilename);Z = standzemissing (readgeoraster(demFilename)， deminfo. missingdatindicator);%删除临时压缩文件删除(demFilename);

demdata=Z；[xdim，zdim]=size（demdata）；xspace=30；x维度以米为单位缩放百分比zspace = 30;%z尺寸的换算单位（米）%将数据重塑为一维数组demdata = demdata (:);

打开模板myworld = vrworld (“vr_template_terrain.x3d”);打开虚拟世界开放（myworld）；%创建一个节点VRTerrain的句柄，该节点将包含DEM数据地形节点=vrnode（myworld，“VRTerrain”);

%创建地形形状的子对象newShape = vrnode (Terrain_node,“孩子”，“Terrain_Shape”，“形状”);为形状创建外观字段新出现=vrnode（新形状，“外观”，“Terrain_Appearance”，“外观”);为外观创建材质字段newMat = vrnode (newAppear,“材料”，“地形材料”，“材料”);为材质字段分配属性newMat。ambientIntensity = 0.25;newMat。diffuseccolor = [0.9 0.6 0.6];newMat。亮度= 0.078125;newMat。speccolor = [0.0955906 0.0955906 0.0955906];%为形状创建几何体字段newEGrid = vrnode (newShape,“几何学”，“DEM_EGrid”，“升降网格”);%为几何域分配属性-使用DEM数据newEGrid.creaseAngle=3.14；newEGrid.xDimension=xdim；newEGrid.zDimension=zdim；newEGrid.xspace=xspace；newEGrid.zspace=zspace；newEGrid.height=demdata；newEGrid.ccw=“真的”；这个设置将使地形表面从两边都可见newEGrid。固体=“假”；

%地形高程用于为图像着色cmap=demcmap（Z，256）；%创建当前目录的纹理子目录%输出参数仅用于在目录%已经存在[~，~]=mkdir(“纹理”);%缩放高度值以使用完整的颜色贴图范围%缩放依赖于此地形从零高度开始的事实Zscaled=Z.*（大小（cmap，1）-1）。/max（Z（：）；%在纹理子目录中将纹理保存为PNG图像%向左旋转图像以匹配X3D模型中所需的图像方向Zscaled中的%元素表示cmap中的索引imwrite（rot90（Zscaled）、cmap、，“纹理/ sanfrancisco_elev.png”);

newTexture = vrnode (newAppear,“纹理”，“地形纹理”，“ImageTexture”); newTexture.url=“纹理/ sanfrancisco_elev.png”；

%本架飞机。X3d从/sl3ddemos到当前目录pt = fileparts ((“vrterrain_simple.m”))；复制文件（完整文件），“b747.x3d”pwd),“f”);

%创建一个名为“Boeing”的新变换节点平面=vrnode（myworld，“波音”，“转变”)；平面_inline=vrnode（平面，“孩子”，“Boeing_Inline”，“内联”);在/sl3ddemos目录中准备了一个简单的Boeing模型plane_inline.url =“b747.x3d”；

ypeak=max（Z（：）；[xmax，zmax]=find（Z==ypeak）；%如果更多顶点具有相同的最大高度，则使用第一个峰值%将矩阵索引转换为X3D中x和z方向上的米xpeak = xspace * (xmax (1) 1);zpeak = zspace * (zmax (1) 1);

飞机。翻译= [xpeak ypeak+200 zpeak];%将飞机的大小缩放20倍，以便%在虚拟场景中可见，无需进行任何额外的缩放飞机。[20 20 20];

％%三元组由从一个顶点开始的3条线（1米长）组成%沿x、y和z方向。线条的颜色如下：% +x -红色% +y -绿色%+z-蓝色％%如果这个三元组被包含在场景的顶层%层次，它表示全局场景坐标。%如果它作为变换节点的子节点包含，则表示局部%场景中该节点的坐标系（方向）。%将三元组添加到场景的层次结构的顶层三和弦= vrnode (myworld,“Triad1”，“三位一体”);%缩放空间坐标轴的大小，使其可见%在虚拟场景中，无需进行任何额外的缩放triad.scale=[xdim*xspace/8分钟（xdim*xspace/8，zdim*zspace/8）zdim*zspace/8]；将三角放置在波音747的中心三元组位置=[xpeak ypeak+200 zpeak]；

myworld.World\u Info.title=“B747飞越旧金山地区”；

拯救（我的世界），“Vrratrain_sanfrancisco.x3d”);

关闭（myworld）；删除（myworld）；

%这是如何在外部查看器中打开X3D文件:%vrview（'terrain.x3d'，'-web'）；%这是如何在内部查看器中打开X3D文件:% vrview(“地形。x3d”、“内部”);%以下是如何在默认查看器中打开X3D文件：createdworld = vrview (“Vrratrain_sanfrancisco.x3d”);设置反锯齿以平滑地形纹理myfig=get（createdworld，“数字”)；设置（myfig，“抗锯齿”，“上”);

清除所有已使用的变量清楚的地形节点ZZscaled提出创建世界demFilename德姆达塔前女友文件名．..id拉特朗混乱myfigmyworldnewAppear纽格里德．..newMatnewShape新织物纳米好啊飞机平面线pt三人一组．..已经xdim最大值xpeakxspaceypeak兹迪姆zmaxzpeakzspace