Tags:
create new tag
view all tags

whereIs Plugin

This is the documentation for the whereIs plugin.

The current version v1.8 available here.

PURPOSE:

Returns a list of OBSIDs where a sky location (ra,dec) can be found in Herschel SPIRE and/or PACS maps.

USAGE:

Command-line:

task = whereIsTask()
result = task(ra=123.456, dec=20.00, doPacs=True, doSpire=True, showIt=False, verbose=False, covMapsDir='/home/user/PhotCoveragePLW/', cropDisplayMap=True)
#
print result

After the plugin is installed, the task will appear in Task -> All. Double click on whereIs to open the GUI of the task. You can click on HELP button to get to this page.

Help A zip with all SPIRE 500 Ám coverage maps is available and can be provided upon request (total 408 MB, 6878 files).

INPUTS:

ra - the input Right Ascention (J2000) in decimal degrees
dec - the input Declintation (J2000) in decimal degrees
doPacs - True/False, if PACS maps should be searched (default is False)
doSpire - True/False, if SPIRE maps should be searched (default is True)
showIt - True/False, if to display the results. Default (False). Warning! it can be very slow as is access the Herschel Science Archive.
verbose - True/False, if to be verbose (default False).
cropDisplayMap - True/False, only if showIt=True, if level3 or leve2-5 maps exists then they are cropped to +/- 100 pixels around ra,dec

OUTPUT:

A dictionary with keys PACS and SPIRE each one with Long1d() arrays with the list of the obsids where the coverage at (ra, dec) is greater than zero

TASK DEFS:

whereIsInPacsPhot
whereIsInSpirePhot

CAVEATS:

Can be very slow if the point is in very large (level-2.5 or level-3 map).

HISTORY:

#    Created: April 2014, version 1.0
#        Ivan Valtchanov, Herschel Science Centre, ESAC, ESA
#    v1.1: Bug fix: 
#       1) the check if the pixel coordinates of th source are within a map had swapped row, column.
#       2) if showIt=True then only a subset of OBSIDs is returned if there were level2_5 or level 3.
#       3) only show one map for all in level2_5 or level3
#       4) reverse sort the input table by the OBSID, so the output is revers sorted for each instrument. 
#
#   v1.2 (11 Apr 2014): PACS parallel mode maps were missing from PACS-maps-polygons-red.csv
#   v1.3 (07 May 2014): introduced try: except: to exit gracefully if maps are not available when checking for coverage
#   v1.4 (30 Oct 2014): minor bug fix, added option to allow reading the coverage maps from a folder instead of HSA
#   v1.5 (10 Dec 2014): minor fix: changing the link to the help URL as the plugin is now public.
#   v1.6 (10 Jun 2015): minor fix: PACS coverage in verbose mode was shown as integer, which in some cases resulted in zero
#   v1.7 (24 Sep 2015): fixed PACS display, only one login to the HAIO (avoid =getObservation= in loops)
#                       improvements for speed and memory, allow for large maps to be cropped for display (=cropDisplayMaps=True=)
# v1.8 (17 Nov 2015): minor fix for displaying PACS level-2 maps #

Example run:

result = task(ra=40.0,dec=-1.5667,doPacs=True, doSpire=True, showIt=True, verbose=True)

The output in HIPE will look like this, note the verbose is on and the maps will also be displayed:

