Hello
I am trying to apply a WCS Database (created by someone else) to other exposures. I am doing:
mscred
for(i=1;i<=9;i+=1) {
imgets("sciima0000"+i//"[0]","RAOFFSET"
raoff=(-1)*real(imgets.value)
imgets("sciima0000"+i//"[0]","DECOFFSE"
decoff=(-1)*real(imgets.value)
print("RAOFFSET = "//raoff//", DECOFFSET = "//decoff) mscsetwcs("sciima0000"+i//".fits","distortion.db",ra="RA",dec="DEC", equinox="EQUINOX",ra_offset=raoff,dec_offset=decoff)
}
I have the following message into the logfile
Image: sciima0001.fits[im1] Database: distortion.db Solution: im1
Cannot find or decode record im1 in database file distortion.db
Reference point: 7:13:24.62 -29:10:12.0 2000.0 (hours degrees)
Image: sciima0001.fits[im2] Database: distortion.db Solution: im2
Cannot find or decode record im2 in database file distortion.db
Reference point: 7:13:24.62 -29:10:12.0 2000.0 (hours degrees)
....
Image: sciima0009.fits[im4] Database: distortion.db Solution: im4
Cannot find or decode record im4 in database file distortion.db
Reference point: 7:13:24.64 -29:10:12.4 2000.0 (hours degrees)
The part of the date base is:
# Thu 17:02:46 06-Jun-2013
begin im1
xrefmean 964.700203125
yrefmean 1141.792984375
lngmean 5.364454383680553
latmean -69.50680681423612
pixsystem physical
coosystem j2000
projection tnx
lngref 5.363651111111106
latref -69.4999333333331
lngunits hours
latunits degrees
xpixref 173.7071643271124
ypixref 2403.391646719334
geometry general
function polynomial
xishift -3.767293338535515
etashift -47.40723615004112
xmag 0.01945569942499501
ymag 0.01971450124499278
xrotation 359.5396135462445
yrotation 0.4689580328811265
wcsxirms 0.002304323244250163
wcsetarms 0.004517713702252843
xirms 0.00230432324425027
etarms 0.004517713702252851
surface1 11
3. 3.
2. 2.
2. 2.
0. 0.
5. 5.
2044. 2044.
5. 5.
2044. 2044.
-3.767293338535515 -47.40723615004112
0.01945507134593546 1.563299085656958E-4
1.613586631961425E-4 0.01971384089214057
surface2 18
3. 3.
4. 4.
4. 4.
2. 2.
5. 5.
2044. 2044.
5. 5.
2044. 2044.
-0.2000441189104647 -0.01304663181010625
2.493979183483538E-4 2.174054284822841E-5
-1.114794181513601E-7 -2.398957087927592E-8
-4.531542177492621E-12 7.310481620132716E-12
3.137509342682636E-4 3.357755871516543E-5
-2.641687769102021E-8 -8.825184536163208E-9
4.322499461312955E-12 2.050018792243589E-12
-1.623821782398624E-7 -2.478630823986458E-8
7.862705238522332E-12 3.538515044761089E-12
8.649900385974039E-12 4.562723271772565E-12
Any idea what I'm doing wrong?
Thanks a lot in advance