;+ ; PROJECT: ; Solar-B/EIS ; ; NAME: ; EIS_ITOOL_PTOOL ; ; PURPOSE: ; Make a widget base for pointing for an appropriate SOHO instrument ; ; EXPLANATION: ; This routine creates a widget interface based on a given "parent" ; base in IMAGE_TOOL for the pointing purpose. Events genereated from ; this widget set is handled by routine PT_PTOOL_EVENT. It requies a ; pointing structure be defined (in the common block included in ; image_tool_com.pro) for a given instrument. This pointing structure ; should include the following tags: ; ; MESSENGER - ID of widget in the caller that triggers a ; timer event in the planning tool to signal the ; completion of pointing; must be a widget that ; does not usually generate any event ; INSTRUME - Code specifying the instrument; e.g., 'C' for CDS ; SCI_SPEC - Science specification ; STD_ID - Study ID ; G_LABEL - Generic label for the pointing; e.g., 'RASTER' ; X_LABEL - Label for X coordinate of pointing; e.g., 'INS_X' ; Y_LABEL - Label for Y coordinate of pointing; e.g., 'INS_Y' ; DATE_OBS - Date/time of beginning of observation, in TAI format ; DO_POINTING- An integer of value 0 or 1 indicating whether pointing ; should be handled at the planning level (i.e., by ; IMAGE_TOOL) ; N_POINTINGS- Number of pointings to be performed by IMAGE_TOOL ; POINTINGS - A structure array (with N_POINTINGS elements) of type ; "DETAIL_POINT" to be handled by IMAGE_TOOL. It has the ; following tags: ; ; POINT_ID - A string scalar for pointing ID ; INS_X - X coordinate of pointing area center in arcs ; INS_Y - Y coordinate of pointing area center in arcs ; WIDTH - Area width (E/W extent) in arcsec ; HEIGHT - Area height (N/S extent) in arcsec ; OFF_LIMB - An interger with value 1 or 0 indicating ; whether or not the pointing area should ; be off limb ; ; N_RASTERS - Number of rasters for each pointing (this is ; irrelevant to the SUMER) ; RASTERS - A structure array (N_RASTERS-element) of type ; "RASTER_POINT" that contains raster size and pointing ; information (this is irrelevant to the SUMER). It has ; the following tags: ; ; POINTING - Pointing handling code; valis ; values are: 1, 0, and -1 ; INS_X - Together with INS_Y, the pointing to use ; when user-supplied values are not ; allowed. Only valid when POINTING=0 ; (absolute) or POINTING=-1 (relative to ; 1st raster). ; INS_Y - ... ; WIDTH - Width (E/W extent) of the raster, in arcs ; HEIGHT - Height (N/S extent) of the raster, in arcs ; ; Note: For the case of CDS, pointings.width, pointings.height, ; pointings.ins_x, and pointings.ins_y should match the first ; raster's rasters.width, rasters.height, rasters.ins_x, and ; rasters.ins_y, respectively. ; ; CALLING SEQUENCE: ; ITOOL_PTOOL, parent, child [,font=font] ; ; INPUTS: ; PARENT - ID of parent widget upon which the pointing widget is built ; ; OPTIONAL INPUTS: ; FONT - Font name to be used for labelling ; ; OUTPUTS: ; CHILD - ID of the base widget being built and leter on remapped ; ; OPTIONAL OUTPUTS: ; None. ; ; KEYWORD PARAMETERS: ; None. ; ; COMMON BLOCKS: ; FOR_POINTING - Internal common block used by this routine and ; ITOOL_PTOOL ; Others - Included in image_tool_com.pro ; ; CATEGORY: ; Planning, pointing ; ; PREVIOUS HISTORY: ; Written January 23, 1995, Liyun Wang, NASA/GSFC ; Separated from the origina code of IMAGE_TOOL ; ; MODIFICATION HISTORY: ; Version 2, Liyun Wang, NASA/GSFC, March 8, 1995 ; Added indicator showing if a selected pointed area is pointed ; Added Warnings if there is any pointing area remain unpointed ; Version 3, Liyun Wang, NASA/GSFC, March 31, 1995 ; Added undo feature ; Version 4, Liyun Wang, NASA/GSFC, April 27, 1995 ; Added the FONT keyword ; Version 5, November 2, 1995, Liyun Wang, NASA/GSFC ; Modified to cope with image icon stacking ; Version 6, November 17, 1995, Liyun Wang, NASA/GSFC ; Added validity check for CDS ; Version 7, December 27, 1995, Liyun Wang, NASA/GSFC ; Added pointing for MK_SOHO ; Version 8, April 17, 1996, Liyun Wang, NASA/GSFC ; Fixed a bug that caused the program crash when switching ; studies with different pointing numbers ; Version 9, August 19, 1996, Liyun Wang, NASA/GSFC ; Corrected an error in plotting rasters with relative pointings ; Version 10, December 18, 1996, Liyun Wang, NASA/GSFC ; Modified such that the time projection is turned off if ; raster is placed off the limb ; Version 11, February 12, 1997, Liyun Wang, NASA/GSFC ; Fixed a bug causing text widgets remaining inactive when ; studies with variable pointings are received ; Version 12, February 18, 1997, Liyun Wang, NASA/GSFC ; Added fields SCI_SPEC and ID in pointing tool ; Version 13, April 1, 1997, Liyun Wang, NASA/GSFC ; Fixed a bug occurred during plotting off_limb rasters ; Version 14, William Thompson, GSFC, 8 April 1998 ; Changed !D.N_COLORS to !D.TABLE_SIZE for 24-bit displays ; Modified, 12-Jan-2006, Zarro (L-3Com/GSFC) - added /NO_QUIT ; keyword ; Modified, 14-August-2006, Zarro (ADNET/GSFC) - made EIS-specific ; ;- PRO itool_pnt_send, event=event ;--------------------------------------------------------------------------- ; Send pointing structure back to planning tool ;--------------------------------------------------------------------------- @image_tool_com ; IF NOT exit_ok THEN BEGIN ; exit_ok = xanswer(['Warning!!!', $ ; 'Not all pointing areas are pointed.',$ ; 'Proceed any way?'], /beep, /supp,/center, group=event.top) ; ENDIF ; IF NOT exit_ok THEN RETURN exit_ok=1 chk=where(pointing_stc.pointings.out_fov eq 1,count) if count ne 0 then begin xack,'At least one raster is pointed outside the EIS FOV' return endif exit_ok=1 TVLCT, r, g, b, /get rgb = [[r], [g], [b]] IF WIDGET_INFO(pointing_stc.messenger, /valid) THEN BEGIN widget_control2, pointing_stc.messenger, set_uvalue=eis_plan_stc(pointing_stc), $ timer=1.0 caller = WIDGET_INFO(pointing_stc.messenger, /parent) IF WIDGET_INFO(caller, /valid) THEN xshow, caller ENDIF END ;--------------------------------------------------------------------------- ; PURPOSE: ; Get two vectors that reflect the raster shape for plotting purpose ; ; CALLING SEQUENCE: ; pt_fov_shape, width, height, ppx, ppy [, initial=initial] ; ; INPUTS: ; WIDTH - Width (E/W extent) of the raster, in arcsecs ; HEIGHT - Height (N/S extent) of the raster, in arcsecs ; ; OPTIONAL INPUTS: ; INITIAL - A 2-element vector that contains the coordinates (in ; arc seconds) of center of the raster. ; ; OUTPUTS: ; PPX - Array of X positions of each point of a pointing area, in arcsec ; PPY - Array of Y positions of each point of a pointing area, in arcsec ; ;--------------------------------------------------------------------------- ; PRO pt_fov_shape, width, height, ppx, ppy, initial=initial IF N_ELEMENTS(initial) EQ 0 THEN initial = [0, 0] ;---------------------------------------------------------------------- ; x0 and y0 are coordinates of the center of a raster ;---------------------------------------------------------------------- x0 = initial(0) & y0 = initial(1) ;---------------------------------------------------------------------- ; Almost all rasters have even number of pixels in size, the plotted raster ; is not exactly centered at the pixel indicated by (x0, y0) ;---------------------------------------------------------------------- hlf_wid = FLOAT(width)/2.0 hlf_hgt = FLOAT(height)/2.0 IF hlf_wid EQ 0.0 AND hlf_hgt EQ 0.0 THEN BEGIN ppx = x0 ppy = y0 ENDIF ELSE BEGIN IF hlf_wid EQ 0.0 THEN BEGIN ppx = [x0, x0] IF (hlf_hgt*2.0 EQ height) THEN yll = y0-hlf_hgt+1.0 ELSE yll = y0-hlf_hgt yur = y0+hlf_hgt ppy = [yll, yur] ENDIF ELSE IF hlf_hgt EQ 0.0 THEN BEGIN ppy = [y0, y0] IF (hlf_wid*2.0 EQ width) THEN xll = x0-hlf_wid+1.0 ELSE xll = x0-hlf_wid xur = x0+hlf_wid ppx = [xll, xur] ENDIF ELSE BEGIN IF (hlf_wid*2.0 EQ width) THEN xll = x0-hlf_wid+1.0 ELSE xll = x0-hlf_wid xur = x0+hlf_wid IF (hlf_hgt*2.0 EQ height) THEN yll = y0-hlf_hgt+1.0 ELSE yll = y0-hlf_hgt yur = y0+hlf_hgt ppx = [xll, xur, xur, xll, xll] ppy = [yll, yll, yur, yur, yll] ENDELSE ENDELSE END ;--------------------------------------------------------------------------- ; To modify the position of a point according to the solar differential ; rotation ; ; Pass in and return out 2-element vector in arcsecond ;--------------------------------------------------------------------------- FUNCTION pt_rotate_point, position @image_tool_com sz = SIZE(position) oned = (sz(0) EQ 1) IF oned EQ 1 THEN temp = TRANSPOSE(position) ELSE temp = position time_diff = (utc2tai(csi.date_obs)-tai_start) temp = rot_xy(temp(*, 0), temp(*, 1), time_diff, $ date=tai2utc(tai_start), off=off) IF off EQ 1 THEN BEGIN ; flash_msg, comment_id, $ ; 'Pointing area was off the limb and is ' + $ ; 'plotted at current study start time.', num=2 ; widget_control2, point_wid.point_tproj, set_button=0 ; time_proj = 0 RETURN, position ENDIF ELSE BEGIN IF oned EQ 1 THEN RETURN, TRANSPOSE(temp) ELSE RETURN, temp ENDELSE END ;--------------------------------------------------------------------------- ; PURPOSE: ; To return a 5x2 array as pointing area shape in device coordinates ; ; EXPLANATION: ; Given the initial posistion (in solar disc coordinates), width and ; height of a pointing area, this routine returns a 2x5 array ; representing the shape of the pointing area in device coordinates that ; can be used for plotting. Time projection and off limb criteria are ; taken into account via the keywords TIME_PROJ and OFF_LIMB. ; ; CALLING SEQUENCE: ; result = pt_fov_dshape(width, height, initial=initial, csi=csi,$ ; time_proj=time_proj, off_limb=off_limb) ; ; INPUTS: ; WIDTH - Width of the pointing area in arcsecs ; HEIGHT - Height of the pointing area in arcsecs ; INITIAL - Initial position of the pointing area in arcsecs ; CSI - Coordinate system info structure ; ; OUTPUTS: ; RESULT - 5x2 integer array representing the pointing area shape in ; device pixels ; ; KEYWORD PARAMETERS: ; TIME_PROJ - Do differential rotating if set ; OFF_LIMB - Do not do differential rotating if set ; ADJUST - If set, the pointing area position will be adjusted so ; that it will be within the display area; Needed when you ; just want to starting pointing ; NOROTATE - Set this keyword not to differentially rotate ; ; COMMON BLOCKS: ; ITOOL_AXES_COM (from itool_plot_axes.pro) ; ;--------------------------------------------------------------------------- ; FUNCTION pt_fov_dshape, width, height, initial=initial, time_proj=time_proj, $ off_limb=off_limb, csi=csi, adjust=adjust, $ norotate=norotate COMMON itool_axes_com, x_range, y_range initial0 = initial IF time_proj and (1-off_limb) and (1-KEYWORD_SET(norotate)) THEN $ initial0 = pt_rotate_point(initial0) ;--------------------------------------------------------------------------- ; If ADJUST is set and initial0 is out of displaying range, make it in ;--------------------------------------------------------------------------- IF KEYWORD_SET(adjust) THEN BEGIN IF initial0(0) LT x_range(0) OR initial0(0) GT x_range(1) THEN $ initial0(0)=0.5*(x_range(0)+x_range(1)) IF initial0(1) LT y_range(0) OR initial0(1) GT y_range(1) THEN $ initial0(1)=0.5*(y_range(0)+y_range(1)) ENDIF pt_fov_shape, width, height, ppx, ppy, initial = initial0 ;---------------------------------------------------------------------- ; convert ppx and ppy into graphic device coordinates ;---------------------------------------------------------------------- RETURN, cnvt_coord([[ppx], [ppy]], from=3, to=1, csi=csi) END ;---------------------------------------------------------------------- ; Update contents on the point widget base ;---------------------------------------------------------------------- ; PRO pt_refresh_base @image_tool_com COMMON for_pointing, pointing_done, saved_pointings, central_pos, $ saved_x, saved_y, saved_pt, index IF N_ELEMENTS(index) EQ 0 THEN index = 0 IF N_ELEMENTS(csi) EQ 0 THEN RETURN ; IF tools(curr_tool).uvalue EQ 'ptool' AND $ ; WIDGET_INFO(tools(curr_tool).base, /valid) THEN $ ; widget_control2, tools(curr_tool).base, sensitive=(csi.flag NE 0) ;--------------------------------------------------------------------------- ; Warn the user if the pointing area is wider than the displayed image ;--------------------------------------------------------------------------- sz = SIZE(image_arr) index = 0 > index < (pointing_stc.n_pointings-1) IF index LT 0 THEN index = 0 active = pointing_stc.do_pointing NE 0 active=1b widget_control2, point_wid.point_go, sensitive=active widget_control2, point_wid.point_x, sensitive=active widget_control2, point_wid.point_y, sensitive=active widget_control2, point_wid.pointed, sensitive=active ; IF NOT active THEN BEGIN ;--------------------------------------------------------------------------- ; No pointing change is allowed or needed ;--------------------------------------------------------------------------- ; widget_control2, point_wid.pointed, set_button=1 ; pointing_done(index) = 1 ; widget_control2, point_wid.raster_plot, $ ; sensitive=(pointing_stc.n_rasters GT 1) ; ENDIF widget_control2, point_wid.point_list, set_list_select=index ; IF WIDGET_INFO(comment_id, /valid) THEN BEGIN ; IF pointing_stc.pointings(index).width GT csi.cdelt1*sz(1) OR $ ; pointing_stc.pointings(index).height GT csi.cdelt2*sz(2) THEN $ ; widget_control2, comment_id, set_value=$ ; 'Warning: Pointing area is wider than the image area!' ; ENDIF widget_control2, point_wid.point_size, set_value=$ STRTRIM(FIX(pointing_stc.pointings(index).width),2)+' X '+$ STRTRIM(FIX(pointing_stc.pointings(index).height),2) widget_control2, point_wid.point_cnt, set_value=num2str(index) widget_control2, point_wid.point_x, set_value=$ num2str(pointing_stc.pointings(index).ins_x, FORMAT='(f10.2)') widget_control2, point_wid.point_y, set_value=$ num2str(pointing_stc.pointings(index).ins_y, FORMAT='(f10.2)') dprint,'EIS x-pointing: ',pointing_stc.pointings(index).ins_x dprint,'EIS y-pointing: ',pointing_stc.pointings(index).ins_y widget_control2, point_wid.sci_spec, set_value=pointing_stc.pointings(index).zone END ;--------------------------------------------------------------------------- ; Make sure that buttons look right ;--------------------------------------------------------------------------- PRO pt_button_update, widgets=point_wid, status=status, n_rasters=n_rasters IF status THEN BEGIN widget_control2, point_wid.pointed, set_button=1 widget_control2, point_wid.point_id, sensitive=1 widget_control2, point_wid.raster_plot, sensitive=(n_rasters GT 1) ENDIF ELSE BEGIN widget_control2, point_wid.point_id, sensitive=0 widget_control2, point_wid.pointed, set_button=0 widget_control2, point_wid.raster_plot, sensitive=0 ENDELSE END ;--------------------------------------------------------------------------- ; Reset the pointing base upon a given pointing structure ;--------------------------------------------------------------------------- PRO pt_fov_reset, pointing_stc, widgets=point_wid COMMON for_pointing IF pointing_stc.n_pointings EQ 1 THEN $ widget_control2, point_wid.point_label, set_value='1 RASTER' $ ELSE $ widget_control2, point_wid.point_label, $ set_value=num2str(pointing_stc.n_pointings)+' RASTERS' IF pointing_stc.n_pointings GT 0 THEN BEGIN pointing_done = INTARR(pointing_stc.n_pointings) FOR i=0, pointing_stc.n_pointings-1 DO BEGIN IF N_ELEMENTS(pointing_list) EQ 0 THEN $ pointing_list = ' '+$ STRTRIM(pointing_stc.pointings(i).point_id, 2) $ ELSE $ pointing_list = [pointing_list,' '+$ STRTRIM(pointing_stc.pointings(i).point_id,2)] ENDFOR ENDIF ELSE BEGIN pointing_list = ' '+'0' pointing_done = INTARR(1) ENDELSE IF pointing_stc.std_id GT 0 THEN $ widget_control2, point_wid.std_id, $ set_value=STRTRIM(pointing_stc.std_id, 2) $ ELSE $ widget_control2, point_wid.std_id, set_value='N/A' widget_control2, point_wid.point_list, set_value=pointing_list widget_control2, point_wid.undo_one, sensitive=0 widget_control2, point_wid.undo_all, sensitive=0 widget_control2, point_wid.raster_plot, sensitive=0 widget_control2, point_wid.point_id, sensitive=0 widget_control2, point_wid.pointed, set_button=0 saved_pointings = pointing_stc.pointings pt_refresh_base END ;------------------------------------------------------------------------------- PRO pt_fov_update @image_tool_com COMMON for_pointing if n_elements(index) eq 0 then index=0 IF (pointing_stc.pointings(index).ins_x NE central_pos(0) OR $ pointing_stc.pointings(index).ins_y NE central_pos(1)) THEN BEGIN pointing_stc.pointings(index).ins_x = central_pos(0) pointing_stc.pointings(index).ins_y = central_pos(1) pointing_done(index) = 1 ENDIF pt_button_update, widgets=point_wid, status=pointing_done(index), $ n_rasters=pointing_stc.n_rasters widget_control2, point_wid.undo_one, sensitive=1 widget_control2, point_wid.undo_all, sensitive=1 pt_refresh_base RETURN END ;--------------------------------------------------------------------------- ; PURPOSE: ; Event handler for events generated by the pointing base widget ; ; CALLING SEQUENCE: ; pt_ptool_event, event, uvalue ; ; INPUTS: ; EVENT - The event structure handled by XMANAGER ; UVALUE - uvalue of the EVENT ; ; COMMON BLOCKS: ; FOR_POINTING - Internal common block used by this routine and ; ITOOL_PTOOL ; Others - Included in image_tool_com.pro ;--------------------------------------------------------------------------- ; PRO pt_ptool_event, event, uvalue @image_tool_com COMMON for_pointing ;--------------------------------------------------------------------------- ; Use the rather crude help system for now ;--------------------------------------------------------------------------- IF help_mode THEN BEGIN xshow_help, help_stc, 'POINTING', tbase = help_wbase RETURN ENDIF tmp = WHERE(pointing_done EQ 1, count) IF count LT pointing_stc.n_pointings THEN exit_ok = 0 ELSE exit_ok = 1 widget_control2, tools(curr_tool).base, get_uvalue = point_wid CASE (uvalue) OF 'send': BEGIN IF itool_getxy_field(event) EQ 0 THEN RETURN itool_pnt_send, event=event END 'POINT_DONE': BEGIN IF NOT exit_ok THEN BEGIN exit_ok = xanswer(['Warning!!!', $ 'Not all pointing values are changed.', $ 'Do you wish to exit Pointing Tool any way?'], $ /beep, group=event.top, /center,/suppress) ENDIF IF exit_ok THEN BEGIN itool_switcher, 0 widget_control2, comment_id, set_value=' ' itool_refresh ENDIF END 'POINT_TPROJ': BEGIN IF time_proj EQ 1 THEN time_proj = 0 ELSE time_proj = 1 itool_refresh END 'ptool_peek': BEGIN hd = icon_stc.header_cur mess = 'Selected image does not have raster information' xcen = fxpar(hd, 'XCEN') IF !err EQ -1 THEN BEGIN xack, mess, group=event.top RETURN ENDIF ycen = fxpar(hd, 'YCEN') IF !err EQ -1 THEN BEGIN xack, mess, group=event.top RETURN ENDIF ixwidth = fxpar(hd, 'IXWIDTH')/csi.cdelt1/csi.ddelt1 IF !err EQ -1 THEN BEGIN xack, mess, group=event.top RETURN ENDIF iywidth = fxpar(hd, 'IYWIDTH')/csi.cdelt2/csi.ddelt2 IF !err EQ -1 THEN BEGIN xack, mess, group=event.top RETURN ENDIF pos = cnvt_coord(xcen, ycen, csi=csi, from=3, to=1, $ date=disp_utc) cx = pos(0, 0) & cy=pos(0, 1) hlf_wid = FIX(ixwidth)/2 hlf_hgt = FIX(iywidth)/2 IF hlf_wid EQ 0 AND hlf_hgt EQ 0 THEN BEGIN xx = cx yy = cy ENDIF ELSE BEGIN IF hlf_wid EQ 0 THEN BEGIN yll = cy-hlf_hgt yur = cy+hlf_hgt xx = [cx, cx] yy = [yll, yur] ENDIF ELSE IF hlf_hgt EQ 0 THEN BEGIN xll = cx-hlf_wid xur = cx+hlf_wid xx = [xll, xur] yy = [cy, cy] ENDIF ELSE BEGIN xll = cx-hlf_wid xur = cx+hlf_wid yll = cy-hlf_hgt yur = cy+hlf_hgt xx = [xll, xur, xur, xll, xll] yy = [yll, yll, yur, yur, yll] ENDELSE ENDELSE PLOTS, xx, yy, col=!d.table_size-1, /dev, thick=1, lines=0 END 'cancel': BEGIN pointing_go = 0 widget_control2, draw_icon, draw_button=1 widget_control2, tools_bt, sensitive=1 if widget_valid(point_wid.ptool_exit) then $ widget_control2, point_wid.ptool_exit, sensitive=1 widget_control2, point_wid.cancel, sensitive=0 widget_control2, point_wid.point_go, set_value='Go' widget_control2, comment_id, set_value='' itool_restore_pix, pix_win END 'POINT_GO': BEGIN IF pointing_go EQ 0 THEN BEGIN if not exist(image_arr) then begin flash_msg, comment_id, $ $ ['Please load in your images via:',$ 'File-> Import Images'],$ num=2 return endif ;--------------------------------------------------------------------------- ; Begin pointing mode ;--------------------------------------------------------------------------- widget_control2, draw_icon, draw_button=0, /hour widget_control2, tools_bt, sensitive=0 if widget_valid(point_wid.ptool_exit) then $ widget_control2, point_wid.ptool_exit, sensitive=0 widget_control2, point_wid.cancel, sensitive=1 widget_control2, point_wid.point_go, set_value='Ok' pointing_go = 1 widget_control2, comment_id, $ set_value=['You are now looking at the raster ',$ 'that you have chosen.', $ 'Press and drag the right mouse button',$ 'to move the raster.',$ 'Click the "Ok" button when done.'] initial0 = [pointing_stc.pointings(index).ins_x, $ pointing_stc.pointings(index).ins_y] pointing_fov = pt_fov_dshape(pointing_stc.pointings(index).width, $ pointing_stc.pointings(index).height, $ initial=initial0, time_proj=time_proj, $ off_limb=pointing_stc.pointings(index).off_limb, $ csi=csi) IF !d.window NE root_win THEN setwindow, root_win PLOTS, pointing_fov(*, 0), pointing_fov(*, 1), /dev ENDIF ELSE BEGIN ;--------------------------------------------------------------------------- ; Get out of pointing mode ;--------------------------------------------------------------------------- pointing_go = 0 widget_control2, draw_icon, draw_button=1 widget_control2, tools_bt, sensitive=1 if widget_valid(point_wid.ptool_exit) then $ widget_control2, point_wid.ptool_exit, sensitive=1 widget_control2, point_wid.cancel, sensitive=0 widget_control2, point_wid.point_go, set_value='Go' widget_control2, comment_id, set_value='' xx = pointing_fov(*, 0) yy = pointing_fov(*, 1) x0 = 0.5*(MIN(xx)+MAX(xx)) y0 = 0.5*(MIN(yy)+MAX(yy)) central_pos = TRANSPOSE(cnvt_coord(x0, y0, csi=csi, from=1, to=3)) ;-- check if raster falls in EIS fov rfov_stc=fov_stc limb_fov=eis_limb_fov(fov_stc,offlimb=offlimb) if time_proj and ~offlimb then $ rfov_stc=eis_drot_fov(fov_stc,tend=csi.date_obs) in_fov=eis_valid_point(central_pos[0],central_pos[1],pointing_stc,rfov_stc,index=index) if not in_fov then begin xack,'Warning: Raster falls outside EIS field-of-view',/center pointing_stc.pointings[index].out_fov=1b endif else pointing_stc.pointings[index].out_fov=0b radial = SQRT(central_pos(0)*central_pos(0)+$ central_pos(1)*central_pos(1)) pointing_stc.pointings[index].off_limb=(radial gt csi.radius) ;IF (radial GT csi.radius) AND time_proj and offlimb THEN BEGIN ; flash_msg, comment_id, $ ; ['Warning: Time projection option is turned off because you '+$ ; 'have placed the raster off the limb.'], num=2 ; time_proj = 0 ; widget_control2, point_wid.point_tproj, set_button=0 ;ENDIF ;---------------------------------------------------------------------- ; Currently I decide that only when the time_proj flag ; is on (with value 1) does the time projection rotate take place ;---------------------------------------------------------------------- IF time_proj AND (1-pointing_stc.pointings(index).off_limb) THEN BEGIN ;---------------------------------------------------------------------- ; Depending on the image observation time and the starting time of ; the raster, we need to rotate the raster pointing forward from the ; image observation time to the actual time at which the study ; begins. If the time lag is greater such that the poiting of the ; raster would be off, the user is warned, and no correction due to ; rotation is made ;---------------------------------------------------------------------- time_diff = (tai_start-utc2tai(csi.date_obs))/86400.0 helio = cnvt_coord(central_pos(0), central_pos(1), csi=csi, $ from=3, to=4) helio(0, 1) = helio(0, 1)+diff_rot(time_diff, helio(0, 0), /synodic) long_diff = 90.0-helio(0, 1) IF (long_diff LE 0.0 OR long_diff GE 180.0) THEN BEGIN ;--------------------------------------------------------------------------- ; Pointing center off the limb ;--------------------------------------------------------------------------- ; flash_msg, comment_id, $ ; ['Warning: The projected raster at the time of '+$ ; 'observation will be off the limb!', $ ; 'Pointing shown is for the time at '+$ ; 'which this displayed image was made.'], num=2 ; widget_control2, point_wid.point_tproj, set_button=0 ; time_proj = 0 do_nothing=1 ENDIF ELSE BEGIN ;---------------------------------------------------------------------- ; When converting rotated pointing back, actual study time ; should be used ;---------------------------------------------------------------------- new_date = tai2utc(tai_start) central_pos = TRANSPOSE(cnvt_coord(helio, csi=csi, from=4, to=3, $ date=new_date)) ENDELSE ENDIF widget_control2, draw_icon, draw_button=1 pt_fov_update ENDELSE END 'UNDO_ONE': BEGIN IF N_ELEMENTS(saved_x) NE 0 THEN BEGIN pointing_stc.pointings(index).ins_x = saved_x pointing_stc.pointings(index).ins_y = saved_y pointing_done(index) = saved_pt widget_control2, point_wid.pointed, set_button = saved_pt ENDIF ELSE BEGIN pointing_stc.pointings(index).ins_x = saved_pointings(index).ins_x pointing_stc.pointings(index).ins_y = saved_pointings(index).ins_y pointing_done(index) = 0 saved_pt = 0 widget_control2, point_wid.pointed, set_button = 0 ENDELSE widget_control2, point_wid.point_id, sensitive = saved_pt widget_control2, point_wid.undo_one, sensitive = 0 widget_control2, point_wid.raster_plot, sensitive = (saved_pt $ AND pointing_stc.n_rasters GT 1) pt_refresh_base END 'UNDO_ALL': BEGIN IF xanswer([' ', 'Are you sure?', ' '], /beep, $ group=event.top, /center) THEN BEGIN pointing_stc.pointings = saved_pointings pointing_done(*) = 0 saved_pt = 0 saved_x = pointing_stc.pointings(index).ins_x saved_y = pointing_stc.pointings(index).ins_y widget_control2, point_wid.undo_all, sensitive = 0 widget_control2, point_wid.undo_one, sensitive = 0 widget_control2, point_wid.pointed, set_button = 0 widget_control2, point_wid.point_id, sensitive = 0 widget_control2, point_wid.raster_plot, sensitive = 0 pt_refresh_base itool_refresh ENDIF END 'POINTED': BEGIN ;--------------------------------------------------------------------------- ; Toggles the current selected pointing area pointed or unpointed ;--------------------------------------------------------------------------- IF pointing_done(index) THEN BEGIN pointing_done(index) = 0 widget_control2, point_wid.point_id, sensitive = 0 widget_control2, point_wid.raster_plot, sensitive = 0 ; widget_control2, point_wid.undo_one, sensitive = 0 ;--------------------------------------------------------------------------- ; Save current pointing ;--------------------------------------------------------------------------- saved_x = pointing_stc.pointings(index).ins_x saved_y = pointing_stc.pointings(index).ins_y ENDIF ELSE BEGIN pointing_done(index) = 1 widget_control2, point_wid.point_id, sensitive = 1 widget_control2, point_wid.raster_plot, $ sensitive = (pointing_stc.n_rasters GT 1) ENDELSE tmp = WHERE(pointing_done EQ 1, count) IF count LT pointing_stc.n_pointings THEN exit_ok = 0 ELSE exit_ok = 1 END 'POINT_ID': BEGIN IF NOT itool_getxy_field(event) THEN RETURN IF pointing_done(index) THEN BEGIN initial = [pointing_stc.pointings(index).ins_x, $ pointing_stc.pointings(index).ins_y] temp = pt_fov_dshape(pointing_stc.pointings(index).width, $ pointing_stc.pointings(index).height, $ initial=initial, time_proj=time_proj, $ off_limb=pointing_stc.pointings(index).off_limb, $ csi=csi) flash_plots, temp(*, 0), temp(*, 1), color=(l_color/2 > 0) ENDIF ELSE BEGIN flash_msg, comment_id, 'That raster has not been ' + $ 'pointed yet.', num = 2 ENDELSE END 'POINT_PLOT': BEGIN itool_point_plot END 'RASTER_PLOT': BEGIN ;---------------------------------------------------------------------- ; Raster plotting only applies to the case where there are more than 1 ; raster in each pointing area. This is irrelevant to the SUMER. ; ; First find the target raster (the first in the list that has pointing ; value being 1, or if there is no such raster, being 0) ;---------------------------------------------------------------------- IF NOT itool_getxy_field(event) THEN RETURN ii = WHERE(pointing_stc.rasters.pointing EQ 1, cnt) IF cnt NE 0 THEN BEGIN initial0 = [pointing_stc.pointings(index).ins_x, $ pointing_stc.pointings(index).ins_y] ENDIF ELSE BEGIN dhelp, 'No raster with deferred point_wid..' ii = (WHERE(pointing_stc.rasters.pointing EQ 0, cnt))(0) IF cnt NE 0 THEN BEGIN initial0 = [pointing_stc.rasters(ii).ins_x, $ pointing_stc.rasters(ii).ins_y] ENDIF ELSE BEGIN PRINT, 'Target raster not found.' RETURN ENDELSE ENDELSE widget_control2, /hour IF !d.window NE root_win THEN setwindow, root_win IF time_proj AND (1-pointing_stc.pointings(index).off_limb) THEN $ initial0 = pt_rotate_point(initial0) itool_restore_pix, pix_win FOR i = 0, pointing_stc.n_rasters-1 DO BEGIN CASE (pointing_stc.rasters(i).pointing) OF -1: BEGIN ;--------------------------------------------------------------------------- ; Raster position is offset from the first one ;--------------------------------------------------------------------------- initial = initial0+[pointing_stc.rasters(i).ins_x,$ pointing_stc.rasters(i).ins_y] pt_fov_shape, pointing_stc.rasters(i).width, $ pointing_stc.rasters(i).height, ppx, ppy, $ initial = initial linestyle = 1 END 0: BEGIN ;--------------------------------------------------------------------------- ; Raster position is determined at the study level ;--------------------------------------------------------------------------- initial = [pointing_stc.rasters(i).ins_x, $ pointing_stc.rasters(i).ins_y] IF time_proj THEN initial = pt_rotate_point(initial) pt_fov_shape, pointing_stc.rasters(i).width, $ pointing_stc.rasters(i).height, ppx, ppy, $ initial = initial linestyle = 2 END 1: BEGIN ;--------------------------------------------------------------------------- ; Raster position is the same as the pointings(index).ins_x ; and pointings(index).ins_y ;--------------------------------------------------------------------------- dprint, 'This is a deferred-pointing raster...' pt_fov_shape, pointing_stc.rasters(i).width, $ pointing_stc.rasters(i).height, ppx, ppy, $ initial = initial0 linestyle = 0 END ENDCASE temp = cnvt_coord([[ppx], [ppy]], csi=csi, from=3, to=1) PLOTS, temp(*, 0), temp(*, 1), /dev, color=l_color, $ linestyle=linestyle ENDFOR itool_copy_to_pix END 'POINT_X': BEGIN IF NOT itool_getxy_field(event) THEN RETURN widget_control2, point_wid.point_y, /input_focus END 'POINT_Y': BEGIN IF NOT itool_getxy_field(event) THEN RETURN widget_control2, point_wid.point_x, /input_focus END 'POINT_LIST': BEGIN IF NOT itool_getxy_field(event) THEN RETURN index = event.index widget_control2, point_wid.undo_one, sensitive = 0 pt_button_update, widgets=point_wid, status=pointing_done(index), $ n_rasters=pointing_stc.n_rasters ;--------------------------------------------------------------------------- ; Save current pointing ;--------------------------------------------------------------------------- saved_x = pointing_stc.pointings(index).ins_x saved_y = pointing_stc.pointings(index).ins_y saved_pt = pointing_done(index) pt_refresh_base END ELSE: ENDCASE RETURN END FUNCTION eis_itool_ptool, parent, font=font,no_quit=no_quit ;--------------------------------------------------------------------------- ; PURPOSE: ; Main program to create the pointing widget ; ; CALLING SEQUENCE: ; child = itool_ptool(parent, [, font=font]) ; ; INPUTS: ; PARENT - ID of the parent base widget on which pointing widget ; is built up ; OUTPUTS: ; CHILD - ID of the pointing base widget ; ; KEYWORD PARAMETERS: ; FONT - Basic font to be used in pointing widget ;--------------------------------------------------------------------------- ; @image_tool_com COMMON for_pointing ;--------------------------------------------------------------------------- ; Check validity of the given STUDY structure ;--------------------------------------------------------------------------- IF datatype(pointing_stc) NE 'STC' THEN BEGIN MESSAGE, 'Invalid parameter type.',/cont RETURN, -1 ENDIF if not have_tag(pointing_stc,'pointings') then begin MESSAGE, 'Invalid pointing structure.',/cont RETURN, -1 ENDIF IF N_ELEMENTS(font) EQ 0 THEN font = 'fixed' ;--------------------------------------------------------------------------- ; Clean some variables in common block ;--------------------------------------------------------------------------- IF N_ELEMENTS(saved_x) NE 0 THEN delvarx, saved_x, saved_y, saved_pt mk_dfont,lfont=lfont,bfont=nfont title = 'Pointing Tool' ;--------------------------------------------------------------------------- ; Since this is a general routine to create an instrument base, it should ; always be unmapped initially. The return widget ID "child" is used to map ; the base ;--------------------------------------------------------------------------- child = WIDGET_BASE(parent, /column) ; temp = WIDGET_LABEL(child, value='') ; temp = WIDGET_LABEL(child, value=title, font=lfont) ; temp = WIDGET_LABEL(child, value='') junk = WIDGET_BASE(child, /row) if keyword_set(no_quit) then ptool_exit=-1 else $ ptool_exit = WIDGET_BUTTON(junk, value='Exit',uvalue='POINT_DONE',font=lfont) temp = WIDGET_BUTTON(junk, value='Plot All Rasters',font=lfont,uvalue='POINT_PLOT') ; point_plot = WIDGET_BUTTON(temp, value='All Pointing Areas', $ ; uvalue='POINT_PLOT',font=lfont) ; raster_plot = WIDGET_BUTTON(temp, uvalue='RASTER_PLOT',font=lfont, $ ; value='All Rasters for Current Pointing') ; ptool_peek = WIDGET_BUTTON(temp, uvalue='ptool_peek', font=lfont,$ ; value='Locate Raster Position') ; widget_control2, ptool_peek, sensitive=0 raster_plot=-1l ptool_send = WIDGET_BUTTON(junk, uvalue='send', font=lfont, $ value='Send to Planning Tool') temp = WIDGET_LABEL(child, value='') ;---------------------------------------------------------------------- ; The pointing button row ;---------------------------------------------------------------------- tmp = WIDGET_BASE(child, /row) point_go = WIDGET_BUTTON(tmp, value='Go', uvalue='POINT_GO',font=lfont) cancel = WIDGET_BUTTON(tmp, value='Cancel', uvalue='cancel',font=lfont) widget_control2, cancel, sensitive=0 undo_one = WIDGET_BUTTON(tmp, value='Undo', uvalue='UNDO_ONE',font=lfont) undo_all = WIDGET_BUTTON(tmp, value='Undo All', uvalue='UNDO_ALL',font=lfont) point_id = WIDGET_BUTTON(tmp, value='Identify Raster', uvalue='POINT_ID',font=lfont) ;---------------------------------------------------------------------- ; Make a toggle button to enable/disable time projection. Default is ; to allow time projection ;---------------------------------------------------------------------- temp = WIDGET_BASE(child, /row) tproj_base = WIDGET_BASE(temp, /nonexclusive, /frame) point_tproj = WIDGET_BUTTON(tproj_base, value='Allow Time Projection', $ uvalue='POINT_TPROJ',font=lfont) IF N_ELEMENTS(time_proj) NE 0 THEN BEGIN IF time_proj THEN widget_control2, point_tproj, set_button=1 $ ELSE widget_control2, point_tproj, set_button=0 ENDIF ELSE BEGIN time_proj = 1 widget_control2, point_tproj, set_button=1 ENDELSE ; tmp = WIDGET_LABEL(temp, value=' ', font=lfont) ; tmp = WIDGET_BASE(temp, /row, /frame) ; tmp0 = WIDGET_LABEL(tmp, value='ID', font=lfont) ; std_id = WIDGET_TEXT(tmp, value='', font=lfont, xsize=5) std_id=-1 tmp0 = WIDGET_BASE(child, /column) tmp2 = WIDGET_BASE(tmp0, /row) lb = WIDGET_LABEL(tmp2, value='Study Start Time:', font=lfont) start_text = WIDGET_TEXT(tmp2, value='', font=lfont, xsize=25) tmp = WIDGET_BASE(tmp0, /row) lb = WIDGET_LABEL(tmp, value= 'Raster Acronym: ', font=lfont) sci_spec = WIDGET_TEXT(tmp, value='', font=lfont, xsize=25) sr_row3 = WIDGET_BASE(child, /row) sr_left = WIDGET_BASE(sr_row3, /column, /frame) point_label = WIDGET_LABEL(sr_left, value='', font=lfont) tmp = WIDGET_LABEL(sr_left, value=pointing_stc.g_label, font=lfont) point_list = WIDGET_LIST(sr_left, ysize=9, xsize=4, uvalue='POINT_LIST', $ font=lfont) tmp = WIDGET_BASE(sr_row3, /column, /frame) ;--------------------------------------------------------------------------- ; Get maximum length of labels ;--------------------------------------------------------------------------- max_len = 8 IF STRLEN(pointing_stc.x_label) GT max_len THEN $ max_len = STRLEN(pointing_stc.x_label) IF STRLEN(pointing_stc.y_label) GT max_len THEN $ max_len = STRLEN(pointing_stc.y_label) tmp1 = WIDGET_BASE(tmp, /row) tmp2 = WIDGET_LABEL(tmp1, font=lfont, value=$ justify('POINTING', MAX=max_len, just='>')) point_cnt = WIDGET_TEXT(tmp1, value='1', font=lfont, xsize=2) ; tmp2 = WIDGET_BASE(tmp1, /row, /nonexclusive) ; pointed = WIDGET_BUTTON(tmp2, value='POINTED', uvalue='POINTED', $ ; font=lfont) pointed=-1 tmp1 = WIDGET_BASE(tmp, /row) tmp2 = WIDGET_LABEL(tmp1, font=lfont, value=$ justify('SIZE', MAX=max_len, just='>')) point_size = WIDGET_TEXT(tmp1, value=' ', font=lfont, $ xsize=16) ;--------------------------------------------------------------------------- ; cw_field can be used here, but it won't accept /input_focus ;--------------------------------------------------------------------------- temp = WIDGET_BASE(tmp,/row) tmp1 = WIDGET_LABEL(temp, value=justify(pointing_stc.x_label, MAX=max_len, $ just='>'), font=lfont) point_x = WIDGET_TEXT(temp, value='0.0', uvalue='POINT_X', $ font=lfont, xsize=8,/editable) tmp1 = WIDGET_LABEL(temp, value='"', font=lfont) temp = WIDGET_BASE(tmp,/row) tmp1 = WIDGET_LABEL(temp, value=justify(pointing_stc.y_label, MAX=max_len, $ just='>'), font=lfont) point_y = WIDGET_TEXT(temp, value='0.0', uvalue='POINT_Y', $ font=lfont, xsize=8,/editable) tmp1 = WIDGET_LABEL(temp, value='"', font=lfont) ; temp = WIDGET_BASE(tmp,/row) ; tmp1 = WIDGET_LABEL(temp, font=lfont, value=$ ; justify('ZONE', MAX=max_len, just='>')) ; zone = WIDGET_TEXT(temp, value=' ', font=lfont, xsize=16) zone=-1 point_wid = {point_x:point_x, point_y:point_y, point_size:point_size, $ point_cnt:point_cnt, point_go:point_go, pointed:pointed, $ undo_one:undo_one, point_label:point_label, $ point_list:point_list, point_tproj:point_tproj, tproj_base:tproj_base,$$ point_id:point_id, undo_all:undo_all, $ raster_plot:raster_plot, zone:zone, cancel:cancel, $ sci_spec:sci_spec, std_id:std_id, ptool_exit:ptool_exit} IF pointing_stc.n_pointings GT 0 THEN $ pointing_done = INTARR(pointing_stc.n_pointings) $ ELSE BEGIN pointing_done = INTARR(1) ENDELSE pt_fov_reset, pointing_stc, widgets=point_wid ;--------------------------------------------------------------------------- ; Save widget IDs and other info into the UVALUE of base CHILD as a structure ;--------------------------------------------------------------------------- widget_control2, child, set_uvalue = point_wid RETURN, child END