Minor bug fix. obj doesn't necessarily have attribute original.

Author: rmjarvis

imsim/lsst_image.py

@@ -1,4 +1,5 @@
 import numpy as np
+import hashlib
 import logging
 import galsim
 from galsim.config import RegisterImageType, GetAllParams, GetSky, AddNoise
@@ -293,11 +294,15 @@ def addNoise(self, image, config, base, image_num, obj_num, current_var, logger)
             camera = get_camera(self.camera_name)
             det_name = base['det_name']
             serial_number = camera[det_name].getSerial()
+            # Note: the regular Python hash function is non-derterministic, which is not good.
+            # Instead we use hashlib.sha256, which is deterministic and convert that to an integer.
+            # https://stackoverflow.com/questions/27954892/deterministic-hashing-in-python-3
+            seed = int(hashlib.sha256(serial_number.encode('UTF-8')).hexdigest(), 16) & 0xFFFFFFFF
             # Only apply fringing to e2v sensors.
             if serial_number[:3] == 'E2V':
                 ccd_fringing = CCD_Fringing(true_center=image.wcs.toWorld(image.true_center),
                                             boresight=self.boresight,
-                                            seed=hash(serial_number), spatial_vary=True)
+                                            seed=seed, spatial_vary=True)
                 ny, nx = sky.array.shape
                 xarr, yarr = np.meshgrid(range(nx), range(ny))
                 logger.info("Apply fringing")

imsim/stamp.py

@@ -293,7 +293,7 @@ def _getGoodPhotImageSize1(self, obj, keep_sb_level, pixel_scale):
                 else:
                     obj_list.append(item)
             obj = galsim.Add(obj_list)
-        elif isinstance(obj.original, galsim.RandomKnots):
+        elif hasattr(obj, 'original') and isinstance(obj.original, galsim.RandomKnots):
             # Handle RandomKnots object directly
             obj = obj.original._profile
 

tests/test_fringing.py

@@ -1,4 +1,5 @@
 import numpy as np
+import hashlib
 import pytest
 from imsim import make_batoid_wcs, CCD_Fringing, get_camera
 import galsim
@@ -18,6 +19,7 @@ def test_fringing():
     camera = get_camera()
     det_name = 'R22_S11'
     serial_num = camera[det_name].getSerial()
+    seed = int(hashlib.sha256(serial_num.encode('UTF-8')).hexdigest(), 16) & 0xFFFFFFFF
 
     xarr, yarr = np.meshgrid(range(4096), range(4004))
 
@@ -41,7 +43,7 @@ def test_fringing():
 
     ccd_fringing = CCD_Fringing(true_center=image.wcs.toWorld(image.true_center),
                                 boresight=world_center,
-                                seed=hash(serial_num), spatial_vary=True)
+                                seed=seed, spatial_vary=True)
     # Test zero value error
     with pytest.raises(ValueError):
         ccd_fringing.calculate_fringe_amplitude(xarr, yarr, amplitude=0)
@@ -91,7 +93,7 @@ def test_fringing():
 
     ccd_fringing_1 = CCD_Fringing(true_center=image.wcs.toWorld(image.true_center),
                                 boresight=world_center,
-                                seed=hash(serial_num), spatial_vary=False)
+                                seed=seed, spatial_vary=False)
     fringe_map1 = ccd_fringing_1.calculate_fringe_amplitude(xarr,yarr)
 
     # Try another random location on the focal plane.
@@ -101,7 +103,7 @@ def test_fringing():
 
     ccd_fringing_2 = CCD_Fringing(true_center=image.wcs.toWorld(image.true_center),
                                 boresight=world_center,
-                                seed=hash(serial_num), spatial_vary=False)
+                                seed=seed, spatial_vary=False)
     fringe_map2 = ccd_fringing_2.calculate_fringe_amplitude(xarr,yarr)
     # Check if the two fringing maps are indentical.
     if np.array_equal(fringe_map1,fringe_map2) != True: