Definition of n_annuli
@@ -794,17 +794,17 @@ def noise_per_annulus(array, separation, fwhm, init_rad=None, wedge=(0, 360),
raise TypeError('Wedge must be a tuple with the initial and final '
'angles')
+ if init_rad is None:
+ init_rad = fwhm
+
init_angle, fin_angle = wedge
centery, centerx = frame_center(array)
- n_annuli = int(np.floor((centery)/separation)) - 1
+ n_annuli = int(np.floor((centery - init_rad)/separation))
noise = []
vector_radd = []
if verbose:
print('{} annuli'.format(n_annuli))
- if init_rad is None:
- init_rad = fwhm
-
if debug:
_, ax = plt.subplots(figsize=(6, 6))
ax.imshow(array, origin='lower', interpolation='nearest',
