def readMar(f): from fabio import marccdimage as mccd print "reading image",f i=mccd.marccdimage() i.read(f) # conversion to float32 is fastest ... (but then in bincount is converted to float anyway) #return i.data return i.data.astype(np.float)
def checkFile(self, filename):
""" check we can read it"""
obj = openimage(filename)
obj2 = marccdimage()
obj2.read(filename)
self.assertEqual(obj.data[10, 10], obj2.data[10, 10])
self.assertEqual(type(obj), type(obj2))
self.assertEqual(abs(obj.data.astype(int) - obj2.data.astype(int)).sum(), 0)
def test_read(self):
""" check we can read MarCCD images"""
for line in TESTIMAGES.split("\n"):
vals = line.split()
name = vals[0]
dim1, dim2 = [int(x) for x in vals[1:3]]
mini, maxi, mean, stddev = [float(x) for x in vals[3:]]
obj = marccdimage()
obj.read(self.fn[name])
self.assertAlmostEqual(mini, obj.getmin(), 2, "getmin")
self.assertAlmostEqual(maxi, obj.getmax(), 2, "getmax")
self.assertAlmostEqual(mean, obj.getmean(), 2, "getmean")
self.assertAlmostEqual(stddev, obj.getstddev(), 2, "getstddev")
self.assertEqual(dim1, obj.dim1, "dim1")
self.assertEqual(dim2, obj.dim2, "dim2")
def test_read(self):
""" check we can read MarCCD images"""
for line in TESTIMAGES.split("\n"):
vals = line.split()
name = vals[0]
dim1, dim2 = [int(x) for x in vals[1:3]]
mini, maxi, mean, stddev = [float(x) for x in vals[3:]]
obj = marccdimage()
obj.read(self.fn[name])
self.assertAlmostEqual(mini, obj.getmin(), 2, "getmin")
self.assertAlmostEqual(maxi, obj.getmax(), 2, "getmax")
self.assertAlmostEqual(mean, obj.getmean(), 2, "getmean")
self.assertAlmostEqual(stddev, obj.getstddev(), 2, "getstddev")
self.assertEqual(dim1, obj.dim1, "dim1")
self.assertEqual(dim2, obj.dim2, "dim2")
