Changeset 35 for trunk/ToBeReviewed/PLOTS/DESSINE/tvplus.pro

Timestamp:

05/02/06 16:44:47 (18 years ago)

Author:

pinsard

Message:

upgrade of PLOTS/DESSINE according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/ : files

File:

• trunk/ToBeReviewed/PLOTS/DESSINE/tvplus.pro (copied) (copied from trunk/PLOTS/DESSINE/tvplus.pro) (2 diffs)

trunk/ToBeReviewed/PLOTS/DESSINE/tvplus.pro

@@ -5 +5 @@
 ; NAME: tvplus
 ;
-; PURPOSE: tvscl ameliore, permet de visualiser un tableau sans interpolation
-;
-; CATEGORY: pour regarder un tableau au point par point, chaque point
-; etant represente par un carre de coefzoom pixels de cote. la couleur
-; du carre, representant sa valeur, est une interpolation lineaire
-; entre le min du tableau de couleur 0 et le max du tableau de couleur
-; !d.n_colors-1 <255
-;
-; CALLING SEQUENCE:tvplus,z2d [,coefzoom]
-;
-; INPUTS: z2d le tableau 2d a visualiser
-;
-;         coefzoom : (optionnel) c''est la taille des cotes des carres
-;         representant un point, mesure en pixel. si cet argument
-;         n''est pas specifie, la taille des carres est calculee pour
-;         que la fenetre soit au plus pres de 95% de la taille de l''ecran
+; PURPOSE: enhanced version of tvscl
+;
+; CATEGORY: quick exploration of 2D arrays 
+;
+; CALLING SEQUENCE: tvplus,z2d [,cellsize]
+;
+; INPUTS: 
+;
+;    z2d: 2D array to visualize
+;
+;    cellsize: (optional) this is the size (in pixel) of the square
+;    representing 1 array element. By default, this size is computed
+;    automatically in order that the size of the plotting window do
+;    not exceed the screen size. If the user specify a large value
+;    of cellsize that forces tvplus to create a window larger than
+;    the screen, a "scrolling window" will be displayed instead of a
+;    regular window. Unfortunately the nice fonctionnalities of tvplus
+;    are not coded for "scrolling window" case...
 ;
 ; KEYWORD PARAMETERS:
@@ -27 +29 @@
 ;       the bar. default is 0.
 ;
-;       C_NAN: c'est le numero de la couleur des points a
-;       !values.f_nan. par defaut c''est !d.n_colors< 256 -1
-;
-;       C_MASK: c'est le numero de la couleur des points a
-;       la valeur specifiee par MASK. par defaut c''est 0
-;
-;       OFFSET: c''est un vecteur de 2 elements donnant la valeur dont
-;       il faut decaler la numerotation des points en x et en y.
-;       Par ex:
-;          tvplus,sst[x1:x2,y1:y2],offest=[x1,y1]
-;
-;       MASK: la valeur qui sert de mask a l''ecran. si abs(mask) lt
-;       1e6 les points egaux a mask seronts en noir, dans le cas
-;       contraire ce seront les points ge que abs(mask)/10. Ceci est
-;       nessaire pour eviter les erreurs d''arrondi. 
-;
-;       MIN et MAX permet de borner z2d a un min ou a un max, les
-;       points inferieurs a min etant de la couleur 0 (noir) et ceux
-;       superieurs a max etant de la couleur (!d.n_colors-1) < 255 (blanc)
+;       C_NAN: The color number that should be used for the NaN values.
+;       default value is !d.n_colors < 255
+;
+;       C_MASK: The color number that should be used for the mask values.
+;       default value is 0
+;
+;       OFFSET: 2 elements vector used by tvplus itself when shoing zoom.
+;       It is used to shift the ranges of xaxis and yaxis.
+;       For example: tvplus,sst[x1:x2,y1:y2],offest=[x1,y1]
+;
+;       MASK: The mask value. Note that if abs(mask) < 1.e6, then the
+;       exact value of the mask is used to find the maskwd point.
+;       if abs(mask) > 1.e6, the test to find the masked value is ge
+;       abs(mask)/10. This is necessary to avoid the rounding errors
+;
+;       MIN and MAX: scalars used to specify the min and max values
+;       of the color bar. default is from 0 to !d.n_colors < 255
 ;
 ;       NCOLORS:  This is the number of colors in the color
 ;       bar. default is !d.n_colors<255 -1
 ;
