Skip to content

Consider polar projections in Raster spatial join bounding box calculation #1130

Description

@paleolimbot

Filing this as an issue so we don't forget!

Extent reprojection requires internal points across the whole domain, not just the boundary. I can demonstrate if desired, my approach is in reproj::reproj_extent

Originally posted by @mdsumner in #1126

This is a great point...in addition to the example in the linked thread, any polar projection when converted to lat/lon falls in this category (e.g., the south pole is directly in the middle of some projections, so even densifying the edges is not going to work, and even including points in the middle is not going to necessarily catch the minimum). Some of the minimum values it covers might also not be coverable by the other projection (e.g., mercator projections with a sigularity at the south pole).

Probably a geography covering is more accurate than a bounding box.

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

      Milestone

      No milestone

      Relationships

      None yet

      Development

      No branches or pull requests

      Issue actions