本章概要
在上一章节了学习了一些图形显示技 术。在这一章节将学习几个新的图形显示技巧，以便让图形显示具有专业的感观效果。
具体来说，将学会：
1. 怎样让鼠标交互作用于图形显示
2. 怎 样从图形显示中删除注释
3. 怎样在图形显示上画“橡皮条”
4. 怎样在图形显示技巧中使用Z图形缓冲区
将光标用于图形显示
数 据可视化显示的原因之一是用户可用不同的方式对数据进行交互式的操作。用户喜欢的一种对数据交互式操作的方式是使用光标去选择或者标注部分数据。这种交互 作用在IDL中用Cursor命令很容易完成。
用LoadData命令装入Time Series数据集，可看到Cursor命令是如何工作的。
IDL> curve = LoadData (1)
输 入下述命令，显示曲线：
IDL> Window, Xsize = 400, Ysize = 400
IDL> LoadCT, 0
IDL> TvLCT, 255,255,0,1
IDL> Plot, curve
Cursor 命令接受两个参数。这些参数必须是记录鼠标键按下时光标位置的变量。Cursor命令要求光标位于当前图形窗口中。（即被!D.Window 系统变量指 向的窗口。）例如，如果输入这个命令，IDL将会等待光标被移动到当前图形窗口（如果输入的是上述命令，就是0号索引窗口）并单击鼠标键。当执行上述动作 后，IDL将光标位置返回到变量xLocation和yLocation中。输入：
IDL> Cursor, xLocation, yLocation
如 果打印出这些变量的值，将发现这些值被赋予的是数据坐标空间。xLocation的数值从0到100，yLocation的数值从0到30。（如果是在图 形边界内点击的鼠标，它们至少是这么多。如果不是在图形边界内点击的鼠标会怎么样？）缺省时，Cursor命令返回数据坐标位置。
IDL> Print, xLocation, yLocation
什 么时候返回的光标位置？
从上面的命令看，似乎鼠标键被按下时返回光标位置，但并非总是这样。事实上，Cursor命令什么时候报告光标的位置是由 Cursor命令的关键字所决定的。这些关键字是：
Change 当光标位置发生改变或用户移动光标时，返回光标位置。
Down 当鼠标 键被按下时，返回光标位置。
NoWait 当Cursor命令执行时，光标位置被立即返回。没有任何延迟或等待鼠标的按键。这个关键字有时用于当 对象正在显示窗口中被移动时的循环中。
UP 不是在鼠标键被按下时，而是放开或释放后返回光标位置。
Wait Cursor命令等待鼠标 键被按下后返回光标的位置。只要鼠标键被按下，此关键字对Cursor命令的作用就类似于用NoWait关键字调用Cursor命令。此关键字是 Cursor命令的缺省状态。
在Cursor命令中，小心使用合适的关键字，特别是在循环过程中使用Cursor命令。用户有时习惯地认为 Cursor命令的缺省属性是只有鼠标键被按下时才返回光标的位置。其实不然，缺省属性只是等待一个单击动作，以后的行为就和NoWait关键字一样。在 循环中这个区别是至关重要的。
哪一个鼠标键和光标共同作用呢？
除了设置光标属性外，有时还想知道哪个鼠标键用于对Cursor命令作出反 应。例如，想要用鼠标左键做某件事，而做另外不同的事情要用鼠标右键 Cursor命令作出的反应。可以检查系统变量!Mouse中的Button字段， 来判断哪一个鼠标键在和Cursor命令共同作用。（老版本的IDL是用系统变量!Err的值来判断的。）这个字段是一个整型位映象。Button这个字 段的有效值及其意义如下：
!Mouse.Button = 0 当前没有按键被使用
!Mouse.Button = 1 左键用于 Cursor命令
!Mouse.Button = 2 中间键用于Cursor命令
!Mouse.Button = 4 右键用于 Cursor命令
用光标标注图形输出
使用Cursor命令的一种方法是允许用户交互地在线画图上放置符号标记。例如，正确无误地输入下列 命令。当输完最后一个回车键后，在当前的图形窗口上单击鼠标五次。五个符号将放置在窗口中。（如果在输入下列代码时出现打字错误，必须从头开始重新输 入。）输入：
IDL> For j = 0, 4 DO BEGIN $
IDL> Cursor, xloc, yloc, /DOWN & $
IDL> Plots, xloc, yloc, Psym = 4, SymSize = 2, Color = 1 & ENDFOR
画 方框
有时可能为了选取图形显示中的某部分，而在它的周围画上方框。这里有些命令可用来选择由Cursor命令产生的方框的对角，画出该方框，并将 图形缩放到该方框坐标范围。首先画图：
IDL> Plot, curve
接着，使用光标选择想画的方框的一角。要确保在当前图形窗 口上点击光标。为确定哪个是当前窗口，并让它不被隐藏，可输入：
IDL> WShow
现在键入第一个Cursor命令。在图形轴的 范围内某处点击：
IDL> Cursor, x1,y1, /DOWN; Select one corner of box.
接 着输入第二个Cursor命令。在图形轴的范围内某处点击：
IDL> Cursor, x2,y2, /DOWN; Select diagonal corner of box.
上 述Cursor命令返回的坐标是数据空间坐标。按如下画方框：
IDL> Plots, [x1,x1,x2,x2,x1],[y1,y2,y2,y1,y1], color = 1
输 出结果应类似于图58中所示，尽管实际的图形上方框取决于在窗口中点击的位置。
为了放大这部分图形，必须保证方框坐标的正确顺序。这是非常必要 的，因为可能先选择的是方框的右下角，然后是左上角，这样x1将大于x2。还可以想象其它的假设。为了适应所有的情况，键入：
图58：在部分数据 周围画上方框的线画图。用Cursor命令选择方框的坐标，用PlotS命令画方框。

