def mri_head_reco_op_4_channel():
"""Reconstruction operator for 4 channel MRI of a head.
This is a T2 weighted TSE scan of a healthy volunteer.
The reconstruction operator is the sum of the modulus of each channel.
See the data source with DOI `10.5281/zenodo.800525`_ or the
`project webpage`_ for further information.
See Also
--------
mri_head_data_4_channel
References
----------
.. _10.5281/zenodo.800525: https://zenodo.org/record/800525
.. _project webpage: http://imsc.uni-graz.at/mobis/internal/\
platform_aktuell.html
"""
# To get the same rotation as in the reference article
space = odl.uniform_discr(min_pt=[-115.2, -115.2],
max_pt=[115.2, 115.2],
shape=[256, 256],
dtype=complex)
trafo = odl.trafos.FourierTransform(space)
return odl.ReductionOperator(odl.ComplexModulus(space) * trafo.inverse, 4)
def mri_knee_reco_op_8_channel():
"""Reconstruction operator for 8 channel MRI of a knee.
This is an SE measurement of the knee of a healthy volunteer.
The reconstruction operator is the sum of the modulus of each channel.
See the data source with DOI `10.5281/zenodo.800529`_ or the
`project webpage`_ for further information.
See Also
--------
mri_knee_data_8_channel
References
----------
.. _10.5281/zenodo.800529: https://zenodo.org/record/800529
.. _project webpage: http://imsc.uni-graz.at/mobis/internal/\
platform_aktuell.html
"""
# To get the same rotation as in the reference article
space = odl.uniform_discr(min_pt=[-74.88, -74.88],
max_pt=[74.88, 74.88],
shape=[192, 192],
dtype=complex)
trafo = odl.trafos.FourierTransform(space)
return odl.ReductionOperator(odl.ComplexModulus(space) * trafo.inverse, 8)