+ // lungs' barycenter
+ typename LabelType::CentroidType lung_centroid = label->GetCentroid();
+
+ // try to guess ideal ellipse axes from the lungs' bounding box and centroid
+ // with some hard-coded "magic" constants...
+#if ITK_VERSION_MAJOR >= 4
+ typename LabelType::RegionType lung_bbox = label->GetBoundingBox();
+#else