IDL> Xmin = Min([x1,x2], Max = xmax)
IDL> Ymin = Min([y1,y2], Max = ymax)
最 后，已经为放大对方框内的数据做好了准备。除了正确地设置数据范围外，还必须设置Style关键字。知道为什么吗？如果不知道，可在不使用这两个关键字的 情况下试试下面的命令。将会发生什么呢？
IDL> Plot, curve, XRange = [xmin, xmax], Yrange = [ymin, ymax], $
Xstyle = 1, Ystyle = 1

在 图像上使用Cursor命令
通常当在处理图像数据时使用Cursor命令，希望用设备坐标而不是数据坐标返回光标位置。这是因为设备坐标和图像中 对应的位置之间通常存在一种简单的关系（大多数是一对一的关系）。为了解如何工作的，可用LoadData命令打开360*360的 World Elevation数据集，键入：
IDL> image = LoadData(7)
显示图像，并装入某些颜色。如 下：
IDL> topColor = !D.Table_Size-1
IDL> LoadCT, 3, Ncolors = !D.Table_Size-1
IDL>TvLCT, 255,255,0, TopColor
IDL> Window, XSize = 360, YSize = 360
IDL> TV, BytScl (image, Top = !D.Table_Size-2)
利 用光标在图像中选择某一特定行和列。注意Cursor和PlotS命令中的Device关键字。这是确保返回的坐标是设备坐标而不是数据坐标。在该位置上 画一个十字线。（确保在输入Cursor命令后，在图像窗口中点击一下。）键入：
IDL> S = Size(image)
IDL> Cursor, col, row, / Device; Click in the window!
IDL> Plots, [col, col], [0,s (2) ], / Device, Color = topColor
IDL> Plots, [0,s (1) ], [row, row ], / Device, Color = topColor
注 意，在图像中某一特定的行和列上获得图像数据是多么的容易。例如，可以轻易地绘制出图像中行和列的数据剖面，键入：
IDL> Window, 1, Xsize = 500, Ysize = 300
IDL> !P.Multi = [0, 2, 1]
IDL> Plot, image [*, row], Title =’ Row Profile’
IDL> Plot, image [col, *], Title =’ Column Profile’
IDL> !P.Multi = 0
IDL> Wset, 0
输 出结果类似于图59所示。
在循环中使用Cursor命令
有时想在循环中使用Cursor命令。例如，当用光标选择图像上的单个像素时，可 能想知道它的像素值。下面是个简单的循环程序，它将一直执行下去，直到单击右键或中键退出。打开文本编辑器，准确无误地输入如下代码。
TopColor = !D.Table_Size-1
LoadCT, 3, Ncolors = !D.Table_Size-1
TvLCT, 255, 255, 0, topColor
TV, BytScl (image, Top =!D.Table_Size-2)
!Mouse.Button = 1
REPEAT BEGIN
Cursor, col, row, /Down, /Device
Print, ‘Pixel Value:’, image[col, row]
ENDREP UNTIL !Mouse.Button NE 1
END

图 59：用Cursor命令在图像中选择的行和列的剖面。

