In [ ]:

    

from __future__ import print_function
import matplotlib.pyplot as plt
%matplotlib inline

import SimpleITK as sitk
print(sitk.Version())
from myshow import myshow
# Download data to work on
from downloaddata import fetch_data as fdata


    

In [ ]:

    

img = sitk.ReadImage(fdata("a_vm1108.png"))
myshow(img)


    

In [ ]:

    

seeds = [[400,500], [1500,800]]
seg = sitk.VectorConfidenceConnected(img, seedList=seeds, numberOfIterations=5, multiplier=7)


    

In [ ]:

    

myshow(sitk.Mask(img, seg==0))


    

In [ ]:

    

seg = sitk.BinaryClosingByReconstruction(seg, 100)


    

In [ ]:

    

myshow(sitk.Mask(img, seg==0))