to Mr. valdes
omizzolo wrote on Nov 06, 2006
I am working on mosaic images obtained with the 90prime ad Kitt Peak. I obtained an astrometric solution for one mosaic and now i would apply this solution to the other mosaic of the same field. I tried to use the iraf task mscimage to do this, using a non-mef image as reference image (the image with the astrometric solution), but after launching i obtained this message:
WCS reference image is 39masckwcs_nomef.fits
Resampling jun200046mask.fits[im1] ...
x_images.e(425) malloc: *** vm_allocate(size=698503168) failed (error code=3)
x_images.e(425) malloc: *** error: can't allocate region
x_images.e(425) malloc: *** set a breakpoint in szone_error to debug
PANIC in `/scisoft/iraf/iraf/bin.macosx/x_images.e': Out of memory
ERROR on line 345: bus error
mscimage (mode=h)
mscimage (mode=h).
The parameters i used in mscimage (from epat mscimage) are the following:
PACKAGE = mscred
TASK = mscimage
input = jun200046mask.fits List of input mosaic exposures
output = 46mask_nomef.fits List of output images
(format = mef) Output format (image|mef)
(pixmask= yes) Create pixel mask?
(verbose= )_.verbose) Verbose output?
# Output WCS parameters
(wcssour= image) Output WCS source (image|parameters)
(referen= 39masckwcs_nomef.fits) Reference image
(ra = INDEF) RA of tangent point (hours)
(dec = INDEF) DEC of tangent point (degrees)
(scale = INDEF) Scale (arcsec/pixel)
(rotatio= INDEF) Rotation of DEC from N to E (degrees)
# Resampling parmeters
(blank = 0.) Blank value
(interpo= linear) Interpolant for data
(minterp= linear) Interpolant for mask
(boundar= reflect) Boundary extension
(constan= 0.) Constant boundary extension value
(fluxcon= no) Preserve flux per unit area?
(ntrim = 8) Edge trim in each extension
(nxblock= INDEF) X dimension of working block size in pixels
(nyblock= INDEF) Y dimension of working block size in pixels
# Geometric mapping parameters
(interac= no) Fit mapping interactively?
(nx = 10) Number of x grid points
(ny = 20) Number of y grid points
(fitgeom= general) Fitting geometry
(xxorder= 4) Order of x fit in x
(xyorder= 4) Order of x fit in y
(xxterms= half) X fit cross terms type
(yxorder= 4) Order of y fit in x
(yyorder= 4) Order of y fit in y
(yxterms= half) Y fit cross terms type
(fd_in = )
(fd_ext = )
(fd_coor= )
(mode = ql)
Can you halp me to solve this problem? Thanks a lot
Alessandro Omizzolo
WCS reference image is 39masckwcs_nomef.fits
Resampling jun200046mask.fits[im1] ...
x_images.e(425) malloc: *** vm_allocate(size=698503168) failed (error code=3)
x_images.e(425) malloc: *** error: can't allocate region
x_images.e(425) malloc: *** set a breakpoint in szone_error to debug
PANIC in `/scisoft/iraf/iraf/bin.macosx/x_images.e': Out of memory
ERROR on line 345: bus error
mscimage (mode=h)
mscimage (mode=h).
The parameters i used in mscimage (from epat mscimage) are the following:
PACKAGE = mscred
TASK = mscimage
input = jun200046mask.fits List of input mosaic exposures
output = 46mask_nomef.fits List of output images
(format = mef) Output format (image|mef)
(pixmask= yes) Create pixel mask?
(verbose= )_.verbose) Verbose output?
# Output WCS parameters
(wcssour= image) Output WCS source (image|parameters)
(referen= 39masckwcs_nomef.fits) Reference image
(ra = INDEF) RA of tangent point (hours)
(dec = INDEF) DEC of tangent point (degrees)
(scale = INDEF) Scale (arcsec/pixel)
(rotatio= INDEF) Rotation of DEC from N to E (degrees)
# Resampling parmeters
(blank = 0.) Blank value
(interpo= linear) Interpolant for data
(minterp= linear) Interpolant for mask
(boundar= reflect) Boundary extension
(constan= 0.) Constant boundary extension value
(fluxcon= no) Preserve flux per unit area?
(ntrim = 8) Edge trim in each extension
(nxblock= INDEF) X dimension of working block size in pixels
(nyblock= INDEF) Y dimension of working block size in pixels
# Geometric mapping parameters
(interac= no) Fit mapping interactively?
(nx = 10) Number of x grid points
(ny = 20) Number of y grid points
(fitgeom= general) Fitting geometry
(xxorder= 4) Order of x fit in x
(xyorder= 4) Order of x fit in y
(xxterms= half) X fit cross terms type
(yxorder= 4) Order of y fit in x
(yyorder= 4) Order of y fit in y
(yxterms= half) Y fit cross terms type
(fd_in = )
(fd_ext = )
(fd_coor= )
(mode = ql)
Can you halp me to solve this problem? Thanks a lot
Alessandro Omizzolo
Francisco Valdes wrote on Nov 06, 2006
Hello Alessandro,
It is not clear what you mean by obtaining a astrometric solution. If this means you used msctpeak on one image then this will produce some solution in a database file. To apply the database file to another image you use mscsetwcs. This creates the same higher order terms but sets the origin of the solution to some keyword, usually RA/DEC, to account for a different telescope pointing. Then you would usually need to adjust this zeropoint with msczero or msccmatch because the telescope pointing is often not precise. Once you do this then you can use mscimage to resample the image.
What is probably happening is that you have your reference with a reasonable WCS but the other image doesn't have the right WCS. In this case the resampling may try to produce something that is way off and then think it is trying to make a very large image.
If this is not the case and then I do not know the source of your error. I you were to send me some image header listings of the reference and extension headers I could see if I see the problem.
Yours,
Frank Valdes
It is not clear what you mean by obtaining a astrometric solution. If this means you used msctpeak on one image then this will produce some solution in a database file. To apply the database file to another image you use mscsetwcs. This creates the same higher order terms but sets the origin of the solution to some keyword, usually RA/DEC, to account for a different telescope pointing. Then you would usually need to adjust this zeropoint with msczero or msccmatch because the telescope pointing is often not precise. Once you do this then you can use mscimage to resample the image.
What is probably happening is that you have your reference with a reasonable WCS but the other image doesn't have the right WCS. In this case the resampling may try to produce something that is way off and then think it is trying to make a very large image.
If this is not the case and then I do not know the source of your error. I you were to send me some image header listings of the reference and extension headers I could see if I see the problem.
Yours,
Frank Valdes
Last post on Nov 06, 2006