Memory leak in ITKSkullStrip

[AngryCoconut2]

Here's my code.

constexpr unsigned int Dimension = 3;
using PixelType = short;

using ImageType = itk::Image<PixelType, Dimension>;
using AtlasType = itk::Image<float, Dimension>;
using VTKImageToAtlasType = itk::VTKImageToImageFilter<AtlasType>;
using VTKImageToImageType = itk::VTKImageToImageFilter<ImageType>;
typedef itk::Image<short, 3> AtlasImageType;
typedef itk::Image<unsigned char, 3> AtlasLabelType;

int main(){
        string atlas_ori_path = "./atlas.nii";
    string atlas_mask_path = "./mask.nii";
    string image_path = "./ori.nii";

    using ImageIOTypeNII = itk::NiftiImageIO;
    using ReaderTypeNII= itk::ImageFileReader< ImageType  >;
    using ReaderTypeNII_atlas = itk::ImageFileReader< AtlasType  >;

    ReaderTypeNII::Pointer imgReader = ReaderTypeNII::New();
    imgReader->SetFileName(image_path);
    ImageIOTypeNII::Pointer oriNiftiImageIORead = ImageIOTypeNII::New();
    imgReader->SetImageIO(oriNiftiImageIORead);
    imgReader->Update();

    ReaderTypeNII_atlas::Pointer atlasReader = ReaderTypeNII_atlas::New();
    atlasReader->SetFileName(atlas_ori_path);
    ImageIOTypeNII::Pointer atlasNiftiImageIORead = ImageIOTypeNII::New();
    atlasReader->SetImageIO(atlasNiftiImageIORead);
    atlasReader->Update();

    ReaderTypeNII_atlas::Pointer maskReader = ReaderTypeNII_atlas::New();
    maskReader->SetFileName(atlas_mask_path);
    ImageIOTypeNII::Pointer maskNiftiImageIORead = ImageIOTypeNII::New();
    maskReader->SetImageIO(maskNiftiImageIORead);
    maskReader->Update();

    AtlasType::Pointer atlas_itk = AtlasType::New();
    AtlasType::Pointer atlas_mask_itk = AtlasType::New();
    ImageType::Pointer image_itk = ImageType::New();
    ImageType::Pointer brain_itk = ImageType::New();
    image_itk = imgReader->GetOutput();
    atlas_itk = atlasReader->GetOutput();
    atlas_mask_itk = maskReader->GetOutput();

        cout << "Start!" << endl;

    {
        using StripTsFilterType = itk::StripTsImageFilter<ImageType, AtlasType, AtlasType>;
        StripTsFilterType::Pointer stripTsFilter = StripTsFilterType::New();
        stripTsFilter->SetInput(image_itk);
        stripTsFilter->SetAtlasImage(atlas_itk);
        stripTsFilter->SetAtlasBrainMask(maskReader->GetOutput());
        stripTsFilter->ReleaseDataFlagOn();
        ImageType::SizeType imagesize = atlas_itk->GetLargestPossibleRegion().GetSize();
        ImageType::SizeType labelsize = atlas_mask_itk->GetLargestPossibleRegion().GetSize();
        //stripTsFilter->Get
        try
        {
            stripTsFilter->Update();
        }
        catch (itk::ExceptionObject &exception)
        {
            std::cerr << "Exception caught ! " << std::endl;
            std::cerr << exception << std::endl;
        }
        //stripTsFilter->ReleaseDataFlagOn();
    }
    cout << "Done!" << endl;
}

---

All contents in braces shall be released after the end of braces.
By that, It means I should get Exactly using memory by "Done!" as same as "Start!".
However there's 500MB memory leaked.

I ran a few cases to test, hope this would help:

• If I close the program, all leaked memory will be released.
• If I put those code into a function(see below), every time I call "func()" there will be 500MB memory leaked.
  => 2 call = 1000MB
  => 3 call = 1500MB and so on.

  void func(ImageType image_itk,AtlasType atlas_itk,AtlasType atlas_mask_itk)
  {
      using StripTsFilterType = itk::StripTsImageFilter<ImageType, AtlasType, AtlasType>;
      StripTsFilterType::Pointer stripTsFilter = StripTsFilterType::New();
      stripTsFilter->SetInput(image_itk);
      stripTsFilter->SetAtlasImage(atlas_itk);
      stripTsFilter->SetAtlasBrainMask(atlas_mask_itk);
      stripTsFilter->ReleaseDataFlagOn();
      try
      {
          stripTsFilter->Update();
      }
      catch (itk::ExceptionObject &exception)
      {
          std::cerr << "Exception caught ! " << std::endl;
          std::cerr << exception << std::endl;
      }
      //stripTsFilter->ReleaseDataFlagOn();
  }

  Is there anything wrong with the code I wrote?

[AngryCoconut2]

What's more, I'm using ITK-5.2.1