本文整理汇总了Python中qubic.QubicAcquisition.get_projection_operator方法的典型用法代码示例。如果您正苦于以下问题:Python QubicAcquisition.get_projection_operator方法的具体用法?Python QubicAcquisition.get_projection_operator怎么用?Python QubicAcquisition.get_projection_operator使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类qubic.QubicAcquisition
的用法示例。
在下文中一共展示了QubicAcquisition.get_projection_operator方法的2个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: test_primary_beam
# 需要导入模块: from qubic import QubicAcquisition [as 别名]
# 或者: from qubic.QubicAcquisition import get_projection_operator [as 别名]
def test_primary_beam():
def primary_beam(theta, phi):
import numpy as np
with settingerr(invalid='ignore'):
return theta <= np.radians(9)
a = QubicAcquisition(150, s, primary_beam=primary_beam)
p = a.get_projection_operator()
assert_equal(p.matrix.ncolmax, 5)
示例2: create_sweeping_pointings
# 需要导入模块: from qubic import QubicAcquisition [as 别名]
# 或者: from qubic.QubicAcquisition import get_projection_operator [as 别名]
delta_az = 15.
angspeed_psi = 0.1
maxpsi = 45.
nsweeps_el = 300
duration = 24 # hours
ts = 20 # seconds
sampling = create_sweeping_pointings(
[racenter, deccenter], duration, ts, angspeed, delta_az, nsweeps_el,
angspeed_psi, maxpsi)
# acquisition model
acq = QubicAcquisition(150, sampling, kind='I', synthbeam_fraction=0.99,
detector_sigma=sigma, detector_fknee=fknee,
detector_fslope=fslope, detector_ncorr=ncorr)
C = acq.get_convolution_peak_operator()
P = acq.get_projection_operator()
H = P * C
# produce the Time-Ordered data
y = H(x0)
# noise
psd = _gaussian_psd_1f(len(acq.sampling), sigma=sigma, fknee=fknee,
fslope=fslope, sampling_frequency=1/ts)
invntt = acq.get_invntt_operator()
noise = acq.get_noise()
# map-making
coverage = P.pT1()
mask = coverage > 10
P = P.restrict(mask, inplace=True)