Я хочу использовать:
st_join(polygon_A, polygon_B, join=st_intersects)
, чтобы получить те же выходные данные, что и в текущей версии, с дополнительным столбцом, который указывает долю области, захваченной с polygon_A на polygon_B? В идеале это будет работать с largest=TRUE (возвращая уникальные записи polygon_A только с одним совпадением из polygon_B) или с largest=FALSE (возвращая все совпадения с polygon_B с пропорциями захваченных областей на polygon_A).
Мой обходной путь должен был использовать st_intersection(polygon_A, polygon_B), вычислить площадь новых фигур (результат пересечения), используя st_area, а затем разделить его на площадь исходной формы polygon_A. Есть ли лучший способ добиться этого? Или, возможно, планирует добавить это как функцию к st_join.
Для записи и во избежание проблемы XY моя цель состоит в том, чтобы определить случаи, когда пропорции близки к 0,50 / 0,50 или 0,70 / 0,30, и более внимательно посмотреть, чтобы определить, является ли largest=TRUE осмысленный. Меня не очень интересуют случаи, когда пропорции 0,99 / 0,01 - largest=TRUE достаточно для этих случаев.