Total PACS maps: 18324
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342223333
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342223332
Total number of PACS OBSIDs where the point may be: 2
Now checking the coverage from HSA
getObservation is retrieving the observation from the HSA 
Coverage at 160 um in pixel column=89.127027, row=158.389185: value=5994 in OBSID 1342223333
Found level3 maps, will only show this, may take time to download...
Display will show the map for the following OBSIDs: [1342223332, 1342223333]
getObservation is retrieving the observation from the HSA 
Coverage at 160 um in pixel column=76.203124, row=174.569990: value=5902 in OBSID 1342223332
Total SPIRE maps: 6878
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342257362
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342248004
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342248003
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342248002
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342246580
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342237553
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342212373
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342201318
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342201317
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342201316
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342201315
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342201314
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342201313
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342201312
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342201311
Total number of SPIRE OBSIDs where the point may be: 15
Now checking the coverage from HSA
getObservation is retrieving the observation from the HSA 
getObservation is retrieving the observation from the HSA 
Coverage at 500 um in pixel column=120.509429, row=150.039807: value=28 in OBSID 1342248004
Found level3 maps, will only show this, may take time to download...
Display will show the map for the following OBSIDs: [1342201311, 1342201312, 1342201313, 1342201314, 1342201315, 1342201316, 1342201317, 1342201318, 1342248004, 1342248003, 1342248002]
getObservation is retrieving the observation from the HSA 
Coverage at 500 um in pixel column=124.898330, row=160.321371: value=41 in OBSID 1342248003
getObservation is retrieving the observation from the HSA 
Coverage at 500 um in pixel column=139.586385, row=147.972361: value=22 in OBSID 1342248002
getObservation is retrieving the observation from the HSA 
getObservation is retrieving the observation from the HSA 
getObservation is retrieving the observation from the HSA 
getObservation is retrieving the observation from the HSA 
Coverage at 500 um in pixel column=34.965369, row=52.693440: value=156 in OBSID 1342201318
getObservation is retrieving the observation from the HSA 
Coverage at 500 um in pixel column=35.038806, row=45.999839: value=155 in OBSID 1342201317
getObservation is retrieving the observation from the HSA 
Coverage at 500 um in pixel column=26.315366, row=55.746053: value=57 in OBSID 1342201316
getObservation is retrieving the observation from the HSA 
Coverage at 500 um in pixel column=32.249454, row=38.095532: value=153 in OBSID 1342201315
getObservation is retrieving the observation from the HSA 
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342201314
But the 500 um coverage is 0
getObservation is retrieving the observation from the HSA 
Coverage at 500 um in pixel column=25.059385, row=38.420789: value=94 in OBSID 1342201313
getObservation is retrieving the observation from the HSA 
(Ra, Dec) point (40.000000,-1.566700) found in OBSID 1342201312
But the 500 um coverage is 0
getObservation is retrieving the observation from the HSA 
Coverage at 500 um in pixel column=16.411225, row=41.472156: value=44 in OBSID 1342201311
HIPE>  print result
{'SPIRE': [1342248004,1342248003,1342248002,1342201318,1342201317,1342201316,1342201315,1342201313,1342201311], 'PACS': [1342223333,1342223332]}
HIPE> 

The input coordinates will be shown in the displayed maps as a blue circle. The maps will be centred on the input position and zoomed in.

Note that if there are level-3 (or level-2.5) maps only this big map will be shown, the display window title have an information if the level of the map.

  • Screen_Shot_2014-05-07_at_17.23.10.png:
    Screen_Shot_2014-05-07_at_17.23.10.png

-- IvanV - 07 May 2014

Topic attachments
I Attachment History Action Size Date Who Comment
PNGpng Screen_Shot_2014-05-07_at_17.23.10.png r1 manage 513.4 K 2014-05-07 - 16:26 IvanV  
Unknown file formatjar whereIs_1.2.jar r1 manage 874.5 K 2014-05-07 - 15:50 IvanV  
Unknown file formatjar whereIs_1.3.jar r2 r1 manage 875.3 K 2014-05-07 - 16:03 IvanV  
Unknown file formatjar whereIs_1.4.jar r1 manage 876.3 K 2014-10-30 - 16:09 IvanV Version 1.4 of the plugin
Unknown file formatjar whereIs_1.5.jar r1 manage 876.4 K 2014-12-10 - 09:40 IvanV Version 1.5 of the whereIs plugin
Unknown file formatjar whereIs_1.6.jar r3 r2 r1 manage 875.8 K 2015-06-10 - 09:05 IvanV Version 1.6 of the plugin
Unknown file formatjar whereIs_1.7.jar r1 manage 876.8 K 2015-09-24 - 14:13 IvanV Version 1.7 of the plugin
Unknown file formatjar whereIs_1.8.jar r1 manage 877.4 K 2015-11-17 - 15:22 IvanV WhereIs plugin v1.8
XMLxml whereIs_registry.xml r3 r2 r1 manage 0.8 K 2015-11-17 - 15:22 IvanV  
Edit | Attach | Watch | Print version | History: r9 < r8 < r7 < r6 < r5 | Backlinks | Raw View | Raw edit | More topic actions
Topic revision: r9 - 2015-11-17 - IvanV
 
This site is powered by the TWiki collaboration platform Powered by Perl