STYLE: Fix several spelling errors

cd ITK
find . -type f | fgrep -v ThirdParty | egrep '\.[hc](xx|pp)?$' | \
  xargs sed -i -r \
    -e 's/(o|O)uput/\1utput/g' \
    -e 's/(u|U)sefull/\1seful/g' \
    -e 's/(a|A)sign/\1ssign/g' \
    -e 's/(s|S)waping/\1wapping/g' \
    -e 's/(m|M)easurment/\1easurement/g' \
    -e 's/(r|R)essource/\1esource/g' \
    -e 's/(l|L)enght/\1ength/g' \
    -e 's/the all the/all the/g' \
    -e 's/[nN]-[dD]imensional/n-dimensional/g' \
    -e 's/([^A-Za-z])a n-dimensional/\1an n-dimensional/g' \
    -e 's/([^A-Za-z])a itk/\1an itk/g' \
    -e 's/([^A-Za-z])a an([^A-Za-z])/\1an\2/g'

Also, change "a" -> "an" where appropriate

Examples/DataRepresentation/Image/Image1.cxx

@@ -91,7 +91,7 @@ main(int, char *[])
   // the image grid starts. These two parameters make it possible to process
   // selected regions.
   //
-  // The \code{Index} is represented by a n-dimensional array where each
+  // The \code{Index} is represented by an n-dimensional array where each
   // component is an integer
   // indicating---in topological image coordinates---the initial pixel
   // of the image.

Examples/Filtering/BilateralImageFilter.cxx

@@ -215,7 +215,7 @@ main(int argc, char * argv[])
 
   //  Software Guide : BeginLatex
   //
-  //  The output of the filter is connected here to a intensity rescaler
+  //  The output of the filter is connected here to an intensity rescaler
   //  filter and then to a writer. Invoking \code{Update()} on the writer
   //  triggers the execution of both filters.
   //

Examples/Filtering/SmoothingRecursiveGaussianImageFilter2.cxx

@@ -18,7 +18,7 @@
 
 //  Software Guide : BeginLatex
 //
-//  Setting up a pipeline of $m$ filters in order to smooth an N-dimensional
+//  Setting up a pipeline of $m$ filters in order to smooth an n-dimensional
 //  image may be a lot of work to do for achieving a simple goal. In order to
 //  avoid this inconvenience, a filter packaging this $m$ filters internally
 //  is available. This filter is the

Examples/IO/VisibleHumanStreamReadWrite.cxx

@@ -72,7 +72,7 @@ main(int argc, char * argv[])
   nameGenerator->SetEndIndex(2878);
   nameGenerator->SetIncrementIndex(1);
 
-  // create a ImageIO for the red channel
+  // create an ImageIO for the red channel
   using ImageIOType = itk::RawImageIO<PixelType, 2>;
   auto rimageio = ImageIOType::New();
   rimageio->SetDimensions(0, 2048);
@@ -82,7 +82,7 @@ main(int argc, char * argv[])
   rimageio->SetHeaderSize(rimageio->GetImageSizeInPixels() * 0);
 
 
-  // create a ImageIO for the green channel
+  // create an ImageIO for the green channel
   auto gimageio = ImageIOType::New();
   gimageio->SetDimensions(0, 2048);
   gimageio->SetDimensions(1, 1216);
@@ -91,7 +91,7 @@ main(int argc, char * argv[])
   gimageio->SetHeaderSize(gimageio->GetImageSizeInPixels() * 1);
 
 
-  // create a ImageIO for the blue channel
+  // create an ImageIO for the blue channel
   auto bimageio = ImageIOType::New();
   bimageio->SetDimensions(0, 2048);
   bimageio->SetDimensions(1, 1216);

Examples/Segmentation/GeodesicActiveContourShapePriorLevelSetImageFilter.cxx

@@ -75,7 +75,7 @@
 // potential image.
 //
 // The \doxygen{FastMarchingImageFilter} creates an initial level set using
-// three user specified seed positions and a initial contour radius. Three
+// three user specified seed positions and an initial contour radius. Three
 // seeds are used in this example to facilitate the segmentation of long
 // narrow objects in a smaller number of iterations. The output of the
 // FastMarchingImageFilter is passed as the input to the

Examples/SpatialObjects/BlobSpatialObject.cxx

@@ -20,7 +20,7 @@
 //
 // \index{itk::BlobSpatialObject}
 //
-// \doxygen{BlobSpatialObject} defines an N-dimensional blob.
+// \doxygen{BlobSpatialObject} defines an n-dimensional blob.
 // This class derives from \doxygen{itkPointBasedSpatialObject}.  A blob
 // is defined as a list of points which compose the object.
 //

