In [1]:
from __future__ import print_function
%matplotlib inline
import matplotlib.pyplot as plt

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

SimpleITK Version: 1.2.4 (ITK 4.13)
Compiled: Nov 12 2019 23:05:07
In [2]:
img = sitk.ReadImage(fdata("a_vm1108.png"))
myshow(img)

Fetching a_vm1108.png
In [3]:
seeds = [[400,500], [1500,800]]
seg = sitk.VectorConfidenceConnected(img, seedList=seeds, numberOfIterations=5, multiplier=7)
In [4]:
myshow(sitk.Mask(img, seg==0))
In [5]:
seg = sitk.BinaryClosingByReconstruction(seg, 100)
In [6]:
myshow(sitk.Mask(img, seg==0))