Script execution
This commit is contained in:
@@ -1,10 +1,17 @@
|
||||
import numpy
|
||||
import time
|
||||
|
||||
from scipy import signal, ndimage
|
||||
|
||||
def img_get_int(fname, thres1, thres2, thres3, thres4, header, width, height, depth, x1,y1,x2,y2, bx1,by1,bx2,by2 , filter_median = False, filter_nsigma = 0):
|
||||
# read actual image file
|
||||
img = numpy.fromfile(fname, dtype=numpy.uint32)
|
||||
img.shape = height, width
|
||||
if filter_nsigma>0:
|
||||
img = ndimage.gaussian_filter(img, filter_nsigma)
|
||||
elif filter_median:
|
||||
#img = signal.medfilt2d(img.astype('d'), kernel_size=3)
|
||||
img = ndimage.median_filter(img, size=3)
|
||||
# signal roi
|
||||
area_I = ( x2 - x1 + 1) * ( y2 - y1 + 1)
|
||||
I_sum = img[y1:y2, x1:x2].sum()
|
||||
@@ -12,7 +19,6 @@ def img_get_int(fname, thres1, thres2, thres3, thres4, header, width, height, d
|
||||
thresh2_count = len(numpy.where(img>thres2)[0])
|
||||
thresh3_count = len(numpy.where(img>thres3)[0])
|
||||
thresh4_count = len(numpy.where(img>thres4)[0])
|
||||
time.sleep(1.0)
|
||||
# background roi
|
||||
I_sum_bgr = img[by1:by2, bx1:bx2].sum()
|
||||
area_bgr= (bx2 - bx1 + 1) * (by2 - by1 + 1)
|
||||
|
||||
Reference in New Issue
Block a user