code: def mri_head_reco_op_32_channel():
      
      # 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, 32)
docs: Reconstruction operator for 32 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.800527`_ or the
      `project webpage`_ for further information.
  
      See Also
      --------
      mri_head_data_32_channel
  
      References
      ----------
      .. _10.5281/zenodo.800529: https://zenodo.org/record/800527
      .. _project webpage: http://imsc.uni-graz.at/mobis/internal/\
  platform_aktuell.html