script/devices/eiger.py

@@ -0,0 +1,76 @@
+from ijutils import *
+from ch.psi.pshell.imaging.Overlays import *
+from ch.psi.pshell.imaging.Utils import *
+import ch.psi.pshell.imaging.Pen as Pen
+import java.awt.Rectangle as Rectangle
+import ch.psi.pshell.imaging.Data as Data
+
+
+###############################################################################
+# ROI Integration
+###############################################################################
+
+def integrate_roi(source, x,y, w, h):
+    if source.data is None:
+        source.update()
+    roi = source.data.getRoi(Rectangle(x,y, w, h))
+    return roi.integrate(False)
+
+
+
+class RoiIntensity(ReadonlyRegisterBase):
+    def __init__(self, name, source, x,y, w, h):
+        ReadonlyRegisterBase.__init__(self, name)
+        self.source=source
+        self.roi = x,y, w, h
+        
+    def doRead(self):
+        x,y, w, h = self.roi
+        return integrate_roi(self.source, x,y, w, h)
+
+###############################################################################
+# Frame integration
+###############################################################################
+
+def get_image(source, roi=None, wait_next=False):
+    if wait_next:
+        source.waitNext(-1)
+    ret = load_image(Utils.grayscale(source.output, Rectangle(roi[0], roi[1], roi[2], roi[3]) if (roi is not None) else None))
+    return ret
+    
+def grab_frames(source, samples, roi=None, wait_next=False, sleep=0):
+    frames = []            
+    for i in range(samples):
+        if sleep>0:
+            time.sleep(sleep)
+        aux = get_image(source, roi, wait_next)        
+        frames.append(aux)        
+    return frames
+    
+def average_frames(source, samples=1, roi=None, wait_next=False, sleep=0, as_float=True):
+    return average_ips(grab_frames(source, samples, roi, wait_next, sleep), as_float)
+
+def integrate_frames(source, samples=1, roi=None, wait_next=False, sleep=0, as_float=True):
+    return integrate_ips(grab_frames(source, samples, roi, wait_next, sleep), as_float)
+       
+
+
+print integrate_roi(image, 10, 5, 20, 10)
+        
+add_device(RoiIntensity("Region1", image, 10, 5, 20, 10), True)
+add_device(RoiIntensity("Region2", image, 10, 5, 40, 20), True)
+
+
+import ch.psi.pshell.data.ProviderCSV as ProviderCSV 
+ProviderCSV.setDefaultItemSeparator(" ")
+tscan((Region1, Region2), 10, 0.1, layout="table", provider = "csv")
+
+ret = grab_frames(image, 10, sleep=0.1)
+av = average_ips(ret, None, True)
+data = get_ip_array(av)
+plot(data)
+save_image(av, "{images}/float.tif","tiff")
+
+av = average_ips(ret, None, False)
+save_image(av, "{images}/int.tif","tiff")
+