#==============================================================================
#
# Project:   SharpImage
# Module:    MorphologicalGradientMagnitude.py
# Language:  IronPython
# Author:    Dan Mueller
# Date:      $Date: 2007-07-13 06:30:23 +1000 (Fri, 13 Jul 2007) $
# Revision:  $Revision: 8 $
#
# Copyright (c) Queensland University of Technology (QUT) 2007.
# All rights reserved.
#
# This software is distributed WITHOUT ANY WARRANTY; without even
# the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
# PURPOSE.  See the above copyright notices for more information.
#
#==============================================================================

# Import the base script class
import ImageToImageScript
from ImageToImageScript import *

# Add reference and import Common library
clr.AddReference("ManagedITK.MorphologicalFilters")
from itk import *

class MorphologicalGradientMagnitudeScript(ImageToImageScriptObject):
    # -------------------------------------------------------------------------
    Name = "MorphologicalGradientMagnitude"
    Help = """Computes the gradient magnitude of an image region at each pixel
              using morphological operators."""
    Parameters = """(itkSize) KernelRadius = the radius of the ball. (itkSize(1,1,.))
                    (itkFlatStructuringElement) Kernel = the structuring element. (Ball)"""
    Kernel = None
    KernelRadius = None
    # -------------------------------------------------------------------------

    def ThreadedDoWork(self):
        """ Perform the main functions of the script. """
        """ NOTE: This function is invoked on a background thread. """
        try:
            self.StartedWork()
            self.WriteInputName()
            filterGrad = itkMorphologicalGradientImageFilter.New( self.Input, self.Output )
            self.AddEventHandlersToProcessObject( filterGrad )
            filterGrad.SetInput( self.Input )
            if (self.KernelRadius == None):
                self.KernelRadius = itkSize( self.Input.Dimension )
                self.KernelRadius.Fill( 1 )
            if (self.Kernel == None):
                self.Kernel = itkFlatStructuringElement.Ball( self.KernelRadius );
            self.WriteLineToConsole( "Kernel=" + self.Kernel.ToString() )
            filterGrad.Kernel = self.Kernel
            filterGrad.UpdateLargestPossibleRegion()
            filterGrad.GetOutput( self.Output )
            self.DisconnectInputAndOutput()
            self.DisposeOfObject( filterGrad )
            self.WriteOutputName()
            self.FinishedWork( True )
        except Exception, ex:
            self.HandleException( ex )
            self.FinishedWork( False )
            self.Finalise( )