STYLE: Remove p->Initialize() calls directly after p = T::New()

`T::New()` already initializes the object it returns.

Using Notepad++, Replace in Files, doing:

    Find what:  (\w+)([ ]+= .+::New\(\);\r\n)[\r\n]* [ ]+\1->Initialize\(\);
    Replace with: $1$2
    Filters: itk*.* !+\test
    Directory: D:\src\ITK\Modules
    [v] Match case
    (*) Regular expression

Manually removed an obsolete comment saying `Clear elements`.

Modules/Core/Transform/include/itkBSplineTransformInitializer.hxx

@@ -85,8 +85,8 @@ BSplineTransformInitializer<TTransform, TImage>::InitializeTransform() const
   using CoordRepType = typename ImagePointType::CoordRepType;
 
   using PointSetType = PointSet<CoordRepType, SpaceDimension>;
-  auto cornerPoints = PointSetType::New();
-  cornerPoints->Initialize();
+  autocornerPoints = PointSetType::New();
+
 
   using PointType = typename PointSetType::PointType;
   using PointIdentifier = typename PointSetType::PointIdentifier;

Modules/Filtering/BiasCorrection/include/itkN4BiasFieldCorrectionImageFilter.hxx

@@ -498,15 +498,15 @@ N4BiasFieldCorrectionImageFilter<TInputImage, TMaskImage, TOutputImage>::UpdateB
 
   const typename ImporterType::OutputImageType * parametricFieldEstimate = importer->GetOutput();
 
-  PointSetPointer fieldPoints = PointSetType::New();
-  fieldPoints->Initialize();
+  PointSetPointerfieldPoints = PointSetType::New();
+
   auto & pointSTLContainer = fieldPoints->GetPoints()->CastToSTLContainer();
   pointSTLContainer.reserve(numberOfIncludedPixels);
   auto & pointDataSTLContainer = fieldPoints->GetPointData()->CastToSTLContainer();
   pointDataSTLContainer.reserve(numberOfIncludedPixels);
 
-  auto weights = BSplineFilterType::WeightsContainerType::New();
-  weights->Initialize();
+  autoweights = BSplineFilterType::WeightsContainerType::New();
+
   auto & weightSTLContainer = weights->CastToSTLContainer();
   weightSTLContainer.reserve(numberOfIncludedPixels);
 

Modules/Filtering/DisplacementField/include/itkDisplacementFieldToBSplineImageFilter.hxx

@@ -121,8 +121,8 @@ DisplacementFieldToBSplineImageFilter<TInputImage, TInputPointSet, TOutputImage>
     itkExceptionMacro("The number of input points does not match the number of weight elements.");
   }
 
-  auto fieldPoints = InputPointSetType::New();
-  fieldPoints->Initialize();
+  autofieldPoints = InputPointSetType::New();
+
 
   auto weights = WeightsContainerType::New();
 

Modules/Filtering/FastMarching/include/itkFastMarchingImageToNodePairContainerAdaptor.hxx

