wassname/scikit-image
1
0
Fork 0
mirror of https://github.com/wassname/scikit-image.git synced 2026-07-08 17:36:35 +08:00
Code Issues Packages Projects Releases Wiki Activity

added test_blob_log

Browse Source
This commit is contained in:
Vighnesh Birodkar
2014-03-10 23:33:33 +05:30
parent a067485dec
commit 47c5f1bda6
2 changed files with 67 additions and 31 deletions
Download Patch File Download Diff File Expand all files Collapse all files
skimage/feature/blob.py
+22 -30
View File
@@ -202,11 +202,11 @@ def blob_dog(image, min_sigma=1, max_sigma=50, sigma_ratio=1.6, threshold=2.0,
return []
def blob_log(image, min_sigma=1, max_sigma=50, num_sigma=10, threshold=.1,
def blob_log(image, min_sigma=1, max_sigma=50, num_sigma=10, threshold=.2,
overlap=.5, log_scale=False):
"""Finds blobs in the given grayscale image.
Blobs are found using the Laplacian of Gaussian (DoG) method[1]_.
Blobs are found using the Laplacian of Gaussian (LoG) method[1]_.
For each blob found, its coordinates and area are returned.
Parameters
@@ -249,33 +249,25 @@ def blob_log(image, min_sigma=1, max_sigma=50, num_sigma=10, threshold=.1,
--------
>>> from skimage import data, feature, exposure
>>> img = data.coins()
>>> img = exposure.equalize_hist(img) # imporves detection
>>> img = exposure.equalize_hist(img) # improves detection
>>> feature.blob_log(img,threshold = .3)
array([[ 107, 333, 6],
[ 107, 337, 25],
[ 108, 329, 6],
[ 113, 323, 6],
[ 114, 322, 6],
[ 121, 273, 1608],
[ 124, 336, 904],
[ 125, 45, 1061],
[ 125, 207, 904],
array([[ 113, 323, 6],
[ 121, 272, 1815],
[ 124, 336, 760],
[ 126, 46, 760],
[ 126, 208, 760],
[ 127, 102, 760],
[ 128, 155, 760],
[ 178, 261, 25],
[ 186, 345, 2268],
[ 193, 276, 1413],
[ 194, 213, 1413],
[ 196, 102, 1061],
[ 197, 43, 904],
[ 198, 155, 904],
[ 198, 255, 56],
[ 214, 282, 25],
[ 260, 174, 1608],
[ 262, 244, 1413],
[ 262, 302, 1413],
[ 266, 114, 1061],
[ 268, 358, 1061]])
[ 128, 154, 760],
[ 185, 344, 1815],
[ 194, 213, 1815],
[ 194, 276, 1815],
[ 197, 44, 760],
[ 198, 103, 760],
[ 198, 155, 760],
[ 260, 174, 1815],
[ 263, 244, 1815],
[ 263, 302, 1815],
[ 266, 115, 760]])
"""
@@ -285,10 +277,10 @@ def blob_log(image, min_sigma=1, max_sigma=50, num_sigma=10, threshold=.1,
image = img_as_float(image)
if log_scale:
sigma_list = np.linspace(min_sigma, max_sigma, num_sigma)
else:
start, stop = log(min_sigma, 10), log(max_sigma, 10)
sigma_list = np.logspace(start, stop)
sigma_list = np.logspace(start, stop, num_sigma)
else:
sigma_list = np.linspace(min_sigma, max_sigma, num_sigma)
gl_images = [-gaussian_laplace(image, s) * s ** 2 for s in sigma_list]
image_cube = np.dstack(gl_images)
skimage/feature/tests/test_blob.py
+45 -1
View File
@@ -1,6 +1,6 @@
import numpy as np
from skimage.draw import circle
from skimage.feature import blob_dog
from skimage.feature import blob_dog, blob_log
import math
@@ -36,3 +36,47 @@ def test_blob_dog():
assert abs(b[0] - 200) <= thresh
assert abs(b[1] - 350) <= thresh
assert abs(radius(b[2]) - 45) <= thresh
def test_blob_log():
img = np.ones((512, 512))
xs, ys = circle(400, 130, 5)
img[xs, ys] = 255
xs, ys = circle(160, 50, 15)
img[xs, ys] = 255
xs, ys = circle(100, 300, 25)
img[xs, ys] = 255
xs, ys = circle(200, 350, 30)
img[xs, ys] = 255
blobs = blob_log(img, min_sigma=5, max_sigma=20, threshold=1)
area = lambda x: x[2]
radius = lambda x: math.sqrt(x / math.pi)
s = sorted(blobs, key=area)
thresh = 3
b = s[0]
assert abs(b[0] - 400) <= thresh
assert abs(b[1] - 130) <= thresh
assert abs(radius(b[2]) - 5) <= thresh
b = s[1]
assert abs(b[0] - 160) <= thresh
assert abs(b[1] - 50) <= thresh
assert abs(radius(b[2]) - 15) <= thresh
b = s[2]
assert abs(b[0] - 100) <= thresh
assert abs(b[1] - 300) <= thresh
assert abs(radius(b[2]) - 25) <= thresh
b = s[3]
assert abs(b[0] - 200) <= thresh
assert abs(b[1] - 350) <= thresh
assert abs(radius(b[2]) - 30) <= thresh