-;       NOINTERP: activer si on ne veut pas que les valeurs du tableau
-;       soient interpoles entre 0 et !d.n_colors-1 < 255 pour obtenir
-;       les correspondances avec les couleurs. Utile qd on travaille
-;       avec des tableaux de type byte
-;
-;       WINDOW: numero de la fenetre ou l''on veut faire la visualisation (permet
-;       d''ouvrir plusieurs fenetres). par defaut on ouvre 'IDL0'
-;
-; OUTPUTS: affichage de parametres selectionnes a la souris, cf. apres
-;
-; COMMON BLOCKS:
+;       NOINTERP: used this keyword if you don't whant that the values
+;       are interpolated between 0 (or min) and !d.n_colors < 255 (or max).
+;       This can be for example usefull when working on byte type arrays.
+;
+;       NOUSEINFOS: activite to supress the printed message explaining
+;       how to use tvplus
+;
+;       WINDOW: Number of the window used to display the array values.
+;       default is window number 0.
+;
+; OUTPUTS: no (except the printed finformations, see below)
+;
+; COMMON BLOCKS: no
 ;
 ; SIDE EFFECTS:
-;       1) qd la fenetre qui doit etre ouverte a une dimension -imposee
-;       par la presence de coefzoom- superieure a 75% de la taille de
-;       l''ecran, c''est une fenetre avec un scrolling bar qui est
-;       ouverte.
-;
-;       2) qd le dessin est realise ds une fenetre normale, le
-;       programme s''arrete apres le trace pour faire des l'analyses au
-;       point par point a l''aide de la souris:
-;          Bouton de gauche: position du curseur et valeur du tableau en ce point
-;          Bouton du milieu: zoom
-;          Bouton de droite: quitter
-;
-; RESTRICTIONS: l''analyse a la souris ne marche pas ds les fenetres
-; avec une scrolling bar...
-;
-; EXAMPLE:  IDL> tvplus,nlec('sst'),max=30
+;
+; use your mouse to scan the array values...
+;     left button  : mouse position and associated array value
+;     middle button: use it twice to define a zoom box
+;     right button : quit
+;
+; RESTRICTIONS: the nice fonctionnalities of tvplus are not coded
+;               for "scrolling window" case...
+;
+; EXAMPLE:  IDL> tvplus, dist(100)
 ;
 ; MODIFICATION HISTORY: Sebastien Masson (smasson@lodyc.jussieu.fr)
 ;                       18/12/98
+; Aug 2005: quick cleaning + english
 ;-
 ;------------------------------------------------------------
 ;------------------------------------------------------------
 ;------------------------------------------------------------
-PRO tvplus, z2d, coefzoom, BOTTOM = bottom, C_MASK = c_mask, C_NAN = c_nan, WINDOW = window $
-            , MIN = min, MAX = max, MASK = mask, OFFSET = offset, PASBLABLA = PASBLABLA $
+PRO tvplus, z2d, cellsize, BOTTOM = bottom, C_MASK = c_mask, C_NAN = c_nan, WINDOW = window $
+            , MIN = min, MAX = max, MASK = mask, OFFSET = offset, NOUSEINFOS = NOUSEINFOS $
             , NCOLORS = ncolors, NOINTERP = nointerp, _EXTRA = ex
 ;
