pub fn estimate_homography_rect_to_img( rect_pts: &[Point2<f32>], img_pts: &[Point2<f32>], ) -> Option<Homography>
Estimate H such that: p_img ~ H * p_rect