Project

General

Profile

Actions

Feature #4705

open

Lacey carbon DoG-based hole finder

Added by Neil Voss almost 8 years ago. Updated over 7 years ago.

Status:
In Test
Priority:
Normal
Assignee:
Category:
-
Target version:
-
Start date:
01/13/2017
Due date:
% Done:

0%

Estimated time:
Deliverable:

Description

I wanted to design a lacey carbon DoG-based hole finder.

basically the workflow would be:

users sets a pixel diameter (e.g. 75 pixels) and sloppiness factor (k = 1.3)

then you'll get two sigma values:

k = 1.3
pixdiam = 75

# find xi (E) function of k
Ek = math.sqrt( (k**2 - 1.0) / (2.0 * k**2 * math.log(k)) )

# convert pixel diameter to sigma1
sigma1 = Ek * pixdiam/2

# find sigmaprime
sigmaprime = sigma1 * math.sqrt(k*k-1.0)

#do the blurring
#first blur
imgarray1 = ndimage.gaussian_filter(imgarray0, sigma=sigma1, mode='reflect')
# double blur of first blur
imgarray2 = ndimage.gaussian_filter(imgarray1, sigma=sigmaprime, mode='reflect')

#subtract
dogmap = imgarray1 - imgarray2

threshold and away you go...

Actions

Also available in: Atom PDF