CARTAvis / carta

To CARTA users, this repo holds the CARTA release packages. Please use this repo to log bugs and feature requests. These will be triaged by the development team and prioritised as necessary in the development cycles.
19 stars 0 forks source link

No valid transform was found to match the images #160

Open UseiCHUU opened 2 years ago

UseiCHUU commented 2 years ago

Describe the bug While spatially matching the images, no valid transform was found.

Here is the header of the image that cannot be matched (statistics skipped):

hst_nicmos_stacked.fits

SIMPLE = T / Standard FITS BITPIX = -32 / Floating point (32 bit) NAXIS = 2 NAXIS1 = 144 NAXIS2 = 145 EXTEND = T BSCALE = 1.000000000000E+00 / PHYSICAL = PIXEL*BSCALE + BZERO BZERO = 0.000000000000E+00 BTYPE = Intensity OBJECT = BUNIT = "COUNTS/S" / Brightness (pixel) unit RADESYS = ICRS LONPOLE = 1.800000000000E+02 LATPOLE = -2.067347468645E+01 PC1_1 = 1 PC2_1 = 0 PC1_2 = 0 PC2_2 = 1 CTYPE1 = RA---TAN CRVAL1 = 2.873259180610E+01 CDELT1 = -2.110108949794E-05 CRPIX1 = 47 CUNIT1 = deg CTYPE2 = DEC--TAN CRVAL2 = 2.873259180610E+01 CDELT2 = 2.093289775459E-05 CRPIX2 = 36.5 CUNIT2 = deg EXTNAME = SCI EXTVER = 1 IRAF-TLM = 18:23:50 (12/06/2009) INHERIT = T EXPNAME = N40W04010 FILLCNT = 0 ERRCNT = 0 PODPSFF = F STDCFFF = F STDCFFP = x5569 SAMPNUM = 25 SAMPTIME = 1.026360376000E+03 DELTATIM = 3.023280000000E-01 ROUTTIME = 5.080827060526E+04 TDFTRANS = 0 ENGQUAL = NORMAL

To Reproduce Steps to reproduce the behavior:

HSC NICMOS image, you can find it in http://hst.esac.esa.int/ehst/#home, source mrc 0152-209

Expected behavior The image can be successfully spatially matched as HST WPFC2 images

Platform info (please complete the following information):

Additional context

Here is the header of the HST WFPC2 image on the same target and can be spatially matched:

hst_wfpc_cutout.fits

SIMPLE = T / Standard FITS BITPIX = -32 / Floating point (32 bit) NAXIS = 2 NAXIS1 = 221 NAXIS2 = 213 EXTEND = T BSCALE = 1.000000000000E+00 / PHYSICAL = PIXEL*BSCALE + BZERO BZERO = 0.000000000000E+00 BTYPE = Intensity OBJECT = BUNIT = "COUNTS/S" / Brightness (pixel) unit RADESYS = ICRS LONPOLE = 1.800000000000E+02 LATPOLE = -2.067459341333E+01 PC1_1 = 1.000000000000E+00 PC2_1 = 0.000000000000E+00 PC1_2 = -0.000000000000E+00 PC2_2 = 1.000000000000E+00 CTYPE1 = RA---TAN CRVAL1 = 2.873259180610E+01 CDELT1 = -1.388888888889E-05 CRPIX1 = 86 CUNIT1 = deg CTYPE2 = DEC--TAN CRVAL2 = -2.067441640823E+01 CDELT2 = 1.388888888889E-05 CRPIX2 = 85 CUNIT2 = deg VAFACTOR = 1.000000000000E+00 MIR_REVR = T ORIENTAT = 0.000000000000E+00 FILLCNT = 0 ERRCNT = 0 FPKTTIME = 5.184858260100E+04 LPKTTIME = 5.184858276800E+04 DETECTOR = 1 DEZERO = 3.188088000000E+02 ZP_CORR = 8.615961000000E-03 CTE_1E2 = 0 CTE_1E3 = 0 CTE_1E4 = 0 EXTNAME = SCI EXPNAME = u59l0201r EXTVER = 1 INHERIT = T ROOTNAME = u59l0201r OCRVAL1 = 2.873250000000E+01 OCRVAL2 = -2.067388888889E+01 OCRPIX2 = 424.5 OCRPIX1 = 420 ONAXIS2 = 800 ONAXIS1 = 800 OCD2_2 = 9.875300000000E-06 OCD2_1 = 7.894720000000E-06 OORIENTA = 3.862738769767E+01 OCTYPE1 = RA---TAN OCD1_1 = -9.879860000000E-06 OCD1_2 = 7.891080000000E-06 OCTYPE2 = DEC--TAN OCX10 = 0.000000000000E+00 OCX11 = 4.553200000000E-02 OCY10 = 4.553017896056E-02 OCY11 = 5.075451992452E-05 O_CRVAL1 = 2.873249025000E+01 O_CRVAL2 = -2.067435460000E+01 G_CRVAL1 = 2.873253825833E+01 G_CRVAL2 = -2.067459341333E+01 G_DRA = 1.728300000000E-01 G_DDEC = -8.597280000000E-01 G_RMSRA = 7.511740000000E-01 G_RMSDEC = 6.995460000000E-01 G_NMATCH = 6 T_CRVAL1 = 2.873249025000E+01 T_CRVAL2 = -2.067435460000E+01 T_DRA = 0.000000000000E+00 T_DDEC = 0.000000000000E+00 T_RMSRA = 0.000000000000E+00 T_RMSDEC = 0.000000000000E+00 T_NMATCH = 0 DATE = 2022-07-04T06:21:46.133563 / Date FITS file was written TIMESYS = UTC / Time system for HDU ORIGIN = casacore-3.4.0

kswang1029 commented 2 years ago

@UseiCHUU The HSC NICMOS image, (eg mrc 0152-209 from http://hst.esac.esa.int/ehst/#home) has non-square pixels. Currently the spatial matching can be applied to images with square pixels only.