text: #   Copyright 2020 AUI, Inc. Washington DC, USA
  #
  #   Licensed under the Apache License, Version 2.0 (the "License");
  #   you may not use this file except in compliance with the License.
  #   You may obtain a copy of the License at
  #
  #       http://www.apache.org/licenses/LICENSE-2.0
  #
  #   Unless required by applicable law or agreed to in writing, software
  #   distributed under the License is distributed on an "AS IS" BASIS,
  #   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
  #   See the License for the specific language governing permissions and
  #   limitations under the License.
  def deconvolve_rotation_measure_clean(img_dataset, deconvolve_parms, storage_parms):
      """
      .. todo::
          This function is not yet implemented
      
      An iterative solver to construct a full-polarization model from an observed image(set) and psf(set).
     
      Sky Model : Per flux component, delta-functions in lambda-square space
      
      Algorithm :
      
      - Step (1) : Transform the cube to lambda-square space
      - Step (2) : Construct a 3D RM-synthesis PSF
      - Step (3) : Run CLEAN based-deconvolution
      - Step (4) : Transform back to frequency space.
     
      Input : Requires an input cube (mfs is a cube with nchan=1)
      
      Output : Cube model image, Error map (Spectral index map)
      
      Returns
      -------
      img_dataset : xarray.core.dataset.Dataset
      """
  
  