地学分析与三维可视化实习一(共11页)

1、地学(dxu)分析与三维可视化实习报告一姓名(xngmng)：方志策班级(bnj)：064131学号：20131003574提交日期：2015.11.5作业(zuy)一、任务(rn wu)一、利用(lyng)对象图形法创建一个三维立方体，并将各顶点设置为不同的颜色。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(IDLgr

2、Model);创建多边形oPoly=OBJ_NEW(IDLgrPolygon)oView.add,oModel&oModel.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

3、,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运行出来的结果如下图：任务(rn wu)二、读取head.dat中的数据(shj)，进行体数据(shj)的显示，并实现切片提取操作。IDL代码如下：pro volume device

4、,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 = 10

5、0,ysize = 57 tv,volume23,*,* window,2,xsize = 80,ysize = 57 tv,volume*,20,* ;显示水平(shupng)切片 window,3,xsize = 80,ysize = 100 tv,volume*,*,20end显示(xinsh)结果如下 任务(rn wu)三、读取elev_t.jpg图像以及elevbin.dat中的DEM数据：利用直接(zhji)图形法显示(xinsh)DEM曲面利用对象图形法创建曲面对象和纹理对象并进行叠加显示代码如下：imageFile=FILEPATH(elev_t.jpg,SUBDIRECTOR

6、Y=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 Data&SHADE_SURF,demoModel=OBJ_NEW(IDLgrModel)oView=OBJ

7、_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;计算归一化显示比例，并在各个方向平移-0.5，从而使图像居中xs=NORM_COORD(xr)&xs0=xs0-0.5ys=NORM_COORD(yr)&ys0=

8、ys0-0.5zs=NORM_COORD(zr)&zs0=zs0-0.5oSurface-SetProperty,XCOORD_CONV=xs,YCOORD_CONV=ys,ZCOORD=zs ;TEXTURE_MAPoSurface-SetProperty,TEXTURE_MAP=oImage,color=255,255,255 oModel-Add,oSurface&oView-Add,oModeloModel-ROTATE,1,0,0,-90&oModel-ROTATE,0,1,0,30oModel-ROTATE,1,0,0,30&oWindow-Draw,oView ;OBJ_DEST

9、ROYXOBJVIEW,oModel,/BLOCK,SCALE=1&OBJ_DESTROY,oView,oImage end 作业(zuy)二、任务(rn wu)一、绘制(huzh)函数z=x*exp-x2+ y2，其中-2x,y2的网格曲面图。IDL代码如下：pro test_surface ;对x,y采样 x = (findgen(41) - 20)/10 y = (findgen(41) - 20)/10 ;对x,y网格化 temp_x = make_array(n_elements(y),value = 1) temp_y = make_array(n_elements(x),valu

10、e = 1) xx = x # temp_x yy = temp_y # y ;计算函数的值 z = xx * exp(-xx2 - yy2) ;绘制曲面 surface,zend 得到结果如图所示：任务(rn wu)二、读取head.dat中的三维动画数组，播放(b fn)该动画，并存储该动画的像素映射图IDL代码(di m)如下：pro test_animate file = filepath(head.dat,subdirectory = examples,data) ;读取三维动画数据head = read_binary(file,data_dims = 80,100,57);xint

11、eranimate,set = 80,100,57,/showload ;初始化动画工具 for i = 0,56 do xinteranimate,frame = i,image = head*,*,i ;将数组加载到动画工具的缓冲区 xinteranimate,50,/keep_pixmaps ;播放该动画，并存储像素映射图End得到结果如下图：任务(rn wu)三、载入avhrr.png中的数据，并加以(jiy)显示将其转换到“Interrupted Goode”投影(tuyng)坐标系下，并用iimage命令显示。将第一问中得到的地图转换到“Mollweide”投影坐标系下，并用iim

12、age命令显示。IDL代码如下：pro test_projection file = filepath(avhrr.png,subdirectory = examples,data) ;读取数据 data = read_png(file,r,g,b) red0 = rebin(rdata,360,180) green0 = rebin(gdata,360,180) blue0 = rebin(bdata,360,180) ;对原始数据进行重采样 iimage,red = red0,green = green0,blue = blue0,dimensions = 500,600,view_gri

13、d = 1,3 ;数据显示 smap = map_proj_init(Interrupted Goode) ;创建Interrupted Goode投影 red1 = map_proj_image(red0,map_structure = smap,mask = mask,uvrange = uvrange,xindex = xindex,yindex = yindex) green1 = map_proj_image(green0,xindex = xindex,yindex = yindex) blue1 = map_proj_image(blue0,xindex = xindex,yin

14、dex = yindex) ;投影转换 iimage,red = red1,green = green1,blue = blue1,alpha = mask*255b,/view_next ;显示转换后的地图 mapstruct = map_proj_init(Mollweide,/gctp) ;创建Mollweide投影 red2 = map_proj_image(red1,uvrange,image_structure = smap,map_structure = mapstruct,mask = mask,xindex = xindex2,yindex = yindex2) green2 = map_proj_image(green1,xindex = xindex2,yindex = yindex2) blue2 = map_proj_image(blue1,xindex = xindex2,yindex = yindex2) ;投影转换 iimage,red = red2,gre