fixed HURA alignment so that center is now true

bpops · Aug 2, 2023 · aac0371 · aac0371
1 parent b88b2f0
commit aac0371
Show file tree

Hide file tree

Showing 2 changed files with 55 additions and 16 deletions.
diff --git a/codedapertures/codedapertures.py b/codedapertures/codedapertures.py
@@ -458,7 +458,9 @@ def __init__(self, rank=4, r=5, radius=5, quiet=False):
         self.l = np.zeros((self.rx,self.ry), dtype=np.int16)
         for i in range(self.rx):
             for j in range(self.ry):
-                self.l[i,j] = (i + r*j) % self.v
+                i_idx = i - self.radius
+                j_idx = j - self.radius
+                self.l[i,j] = (i_idx + self.r*j_idx) % self.v
 
         # calculate mask
         self.mask = np.zeros(self.l.shape, dtype=np.int16)
@@ -656,7 +658,9 @@ def __init__(self, rank=4, radius=5, quiet=False):
         self.l = np.zeros((self.rx,self.ry), dtype=np.int16)
         for i in range(self.rx):
             for j in range(self.ry):
-                self.l[i,j] = (i + self.r*j) % self.v
+                i_idx = i - self.radius
+                j_idx = j - self.radius
+                self.l[i,j] = (i_idx + self.r*j_idx) % self.v
 
         # calculate mask
         self.mask = np.zeros(self.l.shape, dtype=np.int16)