Tags:
create new tag
view all tags

Exporting data products from HIPE

This page lists examples and recipes for exporting HIPE data (images, spectra, cubes and source lists) to external applications.

IDL

  • Images

IDL> im = mrdfits('/path/image.fits',1)
% Compiled module: FXMOVE.
% Compiled module: MRD_HREAD.
% Compiled module: FXPAR.
% Compiled module: GETTOK.
% Compiled module: VALID_NUM.
% Compiled module: MRD_SKIP.
MRDFITS: Image array (2012,2009)  Type=Real*8
% Compiled module: SWAP_ENDIAN_INPLACE.
IDL> tv,im

  • Spectra

IDL> spec = mrdfits('/path/1342187893_spectrum_point_HR_unapod.fits',2)
% Compiled module: MATCH.
% Compiled module: MRD_STRUCT.
MRDFITS: Binary table.  4 columns by  2061 rows.
IDL> help,spec,/struc
** Structure <15e03af4>, 4 tags, length=28, data length=28, refs=1:
   WAVE            DOUBLE           31.200000
   FLUX            DOUBLE           8.2931329
   ERROR           DOUBLE           3.4131544
   MASK            LONG                 0
IDL> plot,spec.wave,spec.flux

  • Cubes

IDL> cube = mrdfits('/path/1342192180_MR_SLW_unapod_cube.fits',2)
MRDFITS: Image array (16,18,374)  Type=Real*8
IDL> help,cube
CUBE            DOUBLE    = Array[16, 18, 374]
IDL> plot,cube[8,8,*]

PACS Projected Cubes: The .fits structure is summarized in Sect. 1.15.4.3 of the Herschel Data Analysis Guide (http://herschel.esac.esa.int/hcss-doc-13.0#dag:Dag.DataIO.Fits)

IDL> FITS_HELP,'path/cubeName.fits'
  XTENSION  EXTNAME         EXTVER EXTLEVEL BITPIX GCOUNT  PCOUNT NAXIS  NAXIS*
0                                            32      0      0     0  
1 IMAGE    image                            -64      1      0     3  39 x 39 x 29
2 IMAGE    coverage                         -64      1      0     3  39 x 39 x 29
3 BINTABLE ImageIndex                         8      1      0     2  12 x 29
4 IMAGE    History                           32      1      0     0  
5 BINTABLE HistoryScript                      8      1      0     2  80 x 7
6 BINTABLE HistoryTasks                       8      1      0     2  35 x 1
7 BINTABLE HistoryParameters                  8      1      0     2  103 x 12
IDL> image = mrdfits('path/cubeName', 'image', hd) ; the header contains the image's WCS
IDL> imageIndex = mrdfits('path/cubeName','ImageIndex')
IDL> wave = imageIndex.depthindex  ;cube's wavescale

CLASS

  • The FITS files produced with the hiClass task in HIPE on HIFI data can be read in CLASS with the following commands

    file out MyHIFISpectra.hifi mul
    fits read MyHIFISpectra.fits
#
# Now you have a CLASS file named MyHIFISpectra.hifi (you can use whatever you want as an 
# extension) you can access like you always do in CLASS:
#
    file in MyHIFISpectra.hifi
    find
    get first
    set unit f i
    device image white
    plot

  • For PACS data or any spectum1d product, you can try the this script (remove .txt extension before you load it in HIPE).

DS9

Simply go to File -> Open for images and cubes.


blog comments powered by Disqus

Edit | Attach | Watch | Print version | History: r12 < r11 < r10 < r9 < r8 | Backlinks | Raw View | Raw edit | More topic actions
Topic revision: r12 - 2015-02-18 - AlvarGarcia
 
This site is powered by the TWiki collaboration platform Powered by Perl