From d89bfb5825e2b080b9e0e683503549729db41c15 Mon Sep 17 00:00:00 2001 From: VChristiaens Date: Mon, 25 Nov 2024 00:11:35 +0100 Subject: [PATCH] PCA-ADI step of 2-stage PCA-SADI is now compatible with a threshold in rotation (delta_rot, source_xy and min_frames_pca now propagated) --- vip_hci/psfsub/pca_fullfr.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/vip_hci/psfsub/pca_fullfr.py b/vip_hci/psfsub/pca_fullfr.py index a39a1c9c..ea0bbb90 100644 --- a/vip_hci/psfsub/pca_fullfr.py +++ b/vip_hci/psfsub/pca_fullfr.py @@ -1289,7 +1289,7 @@ def _adimsdi_doublepca( ann_center = dist(yc, xc, y1, x1) pa_thr = _compute_pa_thresh(ann_center, fwhm, delta_rot) - res_ifs_adi = np.zeros_like(cube) + res_ifs_adi = np.zeros_like(residuals_cube_channels) for frame in range(n): ind = _find_indices_adi(angle_list, frame, pa_thr) @@ -1308,7 +1308,7 @@ def _adimsdi_doublepca( left_eigv=left_eigv, min_frames_pca=min_frames_pca, ) - res_ifs_adi[frame] = res_result[-1] + res_ifs_adi[frame] = res_result[-1].reshape((y_in, x_in)) if verbose: print("De-rotating and combining residuals")