Examples/SpatialObjects/EllipseSpatialObject.cxx

@@ -20,7 +20,7 @@
 //
 // \index{itk::EllipseSpatialObject}
 //
-// \doxygen{EllipseSpatialObject} defines an n-Dimensional ellipse. Like
+// \doxygen{EllipseSpatialObject} defines an n-dimensional ellipse. Like
 // other spatial objects this class derives from
 // \doxygen{SpatialObject}. Let's start by including the appropriate header
 // file.

Examples/SpatialObjects/GaussianSpatialObject.cxx

@@ -21,7 +21,7 @@
 // \index{itk::GaussianSpatialObject}
 //
 // This example shows how to create a \doxygen{GaussianSpatialObject} which
-// defines a Gaussian in a N-dimensional space. This object is particularly
+// defines a Gaussian in an n-dimensional space. This object is particularly
 // useful to query the value at a point in physical space. Let's begin by
 // including the appropriate header file.
 //

Examples/Statistics/ScalarImageMarkovRandomField1.cxx

@@ -321,7 +321,7 @@ main(int argc, char * argv[])
   // classification of any given pixel. Note that despite the fact that we
   // call this a radius, it is actually the half size of an hypercube. That
   // is, the actual region of influence will not be circular but rather an
-  // N-Dimensional box. For example, a neighborhood radius of 2 in a 3D image
+  // n-dimensional box. For example, a neighborhood radius of 2 in a 3D image
   // will result in a clique of size 5x5x5 pixels, and a radius of 1 will
   // result in a clique of size 3x3x3 pixels.
   //

Modules/Core/Common/include/itkArray.h

@@ -107,7 +107,7 @@ class ITK_TEMPLATE_EXPORT Array : public vnl_vector<TValue>
     }
   }
 
