北大《空间探测信息处理技术(IDL)》第11章 体绘制.ppt

体绘制工具 XVOLUME XVOLUME, Vol, [,/BLOCK] [,GROUP=widget_id] $ [,/INTERPOLATE] [,/MODAL] [,RENDERER={0|1}] $ [,/REPLACE] [,SCALE=value] [,/TEST] $[, XSIZE=pixels] [, YSIZE=pixels] 参数 Vol : 三维数组 关键字 /BLOCK : 阻止方式 GROUP=widget_id : 设置组件的组 /INTERPOLATE : 设置三线性插值 /MODAL : 模式窗口 RENDERER={0|1} : 渲染方式 =0 : 由OpenGL渲染，默认方式 =1 : 由IDL软件负责渲染 体绘制工具 XVOLUME XVOLUME, Vol, ... 关键字 /REPLACE : 如果在此之前XVOLUME正在运行，则不重起XVOLUME，而是使用新Vol覆盖旧数据。 SCALE=[x,y,z] : 缩放因子，默认为1/SQRT(3) /TEST : 测试，忽略Vol数据，产生测试数据 XSIZE=pixels : 绘图区域的 X 像数点设置 YSIZE=pixels : 绘图区域的 Y 像数点设置 菜单 File-Export Image... : 输出结果到图形格式文件； File-Export VRML : 输出结果到虚拟现实格式数据； File-Print : 打印输出 体绘制工具 例10.2 example11_Xvolume.pro pro exam11_1 file = FILEPATH(head.dat, $ SUBDIRECTORY = [examples, data]) volume = READ_BINARY(file, DATA_DIMS =[80,100,57]) ; Reduce the size of the original array to ; speed up processing: smallVol = CONGRID(volume, 40, 50, 27) ; Using the INTERPOLATE keyword to smooth the data, ; display the volume using the XVOLUME procedure: XVOLUME, smallVol, /INTERPOLATE end 体数据的断层提取 垂直切片 体数据：Vol[x,y,z] 按照坐标(例如：x)，获取二维数组； image = reform(Vol[1,y,z],dimy,dimz) 按照图像存储。 WRITE_IMAGE, Filename, Format, Data $ [, Red, Green, Blue] [, /APPEND] 体数据的断层提取 例10.3 ; Chapter10DisplaySlicesUpright.pro PRO Chapter10DisplaySlicesUpright file = FILEPATH(head.dat, $ SUBDIRECTORY = [examples, data]) image = READ_BINARY(file, $ DATA_DIMS = [80, 100, 57]) LOADCT,5 DEVICE, DECOMPOSED = 0, RETAIN = 2 WINDOW, 0, XSIZE = 800, YSIZE = 600 FOR i = 0, 56, 1 DO TVSCL, 255b - image [*,*,i], /ORDER, i sliceImg = REFORM(image[40,*,*]) sliceImg = CONGRID(sliceImg, 100, 100) TVSCL, 255b - sliceImg, 47 END 体数据的断层提取 任意方向切片 Result = EXTRACT_SLICE( Vol, Xsize, Ysize, $ Xcenter, Ycenter, Zcenter, Xrot, Yrot, Zrot $ [,ANISOTROPY=[xspacing,yspacing, zspacing]] $ [, OUT_VAL=value] [, /RADIANS] [, /SAMPLE] $ [, VERTICES=variable] ) 或者 Result = EXTRACT_SLICE(