def get_operator(self):
"""
Return the fused observation as an operator.
"""
H_qubic = self.qubic.get_operator()
R_qubic = ReshapeOperator(H_qubic.shapeout, H_qubic.shape[0])
H_planck = self.planck.get_operator()
R_planck = ReshapeOperator(H_planck.shapeout, H_planck.shape[0])
return BlockColumnOperator(
[R_qubic(H_qubic), R_planck(H_planck)], axisout=0)
def get_invntt_operator(self):
"""
Return the inverse covariance matrix of the fused observation
as an operator.
"""
invntt_qubic = self.qubic.get_invntt_operator()
R_qubic = ReshapeOperator(invntt_qubic.shapeout, invntt_qubic.shape[0])
invntt_planck = self.planck.get_invntt_operator()
R_planck = ReshapeOperator(invntt_planck.shapeout,
invntt_planck.shape[0])
return BlockDiagonalOperator([
R_qubic(invntt_qubic(R_qubic.T)),
R_planck(invntt_planck(R_planck.T))
],
axisout=0)
def get_polarizer_operator(self, sampling, scene):
"""
Return operator for the polarizer grid.
When the polarizer is not present a transmission of 1 is assumed
for the detectors on the first focal plane and of 0 for the other.
Otherwise, the signal is split onto the focal planes.
"""
nd = len(self)
nt = len(sampling)
grid = self.detector.quadrant // 4
if scene.kind == 'I':
if self.optics.polarizer:
return HomothetyOperator(1 / 2)
# 1 for the first detector grid and 0 for the second one
return DiagonalOperator(1 - grid,
shapein=(nd, nt),
broadcast='rightward')
if not self.optics.polarizer:
raise NotImplementedError(
'Polarized input is not handled without the polarizer grid.')
z = np.zeros(nd)
data = np.array([z + 0.5, 0.5 - grid, z]).T[:, None, None, :]
return ReshapeOperator((nd, nt, 1), (nd, nt)) * \
DenseBlockDiagonalOperator(data, shapein=(nd, nt, 3))