地学分析与三维可视化实习一..docx

地学分析与三维可视化实习报告一姓名：方志策班级：064131学号：20131003574提交日期：2015.11.5作业一、任务一、利用对象图形法创建一个三维立方体，并将各顶点设置为不同的颜色。IDL代码如下：oWindow=OBJ_NEW(IDLgrWindow,dimension=[400,400]);建立一个图像窗口，初始设置图像大小为400*400oView=OBJ_NEW(IDLgrView,viewPlane_Rect=[-1,-1,3,3],zClip=[3,-3],eye=2) oModel=OBJ_NEW(IDLgrModel);创建多边形oPoly=OBJ_NEW(IDLgrPolygon)oView.add,oModeloModel.add,oPolyverts=[[0,0,0],[1,0,0],[1,1,0],[0,1,0],[0,0,1],[1,0,1],[1,1,1],[0,1,1]] ;设置立方体顶点connect=[4,0,1,2,3,4,0,1,5,4,4,1,2,6,5,4,2,3,7,6,4,3,0,4,7,4,4,5,6,7];将一个立方体各个定点依次标明，依次连接八个面构成一个闭合的立方体oPoly.setproperty,data=verts,polygons=connect,style=1oModel.rotate,[1,0,1],60;设置立体图像显示出来的角度为60度oWindow.draw,oViewvertscolor=[[0,0,0],[1,0,0],[1,1,0],[0,1,0],[0,0,1],[1,0,1],[1,1,0],[0,1,1]]*255 ;用RGB成像依次构造立方体各个部位的不同颜色oPoly.setproperty,vert_color=vertscolor,style=2,shading=1oWindow.draw,oViewend运行出来的结果如下图：任务二、读取head.dat中的数据，进行体数据的显示，并实现切片提取操作。IDL代码如下：pro volume device,decomposed = 1 file = filepath(head.dat,subdirectory = [examples,data]) ;读取体数据到volume数组 volume = read_binary(file,data_dims = [80,100,57]) ;利用XVolume交互显示数据 xvolume,volume,/interpolate slice = extract_slice(volume,40,40,40,50,28,30,0,0) window,0,xsize = 40,ysize = 40 tv,slice ;显示垂直切片 window,1,xsize = 100,ysize = 57 tv,volume[23,*,*] window,2,xsize = 80,ysize = 57 tv,volume[*,20,*] ;显示水平切片 window,3,xsize = 80,ysize = 100 tv,volume[*,*,20]end显示结果如下 任务三、读取elev_t.jpg图像以及elevbin.dat中的DEM数据：利用直接图形法显示DEM曲面利用对象图形法创建曲面对象和纹理对象并进行叠加显示代码如下：imageFile=FILEPATH(elev_t.jpg,SUBDIRECTORY=[examples,data]) ;读取图像文件 READ_JPEG,imageFile,imagedemFile=FILEPATH(elevbin.dat,SUBDIRECTORY=[examples,data]) ;读取DEM数据dem=READ_BINARY(demFile,DATA_DIMS=[64,64])dem=CONGRID(dem,128,128,/INTERP)DEVICE,DECOMPOSED=0,RETAIN=2;TITLEWINDOW,0,TITLE=Elevation DataSHADE_SURF,demoModel=OBJ_NEW(IDLgrModel)oView=OBJ_NEW(IDLgrView)oWindow=OBJ_NEW(IDLgrWindow,RETAIN=2,COLOR_MODEL=0) oSurface=OBJ_NEW(IDLgrSurface,dem,STYLE=2)oImage=OBJ_NEW(IDLgrImage,image,INTERLEAVE=0,/INTERPOLATE) oSurface-GetProperty,XRANGE=xr,YRANGE=yr,ZRANGE=zr;计算归一化显