Calculate volume of largest and minimum ellipse around a segmentation

ral_san · February 14, 2020, 10:21am

I think I managed to solve this issue.

gist.github.com

https://gist.github.com/raluca-san/5d8bd4dcb3278b8e2f6dfa959614f9c8

MinVolEllipse.py

# -*- coding: utf-8 -*-
"""
@author: Raluca Sandu
"""

from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import numpy as np
import numpy.linalg as la

This file has been truncated. show original

What is missing is extracting the surface points from a segmentation mask. That could be extracted using the SimpleITK.LabelContour(img, fullyConnected=False).
Also the spacing information is missing, so the ellipsoid should be adjusted by placing it on the correct grid.
This can be achieved scikit-image.draw.ellipsoid https://scikit-image.org/docs/dev/api/skimage.draw.html#skimage.draw.ellipsoid.

Is it possible to propose adding this functionality to the official SimpleITK library?