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

Changed example to make code compact.

Browse Source
This commit is contained in:
Vighnesh Birodkar
2014-04-27 11:52:00 +05:30
parent 7738e65de8
commit 336f117a01
1 changed files with 22 additions and 34 deletions
Download Patch File Download Diff File Expand all files Collapse all files
doc/examples/plot_blob.py
+22 -34
View File
@@ -44,47 +44,35 @@ from skimage.color import rgb2gray
image = data.hubble_deep_field()[0:500, 0:500]
image_gray = rgb2gray(image)
blobs = blob_log(image_gray, max_sigma=30, num_sigma=10, threshold=.1)
blobs_log = blob_log(image_gray, max_sigma=30, num_sigma=10, threshold=.1)
# Compute areas in the 2nd column.
blobs_log[:, 2] = blobs_log[:, 2] * sqrt(2)
fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
ax.set_title('Laplacian of Gaussian')
plt.imshow(image)
blobs_dog = blob_dog(image_gray, max_sigma=30, threshold=.1)
blobs_dog[:, 2] = blobs_dog[:, 2] * sqrt(2)
for blob in blobs:
y, x = blob[0], blob[1]
r = blob[2] * sqrt(2)
c = plt.Circle((x, y), r, color='#00ff00', lw=2, fill=False)
ax.add_patch(c)
blobs_doh = blob_doh(image_gray, max_sigma=30, threshold=.01)
fig1, ax1 = plt.subplots(1, 1)
ax1.set_title('Laplacian of Gaussian')
blobs = blob_dog(image_gray, max_sigma=30, threshold=.1)
fig2, ax2 = plt.subplots(1, 1)
ax2.set_title('Difference of Gaussian')
fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
ax.set_title('Difference of Gaussian')
plt.imshow(image)
fig3, ax3 = plt.subplots(1, 1)
ax3.set_title('Determinant of Hessian')
for blob in blobs:
y, x = blob[0], blob[1]
r = blob[2] * sqrt(2)
c = plt.Circle((x, y), r, color='#ffff00', lw=2, fill=False)
ax.add_patch(c)
axes = [ax1, ax2, ax3]
blobs_list = [blobs_log, blobs_dog, blobs_doh]
colors = ['yellow', 'lime', 'red']
sequence = zip(axes, blobs_list, colors)
blobs = blob_doh(image_gray, max_sigma=30, threshold=.01)
fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
ax.set_title('Determinant of Hessian')
plt.imshow(image)
for blob in blobs:
y, x = blob[0], blob[1]
r = blob[2]
c = plt.Circle((x, y), r, color='#ff0000', lw=2, fill=False)
ax.add_patch(c)
for ax, blobs, color in sequence:
ax.imshow(image, interpolation='nearest')
for blob in blobs:
y, x, r = blob
c = plt.Circle((x, y), r, color=color, linewidth=2, fill=False)
ax.add_patch(c)
plt.show()