-  /** Set the all the elements of the array to the specified value */
+  /** Set all the elements of the array to the specified value */
   void
   Fill(TValue const & v)
   {

Modules/Core/Common/include/itkNumericTraitsVariableLengthVectorPixel.h

@@ -28,7 +28,7 @@ namespace itk
 {
 /**
  * \brief Define numeric traits for VariableLengthVector.
- * \tparam T Component type of VariableLenghtVector
+ * \tparam T Component type of VariableLengthVector
  *
  * We provide here a generic implementation based on creating types of
  * VariableLengthVector whose components are the types of the NumericTraits from

Modules/Core/Common/include/itkVariableLengthVector.h

@@ -497,7 +497,7 @@ class ITK_TEMPLATE_EXPORT VariableLengthVector
   Self &
   operator=(VariableLengthVectorExpression<TExpr1, TExpr2, TBinaryOp> const & rhs);
 
-  /** Set the all the elements of the array to the specified value.
+  /** Set all the elements of the array to the specified value.
    * \pre This function may be called on empty vectors, it's a no-op.
    */
   void

Modules/Core/Common/include/itkVariableLengthVector.hxx

@@ -314,7 +314,7 @@ VariableLengthVector<TValue>::SetSize(unsigned int sz, TReallocatePolicy realloc
   m_NumElements = sz;
 }
 
-/** Set the all the elements of the array to the specified value */
+/** Set all the elements of the array to the specified value */
 template <typename TValue>
 void
 VariableLengthVector<TValue>::Fill(TValue const & v)

Modules/Core/Common/src/itkLightObject.cxx

@@ -73,7 +73,7 @@ LightObject::InternalClone() const
 }
 
 /**
- * Delete a itk object. This method should always be used to delete an object
+ * Delete an itk object. This method should always be used to delete an object
  * when the new operator was used to create it. Using the C++ delete method
  * will not work with reference counting.
  */

Modules/Core/Common/test/itkCommonTypeTraitsGTest.cxx

@@ -30,7 +30,7 @@
 #include <type_traits>
 
 // info: is_pod (c++20 deprecated) is equivalent to is_trivial && is_standard_layout
-// info: type_traits cannot differentiate if a explicit move constructible/assignmet exist.
+// info: type_traits cannot differentiate if an explicit move constructible/assignmet exist.
 
 TEST(CommonTypeTraits, FixedArrayIsPOD)
 {

Modules/Core/Common/test/itkFilterDispatchTest.cxx

@@ -123,7 +123,7 @@ ExampleImageFilter<TInputImage, TOutputImage>::Update()
 
 
 /**
- * General N-dimensional implementation of example filter.
+ * General n-dimensional implementation of example filter.
  * The Dispatch parameter is not used.  It is just used to control
  * instantiation.
  */

Modules/Core/Common/test/itkTimeProbeTest2.cxx

@@ -92,7 +92,7 @@ itkTimeProbeTest2(int, char *[])
   std::cout << std::endl << "Print a normal report" << std::endl;
   localTimer.Report();
 
-  // Print a expanded report (including nameOfProbe, Iteration, Total, Min, Mean-Min
+  // Print an expanded report (including nameOfProbe, Iteration, Total, Min, Mean-Min
   //                          Mean/Min *100 (%), Mean, Max, Max- Mean, Max/Mean(%),
   //                          Total Diff(:Max - Min) and STD)
   std::cout << std::endl << "Print an expanded report" << std::endl;

Modules/Core/Common/test/itkTimeProbesTest.cxx

@@ -156,13 +156,13 @@ itkTimeProbesTest(int, char *[])
   std::cout << std::endl << "Print a normal report of a specific probe" << std::endl;
   collector.Report("Loop1");
 
-  // Print a expanded report (including nameOfProbe, Iteration, Total, Min, Mean-Min
+  // Print an expanded report (including nameOfProbe, Iteration, Total, Min, Mean-Min
   //                          Mean/Min *100 (%), Mean, Max, Max- Mean, Max/Mean(%),
   //                          Total Diff(:Max - Min) and STD)
   std::cout << std::endl << "Print expanded reports from all probes" << std::endl;
   collector.ExpandedReport();
 
-  // Print a expanded report (including nameOfProbe, Iteration, Total, Min, Mean-Min
+  // Print an expanded report (including nameOfProbe, Iteration, Total, Min, Mean-Min
   //                          Mean/Min *100 (%), Mean, Max, Max- Mean, Max/Mean(%),
   //                          Total Diff(:Max - Min) and STD)
   std::cout << std::endl << "Print an expanded report of a specific probe" << std::endl;
@@ -174,7 +174,7 @@ itkTimeProbesTest(int, char *[])
   std::cout << std::endl << "Print a JSON report of a specific probe" << std::endl;
   collector.JSONReport("o:TransformPhysicalPointToIndex");
 
-  // Test writing to a ostream
+  // Test writing to an ostream
   std::ofstream logfile;
   logfile.open("itkTimeProbesTest.txt");
   collector.Report(logfile);

Modules/Core/ImageAdaptors/test/itkVectorImageTest.cxx

@@ -628,7 +628,7 @@ itkVectorImageTest(int, char * argv[])
   }
 
   {
-    // Now read it as a itk::VectorImage.
+    // Now read it as an itk::VectorImage.
     using VectorImageType = itk::VectorImage<PixelType, Dimension>;
     using ReaderType = itk::ImageFileReader<VectorImageType>;
     auto reader = ReaderType::New();
@@ -672,7 +672,7 @@ itkVectorImageTest(int, char * argv[])
 
 
   {
-    // Now read it as a itk::VectorImage.
+    // Now read it as an itk::VectorImage.
     using VectorImageType = itk::VectorImage<PixelType, Dimension>;
     using ReaderType = itk::ImageFileReader<VectorImageType>;
     auto reader = ReaderType::New();

Modules/Core/ImageFunction/include/itkBSplineInterpolateImageFunction.hxx

@@ -470,7 +470,7 @@ BSplineInterpolateImageFunction<TImageType, TCoordRep, TCoefficientType>::Evalua
   // perform interpolation
   double    interpolated = 0.0;
   IndexType coefficientIndex;
-  // Step through eachpoint in the N-dimensional interpolation cube.
+  // Step through eachpoint in the n-dimensional interpolation cube.
   for (unsigned int p = 0; p < m_MaxNumberInterpolationPoints; ++p)
   {
     double w = 1.0;

Modules/Core/Mesh/include/itkImageToMeshFilter.hxx

@@ -35,7 +35,7 @@ ImageToMeshFilter<TInputImage, TOutputMesh>::ImageToMeshFilter()
 }
 
 /**
- *   Make Ouput
+ *   Make Output
  */
 template <typename TInputImage, typename TOutputMesh>
 DataObject::Pointer ImageToMeshFilter<TInputImage, TOutputMesh>::MakeOutput(DataObjectPointerArraySizeType)

Modules/Core/QuadEdgeMesh/test/itkQuadEdgeMeshEulerOperatorJoinVertexTest.cxx

@@ -442,7 +442,7 @@ itkQuadEdgeMeshEulerOperatorJoinVertexTest(int argc, char * argv[])
     }
   std::cout << "OK" << std::endl;
   // We now consider a boundary edge (here we consider [0, 5] but still
-  // with a an adjacent face.
+  // with an adjacent face.
   //
   //   20 --------- 21 --------- 22 --------- 23 --------- 24
   //    |        __/ |        __/ |        __/ |        __/ |

Modules/Core/SpatialObjects/test/itkImageMaskSpatialObjectGTest.cxx

@@ -195,7 +195,7 @@ TEST(ImageMaskSpatialObject, AxisAlignedBoundingBoxRegionIsRegionOfSinglePixelWh
 
 
 // Tests that the AABB region is equal to the image region when only a single pixel has value 0.
-// (This condition should hold for N-dimensional image regions, with N >= 2, and region size >= 2^N.)
+// (This condition should hold for n-dimensional image regions, with N >= 2, and region size >= 2^N.)
 TEST(ImageMaskSpatialObject, AxisAlignedBoundingBoxRegionIsImageRegionWhenOnlyOnePixelIsZero)
 {
   // Test 2D images:

Modules/Core/Transform/include/itkCompositeTransform.hxx

@@ -844,7 +844,7 @@ CompositeTransform<TParametersValueType, NDimensions>::UpdateTransformParameters
   /* NOTE: We might want to thread this over each sub-transform, if we
    * find we're working with longer lists of sub-transforms that do
    * not implement any threading of their own for UpdateTransformParameters.
-   * Since the plan is for a UpdateTransformParameters functor that is
+   * Since the plan is for an UpdateTransformParameters functor that is
    * user-assignable, we would need a method in the
    * functor to return whether or not it does therading. If all sub-transforms
    * return that they don't thread, we could do each sub-transform in its

Modules/Core/Transform/include/itkMultiTransform.hxx

@@ -284,7 +284,7 @@ MultiTransform<TParametersValueType, NDimensions, NSubDimensions>::UpdateTransfo
   /* NOTE: We might want to thread this over each sub-transform, if we
    * find we're working with longer lists of sub-transforms that do
    * not implement any threading of their own for UpdateTransformParameters.
-   * Since the plan is for a UpdateTransformParameters functor that is
+   * Since the plan is for an UpdateTransformParameters functor that is
    * user-assignable, we would need a method in the
    * functor to return whether or not it does therading. If all sub-transforms
    * return that they don't thread, we could do each sub-transform in its

Modules/Filtering/GPUSmoothing/include/itkGPUDiscreteGaussianImageFilter.hxx

@@ -57,7 +57,7 @@ GPUDiscreteGaussianImageFilter<TInputImage, TOutputImage>::GPUDiscreteGaussianIm
   }
   else
   {
-    itkExceptionMacro("GPUDiscreteGaussianImageFilter only supports n-Dimensional image.");
+    itkExceptionMacro("GPUDiscreteGaussianImageFilter only supports n-dimensional image.");
   }
 }
 

Modules/Filtering/ImageFeature/include/itkBilateralImageFilter.hxx

@@ -123,7 +123,7 @@ template <typename TInputImage, typename TOutputImage>
 void
 BilateralImageFilter<TInputImage, TOutputImage>::BeforeThreadedGenerateData()
 {
-  // Build a small image of the N-dimensional Gaussian used for domain filter
+  // Build a small image of the n-dimensional Gaussian used for domain filter
   //
   // Gaussian image size will be (2*std::ceil(2.5*sigma)+1) x
   // (2*std::ceil(2.5*sigma)+1)

Modules/Filtering/ImageFilterBase/test/itkCastImageFilterTest.cxx

@@ -343,7 +343,7 @@ bool
 TestVectorImageCast2()
 {
   // This function casts a VectorImage<float, 2>
-  // to a Image<Vector<unsigned char, 2>, 2>
+  // to an Image<Vector<unsigned char, 2>, 2>
   std::cout << "Casting from a VectorImage<float, 2> to Image<Vector<unsigned char, 2>, 2> ... ";
 
   using UnsignedCharVectorImageType = itk::Image<itk::Vector<unsigned char, 2>, 2>;

Modules/Filtering/ImageGradient/include/itkGradientImageFilter.hxx

@@ -51,9 +51,9 @@ GradientImageFilter<TInputImage, TOperatorValueType, TOutputValueType, TOutputIm
   delete m_BoundaryCondition;
 }
 
-template <typename TInputImage, typename TOperatorValueType, typename TOuputValue, typename TOuputImage>
+template <typename TInputImage, typename TOperatorValueType, typename TOutputValue, typename TOutputImage>
 void
-GradientImageFilter<TInputImage, TOperatorValueType, TOuputValue, TOuputImage>::OverrideBoundaryCondition(
+GradientImageFilter<TInputImage, TOperatorValueType, TOutputValue, TOutputImage>::OverrideBoundaryCondition(
   ImageBoundaryCondition<TInputImage> * boundaryCondition)
 {
   delete m_BoundaryCondition;

Modules/Filtering/ImageIntensity/test/itkRGBToVectorAdaptImageFilterTest.cxx

@@ -24,7 +24,7 @@
  *  convert an RGBPixel image to an image that has
  *  vector pixel type.
  *
- *  This allows to access an RGB image a an image of vectors.
+ *  This allows access to an RGB image as an image of vectors.
  *
  */
 

Modules/Filtering/ImageStatistics/include/itkSumProjectionImageFilter.h

@@ -46,7 +46,7 @@ namespace itk
 
 namespace Functor
 {
-template <typename TInputPixel, typename TOuputPixel>
+template <typename TInputPixel, typename TOutputPixel>
 class SumAccumulator
 {
 public:
@@ -56,7 +56,7 @@ class SumAccumulator
   inline void
   Initialize()
   {
-    m_Sum = NumericTraits<TOuputPixel>::ZeroValue();
+    m_Sum = NumericTraits<TOutputPixel>::ZeroValue();
   }
 
   inline void
@@ -65,13 +65,13 @@ class SumAccumulator
     m_Sum = m_Sum + input;
   }
 
-  inline TOuputPixel
+  inline TOutputPixel
   GetValue()
   {
     return m_Sum;
   }
 
-  TOuputPixel m_Sum;
+  TOutputPixel m_Sum;
 };
 } // namespace Functor
 

Modules/Filtering/ImageStatistics/test/itkAdaptiveHistogramEqualizationImageFilterTest.cxx

@@ -67,7 +67,7 @@ itkAdaptiveHistogramEqualizationImageFilterTest(int argc, char * argv[])
   ITK_TEST_SET_GET_VALUE(beta, filter->GetBeta());
 
   //
-  //  The output of the filter is connected here to a intensity rescaler filter
+  //  The output of the filter is connected here to an intensity rescaler filter
   //  and then to a writer. Invoking \code{Update()} on the writer triggers the
   //  execution of both filters.
   //

Modules/Filtering/Path/include/itkContourExtractor2DImageFilter.hxx

@@ -362,7 +362,7 @@ ContourExtractor2DImageFilter<TInputImage>::GenerateDataForLabels()
     for (InputPixelType label : allLabels)
     {
       const BoundingBoxType & bbox = labelBoundingBoxes[label];
-      // Compute a extendedRegion that includes one-pixel border on all
+      // Compute an extendedRegion that includes one-pixel border on all
       // sides. However, we don't want the 3-by-3 SquareIterator to be centered
       // in the bottom row or right column of this one-pixel-extended region
       // because then its lower-right 2-by-2 sub-square won't be wholly within

Modules/Filtering/Thresholding/include/itkTriangleThresholdCalculator.hxx

@@ -50,7 +50,7 @@ TriangleThresholdCalculator<THistogram, TOutput>::GenerateData()
 
   // Triangle method needs the maximum and minimum indexes
   // Minimum indexes for this purpose are poorly defined - can't just
-  // take a index with zero entries.
+  // take an index with zero entries.
   double         Mx = itk::NumericTraits<double>::min();
   IndexValueType MxIdx = 0;
 

Modules/IO/BMP/test/itkBMPImageIOTest3.cxx

@@ -105,7 +105,7 @@ itkBMPImageIOTest3(int argc, char * argv[])
     {
       std::cout << "Test failed!" << std::endl;
       std::cout << "An image stored in a lower-left bitmap is different than \
-                   the same image stored in a upper-left bitmap."
+                   the same image stored in an upper-left bitmap."
                 << std::endl;
       return EXIT_FAILURE;
     }

Modules/IO/BMP/test/itkBMPImageIOTest4.cxx

@@ -107,7 +107,7 @@ itkBMPImageIOTest4(int argc, char * argv[])
     {
       std::cout << "Test failed!" << std::endl;
       std::cout << "An image stored in a lower-left bitmap is different than \
-                   the same image stored in a upper-left bitmap."
+                   the same image stored in an upper-left bitmap."
                 << std::endl;
       return EXIT_FAILURE;
     }