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_msc_sim
Commits
29350402
Commit
29350402
authored
May 29, 2024
by
Fang Yuedong
Browse files
Merge remote-tracking branch 'origin/master' into current_stable_for_tests
parents
38e2c452
6c85d28e
Changes
4
Hide whitespace changes
Inline
Side-by-side
CSST主巡天Cycle6模拟数据说明.pdf
deleted
100644 → 0
View file @
38e2c452
File deleted
README.md
View file @
29350402
# CSST主巡天仿真软件
# CSST主巡天仿真软件
## 重要更新或问题修复:
## 重要更新或问题修复:
*
v3.0版本相关内容:(待补充)
*
2024.05.15: 更新至v3.0版本,新功能具体内容见
[
软件说明文档
](
https://kdocs.cn/l/cjyiU0SXGyn2
)
。
*
2023.07.29: 更新至v2.1版本,内容包括:
*
2023.07.29: 更新至v2.1版本,内容包括:
*
加入杂散光模块,config文件中添加开关:enable_straylight_model: True/False
*
加入杂散光模块,config文件中添加开关:enable_straylight_model: True/False
*
调整WCS定义:x:E, y:-N
*
调整WCS定义:x:E, y:-N
...
@@ -23,20 +23,33 @@
...
@@ -23,20 +23,33 @@
## 使用方法及相关说明
## 使用方法及相关说明
*
软件安装和使用方法可参考:
[
软件说明文档
](
https://kdocs.cn/l/cjyiU0SXGyn2
)
。目前还较为粗糙,
持续更新中
*
软件安装和使用方法可参考:
[
软件说明文档
](
https://kdocs.cn/l/cjyiU0SXGyn2
)
(
持续更新中
)
。
<!---
<!---
*
仿真数据说明文档:
[
C6数据说明文档
](
https://www.kdocs.cn/l/cuRmPf71Ly7v
)
*
仿真数据说明文档:
[
C6数据说明文档
](
https://www.kdocs.cn/l/cuRmPf71Ly7v
)
*
在线版使用说明:
[
CSST网页接口使用说明
](
https://kdocs.cn/l/cl81flAXe1YQ
)
*
在线版使用说明:
[
CSST网页接口使用说明
](
https://kdocs.cn/l/cl81flAXe1YQ
)
-->
-->
*
软件安装及问题反馈表:
[
问题反馈表
](
https://f.kdocs.cn/w/XYcBkE63/
)
*
问题反馈请使用Issues,若无国台gitlab账号,请使用:
[
问题反馈表
](
https://f.kdocs.cn/w/XYcBkE63/
)
##
相关数据
##
仿真数据产品说明
*
PSF采样数据下载:
[
百度云盘
](
https://pan.baidu.com/s/1_2gbO6-d7x001_FMqZ1tfw
)
,提取码:2awt
*
Cycle 9 仿真数据产品:
[
CSST主巡天Cycle9数据产品说明
](
https://pan.cstcloud.cn/s/gHVB8p2TDA
)
*
往期已发布数据产品说明:
[
CSST主巡天往期仿真数据产品说明
](
https://pan.cstcloud.cn/s/zRF7mu3ET5c
)
## 相关输入数据
<!--
*
PSF采样数据下载:
[
百度云盘
](
https://pan.baidu.com/s/1_2gbO6-d7x001_FMqZ1tfw
)
,提取码:2awt
*
测试用数据,C6星系SED PCA主成分,及恒星SED模板库下载:
[
百度云盘
](
https://pan.baidu.com/s/15B5H-r39vbIBIXdpin63Wg
)
,提取码:nxsy
*
测试用数据,C6星系SED PCA主成分,及恒星SED模板库下载:
[
百度云盘
](
https://pan.baidu.com/s/15B5H-r39vbIBIXdpin63Wg
)
,提取码:nxsy
*
滤光片及无缝光谱相关数据:
[
百度云盘
](
https://pan.baidu.com/s/1JROggg8IKf27i1DnwTaPAw
)
,提取码:07v5
*
滤光片及无缝光谱相关数据:
[
百度云盘
](
https://pan.baidu.com/s/1JROggg8IKf27i1DnwTaPAw
)
,提取码:07v5 -->
*
数据链接:https://pan.baidu.com/s/1cVHovDs4cMsDWWNkoZDRDA?pwd=y63x
*
PSF数据库:
*
多色成像: CSSTSim_C9/PSFCube/set1_dynamic
*
无缝光谱: CSSTSim_C9/PSFCube/SLS_PSF_PCA_fp
*
C9星表数据:
*
恒星星表:CSSTSim_C9/starcat_C9
*
星系星表:CSSTSim_C9/galaxycat_C9/cat2CSSTSim_bundle-50sqDeg
*
类星体星表:CSSTSim_C9/galaxycat_C9/qsosed
## 相关会议及workshop
## 相关会议及workshop
...
...
catalog/Catalog_example.py
View file @
29350402
...
@@ -23,7 +23,7 @@ class Catalog(CatalogBase):
...
@@ -23,7 +23,7 @@ class Catalog(CatalogBase):
star_SED_path : str
star_SED_path : str
path to the star SED data
path to the star SED data
objs : list
objs : list
a list of
O
bservation
S
im.
M
ock
O
bject (Star, Galaxy, or Quasar)
a list of
o
bservation
_s
im.
m
ock
_o
bject
s
(Star, Galaxy, or Quasar)
NOTE: must have "obj" list when implement your own Catalog
NOTE: must have "obj" list when implement your own Catalog
Methods
Methods
...
@@ -41,10 +41,10 @@ class Catalog(CatalogBase):
...
@@ -41,10 +41,10 @@ class Catalog(CatalogBase):
----------
----------
config : dict
config : dict
configuration dictionary which is parsed from the input YAML file
configuration dictionary which is parsed from the input YAML file
chip:
O
bservation
S
im.
I
nstrument.Chip
chip:
o
bservation
_s
im.
i
nstrument
s
.Chip
an
O
bservation
S
im.
I
nstrument.Chip instance, can be used to identify the band etc.
an
o
bservation
_s
im.
i
nstrument
s
.Chip instance, can be used to identify the band etc.
**kwargs : dict
**kwargs : dict
other needed input parameters (in key-value pairs), please modify corresponding
any
other needed input parameters (in key-value pairs), please modify corresponding
initialization call in "ObservationSim.py" as you need.
initialization call in "ObservationSim.py" as you need.
Returns
Returns
...
@@ -70,10 +70,10 @@ class Catalog(CatalogBase):
...
@@ -70,10 +70,10 @@ class Catalog(CatalogBase):
"""Read in all objects in from the catalog file(s).
"""Read in all objects in from the catalog file(s).
This is a must implemented method which is used to read in all objects, and
This is a must implemented method which is used to read in all objects, and
then convert them to
O
bservation
S
im.
M
ock
O
bject (Star, Galaxy, or Quasar).
then convert them to
o
bservation
_s
im.
m
ock
_o
bject
s
(Star, Galaxy, or Quasar).
Currently,
Currently,
the model of
O
bservation
S
im.
M
ock
O
bject.Star class requires:
the model of
o
bservation
_s
im.
m
ock
_o
bject
s
.Star class requires:
param["star"] : int
param["star"] : int
specify the object type: 0: galaxy, 1: star, 2: quasar
specify the object type: 0: galaxy, 1: star, 2: quasar
param["id"] : int
param["id"] : int
...
@@ -88,7 +88,7 @@ class Catalog(CatalogBase):
...
@@ -88,7 +88,7 @@ class Catalog(CatalogBase):
load_norm_filt(obj) function must be implemented to load the filter
load_norm_filt(obj) function must be implemented to load the filter
throughput of that particular photometric system
throughput of that particular photometric system
the model of
O
bservation
S
im.
M
ock
O
bject.Galaxy class requires:
the model of
o
bservation
_s
im.
m
ock
_o
bject
s
.Galaxy class requires:
param["star"] : int
param["star"] : int
specify the object type: 0: galaxy, 1: star, 2: quasar
specify the object type: 0: galaxy, 1: star, 2: quasar
param["id"] : int
param["id"] : int
...
@@ -119,7 +119,7 @@ class Catalog(CatalogBase):
...
@@ -119,7 +119,7 @@ class Catalog(CatalogBase):
Sersic index for galaxy bulge component
Sersic index for galaxy bulge component
param['g1'], param['g2']: float
param['g1'], param['g2']: float
Reduced weak lensing shear components (valid for shear type: catalog)
Reduced weak lensing shear components (valid for shear type: catalog)
the model of
O
bservation
S
im.
M
ock
O
bject.Galaxy class requires:
the model of
o
bservation
_s
im.
m
ock
_o
bject
s
.Galaxy class requires:
Currently a Quasar is modeled as a point source, just like a Star.
Currently a Quasar is modeled as a point source, just like a Star.
NOTE: To construct an object, according to its type, just call:
NOTE: To construct an object, according to its type, just call:
...
@@ -133,7 +133,7 @@ class Catalog(CatalogBase):
...
@@ -133,7 +133,7 @@ class Catalog(CatalogBase):
Parameters
Parameters
----------
----------
**kwargs : dict
**kwargs : dict
other needed input parameters (in key-value pairs), please modify corresponding
any
other needed input parameters (in key-value pairs), please modify corresponding
initialization call in "ObservationSim.py" as you need.
initialization call in "ObservationSim.py" as you need.
Returns
Returns
...
@@ -161,7 +161,7 @@ class Catalog(CatalogBase):
...
@@ -161,7 +161,7 @@ class Catalog(CatalogBase):
Parameters
Parameters
----------
----------
obj :
O
bservation
S
im.
M
ock
O
bject
obj :
instance of a particular class in o
bservation
_s
im.
m
ock
_o
bject
s
the object to get SED data for
the object to get SED data for
**kwargs : dict
**kwargs : dict
other needed input parameters (in key-value pairs), please modify corresponding
other needed input parameters (in key-value pairs), please modify corresponding
...
@@ -198,7 +198,7 @@ class Catalog(CatalogBase):
...
@@ -198,7 +198,7 @@ class Catalog(CatalogBase):
Parameters
Parameters
----------
----------
obj :
O
bservation
S
im.
M
ock
O
bject
obj :
instance of a particular class in o
bservation
_s
im.
m
ock
_o
bject
s
the object to get thourghput data for
the object to get thourghput data for
Returns
Returns
...
...
observation_sim/ObservationSim.py
View file @
29350402
...
@@ -45,7 +45,7 @@ class Observation(object):
...
@@ -45,7 +45,7 @@ class Observation(object):
chip
.
shutter_img
=
np
.
ones_like
(
chip
.
img
.
array
)
chip
.
shutter_img
=
np
.
ones_like
(
chip
.
img
.
array
)
else
:
else
:
chip
.
shutter_img
=
effects
.
ShutterEffectArr
(
chip
.
shutter_img
=
effects
.
ShutterEffectArr
(
chip
.
img
,
t_shutter
=
1.3
,
dist_bearing
=
735
,
dt
=
1E-3
)
chip
.
img
,
t_exp
=
pointing
.
exp_time
,
t_shutter
=
1.3
,
dist_bearing
=
735
,
dt
=
1E-3
)
chip
.
prnu_img
=
effects
.
PRNU_Img
(
xsize
=
chip
.
npix_x
,
ysize
=
chip
.
npix_y
,
sigma
=
0.01
,
chip
.
prnu_img
=
effects
.
PRNU_Img
(
xsize
=
chip
.
npix_x
,
ysize
=
chip
.
npix_y
,
sigma
=
0.01
,
seed
=
int
(
self
.
config
[
"random_seeds"
][
"seed_prnu"
]
+
chip
.
chipID
))
seed
=
int
(
self
.
config
[
"random_seeds"
][
"seed_prnu"
]
+
chip
.
chipID
))
...
...
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