s = 10
normDensity  = local_normalize(density, size=s)
densityTarget = invert(normDensity)**2

vesselnessImage = apply_frangi_filter(densityTarget)
tVesselnessImage = np.ones_like(vesselnessImage)
tVesselnessImage[vesselnessImage < 0.03] = 0

fig, axs = plt.subplots(1,2,figsize=(20, 7))
densityCMAP = axs[0].imshow(vesselnessImage, origin='lower', aspect='auto', extent=[xMin, xMax, yMin, yMax])
axs[1].imshow(tVesselnessImage, origin='lower', aspect='auto', extent=[xMin, xMax, yMin, yMax])
plt.colorbar(densityCMAP)