Skip to content
GitLab
Projects
Groups
Snippets
/
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in
Toggle navigation
Menu
Open sidebar
csst-sims
csst_mci_sim
Commits
5b292754
Commit
5b292754
authored
Apr 14, 2024
by
Yan Zhaojun
Browse files
test
parent
22e2bedd
Pipeline
#4051
failed with stage
in 0 seconds
Changes
1
Pipelines
1
Hide whitespace changes
Inline
Side-by-side
csst_mci_sim/csst_mci_sim.py
View file @
5b292754
...
@@ -1923,7 +1923,7 @@ class MCIsimulator():
...
@@ -1923,7 +1923,7 @@ class MCIsimulator():
z_sat
=
float
(
self
.
orbit_pars
[
self
.
orbit_exp_num
,
3
])
z_sat
=
float
(
self
.
orbit_pars
[
self
.
orbit_exp_num
,
3
])
self
.
log
.
info
(
'load star catlog 111111111111111'
)
wave0
,
zodi0
=
self
.
zodiacal
(
ra
,
dec
,
self
.
TianCe_day
)
# erg/s/cm^2/A/arcsec^2
wave0
,
zodi0
=
self
.
zodiacal
(
ra
,
dec
,
self
.
TianCe_day
)
# erg/s/cm^2/A/arcsec^2
...
@@ -1931,7 +1931,10 @@ class MCIsimulator():
...
@@ -1931,7 +1931,10 @@ class MCIsimulator():
# EarthShine from straylight
# EarthShine from straylight
sl
=
StrayLight
(
self
.
information
[
'dir_path'
],
jtime
=
time_jd
,
sat
=
np
.
array
([
x_sat
,
y_sat
,
z_sat
]),
sl
=
StrayLight
(
self
.
information
[
'dir_path'
],
jtime
=
time_jd
,
sat
=
np
.
array
([
x_sat
,
y_sat
,
z_sat
]),
radec
=
np
.
array
([(
ra
*
u
.
degree
).
value
,
(
dec
*
u
.
degree
).
value
]))
radec
=
np
.
array
([(
ra
*
u
.
degree
).
value
,
(
dec
*
u
.
degree
).
value
]))
self
.
log
.
info
(
'load star catlog 3333333333333333333333333'
)
earth_e
=
sl
.
caculateEarthShineFilter
(
filter
=
'r'
)
earth_e
=
sl
.
caculateEarthShineFilter
(
filter
=
'r'
)
star_e
=
sl
.
caculateStarLightFilter
(
filter
=
'r'
)
star_e
=
sl
.
caculateStarLightFilter
(
filter
=
'r'
)
angle_earth
=
earth_angle
(
time_jd
,
x_sat
,
y_sat
,
z_sat
,
ra
,
dec
)
angle_earth
=
earth_angle
(
time_jd
,
x_sat
,
y_sat
,
z_sat
,
ra
,
dec
)
...
@@ -1941,6 +1944,8 @@ class MCIsimulator():
...
@@ -1941,6 +1944,8 @@ class MCIsimulator():
earthshine_wave0
,
earthshine_flux0
=
ill2flux
(
earth_e
+
star_e
,
self
.
information
[
'dir_pah'
])
earthshine_wave0
,
earthshine_flux0
=
ill2flux
(
earth_e
+
star_e
,
self
.
information
[
'dir_pah'
])
self
.
log
.
info
(
'load star catlog 444444444'
)
# sample as ifs wavelength
# sample as ifs wavelength
wave_mci
=
np
.
linspace
(
2500
,
11000
,
8501
)
#np.arange(2500, 11000, 1)
wave_mci
=
np
.
linspace
(
2500
,
11000
,
8501
)
#np.arange(2500, 11000, 1)
f1
=
interp1d
(
wave0
,
zodi0
)
f1
=
interp1d
(
wave0
,
zodi0
)
...
@@ -1961,7 +1966,7 @@ class MCIsimulator():
...
@@ -1961,7 +1966,7 @@ class MCIsimulator():
########################################################################################
########################################################################################
self
.
log
.
info
(
'load star catlog 5555555555555'
)
#########################################################################################
#########################################################################################
self
.
cal_sky_noise
()
self
.
cal_sky_noise
()
...
@@ -1984,7 +1989,9 @@ class MCIsimulator():
...
@@ -1984,7 +1989,9 @@ class MCIsimulator():
#######################################################################
#######################################################################
if
self
.
TianceEffect
:
if
self
.
TianceEffect
:
self
.
log
.
info
(
'load star catlog 666666666666666'
)
ra_list
=
self
.
star
[
'ra_gaia'
].
tolist
()
ra_list
=
self
.
star
[
'ra_gaia'
].
tolist
()
dec_list
=
self
.
star
[
'dec_gaia'
].
tolist
()
dec_list
=
self
.
star
[
'dec_gaia'
].
tolist
()
pmra_list
=
self
.
star
[
'pmra_gaia'
].
tolist
()
pmra_list
=
self
.
star
[
'pmra_gaia'
].
tolist
()
...
@@ -2022,7 +2029,7 @@ class MCIsimulator():
...
@@ -2022,7 +2029,7 @@ class MCIsimulator():
st_magz
=
[]
st_magz
=
[]
self
.
log
.
info
(
'load star catlog 777777777'
)
#################### generate star image ##########
#################### generate star image ##########
if
self
.
debug
:
if
self
.
debug
:
...
...
Write
Preview
Supports
Markdown
0%
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment