Dear Co-Users
I want to use "mkpsf" for making a point-spread function at a certain
location of the field of view given in sky-coordinates, but the tool
always creates pointspread functions at other locations. Did you have
similar problems, or do you know where I do something wrong? A consistent
interpretation of the mkpsf output would be that the tool always creates
images in "detector coordinates", although dmlist and ds9 treat them as
"sky-coordinates". If this is true, do you know how to manipulate the
created psf-image to display it in "sky-coordinates" with ds9?
Henric
Example (stars indicate the coordinate mismatch):
mkpsf
input coordinate system (SKY|DET) (SKY): SKY
PSF binning in x direction (0.25:256.0) (INDEF): 0.25
PSF binning in y direction (0.25:256.0) (INDEF): 0.25
PSF size in x direction (2:2048) (INDEF): 20
PSF size in y direction (2:2048) (INDEF): 20
input file (events05.img): events05.img
energy in keV (2.5): 2.5
x (3912): 3912
****
y (4125): 4125
****
angle to rotate image in degrees (INDEF): INDEF
number of pixel points in x or y direction for rotation (1) (9):
PSF library file
(2dpsf/aciss1998-11-052dpsf1N0002.fits):
output file (events05.psf):
CASE: Recursive interpolation...
File events05.psf was created
dmlist events05.psf cols
--------------------------------------------------------------------------------
Columns for Image Block AXAF_2DPSF
--------------------------------------------------------------------------------
ColNo Name Unit Type Range
1 AXAF_2DPSF[20,20] Real4(20x20) -Inf:+Inf
--------------------------------------------------------------------------------
Physical Axis Transforms for Image Block AXAF_2DPSF
--------------------------------------------------------------------------------
Group# Axis#
1 1,2 SKY(X) = (+4283.1768)[pixel] +(+1.0)* ((#1)-(+10.50))
*********
(Y) (+4098.5542) (+1.0) ((#2) (+10.50))
*********
--------------------------------------------------------------------------------
World Coordinate Axis Transforms for Image Block AXAF_2DPSF
--------------------------------------------------------------------------------
Group# Axis#
1 1,2 EQPOS(RA ) = (+72.2527)[deg] +TAN[(-0.000136667)*
(SKY(X)-(+4096.50))]
(DEC) (+45.0071) (+0.000136667) ( (Y)
(+4096.50))
** Henric Krawczynski, email: krawcz@astro.yale.edu, tel: (203) 432 5796
Astronomy Dept.-Yale University-P.O. Box 208101-New Haven, CT 06520-8101 USA
This archive was generated by hypermail 2b29 : Wed May 15 2013 - 01:00:09 EDT