Commit 03ab8634 authored by Jay's avatar Jay
Browse files

Updates for affine and CHANGELOG.

parent 9fa3129a
Loading
Loading
Loading
Loading
+4 −0
Original line number Original line Diff line number Diff line
@@ -38,6 +38,10 @@ release.
### Added
### Added
- Debug logging to `place_points_in_overlap` and `distribute_points_in_geom` to make debugging issues easier.
- Debug logging to `place_points_in_overlap` and `distribute_points_in_geom` to make debugging issues easier.


### Changed
- cluster submission object loading to use [`joinedload('*')`](https://docs.sqlalchemy.org/en/20/orm/queryguide/relationships.html#wildcard-loading-strategies) in order to properly load and expunge measures on point objects. This change is at the query level (and not the mapper level).
- Acceptable affine residual bumped to 1px, better support for KaguyaTC data.

### Fixed
### Fixed
- Error in `find_interesting_feature` that was mis-using the ROI API. This bug was introduced in 1.2.0 when the ROI API updated.
- Error in `find_interesting_feature` that was mis-using the ROI API. This bug was introduced in 1.2.0 when the ROI API updated.


+1 −1
Original line number Original line Diff line number Diff line
@@ -109,7 +109,7 @@ def estimate_affine_from_sensors(reference_image,
    # error greater than or equal to 0.1px.
    # error greater than or equal to 0.1px.
    affine.estimate(np.array(base_gcps), np.array(dst_gcps))
    affine.estimate(np.array(base_gcps), np.array(dst_gcps))
    residuals = affine.residuals(np.array(base_gcps), np.array(dst_gcps))
    residuals = affine.residuals(np.array(base_gcps), np.array(dst_gcps))
    mask = residuals <= 0.1
    mask = residuals <= 1
    if len(np.array(base_gcps)[mask]) < 3:
    if len(np.array(base_gcps)[mask]) < 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.')
        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.')