From 95b45372ba78b7d5f51d398ce93bb68946055a25 Mon Sep 17 00:00:00 2001 From: yan Date: Mon, 15 Apr 2024 14:13:56 +0800 Subject: [PATCH] debug --- csst_mci_sim/csst_mci_sim.py | 68 +++++++++++++++++------------------- 1 file changed, 33 insertions(+), 35 deletions(-) diff --git a/csst_mci_sim/csst_mci_sim.py b/csst_mci_sim/csst_mci_sim.py index a6bbabb..d664334 100644 --- a/csst_mci_sim/csst_mci_sim.py +++ b/csst_mci_sim/csst_mci_sim.py @@ -3554,7 +3554,7 @@ class MCIsimulator(): #### World coordinate system and related parameters ##### hdu_g = fits.ImageHDU(data) - hdu_g.header['WCSAXES']=(np.int16(2) , 'number of World Coordinate System axes') + hdu_g.header['WCSAXES']=(np.int16(2) , 'number of WCS axes') hdu_g.header['CRPIX1'] =( round(float(self.information['CRPIX1']) ,1) , 'x-coordinate of reference pixel') hdu_g.header['CRPIX2'] =( round(float(self.information['CRPIX2']) ,1) , 'y-coordinate of reference pixel') @@ -3914,7 +3914,7 @@ class MCIsimulator(): ############################################################################## ######### Object information ######################################### ofd_g.header['OBJECT'] =(self.information['name_obj'][:30], 'object name') - ofd_g.header['TARGET'] =(self.information['target'][:13], 'target name (hhmmss.s+ddmmss)') + ofd_g.header['TARGET'] =(self.information['target'][:15], 'target name (hhmmss.s+ddmmss)') ####################################################################### ofd_g.header['OBSID'] =(str(obsid) , 'observation ID') ofd_g.header['RA_OBJ'] =(float(self.information['ra_obj']) , 'object RA (deg)') @@ -4060,20 +4060,19 @@ class MCIsimulator(): hdu_g.header['INSTRUME']=('MCI', 'instrument name') hdu_g.header['CHANNEL'] =('C1', 'channel number') hdu_g.header['FILTERNO'] =(filternum, 'filter number') - hdu_g.header['DIFFUSER']=( bool(True) , 'insert diffuser for flat calibration') + hdu_g.header['DIFFUSER']=( bool(True) , 'insert diffuser status for flat calibration') hdu_g.header['FLAMP'] =( np.int16(0) , 'status of flat lamp') - hdu_g.header['MCISTAT'] =( np.int16(0) , 'MCI componetns status parameter') + hdu_g.header['MCISTAT'] =( np.int16(0) , 'MCI components status parameter') hdu_g.header['DATE_OBS'] =(data_time[:21] , 'observation date (yyyy-mm-ddThh:mm:ss.s)') hdu_g.header['OBJECT'] =( 'MCI_obj' , 'object name') - hdu_g.header['TARGET'] =(self.information['target'][:13], 'target name (hhmmss.s+ddmmss)') + hdu_g.header['TARGET'] =(self.information['target'][:15], 'target name (hhmmss.s+ddmmss)') ####################################################################### hdu_g.header['OBSID'] =(str(obsid) , 'observation ID') hdu_g.header['RA_OBJ'] =(float(self.information['ra_obj']) , 'object RA (deg)') hdu_g.header['DEC_OBJ']=(float(self.information['dec_obj']) , 'object Dec (deg)') - - + hdu_g.header['DATE-OBS']=(data_time[:21] , 'observation date (yyyy-mm-ddThh:mm:ss.s)') ##### detector and Filter information ##### hdu_g.header['FILTER'] =(self.filter_g[:6] , 'filter band') @@ -4092,8 +4091,8 @@ class MCIsimulator(): hdu_g.header['DETSIZE'] =(str(self.information['ysize'])+'*'+str( self.information['xsize']), 'detector size') hdu_g.header['DATASECT'] =(str(self.image_g.shape[0])+'*'+str(self.image_g.shape[1]), 'data section') - hdu_g.header['PIXSCAL1']=('0.05' , 'pixel scale for axis 1 (arcsec/pixel)') - hdu_g.header['PIXSCAL2']=('0.05' , 'pixel scale for axis 2 (arcsec/pixel)') + hdu_g.header['PIXSCAL1']=(float(0.05) , 'pixel scale for axis 1 (arcsec/pixel)') + hdu_g.header['PIXSCAL2']=(float(0.05) , 'pixel scale for axis 2 (arcsec/pixel)') hdu_g.header['PIXSIZE1']=(int(10) , 'pixel size for axis 1 (micron)') hdu_g.header['PIXSIZE2']=(int(10) , 'pixel size for axis 2 (micron)') @@ -4102,7 +4101,7 @@ class MCIsimulator(): hdu_g.header['PSCAN1'] =(np.int32(27) , 'horizontal prescan width, per readout channel') - hdu_g.header['PSCAN2'] =(np.int32(0) , 'vertical prescan height, per readout channel') + hdu_g.header['PSCAN2'] =(np.int32(0) , 'vertical prescan width, per readout channel') hdu_g.header['OSCAN1'] =(np.int32(320) , 'horizontal overscan width, per readout channel') hdu_g.header['OSCAN2'] =(np.int32(320) , 'vertical overscan width, per readout channel') @@ -4112,7 +4111,7 @@ class MCIsimulator(): #### World coordinate system information ##### - hdu_g.header['WCSAXES']=(np.int16(2) , 'number of World Coordinate System axes') + hdu_g.header['WCSAXES']=(np.int16(2) , 'number of WCS axes') hdu_g.header['CRPIX1'] =( round(float(self.information['CRPIX1']) ,1) , 'x-coordinate of reference pixel') hdu_g.header['CRPIX2'] =( round(float(self.information['CRPIX2']) ,1) , 'y-coordinate of reference pixel') @@ -4187,7 +4186,7 @@ class MCIsimulator(): #### exposure and shutter information ##### hdu_g.header['EXPSTART']=(np.float64(time2mjd(self.dt)), 'exposure start time (MJD)') hdu_g.header['EXPEND'] =(float(time2mjd(t2)), 'exposure end time (MJD)') - hdu_g.header['EXPTIME'] =(np.float32(self.information['exptime']), 'exposure times (s)') + hdu_g.header['EXPTIME'] =(np.float32(self.information['exptime']), 'exposure time (s)') hdu_g.header['DARKTIME'] =(np.float32(self.information['exptime']), 'dark current time (s)') hdu_g.header['SHTSTAT']=( bool(True) , 'shutter status') @@ -4332,7 +4331,7 @@ class MCIsimulator(): ofd_r.header['FITSSWV']=('mci_sim_0.8.03', 'FITS creating software version') ######### Object information ######################################### ofd_r.header['OBJECT'] =(self.information['name_obj'][:30], 'object name') - ofd_r.header['TARGET'] =(self.information['target'][:13], 'target name (hhmmss.s+ddmmss)') + ofd_r.header['TARGET'] =(self.information['target'][:15], 'target name (hhmmss.s+ddmmss)') ####################################################################### ofd_r.header['OBSID'] =(str(obsid) , 'observation ID') ofd_r.header['RA_OBJ'] =(float(self.information['ra_obj']) , 'object RA (deg)') @@ -4442,20 +4441,19 @@ class MCIsimulator(): hdu_r.header['INSTRUME']=('MCI', 'instrument name') hdu_r.header['CHANNEL'] =('C1', 'channel number') hdu_r.header['FILTERNO'] =(filternum, 'filter number') - hdu_r.header['DIFFUSER']=( bool(True) , 'insert diffuser for flat calibration') + hdu_r.header['DIFFUSER']=( bool(True) , 'insert diffuser status for flat calibration') hdu_r.header['FLAMP'] =( np.int16(0) , 'status of flat lamp') - hdu_r.header['MCISTAT'] =( np.int16(0) , 'MCI componetns status parameter') + hdu_r.header['MCISTAT'] =( np.int16(0) , 'MCI components status parameter') hdu_r.header['DATE_OBS'] =(data_time[:21] , 'observation date (yyyy-mm-ddThh:mm:ss.s)') hdu_r.header['OBJECT'] =( 'MCI_obj' , 'object name') - hdu_r.header['TARGET'] =(self.information['target'][:13], 'target name (hhmmss.s+ddmmss)') + hdu_r.header['TARGET'] =(self.information['target'][:15], 'target name (hhmmss.s+ddmmss)') ####################################################################### hdu_r.header['OBSID'] =(str(obsid) , 'observation ID') hdu_r.header['RA_OBJ'] =(float(self.information['ra_obj']) , 'object RA (deg)') - hdu_r.header['DEC_OBJ']=(float(self.information['dec_obj']) , 'object Dec (deg)') - - + hdu_r.header['DEC_OBJ']=(float(self.information['dec_obj']) , 'object Dec (deg)') + hdu_r.header['DATE-OBS']=(data_time[:21] , 'observation date (yyyy-mm-ddThh:mm:ss.s)') ##### detector and Filter information ##### hdu_r.header['FILTER'] =(self.filter_r[:6] , 'filter band') @@ -4474,8 +4472,8 @@ class MCIsimulator(): hdu_r.header['DETSIZE'] =(str(self.information['ysize'])+'*'+str( self.information['xsize']), 'detector size') hdu_r.header['DATASECT'] =(str(self.image_g.shape[0])+'*'+str(self.image_g.shape[1]), 'data section') - hdu_r.header['PIXSCAL1']=('0.05' , 'pixel scale for axis 1 (arcsec/pixel)') - hdu_r.header['PIXSCAL2']=('0.05' , 'pixel scale for axis 2 (arcsec/pixel)') + hdu_r.header['PIXSCAL1']=(float(0.05) , 'pixel scale for axis 1 (arcsec/pixel)') + hdu_r.header['PIXSCAL2']=(float(0.05) , 'pixel scale for axis 2 (arcsec/pixel)') hdu_r.header['PIXSIZE1']=(int(10) , 'pixel size for axis 1 (micron)') hdu_r.header['PIXSIZE2']=(int(10) , 'pixel size for axis 2 (micron)') @@ -4484,7 +4482,7 @@ class MCIsimulator(): hdu_r.header['PSCAN1'] =(np.int32(27) , 'horizontal prescan width, per readout channel') - hdu_r.header['PSCAN2'] =(np.int32(0) , 'vertical prescan height, per readout channel') + hdu_r.header['PSCAN2'] =(np.int32(0) , 'vertical prescan width, per readout channel') hdu_r.header['OSCAN1'] =(np.int32(320) , 'horizontal overscan width, per readout channel') hdu_r.header['OSCAN2'] =(np.int32(320) , 'vertical overscan width, per readout channel') @@ -4494,7 +4492,7 @@ class MCIsimulator(): #### World coordinate system information ##### - hdu_r.header['WCSAXES']=(np.int16(2) , 'number of World Coordinate System axes') + hdu_r.header['WCSAXES']=(np.int16(2) , 'number of WCS axes') hdu_r.header['CRPIX1'] =( round(float(self.information['CRPIX1']) ,1) , 'x-coordinate of reference pixel') hdu_r.header['CRPIX2'] =( round(float(self.information['CRPIX2']) ,1) , 'y-coordinate of reference pixel') @@ -4569,7 +4567,7 @@ class MCIsimulator(): #### exposure and shutter information ##### hdu_r.header['EXPSTART']=(np.float64(time2mjd(self.dt)), 'exposure start time (MJD)') hdu_r.header['EXPEND'] =(float(time2mjd(t2)), 'exposure end time (MJD)') - hdu_r.header['EXPTIME'] =(np.float32(self.information['exptime']), 'exposure times (s)') + hdu_r.header['EXPTIME'] =(np.float32(self.information['exptime']), 'exposure time (s)') hdu_r.header['DARKTIME'] =(np.float32(self.information['exptime']), 'dark current time (s)') hdu_r.header['SHTSTAT']=( bool(True) , 'shutter status') @@ -4713,7 +4711,7 @@ class MCIsimulator(): ############################################################################## ######### Object information ######################################### ofd_i.header['OBJECT'] =(self.information['name_obj'][:30], 'object name') - ofd_i.header['TARGET'] =(self.information['target'][:13], 'target name (hhmmss.s+ddmmss)') + ofd_i.header['TARGET'] =(self.information['target'][:15], 'target name (hhmmss.s+ddmmss)') ####################################################################### ofd_i.header['OBSID'] =(str(obsid) , 'observation ID') ofd_i.header['RA_OBJ'] =(float(self.information['ra_obj']) , 'object RA (deg)') @@ -4828,19 +4826,19 @@ class MCIsimulator(): hdu_i.header['INSTRUME']=('MCI', 'instrument name') hdu_i.header['CHANNEL'] =('C1', 'channel number') hdu_i.header['FILTERNO'] =(filternum, 'filter number') - hdu_i.header['DIFFUSER']=( bool(True) , 'insert diffuser for flat calibration') + hdu_i.header['DIFFUSER']=( bool(True) , 'insert diffuser status for flat calibration') hdu_i.header['FLAMP'] =( np.int16(0) , 'status of flat lamp') - hdu_i.header['MCISTAT'] =( np.int16(0) , 'MCI componetns status parameter') + hdu_i.header['MCISTAT'] =( np.int16(0) , 'MCI components status parameter') hdu_i.header['DATE_OBS'] =(data_time[:21] , 'observation date (yyyy-mm-ddThh:mm:ss.s)') hdu_i.header['OBJECT'] =( 'MCI_obj' , 'object name') - hdu_i.header['TARGET'] =(self.information['target'][:13], 'target name (hhmmss.s+ddmmss)') + hdu_i.header['TARGET'] =(self.information['target'][:15], 'target name (hhmmss.s+ddmmss)') ####################################################################### hdu_i.header['OBSID'] =(str(obsid) , 'observation ID') hdu_i.header['RA_OBJ'] =(float(self.information['ra_obj']) , 'object RA (deg)') hdu_i.header['DEC_OBJ']=(float(self.information['dec_obj']) , 'object Dec (deg)') - + hdu_i.header['DATE-OBS']=(data_time[:21] , 'observation date (yyyy-mm-ddThh:mm:ss.s)') ##### detector and Filter information ##### @@ -4857,11 +4855,11 @@ class MCIsimulator(): - hdu_i.header['DETSIZE'] =(str(self.information['ysize'])+'*'+str( self.information['xsize']), 'detector size') + hdu_i.header['DETSIZE'] =(str(self.information['ysize'])+'*'+str( self.information['xsize']), 'detector size') hdu_i.header['DATASECT'] =(str(self.image_g.shape[0])+'*'+str(self.image_g.shape[1]), 'data section') - hdu_i.header['PIXSCAL1']=('0.05' , 'pixel scale for axis 1 (arcsec/pixel)') - hdu_i.header['PIXSCAL2']=('0.05' , 'pixel scale for axis 2 (arcsec/pixel)') + hdu_i.header['PIXSCAL1']=(float(0.05) , 'pixel scale for axis 1 (arcsec/pixel)') + hdu_i.header['PIXSCAL2']=(float(0.05) , 'pixel scale for axis 2 (arcsec/pixel)') hdu_i.header['PIXSIZE1']=(int(10) , 'pixel size for axis 1 (micron)') hdu_i.header['PIXSIZE2']=(int(10) , 'pixel size for axis 2 (micron)') @@ -4870,7 +4868,7 @@ class MCIsimulator(): hdu_i.header['PSCAN1'] =(np.int32(27) , 'horizontal prescan width, per readout channel') - hdu_i.header['PSCAN2'] =(np.int32(0) , 'vertical prescan height, per readout channel') + hdu_i.header['PSCAN2'] =(np.int32(0) , 'vertical prescan width, per readout channel') hdu_i.header['OSCAN1'] =(np.int32(320) , 'horizontal overscan width, per readout channel') hdu_i.header['OSCAN2'] =(np.int32(320) , 'vertical overscan width, per readout channel') @@ -4880,7 +4878,7 @@ class MCIsimulator(): #### World coordinate system information ##### - hdu_i.header['WCSAXES']=(np.int16(2) , 'number of World Coordinate System axes') + hdu_i.header['WCSAXES']=(np.int16(2) , 'number of WCS axes') hdu_i.header['CRPIX1'] =( round(float(self.information['CRPIX1']) ,1) , 'x-coordinate of reference pixel') hdu_i.header['CRPIX2'] =( round(float(self.information['CRPIX2']) ,1) , 'y-coordinate of reference pixel') @@ -4955,7 +4953,7 @@ class MCIsimulator(): #### exposure and shutter information ##### hdu_i.header['EXPSTART']=(np.float64(time2mjd(self.dt)), 'exposure start time (MJD)') hdu_i.header['EXPEND'] =(float(time2mjd(t2)), 'exposure end time (MJD)') - hdu_i.header['EXPTIME'] =(np.float32(self.information['exptime']), 'exposure times (s)') + hdu_i.header['EXPTIME'] =(np.float32(self.information['exptime']), 'exposure time (s)') hdu_i.header['DARKTIME'] =(np.float32(self.information['exptime']), 'dark current time (s)') hdu_i.header['SHTSTAT']=( bool(True) , 'shutter status') -- GitLab