[GSAS-II] GSAS-II image import int vs float

Gallington, Leighanne C. Gallington at anl.gov
Wed Mar 16 17:18:24 CDT 2022


Hi Bob and Brian,

We've run into some problems when trying to import and apply gain maps, with values ranging from about 0.97 to 1.03, to Jiahui's data.  After doing some digging, I've found that when the PE images are imported, the values are converted from floats to integers.  As far as I can tell, the conversion is happening in the code snippet below, from file G2img_1TIF.py, in function GetTifData.

 elif sizexy == [2048,2048] or sizexy == [1024,1024] or sizexy == [3072,3072]:
        if IFD[273][2][0] == 8:
            if IFD[258][2][0] == 32:
                tifType = 'PE'
                pixy = [200.,200.]
                File.seek(8)
                G2fil.G2Print ('Read APS PE-detector tiff file: '+filename)
                if dataType == 5:
                    image = np.array(np.frombuffer(File.read(4*Npix),dtype=np.float32),dtype=np.int32)  #fastest
                else:
                    image = np.array(np.frombuffer(File.read(4*Npix),dtype=np.int32),dtype=np.int32)

If I change the np.int32 to np.float32, the gain map can be read in, applied to data, and integrated successfully.  I'm not sure what else I might have broken downstream when I do that, and there's still some funkiness with the slider bars for the intensities.  Can you implement a cleaner fix into the G2img code?

Thanks,
Leighanne
--
Leighanne C. Gallington, Ph.D.
Assistant Chemist
Advanced Photon Source (X-ray Science Division)
Argonne National Laboratory
Argonne, IL  60439
gallington at anl.gov
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.aps.anl.gov/pipermail/gsas-ii/attachments/20220316/97ae1031/attachment.html>


More information about the GSAS-II mailing list