Blame view

src/idl_misc/coyote_for_Dustemwrap/cgdraw_roi.pro 10.4 KB
6db3528a   Jean-Philippe Bernard   adding librairies...
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
; docformat = 'rst'
;
; NAME:
;   cgDraw_ROI
;
; PURPOSE:
;   This procedure draws a region or group of regions, defined by the IDLanROI or
;   IDLanROIGroup objects, on a direct graphics device. The type of ROI drawn is
;   based on the TYPE property of a given IDLanROI object. These can be points,
;   polylines, or filled polygons.
;
;******************************************************************************************;
;                                                                                          ;
;  Copyright (c) 2012, by Fanning Software Consulting, Inc. All rights reserved.           ;
;                                                                                          ;
;  Redistribution and use in source and binary forms, with or without                      ;
;  modification, are permitted provided that the following conditions are met:             ;
;                                                                                          ;
;      * Redistributions of source code must retain the above copyright                    ;
;        notice, this list of conditions and the following disclaimer.                     ;
;      * Redistributions in binary form must reproduce the above copyright                 ;
;        notice, this list of conditions and the following disclaimer in the               ;
;        documentation and/or other materials provided with the distribution.              ;
;      * Neither the name of Fanning Software Consulting, Inc. nor the names of its        ;
;        contributors may be used to endorse or promote products derived from this         ;
;        software without specific prior written permission.                               ;
;                                                                                          ;
;  THIS SOFTWARE IS PROVIDED BY FANNING SOFTWARE CONSULTING, INC. ''AS IS'' AND ANY        ;
;  EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES    ;
;  OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT     ;
;  SHALL FANNING SOFTWARE CONSULTING, INC. BE LIABLE FOR ANY DIRECT, INDIRECT,             ;
;  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED    ;
;  TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS;         ;
;  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND             ;
;  ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT              ;
;  (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS           ;
;  SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.                            ;
;******************************************************************************************;
;
;+
; This procedure draws a region or group of regions, defined by the IDLanROI or
; IDLanROIGroup objects, on a direct graphics device. The type of ROI drawn is
; based on the TYPE property of a given IDLanROI object. These can be points,
; polylines, or filled polygons. This procedure is a Coyote Graphics wrapper for
; the built-in IDL procedure DRAW_ROI.
;
; :Categories:
;    Graphics
;    
; :Params:
;     roi: in, required, type=object
;         The input ROI object. Must be either an IDLanROI object or an IDLanROIGroup
;         object (if you wish to draw multiple ROIs at the same time).       
;
; :Keywords:
;     addcmd: in, optional, type=boolean, default=0
;        Set this keyword to add the command to an cgWindow. Setting this keyword
;        automatically sets the WINDOW keyword, but the command does not erase the
;        graphics window as it would normally.
;     color: in, optional, type=string, default='opposite'
;         The name of the polygon color. Color names are those used with cgColor. 
;         This value can also be a long integer or a byte or short integer index 
;         into the current color table.
;     device: in, optional, type=boolean, default=0
;         Set to indicate the polygon vertices are in device coordinates, rather than data coordinates.
;     linestyle: in, optional, type=integer, default=0
;         Set this keyword to the normal linestyle graphics keyword values. The default is a solid line.
;     normal: in, optional, type=boolean, default=0
;         Set to indicate the polygon vertices are in normalized coordinates, rather than data coordinates.
;     outline: in, optional, type=boolean, default=0
;         Set this keyword to draw the outline of the ROI, rather than the default filled ROI. In this
;         case cgPlotS will be used to draw the ROI, rather than the built-in IDL procedure Draw_ROI.
;     psym: in, optional, type=integer
;        Any normal IDL PSYM values, plus any value supported by the Coyote Library
;        routine cgSYMCAT. An integer between 0 and 46. May also be specified as a
;        symbol names. See `cgSymCat` for details.
;     symsize: in, optional, type=float/vector, default=1.0
;        A scalar or vector of symbol sizes. Default is 1.0. May be a vector of the same 
;        length as X.
;     t3d: in, optional, type=boolean, default=0
;        Set this keyword to use the 3D coordinate system established by !P.T.
;     thick: in, optional, type=float, default=1.0
;     window: in, optional, type=boolean, default=0
;         Set this keyword to add the command to the current cgWindow application.
;     _ref_extra: in, optional
;         Any `DRAW_ROI <http://www.exelisvis.com/docs/DRAW_ROI_Procedure.html>` keyword 
;         not defined here may be used
;     
;          
; :Examples:
;    Download `cgDraw_ROI_Example <http://www.idlcoyote.com/tip_examples/cgdraw_roi_example.pro>'
;    for an example program illustrating how cgDraw_ROI can be used.
;       
; :Author:
;    FANNING SOFTWARE CONSULTING::
;       David W. Fanning 
;       1645 Sheely Drive
;       Fort Collins, CO 80526 USA
;       Phone: 970-221-0438
;       E-mail: david@idlcoyote.com
;       Coyote's Guide to IDL Programming: http://www.idlcoyote.com
;
; :History:
;     Change History::
;        Written, 11 November 2012. David W. Fanning.
;
; :Copyright:
;     Copyright (c) 2012, Fanning Software Consulting, Inc.
;-
PRO cgDraw_ROI, roi, $
    ADDCMD=addcmd, $
    COLOR=color, $
    DEVICE=device, $
    LINESTYLE=linestyle, $
    NORMAL=normal, $
    OUTLINE=outline, $
    PSYM=psym, $
    SYMSIZE=symsize, $
    T3D=t3d, $
    THICK=thick, $
    WINDOW=window, $
     _REF_EXTRA=extra

    Compile_Opt idl2

    Catch, theError
    IF theError NE 0 THEN BEGIN
        Catch, /CANCEL
        void = cgErrorMsg()
        IF N_Elements(currentState) NE 0 THEN cgSetColorState, currentState
        RETURN
    ENDIF

    ; Did user pass a parameter?
    IF (N_Params() EQ 0) THEN BEGIN
        Print, 'USE SYNTAX: cgDraw_ROI, thisROI
        RETURN
    ENDIF
    
    ; Do they want this plot in a resizeable graphics window?
    IF Keyword_Set(addcmd) THEN window = 1
    IF Keyword_Set(window) AND ((!D.Flags AND 256) NE 0) THEN BEGIN
    
        currentWindow = cgQuery(/CURRENT, COUNT=wincnt)
        IF wincnt EQ 0 THEN replaceCmd = 0 ELSE replaceCmd=1
        
        ; If adding a command, have to do this differently.
        IF Keyword_Set(addcmd) THEN BEGIN
           cgWindow, 'cgDraw_ROI', roi, $
              COLOR=color, $
              DEVICE=device, $
              LINESTYLE=linestyle, $
              NORMAL=normal, $
              PSYM=psym, $
              SYMSIZE=symsize, $
              T3D=t3d, $
              THICK=thick, $
              ADDCMD=1, $
              _EXTRA=extra
                
            RETURN
            
        ENDIF ELSE BEGIN
        
           ; Otherwise, we are just replacing the commands in a new or existing window.
           cgWindow, 'cgDraw_ROI', roi, $
              COLOR=color, $
              DEVICE=device, $
              LINESTYLE=linestyle, $
              NORMAL=normal, $
              PSYM=psym, $
              SYMSIZE=symsize, $
              T3D=t3d, $
              THICK=thick, $
              REPLACECMD=replaceCmd, $
              _Extra=extra
            
           RETURN
         
         ENDELSE
    ENDIF
    
  ; Set up PostScript device for working with colors.
    IF !D.Name EQ 'PS' THEN Device, COLOR=1, BITS_PER_PIXEL=8
    
    ; Do this in decomposed color, if possible.
    cgSetColorState, 1, CURRENTSTATE=currentState
    
    ; Need a color?
    thisColor = cgDefaultColor(color, DEFAULT='opposite')
    
    ; Get the current color vectors.
    TVLCT, rr, gg, bb, /Get
    
    ; Do you need a window?
    IF ((!D.Flags AND 256) NE 0) && (!D.Window LT 0) THEN cgDisplay
    
    ; Is the color a string? 
    IF Size(thisColor, /TNAME) EQ 'STRING' THEN thisColor = cgColor(thisColor)
    IF N_Elements(psym) EQ 0 THEN psym = 0
        
    ; Draw the ROI
    IF Keyword_Set(outline) THEN BEGIN
        IF Obj_Isa(roi, 'IDLanROI') THEN BEGIN
            roi -> GetProperty, DATA=data
            x = Reform(data[0,*])
            y = Reform(data[1,*])
            cgPlots, x, y, $
                COLOR=thisColor, $
                DEVICE=device, $
                LINESTYLE=linestyle, $
                NORMAL=normal, $
                PSYM=cgSymCat(psym), $
                SYMSIZE=symsize, $
                T3D=t3d, $
                THICK=thick
        ENDIF ELSE BEGIN
            objs = roi -> Get(/ALL, COUNT=count)
            FOR j=0,count-1 DO BEGIN
            thisROI = objs[j]
            thisROI -> GetProperty, DATA=data
            x = Reform(data[0,*])
            y = Reform(data[1,*])
            cgPlots, x, y, $
                COLOR=thisColor, $
                DEVICE=device, $
                LINESTYLE=linestyle, $
                NORMAL=normal, $
                PSYM=cgSymCat(psym), $
                SYMSIZE=symsize, $
                T3D=t3d, $
                THICK=thick
            ENDFOR
        ENDELSE
    ENDIF ELSE BEGIN
        Draw_ROI, roi, $
          COLOR=thisColor, $
          DEVICE=device, $
          LINESTYLE=linestyle, $
          NORMAL=normal, $
          PSYM=cgSymCat(psym), $
          SYMSIZE=symsize, $
          T3D=t3d, $
          THICK=thick, $
          _Extra=extra
    ENDELSE
    
    ; Clean up.
    cgSetColorState, currentState
    IF !D.Name NE 'Z' THEN TVLCT, rr, gg, bb
   
END