Commit 6904f2fc authored by xin's avatar xin
Browse files

fix pointing bug in wcs for astrometry

parent a93829b0

ObservationSim/ObservationSim.py

+4 −4
Original line number Diff line number Diff line
@@ -258,11 +258,11 @@ class Observation(object): 
                    raise ValueError("Unknown shear input")

                

                # chip_output.logger.info("debug point #4")

                h_ext = generateExtensionHeader(

                header_wcs = generateExtensionHeader(

                    xlen=chip.npix_x,

                    ylen=chip.npix_y,

                    ra=pointing.ra,

                    dec=pointing.dec,

                    ra=ra_cen,

                    dec=dec_cen,

                    pa=pointing.img_pa.deg,

                    gain=chip.gain,

                    readout=chip.read_noise,
@@ -273,7 +273,7 @@ class Observation(object): 
                    col_num=chip.colID,

                    extName='raw')



                pos_img, offset, local_wcs, real_wcs = obj.getPosImg_Offset_WCS(img=chip.img, fdmodel=self.fd_model, chip=chip, verbose=False, img_header=h_ext)

                pos_img, offset, local_wcs, real_wcs = obj.getPosImg_Offset_WCS(img=chip.img, fdmodel=self.fd_model, chip=chip, verbose=False, img_header=header_wcs)

                if pos_img.x == -1 or pos_img.y == -1:

                    # Exclude object which is outside the chip area (after field distortion)

                    # print("obj missed!!")