@@ -125,8 +125,8 @@ FastMarchingImageToNodePairContainerAdaptor<TInput, TOutput, TImage>::SetPointsF
 {
   if (iLabel == Traits::Alive || iLabel == Traits::InitialTrial || iLabel == Traits::Forbidden)
   {
-    NodePairContainerPointer nodes = NodePairContainerType::New();
-    nodes->Initialize();
+    NodePairContainerPointernodes = NodePairContainerType::New();
+
 
     using IteratorType = ImageRegionConstIteratorWithIndex<ImageType>;
     IteratorType it(image, image->GetBufferedRegion());

Modules/Numerics/Optimizersv4/include/itkWindowConvergenceMonitoringFunction.hxx

@@ -93,8 +93,8 @@ WindowConvergenceMonitoringFunction<TScalar>::GetConvergenceValue() const -> Rea
   bspliner->SetNumberOfControlPoints(ncps);
   bspliner->SetNumberOfWorkUnits(1);
 
-  auto energyProfileWindow = EnergyProfileType::New();
-  energyProfileWindow->Initialize();
+  autoenergyProfileWindow = EnergyProfileType::New();
+
 
   for (unsigned int n = 0; n < this->m_WindowSize; ++n)
   {

Modules/Registration/Common/include/itkLandmarkBasedTransformInitializer.hxx

@@ -120,8 +120,8 @@ LandmarkBasedTransformInitializer<TTransform, TFixedImage, TMovingImage>::Intern
   }
 
   // Set a pointSet from the input landmarks.
-  auto pointSet = PointSetType::New();
-  pointSet->Initialize();
+  autopointSet = PointSetType::New();
+
 
   PointsContainerConstIterator fixedIt = this->m_FixedLandmarks.begin();
   PointsContainerConstIterator movingIt = this->m_MovingLandmarks.begin();

Modules/Registration/Metricsv4/include/itkLabeledPointSetToPointSetMetricv4.hxx

@@ -129,8 +129,8 @@ typename LabeledPointSetToPointSetMetricv4<TFixedPointSet, TMovingPointSet, TInt
   LabeledPointSetToPointSetMetricv4<TFixedPointSet, TMovingPointSet, TInternalComputationValueType>::
     GetLabeledFixedPointSet(const LabelType label) const
 {
-  auto fixedPointSet = FixedPointSetType::New();
-  fixedPointSet->Initialize();
+  autofixedPointSet = FixedPointSetType::New();
+
 
   typename FixedPointSetType::PointIdentifier count{};
 
@@ -155,8 +155,8 @@ typename LabeledPointSetToPointSetMetricv4<TFixedPointSet, TMovingPointSet, TInt
   LabeledPointSetToPointSetMetricv4<TFixedPointSet, TMovingPointSet, TInternalComputationValueType>::
     GetLabeledMovingPointSet(const LabelType label) const
 {
-  auto movingPointSet = MovingPointSetType::New();
-  movingPointSet->Initialize();
+  automovingPointSet = MovingPointSetType::New();
+
 
   typename MovingPointSetType::PointIdentifier count{};
 

Modules/Registration/RegistrationMethodsv4/include/itkBSplineSyNImageRegistrationMethod.hxx

@@ -278,8 +278,8 @@ typename BSplineSyNImageRegistrationMethod<TFixedImage, TMovingImage, TOutputTra
       }
     }
 
-    auto gradientPointSet = BSplinePointSetType::New();
-    gradientPointSet->Initialize();
+    autogradientPointSet = BSplinePointSetType::New();
+
 
     if (fixedPointSets[0]->GetNumberOfPoints() > 0)
     {

Modules/Registration/RegistrationMethodsv4/include/itkImageRegistrationMethodv4.hxx

@@ -937,8 +937,8 @@ ImageRegistrationMethodv4<TFixedImage, TMovingImage, TTransform, TVirtualImage,
 
   for (SizeValueType n = 0; n < numberOfLocalMetrics; ++n)
   {
-    auto samplePointSet = MetricSamplePointSetType::New();
-    samplePointSet->Initialize();
+    autosamplePointSet = MetricSamplePointSetType::New();
+
 
     using SamplePointType = typename MetricSamplePointSetType::PointType;
 

Modules/Registration/RegistrationMethodsv4/include/itkTimeVaryingBSplineVelocityFieldImageRegistrationMethod.hxx

@@ -149,14 +149,14 @@ TimeVaryingBSplineVelocityFieldImageRegistrationMethod<TFixedImage,
   this->m_OutputTransform->IntegrateVelocityField();
 
   // Keep track of velocityFieldPointSet from the previous iteration
-  VelocityFieldPointSetPointer velocityFieldPointSetFromPreviousIteration = VelocityFieldPointSetType::New();
-  velocityFieldPointSetFromPreviousIteration->Initialize();
+  VelocityFieldPointSetPointervelocityFieldPointSetFromPreviousIteration = VelocityFieldPointSetType::New();
 
-  VelocityFieldPointSetPointer velocityFieldPointSet = VelocityFieldPointSetType::New();
-  velocityFieldPointSet->Initialize();
 
-  auto velocityFieldWeights = WeightsContainerType::New();
-  velocityFieldWeights->Initialize();
+  VelocityFieldPointSetPointervelocityFieldPointSet = VelocityFieldPointSetType::New();
+
+
+  autovelocityFieldWeights = WeightsContainerType::New();
+
 
   // Monitor the convergence
   using ConvergenceMonitoringType = itk::Function::WindowConvergenceMonitoringFunction<RealType>;

Modules/Segmentation/Classifiers/include/itkBayesianClassifierInitializationImageFilter.hxx

@@ -181,7 +181,6 @@ BayesianClassifierInitializationImageFilter<TInputImage, TProbabilityPrecisionTy
   covarianceEstimatorsContainer->Reserve(m_NumberOfClasses);
 
   m_MembershipFunctionContainer = MembershipFunctionContainerType::New();
-  m_MembershipFunctionContainer->Initialize(); // Clear elements
   for (unsigned int i = 0; i < m_NumberOfClasses; ++i)
   {
     meanEstimatorsContainer->InsertElement(i, new typename GaussianMembershipFunctionType::MeanVectorType(1));