113 lines
3.5 KiB
Python
113 lines
3.5 KiB
Python
###################################################################################################
|
|
# Example of using ImageJ functionalities through ijutils.
|
|
###################################################################################################
|
|
|
|
from ijutils import *
|
|
import java.awt.Color as Color
|
|
|
|
import ch.psi.pshell.imaging.Filter as Filter
|
|
from ch.psi.pshell.imaging.Overlays import *
|
|
import ch.psi.pshell.imaging.Pen as Pen
|
|
|
|
|
|
"""
|
|
def detect_led(ip):
|
|
aux = grayscale(ip, in_place=False)
|
|
threshold(aux,0,50)
|
|
binary_fill_holes(aux)
|
|
return analyse_particles(aux, 10000,50000,
|
|
fill_holes = False, exclude_edges = True,print_table=True,
|
|
output_image = "outlines", minCirc = 0.4, maxCirc = 1.0)
|
|
"""
|
|
|
|
roi_center = (2322, 1025)
|
|
roi_radius = 1000
|
|
|
|
|
|
|
|
"""
|
|
def in_roi(x,y):
|
|
global roi_center, roi_radius
|
|
return math.hypot(x-roi_center[0], y-roi_center[1]) < roi_radius
|
|
|
|
def detect_led(ip):
|
|
aux = grayscale(ip, in_place=False)
|
|
invert(aux)
|
|
#subtract_background(aux)
|
|
auto_threshold(aux)
|
|
binary_open(aux)
|
|
(results,output) = analyse_particles(aux, 800,4000,
|
|
fill_holes = True, exclude_edges = True,print_table=False,
|
|
output_image = "outlines", minCirc = 0.4
|
|
, maxCirc = 1.0)
|
|
r=results
|
|
print "\n"
|
|
print r.getColumnHeadings()
|
|
for row in range (r.counter):
|
|
if in_roi(r.getValue("XM",row), r.getValue("YM",row)):
|
|
print r.getRowAsString(row)
|
|
else:
|
|
#print >> sys.stderr, 'Invalid location:' + str(r.getValue("XM", row)) + ", " + str( r.getValue("YM", row))
|
|
print >> sys.stderr, r.getRowAsString(row)
|
|
return (results,output)
|
|
|
|
class MyFilter(Filter):
|
|
def process(self, image, data):
|
|
ip = load_image(image)
|
|
(results,output) = detect_led(ip)
|
|
invert(output)
|
|
op_image(ip, output, "xor")
|
|
return ip.getBufferedImage()
|
|
|
|
"""
|
|
|
|
|
|
|
|
def in_roi(x,y):
|
|
return math.hypot(x-roi_radius, y-roi_radius) < roi_radius
|
|
|
|
|
|
def detect_led(ip):
|
|
global roi_center, roi_radius
|
|
aux = sub_image(ip, roi_center[0] - roi_radius, roi_center[1] - roi_radius, 2* roi_radius, 2*roi_radius)
|
|
grayscale(aux)
|
|
invert(aux)
|
|
#subtract_background(aux)
|
|
auto_threshold(aux)
|
|
binary_open(aux)
|
|
(results,output) = analyse_particles(aux, 800,4000,
|
|
fill_holes = True, exclude_edges = True,print_table=False,
|
|
output_image = "outlines", minCirc = 0.4
|
|
, maxCirc = 1.0)
|
|
r=results
|
|
|
|
print "\n"
|
|
print r.getColumnHeadings()
|
|
for row in range (r.counter):
|
|
if in_roi(r.getValue("XM",row), r.getValue("YM",row)):
|
|
print r.getRowAsString(row)
|
|
else:
|
|
#print >> sys.stderr, 'Invalid location:' + str(r.getValue("XM", row)) + ", " + str( r.getValue("YM", row))
|
|
print >> sys.stderr, r.getRowAsString(row)
|
|
|
|
return (results,output)
|
|
|
|
|
|
class MyFilter(Filter):
|
|
def process(self, image, data):
|
|
global roi_center, roi_radius
|
|
ip = load_image(image)
|
|
(results,output) = detect_led(ip)
|
|
invert(output)
|
|
#set_lut(output, outline_lut1[0], outline_lut1[1], outline_lut1[2])
|
|
#op_image(ip, output, "and")
|
|
output = pad_image(output, roi_center[0] - roi_radius, 0, roi_center[1] - roi_radius, 0)
|
|
|
|
op_image(ip, output, "xor")
|
|
return ip.getBufferedImage()
|
|
#return output.getBufferedImage()
|
|
|
|
#Setting the filter to a source
|
|
img.setFilter(MyFilter())
|
|
|