Commit f71d73b9 authored by Jay's avatar Jay
Browse files

Logging residuals

parent 03ab8634
Loading
Loading
Loading
Loading
+1 −0
Original line number Original line Diff line number Diff line
@@ -111,6 +111,7 @@ def estimate_affine_from_sensors(reference_image,
    residuals = affine.residuals(np.array(base_gcps), np.array(dst_gcps))
    residuals = affine.residuals(np.array(base_gcps), np.array(dst_gcps))
    mask = residuals <= 1
    mask = residuals <= 1
    if len(np.array(base_gcps)[mask]) < 3:
    if len(np.array(base_gcps)[mask]) < 3:
        logging.info(f'Affine residuals: {residuals}.')
        raise ValueError(f'Unable to find enough points to compute an affine transformation. Found {len(np.array(dst_gcps)[mask])} points, but need at least 3.')
        raise ValueError(f'Unable to find enough points to compute an affine transformation. Found {len(np.array(dst_gcps)[mask])} points, but need at least 3.')


    affine.estimate(np.array(base_gcps)[mask], np.array(dst_gcps)[mask])
    affine.estimate(np.array(base_gcps)[mask], np.array(dst_gcps)[mask])