morphology.hpp

Go to the documentation of this file.

/* Copyright (c) 2023 Otto Link. Distributed under the terms of the GNU General
   Public License. The full license is in the file LICENSE, distributed with
   this software. */
 
#pragma once
 
#include "highmap/algebra.hpp"
#include "highmap/array.hpp"
 
namespace hmap
{
 
enum DistanceTransformType : int
{
  DT_EXACT,     
  DT_APPROX,    
  DT_MANHATTAN, 
  DT_JFA,       
};
 
Array area_remove(const Array &array,
                  float        threshold_size,
                  float        background_value = 0.f,
                  float        fill_value = 0.f);
 
Array border(const Array &array, int ir);
 
Array closing(const Array &array, int ir);
 
Array closing_by_reconstruction(const Array &array,
                                int          ir,
                                float        k_smooth_max = 0.f);
 
Array contour_smoothing(const Array &array,
                        int          ir,
                        float        transition_ratio = 0.1f);
 
Array dilation(const Array &array, int ir);
 
Array dilation_expand_border_only(const Array &array, int ir);
 
Array dilation_expand_min_value_border_only(const Array &array);
 
Array distance_transform(const Array &array,
                         bool         return_squared_distance = false);
 
Array distance_transform_approx(const Array &array,
                                bool         return_squared_distance = false);
 
Array distance_transform_manhattan(const Array &array,
                                   bool return_squared_distance = false);
 
Array distance_transform_with_closest(const Array     &array,
                                      Mat<glm::ivec2> &closest,
                                      bool return_squared_distance = false);
 
Array erosion(const Array &array, int ir);
 
void flood_fill(Array &array,
                int    i,
                int    j,
                float  fill_value = 1.f,
                float  background_value = 0.f);
 
Array morphological_black_hat(const Array &array, int ir);
 
Array morphological_gradient(const Array &array, int ir);
 
Array morphological_laplacian(const Array &array, int ir);
 
Array morphological_top_hat(const Array &array, int ir);
 
Array opening(const Array &array, int ir);
 
Array opening_by_reconstruction(const Array &array,
                                int          ir,
                                float        k_smooth_min = 0.f);
 
Array reconstruction_by_dilation(const Array &marker,
                                 const Array &mask,
                                 int          ir,
                                 float        k_smooth_min = 0.f);
 
Array reconstruction_by_erosion(const Array &marker,
                                const Array &mask,
                                int          ir,
                                float        k_smooth_max = 0.f);
 
Array relative_distance_from_skeleton(const Array &array,
                                      int          ir_search,
                                      bool         zero_at_borders = true,
                                      int          ir_erosion = 1);
 
Array relative_distance_from_skeleton(const Array &array,
                                      const Array &skeleton,
                                      int          ir_search,
                                      int ir_erosion = 1); 
 
Array skeleton(const Array &array,
               const Array &skeleton,
               bool         zero_at_borders = true);
 
Array remove_endpoints(const Array &array,
                       int          iterations = 1,
                       float        background_value = 0.f);
 
Array signed_curvature_from_distance(const Array &array, int prefilter_ir = 0);
 
Array signed_distance_transform(const Array &array, int prefilter_ir = 0);
 
Array skeleton(const Array &array, bool zero_at_borders = true);
 
} // namespace hmap
 
namespace hmap::gpu
{
 
Array border(const Array &array, int ir);
 
Array closing(const Array &array, int ir);
 
Array closing_by_reconstruction(const Array &array,
                                int          ir,
                                float        k_smooth_max = 0.f);
 
Array contour_smoothing(const Array &array,
                        int          ir,
                        float        transition_ratio = 0.1f);
 
Array dilation(const Array &array, int ir);
 
Array dilation_expand_border_only(const Array &array, int ir);
 
Array distance_transform_jfa(const Array &array,
                             bool         return_squared_distance = false);
 
Array erosion(const Array &array, int ir);
 
Array morphological_black_hat(const Array &array, int ir);
 
Array morphological_gradient(const Array &array, int ir);
 
Array morphological_laplacian(const Array &array, int ir);
 
Array morphological_top_hat(const Array &array, int ir);
 
Array opening(const Array &array, int ir);
 
Array opening_by_reconstruction(const Array &array,
                                int          ir,
                                float        k_smooth_min = 0.f);
 
Array reconstruction_by_dilation(const Array &marker,
                                 const Array &mask,
                                 int          ir,
                                 float        k_smooth_min = 0.f);
 
Array reconstruction_by_erosion(const Array &marker,
                                const Array &mask,
                                int          ir,
                                float        k_smooth_max = 0.f);
 
Array relative_distance_from_skeleton(const Array &array,
                                      int          ir_search,
                                      bool         zero_at_borders = true,
                                      int          ir_erosion = 1);
 
Array relative_distance_from_skeleton(const Array &array,
                                      const Array &skeleton,
                                      int          ir_search,
                                      int          ir_erosion = 1);
 
Array signed_curvature_from_distance(const Array &array, int prefilter_ir = 0);
 
Array signed_distance_transform(const Array &array, int prefilter_ir = 0);
 
Array skeleton(const Array &array, bool zero_at_borders = true);
 
} // namespace hmap::gpu
 
// ==========================================================================
//  Wrapper
// ==========================================================================
 
namespace hmap
{
 
// clang-format off
enum MorphologyOperation : int
{
    MO_BORDER,    
    MO_CLOSING,   
    MO_DILATION,  
    MO_EROSION,   
    MO_OPENING,   
    MO_BLACK_HAT, 
    MO_TOP_HAT,   
    // features).
    MO_GRADIENT,  
    // magnitude).
    MO_LAPLACIAN, 
    // valleys.
    MO_CLOSING_BY_RECONSTRUCTION,
    MO_OPENING_BY_RECONSTRUCTION,
};
// clang-format on
 
Array morphological_operators(const Array        &array,
                              int                 ir,
                              MorphologyOperation operation);
 
} // namespace hmap
 
namespace hmap::gpu
{
 
Array morphological_operators(const Array        &array,
                              int                 ir,
                              MorphologyOperation operation);
 
} // namespace hmap::gpu