-   tab = float(z2d)
-;------------------------------------------------------------
-; on verifie la taille de tab
-;------------------------------------------------------------
-   if (size(tab))[0] NE 2 then begin
-      ras = report('le tableau doit etre 2D')
-      return
-   endif
-;------------------------------------------------------------
-; on definit au besoin ncoulmax, bottom, topcoul et ncolors
-;------------------------------------------------------------
-   ncoulmax = !d.n_colors < 256
-   IF N_ELEMENTS(bottom) EQ 0 THEN bottom = 0
-   if NOT keyword_set(ncolors) then ncolors = ncoulmax-bottom
-   topcoul = (bottom+ncolors-1) < (ncoulmax-1)
-;------------------------------------------------------------
-; on fait le menage
-;------------------------------------------------------------
-   xenvsauve = !x & yenvsauve = !y & penvsauve = !p
-   reinitplt, /z, /invert
-;------------------------------------------------------------
-; le tableau contient des valeurs a !values.f_nan ???
-;------------------------------------------------------------
-   nan = total(finite(z2d,/nan)) < 1
-   if keyword_set(nan) then begin
-      nanindex = where(finite(z2d,/nan) EQ 1)
-      tab[nanindex] = min(tab, /nan)
-   endif
-;------------------------------------------------------------
-; on calcule coefzoom ou bien on verifit qu''il n''est pas trop grand
-;------------------------------------------------------------
-   dimensions = GET_SCREEN_SIZE()
-   if n_elements(coefzoom) EQ 0 then BEGIN 
-      coefzoom = min(floor(dimensions/(size(z2d))[1: 2]*.85))
-   ENDIF ELSE $ 
-; il faudra ouvrir une fenetre avec une scrolling barre
-   if coefzoom GE min(floor(dimensions/(size(z2d))[1: 2]*.85)) then scrolling = 1
-   if coefzoom LT 1 then begin
-      coefzoom = 1
-      scrolling = 1
-   endif
-;------------------------------------------------------------
-; on met les valeurs masques au max (sur les points non masques) du tableau
-;------------------------------------------------------------
-   if n_elements(mask) then BEGIN 
-      if abs(mask) LT 1e6 then BEGIN 
-         terre = where(tab EQ mask)
-         if terre[0] NE -1 then tab(terre) = min(tab(where(tab NE mask)))
-      ENDIF ELSE BEGIN 
-         terre = where(abs(tab) GE abs(mask)/10.)
-         if terre[0] NE -1 then tab(terre) = min(tab(where(abs(tab) LT abs(mask)/10.)))
-      ENDELSE
-   ENDIF ELSE terre = -1
-;------------------------------------------------------------
-; on borne le tableau a min et max s''il sont declares
-;------------------------------------------------------------
-   if keyword_set(min) then BEGIN 
-      tab = min > tab
-      vraimin = min
-   ENDIF ELSE vraimin = min(tab)
-   if keyword_set(max) then BEGIN 
-      tab = tab < max
-      vraimax = max
-   ENDIF ELSE vraimax = max(tab)
-;------------------------------------------------------------
-; on bidoulle tab 
-;------------------------------------------------------------
-   if NOT keyword_set(nointerp) then BEGIN
-; on interpole tab entre bottom et bottom+ncolors-1
-      m=1.*(ncolors-1)/(vraimax-vraimin)
-      p=bottom-1.*vraimin*m
-      tab=round(1.*m*tab+p) 
-   endif
-; on met les point nan a la couleur numero c_nan
-   if keyword_set(nan) then begin
-      if keyword_set(c_nan) THEN tab[nanindex] = c_nan <  (ncoulmax -1) $
-      ELSE tab[nanindex] = topcoul
-   endif
-; on met les point terre a la couleur numero c_mask
-   if keyword_set(c_mask) AND terre[0] NE -1 THEN $
-    tab[terre] = c_mask < (ncoulmax -1)
-; on passe le tableau en byte
-   tab = byte(tab)
-; on augmente la taille de tab pour pouvoir faire le tv
-   tab = congridseb(tab, (size(tab))[1]*coefzoom,(size(tab))[2]*coefzoom)
-;------------------------------------------------------------
-; on ouvre une fenetre a la bonne taille et on dessine
-;------------------------------------------------------------
-   nx = (size(tab))[1]
-   ny = (size(tab))[2]
-; marges en pixel
-   rapportaspect = 1.*nx/ny
-   if rapportaspect GE 1 THEN margepix = 1.*[25, 25, 75, 25] ELSE margepix = 1.*[25, 100, 25, 25]
-;
-   if n_elements(scrolling) EQ 0 then BEGIN ; on ouvre une fenetre classique
-      if NOT keyword_set(window) then window = 0
-      window,window, xsize = nx+margepix[0]+margepix[1] $
-       , ysize = ny+margepix[2]+margepix[3] ; ouverture de la fenetre
-; qd on utilise des couleurs codees sur 24 bit, je n''arrive pas a stipuler la couleur du
-; fond d''une fenetre a l'aide de !p.background, je suis oblige de faire cette bidouille!!!
-      if !d.n_colors gt 256 then BEGIN ; alors  les couleurs sont codees sur 24 bit
-         device, decomposed=1
-         !p.background='ffffff'x
-         plot,[0],[0], xstyle = 5, ystyle = 5
-         device, decomposed=0
-      endif
-      tv, tab,margepix[0], margepix[2], _EXTRA = ex ; dessin de tab interpole entre 0 et !d.n_colors-1 < 255
-;
-; mise en place des axes et du cadre
-;
-; on recupere en coordonnees normalises la position du tv et donc
-; cette du cadre
-      poscadre = fltarr(4)
-      poscadre[0] = margepix[0]/!d.x_size
-      poscadre[2] = 1.-margepix[1]/!d.x_size
-      poscadre[1] = margepix[2]/!d.y_size
-      poscadre[3] = 1-margepix[3]/!d.y_size
-;    on trace le cadre a l''aide d''un plot
-      if NOT keyword_set(offset) then offset = [0, 0]
-      plot, [0], [0],/noerase, position = poscadre, color = 0 $
-       , xstyle = 1, ystyle = 1, xticklen = 1, yticklen = 1 $
-       , xrange = 1.*[0, nx]/coefzoom-.5+offset[0], yrange = 1.*[0, ny]/coefzoom-.5+offset[1]
-;
-; mise en place de la colorbar
-;
-      if rapportaspect ge 1 then $
-       posbar = [25./!d.x_size, 25./!d.y_size, 1-25./!d.x_size, 50./!d.y_size] $
-      ELSE posbar = [1-75./!d.x_size, 25./!d.y_size, 1-50./!d.x_size, 1-25./!d.y_size]
-   if keyword_set(nointerp) then BEGIN & vraimin = 0 & vraimax = ncoulmax & endif
-      colorbar, min = vraimin, max = vraimax, division = 10, cb_color = 0, position = posbar $
-       , vertical = rapportaspect lt 1, /right, BOTTOM = bottom, NCOLORS = ncolors, _EXTRA = ex
-   ENDIF ELSE BEGIN 
-;------------------------------------------------------------
-; dessin ds une fenetre avec une scrolling barre
-;------------------------------------------------------------
-; qd on utilise des couleurs codees sur 24 bit, je n''arrive pas a stipuler la couleur du
-; fond d''une fenetre a l'aide de !p.background, je suis oblige de faire cette bidouille!!!
-      if !d.n_colors gt 256 then begin
-         window, /pixmap
-         device, decomposed=1
-         !p.background='ffffff'x
-         plot,[0],[0]
-         device, decomposed=0
-      endif
-      slide_image, tab $        ; on le dessine ds une fenetre avec une scrolling bar
-       , xsize = nx, ysize = ny $
-       ,xvisible=round(.7*dimensions[0]) < (size(tab))[1] $
-       ,yvisible=round(.7*dimensions[1]) < (size(tab))[2],/register,congrid=0,show_full=0
-      return
-   ENDELSE 
-;------------------------------------------------------------
-; on prend des infos a la souris
-;------------------------------------------------------------
-   if NOT keyword_set(pasblabla) then begin
-      print, 'Bouton de gauche: position du curseur et valeur du tableau en ce point'
-      print, 'Bouton du milieu: zoom'
-      print, 'Bouton de droite: quitter'
-   endif
-   cursor,x,y,/device, /down
-   while (!mouse.button ne 4) do BEGIN
-      x = x-margepix[0] & x = 0 > floor(x/coefzoom) < ((size(tab))[1]-1)
-      y = y-margepix[2] & y = 0 > floor(y/coefzoom) < ((size(tab))[2]-1)
-      case !mouse.button of
-         1:BEGIN                ; saisie d''une valeur
-            if x LE nx/coefzoom AND y LE ny/coefzoom then begin
-               print,'position: x='+strtrim(x+offset[0],1)+', y= '+strtrim(y+offset[1],1) $
-                +' valeur: '+strtrim(float(z2d[x,y]), 1)
-            ENDIF
-            cursor,x,y,/device, /down
-         END
-         2:BEGIN                ; zoom
-            cursor,x2,y2,/device, /down
-            x2 = x2-margepix[0] & x2 =  0 > floor(x2/coefzoom) < ((size(tab))[1]-1)
-            y2 = y2-margepix[2] & y2 =  0 > floor(y2/coefzoom) < ((size(tab))[2]-1)
-            x =  [x, x2] & x =  x[sort(x)]
-            y =  [y, y2] & y =  y[sort(y)]
-            tvplus, z2d[x[0]:x[1],y[0]:y[1] ], WINDOW = window, MIN = min, MAX = max $
-             , MASK = mask, C_MASK = c_mask, C_NAN = c_nan, /PASBLABLA, OFFSET = [x[0],y[0]] $
-             , NCOLORS = ncolors, NOINTERP = nointerp, BOTTOM = bottom, _EXTRA = ex
-            return
-         END
-      endcase
-   ENDWHILE
-;------------------------------------------------------------
-   !x = xenvsauve & !y = yenvsauve & !p = penvsauve
-;------------------------------------------------------------
-   return
+  IF n_elements(z2d) EQ 0 THEN return
+  arr = reform(float(z2d))
+;------------------------------------------------------------
+; check the size of the input array 
+;------------------------------------------------------------
+  if (size(arr))[0] NE 2 then begin
+    ras = report('Input array must have only 2 dimensions and not '+ strtrim(size(arr, /n_dimensions), 1))
+    return
+  endif
+;------------------------------------------------------------
+; def of ncolmax, bottom, topcol et ncolors
+;------------------------------------------------------------
+  ncolmax = !d.n_colors < 256
+  IF N_ELEMENTS(bottom) EQ 0 THEN bottom = 0
+  if NOT keyword_set(ncolors) then ncolors = ncolmax-bottom
+  topcol = (bottom+ncolors-1) < (ncolmax-1)
+;------------------------------------------------------------
+; get default values of !x, !y, !z
+;------------------------------------------------------------
+;   xenvsauve = !x & yenvsauve = !y & penvsauve = !p
+  reinitplt, /z, /invert
+;------------------------------------------------------------
+; Do we have NaN values in arr???
+;------------------------------------------------------------
+  nan = total(finite(z2d, /nan)) < 1
+  if keyword_set(nan) then begin
+    nanindex = where(finite(z2d, /nan) EQ 1)
+    arr[nanindex] = min(arr, /nan)
+  endif
+;------------------------------------------------------------
+; Compute the size (in pixel) of the square representing 1
+; point of the input array 
+;------------------------------------------------------------
+  dimensions = GET_SCREEN_SIZE()
+  if n_elements(cellsize) EQ 0 then BEGIN 
+    cellsize = min(floor(dimensions/(size(z2d))[1: 2]*.75))
+  ENDIF ELSE $ 
+; we need to use a scrolling bar window
+  if cellsize GE min(floor(dimensions/(size(z2d))[1: 2]*.75)) then scrolling = 1
+  if cellsize LT 1 then begin
+    cellsize = 1
+    scrolling = 1
+  endif
+;------------------------------------------------------------
+; Change the value of the masked value for the min of the non-masked values
+;------------------------------------------------------------
+  if n_elements(mask) then BEGIN 
+    if abs(mask) LT 1e6 then BEGIN 
+      masked = where(arr EQ mask)
+      if masked[0] NE -1 then arr(masked) = min(arr(where(arr NE mask)))
+    ENDIF ELSE BEGIN 
+      masked = where(abs(arr) GE abs(mask)/10.)
+      if masked[0] NE -1 then arr(masked) = min(arr(where(abs(arr) LT abs(mask)/10.)))
+    ENDELSE
+  ENDIF ELSE masked = -1
+;------------------------------------------------------------
+; apply min/max keywords
+;------------------------------------------------------------
+  if n_elements(min) NE 0 then BEGIN 
+    arr = min > arr
+    truemin = min
+  ENDIF ELSE truemin = min(arr)
+  if n_elements(max) NE 0 then BEGIN 
+    arr = arr < max
+    truemax = max
+  ENDIF ELSE truemax = max(arr)
+;
+  IF truemin EQ truemax THEN BEGIN
+    dummy = report('constant value everywhere: '+ strtrim(truemin, 1))
+    return
+  ENDIF
+;------------------------------------------------------------
+; apply other keywords (nointerp, c_nan, c_mask) 
+;------------------------------------------------------------
+  if NOT keyword_set(nointerp) then BEGIN
+; interpolation between bottom and bottom+ncolors-1
+    m = 1.*(ncolors-1)/(truemax-truemin)
+    p = bottom-1.*truemin*m
+    arr = round(1.*m*arr+p) 
+  endif
+; set c_nan for NaN values 
+  if keyword_set(nan) then begin
+    if n_elements(c_nan) NE 0 THEN arr[nanindex] = c_nan <  (ncolmax -1) $
+    ELSE arr[nanindex] = topcol
+  endif
+; c_mask for masked values
+  if n_elements(c_mask) NE 0 AND masked[0] NE -1 THEN $
+    arr[masked] = c_mask < (ncolmax -1)
+; use byte type to save memory
+  arr = byte(arr)
+; increase the size of the array in order to be displayed
+; with the suitable size 
+  arr = congridseb(arr, (size(arr))[1]*cellsize, (size(arr))[2]*cellsize)
+;------------------------------------------------------------
+; open a window with the correct size
+;------------------------------------------------------------
+  nx = (size(arr))[1]
+  ny = (size(arr))[2]
+; margin size (in pixel)
+  xyaspect = 1.*nx/ny
+  if xyaspect GE 1 THEN marginpix = 1.*[25, 25, 75, 25] ELSE marginpix = 1.*[25, 100, 25, 25]
+;
+  if n_elements(scrolling) EQ 0 then BEGIN ; open the regular window
+    if NOT keyword_set(window) then window = 0
+    window, window, xsize = nx+marginpix[0]+marginpix[1] $
+            , ysize = ny+marginpix[2]+marginpix[3] 
+; for 24 bits colors, make sure thate the background color is the good one...
+    if !d.n_colors gt 256 then BEGIN
+      device, decomposed = 1
+      !p.background = 'ffffff'x
+      plot, [0], [0], xstyle = 5, ystyle = 5
+      device, decomposed = 0
+    endif
+    tv, arr, marginpix[0], marginpix[2], _EXTRA = ex
+;
+; axis and plot frame 
+;
+; get the normalized position of the tv (we just done above)
+; to know where the frame should be drawn
+    poscadre = fltarr(4)
+    poscadre[0] = marginpix[0]/!d.x_size
+    poscadre[2] = 1.-marginpix[1]/!d.x_size
+    poscadre[1] = marginpix[2]/!d.y_size
+    poscadre[3] = 1-marginpix[3]/!d.y_size
+; Use plot to draw the frame
+    if NOT keyword_set(offset) then offset = [0, 0]
+    !p.position = poscadre
+    plot, [0], [0], /nodata, /noerase, position = poscadre, color = 0 $
+          , xstyle = 1, ystyle = 1, xticklen = 1, yticklen = 1 $
+          , xrange = 1.*[0, nx]/cellsize-.5+offset[0], yrange = 1.*[0, ny]/cellsize-.5+offset[1]
+    xenvsauve = !x & yenvsauve = !y & penvsauve = !p
+;
+; draw the colorbar 
+;
+    IF truemin ne truemax THEN BEGIN
+      if xyaspect ge 1 then $
+        posbar = [25./!d.x_size, 25./!d.y_size, 1-25./!d.x_size, 50./!d.y_size] $
+        ELSE posbar = [1-75./!d.x_size, 25./!d.y_size, 1-50./!d.x_size, 1-25./!d.y_size]
+    if keyword_set(nointerp) then BEGIN & truemin = 0 & truemax = ncolmax & endif
+      colorbar, min = truemin, max = truemax, division = 10, cb_color = 0, position = posbar $
+                , vertical = xyaspect lt 1, /right, BOTTOM = bottom, NCOLORS = ncolors, _EXTRA = ex
+    ENDIF
+;      !p.position = poscadre
+  ENDIF ELSE BEGIN 
+;------------------------------------------------------------
+; scrolling bar window case...
+;------------------------------------------------------------
+; for 24 bits colors, make sure thate the background color is the good one...
+    if !d.n_colors gt 256 then begin
+      window, /pixmap
+      device, decomposed = 1
+      !p.background = 'ffffff'x
+      plot, [0], [0]
+      device, decomposed = 0
+    endif
+    slide_image, arr $          ; on le dessine ds une fenetre avec une scrolling bar
+                 , xsize = nx, ysize = ny $
+                 , xvisible = round(.7*dimensions[0]) < (size(arr))[1] $
+                 , yvisible = round(.7*dimensions[1]) < (size(arr))[2], /register, congrid = 0, show_full = 0
+    return
+  ENDELSE 
+;------------------------------------------------------------
+; Use the mouse to get nice functionalities
+;------------------------------------------------------------
+; format to print the mouse position
+  CASE 1 OF
+    nx LT 10:fmt1 = '(i1)'
+    nx LT 100:fmt1 = '(i2)'
+    nx LT 1000:fmt1 = '(i3)'
+    nx LT 10000:fmt1 = '(i4)'
+    ELSE:fmt1 = ''
+  ENDCASE
+  CASE 1 OF
+    ny LT 10:fmt2 = '(i1)'
+    ny LT 100:fmt2 = '(i2)'
+    ny LT 1000:fmt2 = '(i3)'
+    ny LT 10000:fmt2 = '(i4)'
+    ELSE:fmt2 = ''
+  ENDCASE
+;
+  if NOT keyword_set(nouseinfos) then begin
+    print, 'left button  : mouse position and associated array value'
+    print, 'middle button: use it twice to define a zoom box'
+    print, 'right button : quit'
+  endif
+  cursor, x, y, /device, /down
+  x = x-marginpix[0] & x = 0 > floor(x/cellsize) < ((size(arr))[1]/cellsize-1)
+  y = y-marginpix[2] & y = 0 > floor(y/cellsize) < ((size(arr))[2]/cellsize-1)
+  while (!mouse.button ne 4) do BEGIN
+    case !mouse.button of
+      0:return
+      1:BEGIN                   ; get value
+        if x LE nx/cellsize AND y LE ny/cellsize then begin
+          print, '(x, y) = (' + string(x+offset[0], format = fmt1) $
+                 + ', ' + string(y+offset[1], format = fmt2) $
+                 + '), value = '+strtrim(float((reform(z2d))[x, y]), 1)
+        ENDIF
+        cursor, x, y, /device, /down
+        x = x-marginpix[0] & x = 0 > floor(x/cellsize) < ((size(arr))[1]/cellsize-1)
+        y = y-marginpix[2] & y = 0 > floor(y/cellsize) < ((size(arr))[2]/cellsize-1)
+      END
+      2:BEGIN                   ; zoom
+        cursor, x2, y2, /device, /down
+        x2 = x2-marginpix[0] & x2 =  0 > floor(x2/cellsize) < ((size(arr))[1]/cellsize-1)
+        y2 = y2-marginpix[2] & y2 =  0 > floor(y2/cellsize) < ((size(arr))[2]/cellsize-1)
+        x =  [x, x2] & x =  x[sort(x)]
+        y =  [y, y2] & y =  y[sort(y)]
+        IF keyword_set(OFFSET) THEN offset = [x[0], y[0]]+offset ELSE offset = [x[0], y[0]]
+        tvplus, z2d[x[0]:x[1], y[0]:y[1] ], WINDOW = window, MIN = min, MAX = max $
+                , MASK = mask, C_MASK = c_mask, C_NAN = c_nan, /NOUSEINFOS, OFFSET = OFFSET $
+                , NCOLORS = ncolors, NOINTERP = nointerp, BOTTOM = bottom, _EXTRA = ex
+        return
+      END
+      ELSE:
+    endcase
+  ENDWHILE
+;------------------------------------------------------------
+  !x = xenvsauve & !y = yenvsauve & !p = penvsauve
+  !x.range = 1.*[0, nx]/cellsize-.5+offset[0]
+  !y.range = 1.*[0, ny]/cellsize-.5+offset[1]
+;------------------------------------------------------------
+  return
 end