imalign
douglas brenner wrote on Jan 07, 2009
This works:
imalign input=@fits.txt reference = 61cyg001.fits coords= ctr.txt output = 61cyg*[0,overwrite] boxsize = 25 bigbox = 50
This doesn't
imalign input=@fits.txt reference = 61cyg001.fits coords= ctr.txt output = 61cyg*[0,overwrite] boxsize = 25 bigbox = 50 niterate = 5
ERROR: Number of input and output images not the same
!!!!
imalign input=@fits.txt reference = 61cyg001.fits coords= ctr.txt output = 61cyg*[0,overwrite] boxsize = 25 bigbox = 50
This doesn't
imalign input=@fits.txt reference = 61cyg001.fits coords= ctr.txt output = 61cyg*[0,overwrite] boxsize = 25 bigbox = 50 niterate = 5
ERROR: Number of input and output images not the same
!!!!
Mike Fitzpatrick wrote on Jan 07, 2009
Try turning on the 'verbose' parameter.
This will print out the result of an initial IMCENTROID call that creates the file list for shifting the images. One problem is that 'bigbox' should be an odd-numbered value, but perhaps by turning on iterate there is some other failure in the task and so the shift file list is invalid. Either way, I expect you'll see some sort of error message.
-Mike
This will print out the result of an initial IMCENTROID call that creates the file list for shifting the images. One problem is that 'bigbox' should be an odd-numbered value, but perhaps by turning on iterate there is some other failure in the task and so the shift file list is invalid. Either way, I expect you'll see some sort of error message.
-Mike
douglas brenner wrote on Jan 07, 2009
Same error, no output
anyway to send an attachment?
anyway to send an attachment?
Mike Fitzpatrick wrote on Jan 07, 2009
Other than cut-n-paste, attachments are difficult here.
If you turned on verbose and didn't see anything before the error message, then it means the IMCENTROID task didn't produce any output at all. You can try running this on your own, e.g. from the imalign script the command would ne something like
Try this both with and without setting the niterate and report any differences in the output.
Otherwise, you can type 'd_trace' before running IMALIGN to trace the script execution, give the command again at the end to toggle it off and post the last screenfull of output (my guess is IMSHIFT will be the task that fails).
-Mike
If you turned on verbose and didn't see anything before the error message, then it means the IMCENTROID task didn't produce any output at all. You can try running this on your own, e.g. from the imalign script the command would ne something like
cl> imcentroid @fits.txt 61cyg001.fits ctr.txt shifts=""
boxsize=25 bigbox=50 niterate=5 verbose+
Try this both with and without setting the niterate and report any differences in the output.
Otherwise, you can type 'd_trace' before running IMALIGN to trace the script execution, give the command again at the end to toggle it off and post the last screenfull of output (my guess is IMSHIFT will be the task that fails).
-Mike
douglas brenner wrote on Jan 07, 2009
First this
#Coords Image X-center Err Y-center Err Num
61cyg002.fits 146.549 (0.496) 118.908 (0.713) 1
61cyg003.fits 146.194 (0.395) 119.504 (0.564) 1
61cyg004.fits 145.802 (0.339) 119.762 (0.494) 1
61cyg005.fits 145.467 (0.314) 119.957 (0.454) 1
Warning: failed to converge near (138,126)
Warning: no sources centered in 61cyg006.fits
61cyg007.fits 144.179 (0.332) 120.735 (0.450) 1
...
Then this
#Refcoords Reference X-center Err Y-center Err Num
61cyg001.fits 147.213 (0.682) 118.613 (0.945) 1
#Shifts Image X-shift Err Y-shift Err N Internal
61cyg002.fits 0.664 (0.843) -0.295 (1.184) 1 (INDEF,INDEF)
61cyg003.fits 1.019 (0.788) -0.891 (1.101) 1 (INDEF,INDEF)
61cyg004.fits 1.411 (0.761) -1.149 (1.066) 1 (INDEF,INDEF)
61cyg005.fits 1.746 (0.750) -1.345 (1.048) 1 (INDEF,INDEF)
...
Then this
#Trim_Section = [16:250,1:240]
# Shifting images:
Then this
ERROR: Number of input and output images not the same.
"constant= constant)"
line 64: immatch$imalign.cl
#Coords Image X-center Err Y-center Err Num
61cyg002.fits 146.549 (0.496) 118.908 (0.713) 1
61cyg003.fits 146.194 (0.395) 119.504 (0.564) 1
61cyg004.fits 145.802 (0.339) 119.762 (0.494) 1
61cyg005.fits 145.467 (0.314) 119.957 (0.454) 1
Warning: failed to converge near (138,126)
Warning: no sources centered in 61cyg006.fits
61cyg007.fits 144.179 (0.332) 120.735 (0.450) 1
...
Then this
#Refcoords Reference X-center Err Y-center Err Num
61cyg001.fits 147.213 (0.682) 118.613 (0.945) 1
#Shifts Image X-shift Err Y-shift Err N Internal
61cyg002.fits 0.664 (0.843) -0.295 (1.184) 1 (INDEF,INDEF)
61cyg003.fits 1.019 (0.788) -0.891 (1.101) 1 (INDEF,INDEF)
61cyg004.fits 1.411 (0.761) -1.149 (1.066) 1 (INDEF,INDEF)
61cyg005.fits 1.746 (0.750) -1.345 (1.048) 1 (INDEF,INDEF)
...
Then this
#Trim_Section = [16:250,1:240]
# Shifting images:
Then this
ERROR: Number of input and output images not the same.
"constant= constant)"
line 64: immatch$imalign.cl
Francisco Valdes wrote on Jan 07, 2009
Hi,
Did you ever get this resolved? The script is not very complicated. I would suggest you copy the script, add some debugging statements at appropriate places and run it. To make your own version of the task:
1. copy immatch$imalign.cl immatch$imalign.par home$
2. load the immatch task
3. redefine imalign=home$imalign.cl
This shows putting the copy in you IRAF home. You could put it anywhere you like as long as the redefine has the appropriate path.
If you'd rather send me the two list files you can to valdes@noao.edu. Note in the script that one of the first things to happen is to turn your "output" specification into a file with the SECTIONS command. You might try that manually to make sure it expands to what you expect:
Yours,
Frank Valdes
Did you ever get this resolved? The script is not very complicated. I would suggest you copy the script, add some debugging statements at appropriate places and run it. To make your own version of the task:
1. copy immatch$imalign.cl immatch$imalign.par home$
2. load the immatch task
3. redefine imalign=home$imalign.cl
This shows putting the copy in you IRAF home. You could put it anywhere you like as long as the redefine has the appropriate path.
If you'd rather send me the two list files you can to valdes@noao.edu. Note in the script that one of the first things to happen is to turn your "output" specification into a file with the SECTIONS command. You might try that manually to make sure it expands to what you expect:
sections 61cyg*[0,overwrite] > listYours,
Frank Valdes
Last post on Jan 07, 2009