1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
|
# Copyright (c) 2016 Zhixian MA <zxma_sjtu@qq.com>
# MIT license
"""
Basic parameters to be used in point sources simulation.
References
----------
[1] Angular diameter distance,
https://en.wikipedia.org/wiki/Angular_diameter_distance
[2] FlatLambdaCDM,
{http://docs.astropy.org/en/stable/api/astropy.cosmology.
FlatLambdaCDM.html#astropy.cosmology.FlatLambdaCDM}
"""
from astropy.cosmology import FlatLambdaCDM
class PixelParams():
"""
A class to transform cosmology distance to angles or pixels.
Parameters
------------
H0: float
The hubble constant at z = 0, whose unit is km/s/Mpc
Om0: float
The total matter density.
ang_res: float
Angular resolution, i.e. degree per pixel. (May be useless)
ang_total: list
Total angles of the simulated sky region,whose unit is degree (\deg)
z : float
Redshift
scale: float
The real object scale.
Example
------------
>>> PixelParams = PixelParams(img_size=(1024,1024),ang_total=(5,5))
"""
# Hubble constant at z = 0
H0 = 71.0
# Omega0, total matter density
Om0 = 0.27
# Redshift
z = 0.0
# Cosmology calculator
cosmo = 0.0
# angular diameter distance, [Mpc]
dA = 0.0
# angular resolution
ang_res = 0.0
def __init__(self, z=0.0):
self.z = z
self.cosmo = FlatLambdaCDM(H0=self.H0, Om0=self.Om0)
# angular diameter distance, [Mpc]
self.dA = self.cosmo.angular_diameter_distance(self.z).value
def get_angle(self, scale=1.0):
"""
Input real object scale, and output the respect observed
angle, and pixels.
"""
ang = scale / self.dA # [rac]
return ang
def get_scale(self, ang=1.0):
"""
Input real observed scale, and output the respect
real object scale, and pixels.
"""
scale = ang * self.dA # [Mpc]
return scale
|
