ENH: Ingest ITKMeshNoise into Modules/Filtering (closes #5174)

Modules/Filtering/MeshNoise/CMakeLists.txt

@@ -0,0 +1 @@
+itk_module_impl()

Modules/Filtering/MeshNoise/README.md

@@ -0,0 +1,19 @@
+# MeshNoise
+
+In-tree ITK module providing additive-noise filters for `itk::Mesh`
+and `itk::QuadEdgeMesh` point coordinates. The flagship class is
+`itk::AdditiveGaussianNoiseMeshFilter`, which perturbs every mesh
+point by an independent Gaussian-distributed offset with
+configurable mean and standard deviation. Useful for testing
+mesh-based registration, surface reconstruction, and shape-analysis
+pipelines under controlled perturbation.
+
+## Origin
+
+Ingested from the standalone remote module
+[**InsightSoftwareConsortium/ITKMeshNoise**](https://github.com/InsightSoftwareConsortium/ITKMeshNoise)
+on 2026-04-28. The upstream repository will be archived read-only
+after this PR merges; it remains reachable at the URL above.
+
+This ingest also resolves the long-open compile error reported in
+[InsightSoftwareConsortium/ITK#5174](https://github.com/InsightSoftwareConsortium/ITK/issues/5174).

Modules/Filtering/MeshNoise/include/itkAdditiveGaussianNoiseMeshFilter.h

@@ -0,0 +1,101 @@
+/*=========================================================================
+ *
+ *  Copyright NumFOCUS
+ *
+ *  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
+ *
+ *         https://www.apache.org/licenses/LICENSE-2.0.txt
+ *
+ *  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.
+ *
+ *=========================================================================*/
+#ifndef itkAdditiveGaussianNoiseMeshFilter_h
+#define itkAdditiveGaussianNoiseMeshFilter_h
+
+#include <itkMeshToMeshFilter.h>
+
+namespace itk
+{
+/** \class AdditiveGaussianNoiseMeshFilter
+ * \brief Add Gaussian noise to the points defining an itkMesh.
+ *
+ * \author Davis Vigneault
+ *
+ * This class adds Gaussian noise with a specified mean and standard
+ * deviation to the coordinates of the points defining an itkMesh.
+ * Mesh topology and other data is passed unaltered.  This may be
+ * useful in testing the robustness of an algorithm to small changes
+ * in the input mesh, augmenting datasets for machine learning, and
+ * counteracting deleterious effects which highly regular regions
+ * of a mesh may occassionally have on mesh processing.
+ *
+ * \ingroup MeshNoise
+ */
+template <typename TInput, typename TOutput = TInput>
+class AdditiveGaussianNoiseMeshFilter : public MeshToMeshFilter<TInput, TOutput>
+{
+
+public:
+  ITK_DISALLOW_COPY_AND_MOVE(AdditiveGaussianNoiseMeshFilter);
+
+  /** Standard class type alias. */
+  using Self = AdditiveGaussianNoiseMeshFilter;
+  using Superclass = MeshToMeshFilter<TInput, TOutput>;
+  using Pointer = SmartPointer<Self>;
+  using ConstPointer = SmartPointer<const Self>;
+
+  using InputMeshType = TInput;
+  using OutputMeshType = TOutput;
+  using InputMeshPointer = typename InputMeshType::Pointer;
+  using OutputMeshPointer = typename OutputMeshType::Pointer;
+
+  /** Type for representing coordinates. */
+  using CoordinateType = typename TInput::CoordinateType;
+
+  /** Method for creation through the object factory. */
+  itkNewMacro(Self);
+
+  /** Run-time type information (and related methods). */
+  itkOverrideGetNameOfClassMacro(AdditiveGaussianNoiseMeshFilter);
+
+  /** Mean of noise. */
+  itkGetConstMacro(Mean, CoordinateType);
+  itkSetMacro(Mean, CoordinateType);
+
+  /** Variance of noise. */
+  itkGetConstMacro(Sigma, CoordinateType);
+  itkSetMacro(Sigma, CoordinateType);
+
+  /** Initialization seed. */
+  itkGetConstMacro(Seed, int);
+  itkSetMacro(Seed, int);
+
+protected:
+  AdditiveGaussianNoiseMeshFilter();
+  ~AdditiveGaussianNoiseMeshFilter() override = default;
+
+  void
+  PrintSelf(std::ostream & os, Indent indent) const override;
+
+  /** Generate Requested Data */
+  void
+  GenerateData() override;
+
+  CoordinateType m_Mean;
+  CoordinateType m_Sigma;
+  int            m_Seed;
+};
+
+} // namespace itk
+
+#ifndef ITK_MANUAL_INSTANTIATION
+#  include "itkAdditiveGaussianNoiseMeshFilter.hxx"
+#endif
+
+#endif

Modules/Filtering/MeshNoise/include/itkAdditiveGaussianNoiseMeshFilter.hxx

@@ -0,0 +1,107 @@
+/*=========================================================================
+ *
+ *  Copyright NumFOCUS
+ *
+ *  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
+ *
+ *         https://www.apache.org/licenses/LICENSE-2.0.txt
+ *
+ *  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.
+ *
+ *=========================================================================*/
+#ifndef itkAdditiveGaussianNoiseMeshFilter_hxx
+#define itkAdditiveGaussianNoiseMeshFilter_hxx
+
+#include "itkNormalVariateGenerator.h"
+
+namespace itk
+{
+
+template <typename TInputMesh, typename TOutputMesh>
+AdditiveGaussianNoiseMeshFilter<TInputMesh, TOutputMesh>::AdditiveGaussianNoiseMeshFilter()
+{
+  this->m_Mean = 0.0;
+  this->m_Sigma = 1.0;
+  this->m_Seed = 0;
+}
+
+template <typename TInputMesh, typename TOutputMesh>
+void
+AdditiveGaussianNoiseMeshFilter<TInputMesh, TOutputMesh>::PrintSelf(std::ostream & os, Indent indent) const
+{
+  Superclass::PrintSelf(os, indent);
+  os << indent << "Mean: " << this->m_Mean << std::endl;
+  os << indent << "Sigma: " << this->m_Sigma << std::endl;
+}
+
+template <typename TInputMesh, typename TOutputMesh>
+void
+AdditiveGaussianNoiseMeshFilter<TInputMesh, TOutputMesh>::GenerateData()
+{
+  using InputPointsContainer = typename TInputMesh::PointsContainer;
+  using OutputPointsContainer = typename TOutputMesh::PointsContainer;
+
+  using InputPointsContainerConstPointer = typename TInputMesh::PointsContainerConstPointer;
+  using OutputPointsContainerPointer = typename TOutputMesh::PointsContainerPointer;
+
+  const InputMeshType * inputMesh = this->GetInput();
+  OutputMeshPointer     outputMesh = this->GetOutput();
+
+  if (!inputMesh)
+  {
+    itkExceptionMacro(<< "Missing Input Mesh");
+  }
+
+  if (!outputMesh)
+  {
+    itkExceptionMacro(<< "Missing Output Mesh");
+  }
+
+  outputMesh->SetBufferedRegion(outputMesh->GetRequestedRegion());
+
+  InputPointsContainerConstPointer inPoints = inputMesh->GetPoints();
+  OutputPointsContainerPointer     outPoints = outputMesh->GetPoints();
+
+  outPoints->Reserve(inputMesh->GetNumberOfPoints());
+  outPoints->Squeeze();
+
+  typename InputPointsContainer::ConstIterator inputPoint = inPoints->Begin();
+  typename OutputPointsContainer::Iterator     outputPoint = outPoints->Begin();
+
+  unsigned int maxDimension = TInputMesh::MaxTopologicalDimension;
+
+  using GeneratorType = itk::Statistics::NormalVariateGenerator;
+  GeneratorType::Pointer generator = GeneratorType::New();
+  generator->Initialize(this->m_Seed);
+
+  while (inputPoint != inPoints->End())
+  {
+    for (unsigned int dim = 0; dim < maxDimension; ++dim)
+    {
+      outputPoint.Value()[dim] = inputPoint.Value()[dim] + generator->GetVariate() * this->m_Sigma + this->m_Mean;
+    }
+    ++inputPoint;
+    ++outputPoint;
+  }
+
+  // Create duplicate references to the rest of data on the mesh
+  this->CopyInputMeshToOutputMeshPointData();
+  this->CopyInputMeshToOutputMeshCellLinks();
+  this->CopyInputMeshToOutputMeshCells();
+  this->CopyInputMeshToOutputMeshCellData();
+
+  for (unsigned int dim = 0; dim < maxDimension; ++dim)
+  {
+    outputMesh->SetBoundaryAssignments(dim, inputMesh->GetBoundaryAssignments(dim));
+  }
+}
+
+} // namespace itk
+
+#endif

Modules/Filtering/MeshNoise/include/itkAdditiveGaussianNoiseQuadEdgeMeshFilter.h

@@ -0,0 +1,95 @@
+/*=========================================================================
+ *
+ *  Copyright NumFOCUS
+ *
+ *  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
+ *
+ *         https://www.apache.org/licenses/LICENSE-2.0.txt
+ *
+ *  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.
+ *
+ *=========================================================================*/
+#ifndef itkAdditiveGaussianNoiseQuadEdgeMeshFilter_h
+#define itkAdditiveGaussianNoiseQuadEdgeMeshFilter_h
+
+#include "itkQuadEdgeMeshToQuadEdgeMeshFilter.h"
+
+namespace itk
+{
+/** \class AdditiveGaussianNoiseQuadEdgeMeshFilter
+ * \brief Add Gaussian noise to the points defining an itkQuadEdgeMesh.
+ *
+ * \author Davis Vigneault
+ *
+ * This class adds Gaussian noise with a specified mean and standard
+ * deviation to the coordinates of the points defining an itkMesh.
+ * Mesh topology and other data is passed unaltered.  This may be
+ * useful in testing the robustness of an algorithm to small changes
+ * in the input mesh, augmenting datasets for machine learning, and
+ * counteracting deleterious effects which highly regular regions
+ * of a mesh may occassionally have on mesh processing.
+ *
+ * \ingroup MeshNoise
+ */
+template <typename TInputMesh, typename TOutputMesh = TInputMesh>
+class AdditiveGaussianNoiseQuadEdgeMeshFilter : public QuadEdgeMeshToQuadEdgeMeshFilter<TInputMesh, TOutputMesh>
+{
+public:
+  ITK_DISALLOW_COPY_AND_MOVE(AdditiveGaussianNoiseQuadEdgeMeshFilter);
+
+  /** Standard class type alias. */
+  using Self = AdditiveGaussianNoiseQuadEdgeMeshFilter;
+  using Superclass = QuadEdgeMeshToQuadEdgeMeshFilter<TInputMesh, TOutputMesh>;
+  using Pointer = SmartPointer<Self>;
+  using ConstPointer = SmartPointer<const Self>;
+
+  /** Type for representing coordinates. */
+  using CoordinateType = typename TInputMesh::CoordinateType;
+
+  /** Method for creation through the object factory. */
+  itkNewMacro(Self);
+
+  /** Run-time type information (and related methods). */
+  itkOverrideGetNameOfClassMacro(AdditiveGaussianNoiseQuadEdgeMeshFilter);
+
+  /** Mean of noise. */
+  itkGetConstMacro(Mean, CoordinateType);
+  itkSetMacro(Mean, CoordinateType);
+
+  /** Variance of noise. */
+  itkGetConstMacro(Sigma, CoordinateType);
+  itkSetMacro(Sigma, CoordinateType);
+
+  /** Initialization seed. */
+  itkGetConstMacro(Seed, int);
+  itkSetMacro(Seed, int);
+
+protected:
+  AdditiveGaussianNoiseQuadEdgeMeshFilter();
+  ~AdditiveGaussianNoiseQuadEdgeMeshFilter() override = default;
+
+  void
+  PrintSelf(std::ostream & os, Indent indent) const override;
+
+  /** Generate Requested Data */
+  void
+  GenerateData() override;
+
+  CoordinateType m_Mean;
+  CoordinateType m_Sigma;
+  int            m_Seed;
+};
+
+} // end namespace itk
+
+#ifndef ITK_MANUAL_INSTANTIATION
+#  include "itkAdditiveGaussianNoiseQuadEdgeMeshFilter.hxx"
+#endif
+
+#endif

Modules/Filtering/MeshNoise/include/itkAdditiveGaussianNoiseQuadEdgeMeshFilter.hxx

@@ -0,0 +1,69 @@
+/*=========================================================================
+ *
+ *  Copyright NumFOCUS
+ *
+ *  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
+ *
+ *         https://www.apache.org/licenses/LICENSE-2.0.txt
+ *
+ *  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.
+ *
+ *=========================================================================*/
+#ifndef itkAdditiveGaussianNoiseQuadEdgeMeshFilter_hxx
+#define itkAdditiveGaussianNoiseQuadEdgeMeshFilter_hxx
+
+#include "itkNormalVariateGenerator.h"
+
+namespace itk
+{
+
+template <typename TInputMesh, typename TOutputMesh>
+AdditiveGaussianNoiseQuadEdgeMeshFilter<TInputMesh, TOutputMesh>::AdditiveGaussianNoiseQuadEdgeMeshFilter()
+{
+  this->m_Mean = 0.0;
+  this->m_Sigma = 1.0;
+  this->m_Seed = 0;
+}
+
+template <typename TInputMesh, typename TOutputMesh>
+void
+AdditiveGaussianNoiseQuadEdgeMeshFilter<TInputMesh, TOutputMesh>::PrintSelf(std::ostream & os, Indent indent) const
+{
+  Superclass::PrintSelf(os, indent);
+  os << indent << "Mean: " << this->m_Mean << std::endl;
+  os << indent << "Sigma: " << this->m_Sigma << std::endl;
+}
+
+template <typename TInputMesh, typename TOutputMesh>
+void
+AdditiveGaussianNoiseQuadEdgeMeshFilter<TInputMesh, TOutputMesh>::GenerateData()
+{
+  typename TInputMesh::ConstPointer inputMesh = this->GetInput();
+  typename TOutputMesh::Pointer     outputMesh = this->GetOutput();
+
+  this->CopyInputMeshToOutputMesh();
+
+  typename TOutputMesh::PointsContainer::Iterator it = outputMesh->GetPoints()->Begin();
+
+  using GeneratorType = itk::Statistics::NormalVariateGenerator;
+  GeneratorType::Pointer generator = GeneratorType::New();
+  generator->Initialize(this->m_Seed);
+
+  while (it != outputMesh->GetPoints()->End())
+  {
+    for (unsigned int d = 0; d < TOutputMesh::MeshTraits::PointDimension; ++d)
+    {
+      it.Value()[d] += (generator->GetVariate() * this->m_Sigma + this->m_Mean);
+    }
+    ++it;
+  }
+}
+} // end namespace itk
+
+#endif