Skip to content
GitLab
Commit 0232a1e0 authored by Zhang Xin's avatar Zhang Xin
Browse files

update

parent 0dceb7a5
Hide whitespace changes
Inline Side-by-side
genSurveyDFSByPointing.py
View file @ 0232a1e0
...@@ -2,7 +2,7 @@ ...@@ -2,7 +2,7 @@
Author: Zhang Xin zhangx@bao.ac.cn Author: Zhang Xin zhangx@bao.ac.cn
Date: 2025-04-21 09:53:28 Date: 2025-04-21 09:53:28
LastEditors: Zhang Xin zhangx@bao.ac.cn LastEditors: Zhang Xin zhangx@bao.ac.cn
LastEditTime: 2025-05-20 13:06:28 LastEditTime: 2025-05-27 15:55:30
FilePath: /CSST_Survey/Users/zhangxin/Work/SurveyPlan/point/pointing_c9/genSurveyDFSByPointing.py FilePath: /CSST_Survey/Users/zhangxin/Work/SurveyPlan/point/pointing_c9/genSurveyDFSByPointing.py
Description: 这是默认设置,请设置`customMade`, 打开koroFileHeader查看配置 进行设置: https://github.com/OBKoro1/koro1FileHeader/wiki/%E9%85%8D%E7%BD%AE Description: 这是默认设置,请设置`customMade`, 打开koroFileHeader查看配置 进行设置: https://github.com/OBKoro1/koro1FileHeader/wiki/%E9%85%8D%E7%BD%AE
''' '''
...@@ -55,74 +55,106 @@ def dataInit(): ...@@ -55,74 +55,106 @@ def dataInit():
"opid": "101000001593", "opid": "101000001593",
"obs_id": "10100001593", "obs_id": "10100001593",
"instrument": "MSC", "instrument": "MSC",
"project_id": "None", "obs_group": "None",
"file_type": "SCI", "proposal_id": "None",
# "project_id": "None",
# "file_type": "SCI",
"obs_type": "WIDE", "obs_type": "WIDE",
"object": "m31", "object": "m31",
"mode": "OM1", "mode": "OM1",
"dataset": "None",
"starttime": "2031-04-23 16:40:19", "starttime": "2031-04-23 16:40:19",
"endtime": "2031-04-23 16:43:37", "endtime": "2031-04-23 16:43:37",
"duration": 150.0, # "duration": 150.0,
"exptime": 150.0,
"params": { "params": {
"target_ra": 146.4567, "t0": "2031-04-23 16:41:04",
"target_dec": 54.0354,
"T1": 150.0, "T1": 150.0,
"t0": 190.0, "bias_voltage": 0,
"T2": 22.0, "r_rate": 0,
"T4": 20, "h_rate": 0,
"exposure_start_time": "2031-04-23 16:41:04", "expstart": "2031-04-23 16:41:04",
"starttime": "2031-04-23 16:40:19", "nframes": 1,
"exposure_duration": 150.0, "cover": 0,
"endtime": "2031-04-23 16:43:37" "gainlvl": "",
"shutter": "F",
"submode": "",
"ledflag": "",
"ledstat": "",
"ledt01": 0.0,
"ledt02": 0.0,
"ledt03": 0.0,
"ledt04": 0.0,
"ledt05": 0.0,
"ledt06": 0.0,
"ledt07": 0.0,
"ledt08": 0.0,
"ledt09": 0.0,
"ledt10": 0.0,
"ledt11": 0.0,
"ledt12": 0.0,
"ledt13": 0.0,
"ledt14": 0.0,
"ra_obj": 146.4567,
"dec_obj": 54.0354,
"pa_obj": 25.0,
"guide_window": 1,
"guide_duration": 10,
"guide_pos": [1.0,2.0],
}, },
"simresult": { "simresult": {
"jd": 2462979.8618552,
"decH": 54.0354, "jd": 2459766.56586422,
"raH": 146.4567, "dec": 18.9566929133858,
"sat_x": -4006.7755, "ra": 55.716621253406,
"sat_y": 5432.4733, "sat_x": 2071.4244979379,
"sat_z": 321.34, "sat_y": 5303.28628499587,
"sun_x": 126733141.1541, "sat_z": 3639.38736139007,
"sun_y": 74307866.5619, "sun_x": -36018853.1282829,
"sun_z": 32209476.4337, "sun_y": 135579783.188155,
"moon_x": 245820.5699, "sun_z": 58773491.2520049,
"moon_y": 284285.6117, "moon_x": -388299.686481602,
"moon_z": 127374.2175, "moon_y": -18964.2980441153,
"deepArea": 0.0, "moon_z": 20511.371245116,
"veryDeepArea": 0.0, "deep_area": 0.0,
"very_deep_area": 0.0,
"ca": -1.0, "ca": -1.0,
"cb": "0", "cb": 1.0,
"exposure": 150.0, "exptime": 150.0,
"turnAngle": 0.95482, "turn_angle": 1.9229313453104,
"isInSun": "1", "is_in_sun": 1,
"cmg": 0.6479, "cmg": 0.0,
"charge": 97200.0, "charge": 97200.0,
"betweenAngle": 19.1178, "between_angle": 16.2831689171475,
"ssaSec": 0.0, "saa_sec": 0.0,
"skyid": "523664", "skyid": 309072.0,
"boardAngle": -0.6703, "board_angle": 25.0,
"quad_0": 1.0, "quad_0": 0.999859206659046,
"quad_i": 0.004, "quad_i": 0.00671061788580071,
"quad_j": 0.0067, "quad_j": -0.00768959137244524,
"quad_k": 0.003, "quad_k": -0.0133193208638005,
"sky_type": "1", "sky_type": 1.0,
"id": 1220, "id": 0,
"decCenter": 60.6576587868062, "sat_vx": -7.104108,
"raCenter": 182.834345163608, "sat_vy": 2.912855,
"decHigh": 60.7576587868062, "sat_vz": 0.332436,
"decLow": 60.5576587868062, "fov_lt_dec": 20.4566929133858,
"topLeftRa": 182.619403486662, "fov_lt_ra": 47.6784104513669,
"topRightRa": 183.049286840554, "fov_lb_dec": 17.4566929133858,
"bottomLeftRa": 182.620734042365, "fov_lb_ra": 47.6784104513669,
"bottomRightRa": 183.047956284851 "fov_rb_dec": 17.4566929133858,
"fov_rb_ra": 47.9427837789041,
"fov_rt_dec": 20.4566929133858,
"fov_rt_ra": 47.9427837789041
} }
} }
return data return data
if __name__ == "__main__": if __name__ == "__main__":
if len(sys.argv) < 3: if len(sys.argv) < 4:
print('input_file_name output_file_name') print('input_file_name output_file_name datasec')
fn1 = 'E17.5_b17.5_beta_11.6_opt_transtime_1_CMG_1_dp_2_0.25_da_10_Texp_1.5_DEC60_500_0.1_800_1000_+5deg.obsid.csv' fn1 = 'E17.5_b17.5_beta_11.6_opt_transtime_1_CMG_1_dp_2_0.25_da_10_Texp_1.5_DEC60_500_0.1_800_1000_+5deg.obsid.csv'
# fn2 = 'pointing_50_5_n.dat' # fn2 = 'pointing_50_5_n.dat'
...@@ -146,42 +178,35 @@ if __name__ == "__main__": ...@@ -146,42 +178,35 @@ if __name__ == "__main__":
data['obs_id'] = d1['obs_id'] data['obs_id'] = d1['obs_id']
data['opid'] = (str(d1['obs_id'])[0:3]+'0'+str(d1['obs_id'])[3:]) data['opid'] = (str(d1['obs_id'])[0:3]+'0'+str(d1['obs_id'])[3:])
data['instrument'] = 'MSC' data['instrument'] = 'MSC'
data['file_type'] = 'SCI' # data['file_type'] = 'SCI'
if d1['sky_type'] == 1: if d1['sky_type'] == 1:
data['obs_type'] = 'WIDE' data['obs_type'] = 'WIDE'
else: else:
data['deep'] = 'WIDE' data['obs_type'] = 'DEEP'
data['object'] = d1['skyid'] data['object'] = d1['skyid']
data['dataset'] = sys.argv[3]
data['starttime'] = JDTimeToBeijingTime( data['starttime'] = JDTimeToBeijingTime(
d1['jd'])[1].isoformat() d1['jd'])[1].isoformat()
data['endtime'] = JDTimeToBeijingTime( data['endtime'] = JDTimeToBeijingTime(
d1['jd']+d1['exposure']/86400.)[1].isoformat() d1['jd']+d1['exposure']/86400.)[1].isoformat()
data['duration'] = d1['exposure'] data['exptime'] = d1['exposure']
# param # param
data['params']['target_ra'] = d1['ra_2'] data['params']['ra_obj'] = d1['ra_2']
data['params']['target_dec'] = d1['dec_2'] data['params']['dec_obj'] = d1['dec_2']
data['params']['pa_obj'] = d1['pos_angle']
data['params']['T1'] = d1['exposure'] data['params']['T1'] = d1['exposure']
data['params']['t0'] = data['params']['T1']+40 data['params']['t0'] = JDTimeToBeijingTime(
data['params']['T2'] = calculateTransTime(
d1['turnAngle'])-20-1.5*2 # 快门打开和关闭的时间，分别是1.5s
data['params']['T4'] = 20
# exposure_time1, exposure_time = JDTimeToBeijingTime(d1['jd'])
data['params']['exposure_start_time'] = JDTimeToBeijingTime(
d1['jd'])[1].isoformat() d1['jd'])[1].isoformat()
data['params']['starttime'] = JDTimeToBeijingTime( data['params']['expstart'] = data['params']['t0']
d1['jd']-(1+1+1 + data['params']['T2'] + data['params']['T4'])/86400.)[1].isoformat() data['params']['cover'] = 1
# exposure_time - \ data['params']['shutter'] = "T"
# timedelta(seconds=1+1+1 + data['params']['T2'] + data['params']['T4'])
data['params']['exposure_duration'] = d1['exposure']
data['params']['endtime'] = JDTimeToBeijingTime(
d1['jd']+(d1['exposure']+3)/86400.)[1].isoformat() # 曝光完了，3s关快门，总控配电1s不计
# simresult # simresult
data['simresult']['jd'] = d1['jd'] data['simresult']['jd'] = d1['jd']
data['simresult']['decH'] = d1['decH'] data['simresult']['dec'] = d1['decH']
data['simresult']['raH'] = d1['raH'] data['simresult']['ra'] = d1['raH']
data['simresult']['sat_x'] = d1['sat_x_1'] data['simresult']['sat_x'] = d1['sat_x_1']
data['simresult']['sat_y'] = d1['sat_y_1'] data['simresult']['sat_y'] = d1['sat_y_1']
data['simresult']['sat_z'] = d1['sat_z_1'] data['simresult']['sat_z'] = d1['sat_z_1']
...@@ -191,29 +216,35 @@ if __name__ == "__main__": ...@@ -191,29 +216,35 @@ if __name__ == "__main__":
data['simresult']['moon_x'] = d1['moon_x_1'] data['simresult']['moon_x'] = d1['moon_x_1']
data['simresult']['moon_y'] = d1['moon_y_1'] data['simresult']['moon_y'] = d1['moon_y_1']
data['simresult']['moon_z'] = d1['moon_z_1'] data['simresult']['moon_z'] = d1['moon_z_1']
data['simresult']['deepArea'] = d1['deepArea'] data['simresult']['deep_area'] = d1['deepArea']
data['simresult']['veryDeepArea'] = d1['veryDeepArea'] data['simresult']['very_deep_area'] = d1['veryDeepArea']
data['simresult']['ca'] = d1['ca'] data['simresult']['ca'] = d1['ca']
data['simresult']['cb'] = d1['cb'] data['simresult']['cb'] = d1['cb']
data['simresult']['exposure'] = d1['exposure'] data['simresult']['exptime'] = d1['exposure']
data['simresult']['turnAngle'] = d1['turnAngle'] data['simresult']['turn_angle'] = d1['turnAngle']
data['simresult']['isInSun'] = d1['isInSun'] data['simresult']['is_in_sun'] = d1['isInSun']
data['simresult']['cmg'] = d1['cmg'] data['simresult']['cmg'] = d1['cmg']
data['simresult']['charge'] = d1['charge'] data['simresult']['charge'] = d1['charge']
data['simresult']['skyid'] = d1['skyid'] data['simresult']['skyid'] = d1['skyid']
data['simresult']['boardAngle'] = d1['boardAngle'] data['simresult']['board_angle'] = d1['boardAngle']
data['simresult']['quad_0'] = d1['quad_0'] data['simresult']['quad_0'] = d1['quad_0']
data['simresult']['quad_i'] = d1['quad_i'] data['simresult']['quad_i'] = d1['quad_i']
data['simresult']['quad_j'] = d1['quad_j'] data['simresult']['quad_j'] = d1['quad_j']
data['simresult']['quad_k'] = d1['quad_k'] data['simresult']['quad_k'] = d1['quad_k']
data['simresult']['sky_type'] = d1['sky_type'] data['simresult']['sky_type'] = d1['sky_type']
data['simresult']['id'] = d1['id'] data['simresult']['id'] = d1['id']
data['simresult']['decCenter'] = d1['dec_1'] data['simresult']['sat_vx'] = d1['sat_vx']
data['simresult']['raCenter'] = d1['ra_1'] data['simresult']['sat_vy'] = d1['sat_vy']
data['simresult']['topLeftRa'] = 0. data['simresult']['sat_vz'] = d1['sat_vz']
data['simresult']['topRightRa'] = 0.
data['simresult']['bottomLeftRa'] = 0. data['simresult']['fov_lt_dec'] = 0.
data['simresult']['bottomRightRa'] = 0. data['simresult']['fov_lt_ra'] = 0.
data['simresult']['fov_lb_dec'] = 0.
data['simresult']['fov_lb_ra'] = 0.
data['simresult']['fov_rb_dec'] = 0.
data['simresult']['fov_rb_ra'] = 0.
data['simresult']['fov_rt_dec'] = 0.
data['simresult']['fov_rt_ra'] = 0.
all_data.append(data) all_data.append(data)
# break # break
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment