【SimpleITK教程】SimpleITK获取医学图像ROI的外接矩形框bounding box

代码

image = sitk.ReadImage('xxx.nrrd')
roi = sitk.ReadImage('xxx.nrrd')
filter = sitk.LabelStatisticsImageFilter()
filter.Execute(image, roi)
bounding_box = filter.GetBoundingBox(1) ### 1指的是要获取boundingbox的label值
# bounding_box为六元数组
# 数值为 (sagittal_min, sagittal_max, coronal_min, coronal_max, axial_min, axial_max)

Tips

  • 考虑到Execute函数必须传递image,而bounding_box又只与roi相关,因此偷懒的写法可以是filter.Execute(roi, roi)

版权声明:本文为JianJuly原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接和本声明。