保存到loop1.pro文件中（此文件已经存在于下载的本书配套程序之中）。编 译，并运行这个小主程序，输入：
IDL> .RUN loop1
移动光标进入图像窗口，开始单击左键。图像像素值就会出现在日志窗 口中，直到点击其它键，而不是左键。
如果在Cursor命令中使用除了Down之外的其它关键字会发生什么呢？实验一下，找出答案。
从显 示中删除注释
当使用光标在图形显示上按照刚才使用的方法来添加注释时，也许会问：“怎样才能删除刚放置在那儿的注释呢？”有两种较好的方法删除注 释。称之为异或法和设备拷贝法。笔者认为在两者之中设备拷贝法更能给出专业的感观效果。两种方法都列举出来，但重点将集中在设备拷贝法上。
删除注 释的异或法
异或法是在图形函数的基础上起作用的。图形函数是两个数的位操作。这两个数分别与已经显示出来的像素（称作所谓的目标像素）以及希望放 置在同一位置的像素（称作所谓的源像素）相关联。
通常，IDL使用的图形函数称作源。在这种图形函数里，IDL忽略了目标像素的值，仅仅在该像素 位置上放置源像素的值。但如果这种图形函数变成XOR（异或法），IDL将目标像素和源像素进行逐位比较。这会产生反向目标像素的效果。换句话说，如果像 素的二进制表示为01100101，那么执行XOR命令后，像素的二进制表示为10011010。
（实际的XOR过程远比这复杂，因为只有IDL 在颜色索引表中的邻近位置上有256种颜色时，它才按这种方式运行，而这是一种少见的情形。大多数人只是认为XOR法是用“相反的”颜色画，并保留原来 的。在实际的XOR法中，可能预测将会使用哪种颜色来画，但在大多数情形下并不是这样。这就是为什么绝大多数 IDL专业程序员宁愿采用设备拷贝法。）
在 任何时侯，图形函数的作用效果都是由Device命令和Set_Graphics_Function关键字设置的。源模式下图形函数为3。XOR模式下图 形函数为6。此时IDL处于系统缺省的源模式下。当处于这种模式时，重新显示图像窗口中的图像。输入：
IDL> TV, BytScl (image, Top =!D.Table_Size-2)
现 在，选择XOR模式：
IDL> Device, Set_Graphics_Function = 6
在图像中画一个方框：
IDL> Plots, [0.2, 0.2, 0.8, 0.8, 0.2], Color = topColor,$
[0.2, 0.8, 0.8, 0.2, 0.2], /Normal
注 意，方框线的颜色不是预测的黄色。取而代之是多彩的，尽管它显得也很合理。在这个模式下线下的像素已经被翻转了。
要删除方框，只须将底下的像素值 翻转回它们的原值。再次用PlotS命令很容易完成。
IDL> Plots, [0.2, 0.2, 0.8, 0.8, 0.2], Color = topColor,$
[0.2, 0.8, 0.8, 0.2, 0.2], /Normal
可 以反复执行上面的命令，让方框随心所欲地出现或消失。在继续下面的练习前，确保将图形函数设回到源模式。键入：
IDL> Device, Set_Graphics_Function = 3
在 IDL程序中可容易地利用图形函数。例如，打开以前写的loop1.pro主程序，按如下修改它。这个程序将在图像窗口中每次点击的位置上画一个大的十字 线。以loop2.pro命名保存程序。键入：
topColor = !D.Table_Size-1
LoadCT, 3, NColors=!D.Table_Size-1
TvLCT, 255, 255, 0, topColor
TV, BytScl(image, Top=!D.Table_Size-2)
!Mouse.Button = 1
; Go into XOR mode.
Device, Set_Graphics_Function=6
; Get initial cursor location. Draw cross-hair.
Cursor, col, row, /Device, /Down
PlotS, [col,col], [0,360], /Device, Color=topColor
PlotS, [0,360], [row,row], /Device, Color=topColor
Print, ‘Pixel Value: ‘, image[col, row]
; Loop.
REPEAT BEGIN
; Get new cursor location.
Cursor, colnew, rownew, /Down, /Device
; Erase old cross-hair.

PlotS, [col,col], [0,360], /Device, Color=topColor
PlotS, [0,360], [row,row], /Device, Color=topColor
Print, ‘Pixel Value: ‘, image(colnew, rownew)
; Draw new cross-hair.
PlotS, [colnew,colnew], [0,360], /Device, Color=topColor
PlotS, [0,360], [rownew,rownew], /Device, Color=topColor
; Update coordinates.
col = colnew
row = rownew
ENDREP UNTIL !Mouse.Button NE 1
;Erase the final cross-hair.
PlotS, [col,col], [0,360], /Device, Color=topColor
PlotS, [0,360], [row,row], /Device, Color=topColor
; Restore normal graphics function.
Device, Set_Graphics_Function=3
END

