JuliaCon 2026

Fast geospatial lookups across projections using SphericalSpatialTrees.jl
2026-08-12 , Muschel — N3

The SpatialTreeInterface defined in GeometryOps.jl provides an efficient way for geometrical queries of polygons that share a crs by relying on search trees whose branches are characterized by rectangular bounding boxes. However, when working across different projections, for example when trying to identify intersecting polygons from different projections on the sphere, rectangles in one projection do not translate into rectangles in another projection, so the tree traversal will not be accurate.

To solve this problem, we extended the SpatialTreeInterface in SphericalSpatialTrees.jl by replacing rectangular bounding boxes with SphericalCaps, using these to characterize the extent of all branches in a spatial tree. This enables users to do combined tree searches across polygons based on different projections, as is common e.g. in Discrete Global Grid System. The presentation will demonstrate the basic concept of the SphericalSpatialTrees.jl package as well as a few downstream applications.

I order to solve the problem of

I am a physicist by training and am currently studying Global Biogeochemical Cycles in the Earth System using Remote Sensing, Meteorological and other data sets based at the Max-Planck-Institute for Biogeochemistry, Jena, Germany.
My first commit to my first Julia package dates back to the year 2012 and since then I have authored and contributed to packages in the Julia Geodata and processing ecosystem, examples are NetCDF.jl, Zarr.jl, DiskArrays.jl, YAXArrays.jl EarthDataLab.jl and others. Some may know me under my github tag @meggart

This speaker also appears in: