User's guide for using INAOE package
Divakara Mayya (ydm@inaoep.mx)
cl> inaoe
+-------------------- INAOE IRAF Package -------------------+
| Version 1.0, APR 24, 2016 |
| Tested with IRAF v2.16 on Fedora |
| Developed by Divakara Mayya at INAOE, Puebla, Mexico |
| |
| Type help before using it for the first time |
| |
| |
|Please contact ydm@inaoep.mx if you have questions/comments|
+-----------------------------------------------------------+
camila. canica. imanal. misc. osiris. spanal.
inaoe> help
camila - Reduction program for CAMILA images [up]
canica - Reduction package for CANICA images [up]
imanal - Utility tasks for analysing images [up]
spanal - Utility tasks for analysing spectra [up]
osiris - Tasks for analysing GTC/OSIRIS imaging data [up]
misc - Miscellaneous utility tasks [up]
inaoe> help spanal
inaoe.spanal:
spanal - Spectroscopic Analysis package [up]
apinputs - Script to create apall database file [up]
bncflat - Prepare a normalized spectroscopic flat [up]
bnchedit - Edit the headers of PIMIS-written files [up]
findspe - Find whether a catalog object is in a GTC slit? [up]
fluxcalc - Measure emission line fluxes in 1-d spectra [up]
omsperms - obtain spectral rms using splot [up]
omspestat- obtain spectral statistics using splot [up]
plotfunc - Prepare a ds9 reg file of apall/trace trajectory [up]
wrnebsub - subtract nebular spectrum from WR spectrum [up]
wr_model - Template fitting to WR spectrum [up]
inaoe.imanal:
#General image analysis tasks
contsub - Subtract the continuum from emission line images [up]
wcsmosaic- Tile the images using WCS data [up]
xytoreg - Convert an XY list to ds9 regfile [up]
hstwcs - Improve astrometry on HST images using catalog star [up]
starcen - Obtain centroids of stars given approximate position [up]
lfosc - Transform LFOSC images to WCS using geomap/geotran [up]
lfosc_wcs- Transform LFOSC images to WCS using geomap/geotran [up]
sedphot - Construct SED using GALEX-SDSS-2MASS images [up]
varphot - Construct light curves for LFOSC data [up]
sdss_to_ubvri - Convert SDSS catalog to ds9 regfile [up]
imscale - Find the image scale from the headers [up]
#Asymmetry parameters
asym - Find Rotational asymmetry Parameter [up]
asym_par - Find asymmetry in model subtracted residual images [up]
nirasym - Calculate asymmetry parameter [up]
Tasks illustrated in this tutorial
- hstwcs
- plotfunc
- xytoreg
- fluxcalc
Using hstwcs
Display on ds9 the image
Using ds9 tool Analysis --> Catalogs ---> 2Mass, load the 2-Mass catalog
Identify whether stars in the center of the image are close (within 10 pixels) to the 2MASS positions.
If not, in the 2MASS catalog tool, File--->Copy to Regions
Note the RA-DEC catalog position of a star nearest to the center.
Note also the XY position of the same star on the image (by visual identification).
eg. 152.75919 -4.69712 (2MASS RADEC)
1030.2345 1277.922
Also save the 2Mass star catalog using File--->Export--->Starbase as SextansA_2m.cat
Use canwcsfix to fix the center coordinates
inaoe> imcopy gtc14aob009_object_849 gtc14aob009_object_849_ref2m
inaoe> canwcsfix gtc14aob009_object_849_ref2m 152.75919 -4.69712 1030.2345 1277.922
epar ccmap to change the xxorder=3 and xyorder=3 yxorder=3 yyorder=3
inaoe> hstwcs gtc14aob009_object_849_ref2m SextansA_2m.cat ccmap+ update+
65 1959.71 1456.952 1947.437 1447.213 152.771869 -4.631445
-------------------------------------------------
XY and RADEC of 41 common stars are stored in gtc14aob009_object_849_ref2m_starcen.coo
-------------------------------------------------
Refsystem: j2000 Coordinates: equatorial FK5
Equinox: J2000.000 Epoch: J2000.00000000 MJD: 51544.50000
Insystem: j2000 Coordinates: equatorial FK5
Equinox: J2000.000 Epoch: J2000.00000000 MJD: 51544.50000
Coords File: gtc14aob009_object_849_ref2m_starcen.coo Image: gtc14aob009_object_849_ref2m
Database: gtc14aob009_object_849_ref2m_starcen.db Solution: gtc14aob009_object_849_ref2m
Ra/Dec or Long/Lat wcs rms: 0.228 0.256 (arcsec arcsec)
Using plotfunc
First extract a spectrum using apall command
inaoe> apall gtc14aob009_object_12_wl_tlt interac+ extras+ resi- lower=-20 upper=20 b_sampl="-50:-40,40:50" t_funct=legendre t_order=5 backgro=median line=1000
onedspec> plotfunc gtc14aob009_object_12_wl_tlt recno=0
The function stored in XY table = plotfunc.dat and ds9 region polygon file = plotfunc.reg
open the spectral image in ds9 and select Regions/load to load plotfunc.reg
To trace a region without continuum:
onedspec> imcopy gtc14aob009_object_12_wl_tlt gtc14aob009_object_12_wl_tlt_ref
onedspec> apall gtc14aob009_object_12_wl_tlt_ref interac+ extras+ resi- lower=-20 upper=20 b_sampl="-50:-40,40:50" t_funct=legendre t_order=5 backgro=median line=1000
onedspec> imcopy gtc14aob009_object_12_wl_tlt gtc14aob009_object_12_wl_tlt_slit17
onedspec> apall gtc14aob009_object_12_wl_tlt_slit17 ref=gtc14aob009_object_12_wl_tlt inter+ trace-
... during aperture edit use the colon command :center 1356 to recenter the trace at the desired pixel.
Using xytoreg
edit SextansA_2m.cat to comment header lines
columns 1 and 2 have Ra and dec respectively.
xytoreg SextansA_2m.cat regfile=SextansA_2m.reg xc=1 yc=2 idcol=0 cirrad=2 color=magenta xyunit=degree sizunit=arcsec
load the SextansA_2m.reg on ds9 using the Region tool.
Using fluxcalc