保 存到文件loop2.pro中。（此文件已经存在于下载的本书配套程序之中。）编译，并执行这个主程序，键入：
IDL> .RUN loop2
将 光标放在图像窗口中，点击左键数次。在每个光标位置上应该有一个十字线。按右键或中间键终止程序。
删除注释的设备拷贝法
设备拷贝法利用像 素映射窗口来删除已显示屏幕上的注释。像素映射窗口和其它IDL图形窗口完全一样，除了它不显示在显示器上外。事实上，它存在于显示设备的视频随机存储器 中。换句话说，它存在于存储器中。从其它任何方面来说，它就象一个正常的IDL图形窗口一样：可用Window命令创建，用WSet命令激活，用 WDelete命令删除，等等。在像素映射窗口内画图和在IDL正常的图形窗口内画图是完全一样的（例如，用Plot、Surface、TV和其它图形输 出命令）。
设备拷贝技术就是从一个窗口（称为源窗口）拷贝一个矩形区，然后将矩形区传入另一个窗口（称为目的窗口）。源窗口和目的窗口有时可以是 同一个窗口，这点等会会看到。图60是设备拷贝技术的图解。
实际的拷贝是通过Device命令和Copy关键字（设备拷贝技术的名字由此而来）完 成。命令的一般形式如下：
Device, Copy = [sx, sy, col, row, dx, dy, sourceWindowID]
在 这个命令中，Copy关键字的组成元素是：
sx, sy 源窗口矩形区左下角的设备坐标（矩形区是从源窗口拷贝的）
col 从源窗口中拷 贝的列数。这是矩形区的宽度。
row 从源窗口中拷贝的行数。这是矩形区的高度。
dx, dy 目标窗口中矩形区左下角的设备坐标。（目 标窗口是要将矩形区拷贝到的窗口。目标窗口总是当前图形窗口。）
sourceWindowID 源窗口的窗口索引号。矩形区从此窗口拷贝到当前图 形窗口（由!D.Window系统变量指定）。源窗口可以是当前图形窗口，但多数情况下它是一个非当前图形窗口的窗口。它通常是一个像素映射窗口。
要 看其是如何工作的，先创建一个像素映射窗口，在里面显示一幅图像。用带Pixmap关键字的Window命令来创建像素映射窗口，键入：
IDL> Window, 1, /Pixmap, XSize = 360, YSize = 360
IDL> TV, BytScl (image, Top =!D.Table_Size-2)
图 60：设备拷贝技术是将源窗口的一个矩形区拷贝到目标窗口的某个位置。实际上可以拷贝整个窗口，而且源窗口和目标窗口可以是同一个窗口。

注 意，当输入这些命令时没有任何可见的线索表明发生了什么。这是因为像素映射窗口仅存在于视频随机存储器中，而不是在显示屏上。为证明这个窗口内存有的内 容，可打开第三个窗口，并将像素映射窗口的内容拷贝到第三个窗口中。如果第三个窗口看起来象原来的图像窗口，那么已经输入的命令是正确的。键入：
IDL> Window, 2, Xsize = 360, Ysize = 360
IDL> Device, Copy =[0, 0, 360, 360, 0, 0, 1]
注 意，已将像素映射窗口的全部内容拷贝到这个新的窗口中。这种操作除了在速度上快几个数量级之外，其它方面与在此新窗口中重新显示图像类似了。将像素映射窗 口中的全部内容拷贝到了新的显示窗口是一种常见的做法，即使只是修改部分的显示窗口。
删除最后创建的两个窗口（包括像素映射窗口），如下：
IDL> Wdelete, 1, 2
当 对像素映射窗口操作完成后，删除它们是很重要的。它们会占用存储空间，这些存储空间可用于其它用途。某些窗口管理器为这些像素映射窗口分配一个定量的存储 空间。如果像素映射窗口超过了视频存储器的容量，可以使用虚拟内存。X终端给像素映射窗口留有非常小的存储空间。
要了解设备拷贝技术在实际中是如 何运用的，修改主程序以前写的主程序文件loop2.pro。可以拷贝到另一个文件，以loop3.pro命名。修改代码如下：
TopColor = !D.Table_Size-1
LoadCT, 3, Ncolors = !D.Table_Size-1
TvLCT, 255, 255, 0, topColor
TV, BytScl (image, Top =!D.Table_Size-2)
!Mouse.Button = 1
; Create a pixmap window and display image in it.
Window, 1, /Pixmap, Xsize = 360, Ysize = 360
TV, BytScl (image, Top =!D.Table_Size-2)
; Make the display window the current graphics window.
Wset, 0

; Get initial cursor location. Draw cross-hair.
Cursor, col, row, /Device, /Down
Plots, [col, col], [0, 360 ], / Device, Color = topColor
Plots, [0, 360 ], [row, row ], / Device, Color = topColor
Print, ‘Pixel Value:’, image[col, row]
; Loop.
REPEAT BEGIN

; Get new cursor location.
Cursor, colnew, rownew, /Device, /Down
; Erase old cross-hair
Device, Copy = [0, 0, 360, 360, 0, 0, 1]
Print, ‘Pixel Value:’, image[colnew, rownew]

