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