UVtools/UVtools.Core/Operations/OperationMorph.cs

/*
 *                     GNU AFFERO GENERAL PUBLIC LICENSE
 *                       Version 3, 19 November 2007
 *  Copyright (C) 2007 Free Software Foundation, Inc. <https://fsf.org/>
 *  Everyone is permitted to copy and distribute verbatim copies
 *  of this license document, but changing it is not allowed.
 */

using Emgu.CV;
using Emgu.CV.CvEnum;
using System;
using System.ComponentModel;
using System.Threading.Tasks;
using UVtools.Core.FileFormats;
using UVtools.Core.Objects;

namespace UVtools.Core.Operations;


public sealed class OperationMorph : Operation
{
    #region Enums
    public enum MorphOperations
    {
        [Description("Erode: Contracts the boundaries within the object")]
        Erode = MorphOp.Erode,

        [Description("Dilate: Expands the boundaries within the object")]
        Dilate = MorphOp.Dilate,

        [Description("Noise removal: Removes small isolated pixels (Erode -> Dilate)")]
        Open = MorphOp.Open,

        [Description("Gap closing: Closes small holes inside the objects (Dilate -> Erode)")]
        Close = MorphOp.Close,
            
        [Description("Gradient: Removes the interior areas of objects and expand the boundaries by half (Dilate - Erode)")]
        Gradient = MorphOp.Gradient,

        // White top-hat transform: This is the difference between the image and its opening.
        [Description("White tophat: Removes small isolated pixels and only return its affected pixels (Image - Noise removal)")]
        WhiteTopHat = MorphOp.Tophat,

        // Black top-hat transform: Difference between the closing and the input image.
        [Description("Black tophat: Closes small holes inside the objects and only return its affected pixels (Gap closing - Image)")]
        BlackTopHat = MorphOp.Blackhat,

        [Description("Hit or miss: Finds pixels in a given kernel pattern")]
        HitMiss = MorphOp.HitMiss,

        [Description("Offset crop: Like erode but discards the outer pixels")]
        OffsetCrop,

        //[Description("Isolate features: Isolates thin features and discards other pixels (Opening - )")]
        //IsolateFeatures,
    }
    #endregion

    #region Members
    private MorphOperations _morphOperation = MorphOperations.Erode;
    private uint _iterationsStart = 1;
    private uint _iterationsEnd = 1;
    private bool _chamfer;
    #endregion

    #region Overrides
    public override string IconClass => "fa-solid fa-dharmachakra";
    public override string Title => "Morph";
    public override string Description =>
        $"Morph Model - " +
        $"Various operations that can be used to change the physical structure of the model or individual layers.";
    public override string ConfirmationText =>
        $"morph model layers {LayerIndexStart} through {LayerIndexEnd}?";

    public override string ProgressTitle =>
        $"Morphing layers {LayerIndexStart} through {LayerIndexEnd}";

    public override string ProgressAction => "Morphed layers";

    #endregion

    #region Properties

    public MorphOp MorphOperationOpenCV
    {
        get
        {
            return _morphOperation switch
            {
                MorphOperations.OffsetCrop => MorphOp.Erode,
                _ => (MorphOp) _morphOperation
            };
        }
    }

    public MorphOperations MorphOperation
    {
        get => _morphOperation;
        set => RaiseAndSetIfChanged(ref _morphOperation, value);
    }
        
    public uint Iterations
    {
        get => IterationsStart;
        set => IterationsStart = IterationsEnd = value;
    }

    public uint IterationsStart
    {
        get => _iterationsStart;
        set => RaiseAndSetIfChanged(ref _iterationsStart, value);
    }

    public uint IterationsEnd
    {
        get => _iterationsEnd;
        set => RaiseAndSetIfChanged(ref _iterationsEnd, value);
    }

    public bool Chamfer
    {
        get => _chamfer;
        set => RaiseAndSetIfChanged(ref _chamfer, value);
    }

    public KernelConfiguration Kernel { get; set; } = new();

    public override string ToString()
    {
        var result = $"[{_morphOperation}] [Iterations: {_iterationsStart}/{_iterationsEnd}] [Chamfer: {_chamfer}]" + LayerRangeString;
        if (!string.IsNullOrEmpty(ProfileName)) result = $"{ProfileName}: {result}";
        return result;
    }

    #endregion

    #region Constructor

    public OperationMorph() { }

    public OperationMorph(FileFormat slicerFile) : base(slicerFile) { }

    #endregion

    #region Equality

    private bool Equals(OperationMorph other)
    {
        return _morphOperation == other._morphOperation && _iterationsStart == other._iterationsStart && _iterationsEnd == other._iterationsEnd && _chamfer == other._chamfer;
    }

    public override bool Equals(object? obj)
    {
        return ReferenceEquals(this, obj) || obj is OperationMorph other && Equals(other);
    }

    public override int GetHashCode()
    {
        return HashCode.Combine((int)_morphOperation, _iterationsStart, _iterationsEnd, _chamfer);
    }

    #endregion

    #region Methods

    protected override bool ExecuteInternally(OperationProgress progress)
    {
        var isFade = Chamfer;
        FileFormat.MutateGetVarsIterationChamfer(
            LayerIndexStart,
            LayerIndexEnd,
            (int)IterationsStart,
            (int)IterationsEnd,
            ref isFade,
            out var iterationSteps,
            out var maxIteration
        );

        Parallel.For(LayerIndexStart, LayerIndexEnd + 1, CoreSettings.GetParallelOptions(progress), layerIndex =>
        {
            int iterations = FileFormat.MutateGetIterationVar(isFade, (int)IterationsStart, (int)IterationsEnd, iterationSteps, maxIteration, LayerIndexStart, (uint)layerIndex);

            using var mat = SlicerFile[layerIndex].LayerMat;
            Execute(mat, iterations);
            SlicerFile[layerIndex].LayerMat = mat;

            progress.LockAndIncrement();
        });

        return !progress.Token.IsCancellationRequested;
    }

    public override bool Execute(Mat mat, params object[]? arguments)
    {
        int iterations = (int) _iterationsStart;
        if (arguments is not null && arguments.Length >= 1)
        {
            iterations = (int) arguments[0];
        }
            
        using var original = mat.Clone();
        var target = GetRoiOrDefault(mat);

        /*if (CoreSettings.CanUseCuda)
        {
            var gpuMat = target.ToGpuMat();
            using var morph = new CudaMorphologyFilter(MorphOperationOpenCV, target.Depth, target.NumberOfChannels, Kernel.Matrix, Kernel.Anchor, iterations);
            morph.Apply(gpuMat, gpuMat);
            gpuMat.Download(target);
        }
        else
        {*/

        var kernel = Kernel.GetKernel(ref iterations);
        CvInvoke.MorphologyEx(target, target, MorphOperationOpenCV, kernel, Kernel.Anchor, iterations, BorderType.Reflect101, default);

        if (_morphOperation == MorphOperations.OffsetCrop)
        {
            var originalRoi = GetRoiOrDefault(original);
            originalRoi.CopyTo(target, target);
        }
        /*else if (_morphOperation == MorphOperations.IsolateFeatures)
        {
            var originalRoi = GetRoiOrDefault(original);
            CvInvoke.Subtract(originalRoi, target, target);
        }*/
        //}

            
        ApplyMask(original, target);
        return true;
    }

    #endregion
}