; Draw new cross-hair.
Plots, [colnew, colnew], [0, 360 ], / Device, Color = topColor
Plots, [0, 360 ], [rownew, rownew ], / Device, Color = topColor
ENDREP UNTIL !Mouse.Button NE 1
; Erase the final cross-hair.
Device, Copy = [0, 0, 360, 360, 0, 0, 1]
END
以 loop3.pro命名存档。（此文件已经存在于下载的本书配套程序之中）编辑，并运行主程序，输入：
IDL> .RUN loop3
放 置光标于图像窗口中，点击左键数次。单击右或中间键终止程序。注意十字线是用黄色绘制。
在继续下一次练习之前，删除像素映射窗口。输入：
IDL> Wdelete, 1
画 一个橡皮筋方框
设备拷贝技术非常适用于在屏幕上画橡皮条选择方框和其它形状。（橡皮条方框是指一角固定，另一角随着光标的变化而变化的方框）。事 实上，修改程序 loop3.pro可以很容易实现。将loop3.pro拷贝到文件rubberbox.pro。（loop3.pro文件已经存在于下 载的本书配套程序之中。）作如下修改，看看创建一个橡皮条方框是多么容易。
topColor = !D.Table_Size-1
LoadCT, 3, NColors=!D.Table_Size-1
TvLCT, 255, 255, 0, topColor
TV, BytScl(image, Top=!D.Table_Size-2)
!Mouse.Button = 1
; Create a pixmap window and display image in it.
Window, 1, /Pixmap, XSize=360, YSize=360
TV, BytScl(image, Top=!D.Table_Size-2)
;Make the display window the current graphics window.
WSet, 0
; Get initial cursor location (static corner of box).

Cursor, sx, sy, /Device, /Down
; Loop.
REPEAT BEGIN
; Get new cursor location (dynamic corner of box).
Cursor, dx, dy, /Wait, /Device
; Erase the old box.
Device, Copy=[0, 0, 360, 360, 0, 0, 1]
; Draw the new box.
PlotS, [sx,sx,dx,dx,sx], [sy,dy,dy,sy,sy], /Device, $
Color=topColor
ENDREP UNTIL !Mouse.Button NE 1
;Erase the final box.
Device, Copy=[0, 0, 360, 360, 0, 0, 1]
END

运 行程序，输入：
IDL> .RUN rubberbox
在继续下一次练习之前，删除像素映射窗口。输入：
IDL> Wdelete,1
图 形窗口的滚动
设备拷贝技术的另一个有效应用是实现图形窗口的滚动。在这个示例中将运用设备拷贝技术让图形显示窗口中的图像滚动起来。图像从左至右 每次滚动四列。使用的算法如下：（1）将窗口右边最后四列拷贝到一个仅有四列宽和360行高的小像素映射窗口。（2）将整个显示窗口的内容（减去已经拷贝 的四列）在同一窗口（也就是说，源窗口和目标窗口是同一个窗口）内向右边移动四列。（3）将像素映射窗口的内容拷贝到显示窗口左边头四列。打开文本编辑 器，输入下列命令。以 scroll.pro命名（此文件已经存在于下载的本书配套程序之中）。
; Open a pixmap window 4 columns wide.
Window, 1, /Pixmap, XSize=4, YSize=360
FOR j=0,360/4 DO BEGIN
; Copy four columns on right of display into pixmap.
Device, Copy=[356, 0, 4, 360, 0, 0, 0]
; Make the display window the active window.
WSet, 0
; Move window contents over 4 columns.
Device, Copy=[0, 0, 356, 360, 4, 0, 0]
; Copy pixmap contents into display window on left.
Device, Copy=[0, 0, 4, 360, 0, 0, 1]
ENDFOR
END
运 行程序，输入:
IDL> .RUN scroll
程序每次滚动一次。再次运行程序，输入：
IDL> .Go
能 修改程序让它一直滚动，直到让它停止吗？
在继续下一次练习之前删除像素映射窗口。输入：
IDL> Wdelete,1

Z 图形缓冲区中的图形显示技巧
可以把IDL中的Z图形缓冲区想象成一个三维盒子，3D对象可以被堆积在里面而不用关心它们的“实体”性质。这个盒子 能用16位深度缓冲区来记录每个对象的深度。盒子的一边是投影平面。可以想象光线通过投影平面的每一个像素最终遇到盒子内实体对象。光线遇到的像素值就是 投影到投影平面的值。用这种方法，Z 图形缓冲区就可以处理曲面和线条的自动消隐。图61就是此概念的图解说明。 
图61：Z图形缓冲区可以被想 象作一个能保留深度信息的3D盒子。光线射到Z图形缓冲区中的物体时，它们的像素值反过来被投影到投影平面上。

其概念是，当物体被装入三 维盒子里时，就可以得到投影平面的一幅快照或一幅图像。这是盒子中三维物体的二维投影。处在其它物体后面的物体就不被显示。（这种属性可以在一些IDL图 形命令中用Transparent关键字来改变，以后将会看到。）快照实际上就是利用TVRD命令对投影平面的屏幕转储。
Z图形缓冲区的实现
在 IDL中，Z图形缓冲区是用软件方式作为另一中图形输出设备实现的，类似于PostScript设备或X窗口，Win或Mac设备。因此，要在Z图形缓冲 区中画图形，必须用Sep_Plot命令将其变成当前的图形输出设备。和其它图形输出设备一样，Z图形缓冲区可用Device命令和适当的关键字来配置。
常 用于Z图形缓冲区的两个关键字是Set_Colors和Set_Resolution。这两个关键字定义如下：
Set_Colors Z图形缓冲 区中的颜色数。在缺省情况下，Z图形缓冲区使用256种颜色。在IDL运行中，这是个非常少的颜色数目。如果希望Z图形缓冲区的输出具有和显示设备相同的 颜色数目，就须设置这个关键字。
Set_Resolution Z图形缓冲区投影平面通常设置为640像素宽和480　　　　　　　　　　　　像 素高。如果要在图形窗口显示Z图形缓冲区的输出结果，就应该将Z图形缓冲区的大小设置成图形窗口的大小。例如：
DEVICE, Set_Resolution=[400,400]
一 个Z图形缓冲区实例：两个曲面
要了解Z图形缓冲区是怎样工作的，先按如下创建两个名为peak和saddle的物体。（完成此例的命令可以在下载 的本书配套程序中的文件twosurf.pro中找到。）
IDL>peak=shift(dist(20,16),10,8)
IDL>peak=Exp(-(peak/5)^2)
IDL>saddle=shift(peak,6,0)+shift(peak,-6,0)/2B

