// try to guess ideal ellipse axes from the lungs' bounding box
+#if ITK_VERSION_MAJOR >= 4
+ typename LabelType::RegionType lung_bbox = label->GetBoundingBox();
+#else
typename LabelType::RegionType lung_bbox = label->GetRegion();
+#endif
axes[0] = 0.95*lung_bbox.GetSize()[0]*spacing[0]/2;
axes[1] = 0.3*lung_bbox.GetSize()[1]*spacing[1]/2;
axes[2] = 0.95*lung_bbox.GetSize()[2]*spacing[2]/2;