ITK 3D图像逐片做Mask再合成3D图像

最新推荐文章于 2022-12-09 12:08:24 发布

sharemon

最新推荐文章于 2022-12-09 12:08:24 发布

阅读量1.3k

点赞数 1

本文链接：https://blog.csdn.net/sharemon/article/details/98744864

版权

itk中做mask的方法如下：

using MaskFilterType = itk::MaskImageFilter< InputImageType, MaskImageType, OutputImageType >;
MaskFilterType::Pointer maskFilter = MaskFilterType::New();
maskFilter->SetInput(reader->GetOutput());
maskFilter->SetMaskImage(mask);
maskFilter->update();

对于3D图像切片再还原，有两种思路：

用SliceBySliceImageFilter
先用ExtractImageFilter切片，然后做Mask，再用JoinImageFilter合并

第一种方法可以参考：itk :: SliceBySliceImageFilter
我在用这个filter的时候，出来的结果一直不对。估计是没有处理好多输入的问题。
第二种方法成功：

for (int i = 0; i < size; i++)
{
	...
	using MaskFilterType = itk::MaskImageFilter< ImageTypeSlice, ImageTypeSlice, ImageTypeSlice >;
	MaskFilterType::Pointer maskFilter = MaskFilterType::New();
	maskFilter->SetInput(imageSlice);
	maskFilter->SetMaskImage(maskSlice);
	maskFilter->Update();

	joinFilter->SetInput(i, maskFilter->GetOutput());
	joinFilter->Update();
}
	
// 最后set一下spacing 和 origin
joinFilter->SetSpacing(0.625);
joinFilter->Update();