要 在Z图形缓冲区中组合两个3D物体，首先看看这两个物体的各自形状。可以让它们在两个窗口中以不同的颜色表显示。首先，在颜色查询表中的不同部位装入蓝色 和红色色谱表：
IDL>colors=!D.Table_Size/2
IDL>LoadCT, 1, ncolors=colors
IDL>LoadCT, 3, ncolors=colors, Bottom=colors-1
　 创建一个窗口，显示第一个物体的阴影曲面图。要注意的是，Set_shading 命令是用来将阴影处理的值限制在颜色查询表中特定部位。键入：
IDL>window, 1, xsize=300, ysize=300
IDL>set_shading, value=[0,coloe-1]
IDL>shade_surf, peak, zrange=[0.0,1.2]
　 　接着将第二物体显示在它自己的显示窗口中。用颜色查询表的不同部位作为阴影处理参数。键入：
IDL>window, 2, xsize=300,ysize=300
IDL>set_shading, value=[colors, 2*colors-1]
IDL>shade_surf, saddle, zrange=[0.0,1.2]
使 Z图形缓冲区成为当前设备
为了在Z图形缓冲区中组合两个物体，必须使Z图形缓冲区成为当前的图形显示设备。这可用Set_Plot命令来实现。 Copy这个关键字可将当前的颜色表复制到Z缓冲区中。保存当前图形显示设备的名字，以便能方便地返回。键入：
IDL>thisDevice=!D.Name
IDL>Set_Plot,’z’,/Copy
配 置Z图形缓冲区
接下来，必须将Z图形设备配置成具体的规格。在这个例子中，要限制颜色的数目，还要使缓冲区的分辨率与当前图形显示窗口尺寸相等。 键入：
IDL>Device,Set_colors=2*colors, Set_Resolution=[300,300]

将 物体装入到Z图形缓冲区中
　　现在，将两个物体装入到Z图形缓冲区中。要注意的是，键入这些命令时，将看不到任何事发生。因为输出已进入内存中的 Z图形缓冲区里面，而不是显示设备。
IDL>Set_shading,values=[0,color-1]
IDL>Shade_surf,peak,zrange=[0.0,1.2]
IDL>Set_shading,values=[colors,2*colors-1]
IDL>shade_surf,saddle,zrange=[0.0,1.2],/noerase
对 投影面进行拍照
接着，对投影表面进行快照。可通过TVRD命令实现。
IDL>picture=TVRD ()
在显示设备上 显示结果
最后，返回到显示设备。打开一个新的窗口来显示“图像”，如下：
IDL>Set_Plot, thisdevice
IDL>Window, 3, xsize=300,ysize=300
IDL>TV, picture

输 出结果应如图62所示。

