itk中的花式数据切割（五）

最新推荐文章于 2023-11-09 16:51:49 发布

winston600

最新推荐文章于 2023-11-09 16:51:49 发布

阅读量1.5k

点赞数

CC 4.0 BY-SA版权

分类专栏： ITK医学图像处理

本文链接：https://blog.youkuaiyun.com/rabbitbride/article/details/72844167

ITK医学图像处理专栏收录该内容

94 篇文章

订阅专栏

本文介绍了ITK库中itkCropImageFilter类的基本使用方法，并通过示例代码展示了如何进行图像裁剪。该方法通过对原始图像进行对称裁剪，实现指定区域的图像提取。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

今天继续研究数据切割，也是对之前四篇切割内容的补充，内容不多，step by step。

itk中有一个类叫做itkCropImageFilter，是itkExtractImageFilter的派生方法。

先上代码，看一下基本用法：

ImageType::SizeType   extractSize = {{8, 12}};  
extractSize[0] = 1; 
extractSize[1] = 1;

typedef itk::CropImageFilter <ImageType, ImageType> CropImageFilterType;
CropImageFilterType::Pointer cropFilter  = CropImageFilterType::New();
cropFilter->SetInput(input_data);
cropFilter->SetBoundaryCropSize(extractSize); //和下面两句重复
//cropFilter->SetUpperBoundaryCropSize(extractSize);
//cropFilter->SetLowerBoundaryCropSize(extractSize);
cropFilter->UpdateLargestPossibleRegion();
cropFilter->Update();

接下来咱们一起看源码是怎么写的，这个类超简单，就一个函数，但我看了还是诧异，+_+：

template <class TInputImage, class TOutputImage>// itk的模版
void
CropImageFilter<TInputImage, TOutputImage>
::GenerateOutputInformation()//itk的数据生成都用的这个函数名
{
  const TInputImage * inputPtr = this->GetInput();
  if( !inputPtr )//防爆安全代码
    {
    return;
    }
  
  // Compute the new region size.
  OutputImageRegionType croppedRegion;
  SizeType   sz;
  OutputImageIndexType   idx;

  InputImageSizeType input_sz =
    inputPtr->GetLargestPossibleRegion().GetSize();
  InputImageIndexType input_idx =
    inputPtr->GetLargestPossibleRegion().GetIndex();

  idx = input_idx + m_LowerBoundaryCropSize; 
  sz  = input_sz  - (m_UpperBoundaryCropSize + m_LowerBoundaryCropSize); 

  croppedRegion.SetSize(sz);
  croppedRegion.SetIndex(idx);

  // Set extraction region in the superclass.
  this->SetExtractionRegion(croppedRegion);//父类的函数生效

  // 
  Superclass::GenerateOutputInformation();
}

看下面两行！@_@，我有点晕了，亲测后，原来是按照给定的size对原始图像做对称切割！
idx = input_idx + m_LowerBoundaryCropSize;
sz = input_sz - (m_UpperBoundaryCropSize + m_LowerBoundaryCropSize);
测试结果：（貌似这是第一个与文章内容相关的图）