图62：Z图形缓冲区可以通过曲面的自动消隐来组合3D物体。
Z图形缓冲区的一些奇怪特点
仔细查 看窗口1和3的输出。特别注意观察那些坐标轴的标记。可以注意到窗口3（也就是来自Z图形缓冲区的输出）的轴的标记要稍微大些。由于某种原因，Z图形缓冲 区使用了一种缺省的不同于IDL在显示窗口中显示图形时所用的字符尺寸。
　　这种简单的事实，如果不意识到的话，会导致在Z图形缓冲区中建立一个 3D坐标空间以及组合IDL图形命令时，造成大量的时间浪费。这主要是因为图形边缘的建立是基于缺省的字符尺寸，并且图形边缘在显示设备上和在Z图形缓冲 区中是有差异的。它们差不多是一样的。但就是 “差不多”将使读者焦头烂额。
　　一种笨规则非常有用，就是，当要在Z图形缓冲区中画图时，总是设 置!P.Charsize系统变量。例如：
IDL>!P.Charsize=1.0
　　这里只是提供一个例子，查看图62中的图 示。这幅图上的轴不是在Z图形缓冲区中创建的，因为那时它们将以屏幕分辨率被着色处理（也就是说,作为一幅图像）而且可用PostScript的分辨率对 它们进行着色处理。如果!P.Charsize关键字没有在对阴影曲面进行着色处理和随后坐标轴的添加之前被设置的话，就不可能在最后的输出中使坐标轴线 对应在正确的位置上。
用Z图形缓冲区使图像变形
使用Z图形缓冲区另一种更强大的技术是用于三维数据的切片显示。这是可行的，因为Z图形缓 冲区具有将图像变形到一个多边形平面上的能力。要了解是如何工作的，先用LoadData命令打开80*100*57 的 3D MRI Head Scan数据集。如下：
IDL>head=LoadData(8)
也许希望能在输入代码时开一个日志文件 去截获这些命令，因为这些命令数量很多，而且必须正确无误地找到它们。日志文件将使读者很轻松地改变并重新运行这些命令。（这种日志文件已经创建好了，可 以下载的本书配套程序中找到warping.pro文件。）
IDL>Journal, ’warping.pro’
一般来说，在 IDL中用Size命令可以获得变量的维数和大小。为定义正确的图像平面，需要知道三维的尺寸。在这个例子中，“尺寸值”将比真实的维数值小1，因为要用 这个数值作为数组的索引号。IDL的索引是号从零开始的。
IDL>s=Size(head)
IDL>xs=s[1]-1
IDL>ys=s[2]-1
IDL>zs=s[3]-1
假 若想在这个数据集的中心显示三个正交切片，也就是说通过三维点（40,50,27）。可以用下述方法定义这些点：
IDL>xpt=40
IDL>ypt=50
IDL>zpt=27
接 着，可以构建用来描绘这三个图像切片或平面的各个多边形。此例子中，每个多边形将是有四个点（矩形的四个角）的简单矩形。矩形的每点都要用一个 （x,y,z）三位值描述。换句话说，每个平面都将是3*4的多边形。键入：
IDL>xplane=[ [xpt,0,0],[xpt,0,zs],[xpt,ys,zs],$
[xpt,ys,0] ]
IDL>yplane=[ [0,ypt,0],[0,ypt,zs],[xs,ypt,zs],$
[xs,ypt,0] ]
IDL>zplane=[ [0,0,zpt],[xs,0,zpt],[xs,ys,zpt],$
[0,ys,zpt]]
下 一步是获得与每个图像平面相对应的图像数据。这在IDL中用数组下标很容易做到。
IDL>ximage=head[xpt,*,*]
IDL>yimage=head[*,ypt,*]
IDL>zimage=head[*,*,zpt]
要 注意的是这些图形都是3D图像（其中有一维是1）。所关心的是与每个图像平面相关的2D图像，因此必须用Reform命令将这些3D图像转换成为2D图像 格式。在此例中，Reform命令将图像重新格式化成80*100*1的图像。当最后的一维为1时，IDL将会舍弃它。这里的结果是一个80*100的图 像。
IDL>ximage=reform(ximage)
IDL>yimage=reform(yimage)
IDL>zimage=reform(zimage)
要 正确显示这些图像，需要确保图像数据已正确地缩放到与显示颜色数相匹配的范围内。相对于整个数据集的范围来调整数据是非常重要的。调整数据并装入颜色：
IDL>mindata=min(head,max=maxdata)
IDL>topcolor=!D.table_size-2
IDL>Loadct, 5, ncolors=!D.table_size-1
IDL>Tvlct, 255, 255,255,topcolor+1
IDL>ximage=bytscl(ximage,top=topcolor,max=maxdata,$
min=mindata)
IDL>yimage=byscl(yimage,top=topcolor,max=maxdata,$
min=mindata)
IDL>zimage=byscl(zimage,top=topcolor,max=maxdata,$
min=mindata)
下 一步要准备建立Z图形缓冲区。 Erase命令将删去先前留在缓冲区内的任何内容。此例子中，用白色擦除，以使显示更加清晰。
IDL>thisdevice=!d.name
IDL>Set_Plot,’z’
IDL>device, Set_colors=topcolor, set_resolution=[400,400]
IDL>erase, color=topcolor+1
用 Scale3命令创建3D坐标空间。这里的坐标轴是用每个维度的大小标记的。
IDL>Scale3, XRange=[0,xs], YRange=[0,ys], ZRange=[0,zs]
最 后将要在Z图形缓冲区中对这些切片进行着色处理。为此需要用Polyfill命令。Pattern关键字将用于设置想显示的图像切 片。 Image_Coord关键字包含一系列与多边形的每个顶点相关联的图像坐标。Image_Interp用于当图像被包裹到多边形上时，指定用双线 性插值，而不是在最近邻域内重采样。T3D关键字通过将三维变换矩阵应用到最后输出中，从而确保多边形出现在3D空间内。键入：
IDL>polyfill,xplane,/t3d,pattern=ximage,/image_interp,$
image_coord=[ [0,0],[0,zs],[ys,zs],[ys,0] ]
IDL>polyfill,yplane,/t3d,pattern=yimage,/image_interp,$
mage_coord=[ [0,0],[0,zs],[xs,zs],[xs,0] ]
IDL>polyfill,zplane,/t3d,pattern=zimage,/image_interp,$
image_coord=[ [0,0],[xs,0],[xs,ys],[0,ys] ]
最 后，将投影平面快速拍照，并显示结果。如下：
IDL>picture=tvrd()
IDL>set_Plot,thisdevice
IDL>window,xsize=400,ysize=400
IDL>tv,picture
如 果已打开了日志文件，现在将其关闭：
IDL>Journal
输出结果应类似于图63所示。如果不是，可用文本编辑器修改日志文件中 的代码来解决问题。要重新运行代码，保存文件并键入：
IDL>@warping

图63：在Z图形缓冲区将图像数据变形到平 面中的一个实例
Z图形缓冲区中的透明效果
注意图63中的每个切片的外部边缘有许多黑块。 这不是图像的部分，而是黑色背景噪音 。
Z 图形缓冲区的一个优点是可以在其中运用透明效果。例如，如果将命令PolyFill中的Transparent关键字设置为大约20或25，那么图像中低 于这个值的像素将是透明的。可以键入试试看看什么效果。（也许想启动另一个日志文件。如果已经关闭了上次的日志文件，应给这次赋一个新名。不幸的是，还不 能添加日志文件。给新的日志文件命名为transparent.pro。这个日志文件的备份可以在下载的本书配套程序中找到。）
IDL>journal, ’transparent’
IDL>Set_plot, ’z’
IDL>erase, color=topcolor+1
IDL>polyfill, xplane,/t3d,pattern=ximage,/image_interp,$
Image_coord=[[0,0],[0,zs],[ys,zs],[xs,0]],$
Transparent=25
IDL>polyfill, yplane,/t3d, pattern=yimage,/image_interp,$
Image_coord=[[0,0], [0,zs],[xs,zs],[xs,0]],$
Transparent=25
IDL> polyfill, zplane,/t3d, pattern=zimage,/image_interp,$
Image_coord=[[0,0], [xs,0],[xs,ys],[0,ys]],$
Transparent=25
IDL>picture=tvrd()
IDL>set_plot,thisdevice
IDL>window,/free,xsize=400,ysize=400
IDL>erase,color=topcolor+1
IDL>tv,picture
IDL>journal
如 果输入出了差错并且最后的输出图像看起来不正确，可对日志文件中的错误进行更改，然后重新启动日志文件，键入：
IDL>@transparent 　　
将Z图形缓冲区效果与体数据着色相结合
Z图形缓冲区效果经常和体数据的着色技术相结合，从而使数据的可视化更生动。例如，假设想创建 这套数据集的等值面（一个等值面是一个在任何地方都有相同值的曲面。它就象数据的一个三维等值线图。）启动一个名为isosurface.pro日志文 件。（此日志文件的备份可以在下载的本书配套程序中找到。）
IDL>Journal,’isosurface.pro’
要创建一个 等值面，首先用Shade_Volume命令创建一组描述曲面的顶点和多边形。在下面的命令中，设置Low关键字以便所有比此等值面值大的值将被此等值面 包围。变量vertices和polygons都是输出变量。它们将用在随后的PolyShade命令中来对曲面着色。看一下值为50的头部数据的直方 图，是一个很适合等值级别。键入：
IDL>plot, histogram(head), Max_value=5000
IDL>shade_volume, head, 50,vertices, polygons , /low
等 值面是用PolyShade命令来进行着色处理的。必须确保使用前面建立的3D变换，键入：
IDL>scale3, xrange=[0,xs],yrange=[0,ys],zrange=[0,zs]
IDL>isosurface=polyshade(vertices,polygons,/t3d)
IDL>loadct, 0, ncolors=topcolor+1
IDL>tv,isosurface
现 在，将这个等值面和穿过数据的Ｚ切片组合起来，这个Z切片就是前面在Z图形缓冲区中获得的。注意现在在Ｚ方向上截掉头部数据。键入：
IDL>shade_volume,head(*,*,0:zpt),50,vertices,$
polygons,/low
IDL>isosurface=polyshade(vertices,polygons,/t3d)
IDL>isosurface(where(isosurface eq 0))=topcolor+1
IDL>tvlct, 70, 70, 70, topcolor+1
IDL>set_plot,’z’, /copy
IDL>tv,isosurface
IDL> scale3, xrange=[0,xs],yrange=[0,ys],zrange=[0,zs]
IDL>polyfill, zplane,/t3d,pattern=zimage,/image_interp,$
Image_coord=[[0,0],[xs,0],[xs,ys],[0,ys]],$
Transparent=25
IDL>picture=tvrd()
IDL>set_plot,thisdevice
IDL>tv,picture
输 出结果应类似于图６４所示。如果不是，修改日志文件并重新运行该文件，如下：
IDL>@isosurface

图64：等值 面和数据切片在Z图形缓冲区中的组合。
本 文引用地址： http://www.sciencenet.cn/m/user_content.aspx?id=303593