WIP: geometry lookup

Project.toml

@@ -16,6 +16,7 @@ FillArrays = "1a297f60-69ca-5386-bcde-b61e274b549b"
 Flatten = "4c728ea3-d9ee-5c9a-9642-b6f7d7dc04fa"
 GeoFormatTypes = "68eda718-8dee-11e9-39e7-89f7f65f511f"
 GeoInterface = "cf35fbd7-0cd7-5166-be24-54bfbe79505f"
+GeometryOps = "3251bfac-6a57-4b6d-aa61-ac1fef2975ab"
 GeometryOpsCore = "05efe853-fabf-41c8-927e-7063c8b9f013"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 Missings = "e1d29d7a-bbdc-5cf2-9ac0-f12de2c33e28"
@@ -25,6 +26,7 @@ ProgressMeter = "92933f4c-e287-5a05-a399-4b506db050ca"
 RecipesBase = "3cdcf5f2-1ef4-517c-9805-6587b60abb01"
 Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 Setfield = "efcf1570-3423-57d1-acb7-fd33fddbac46"
+SortTileRecursiveTree = "746ee33f-1797-42c2-866d-db2fce69d14d"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 
 [weakdeps]
@@ -69,6 +71,7 @@ GeoDataFrames = "0.3"
 GeoFormatTypes = "0.4"
 GeoInterface = "1.0"
 GeometryBasics = "0.4"
+GeometryOps = "0.1.19"
 GeometryOpsCore = "0.1.1"
 Makie = "0.20, 0.21, 0.22"
 Missings = "0.4, 1"
@@ -83,6 +86,7 @@ Reexport = "0.2, 1.0"
 SafeTestsets = "0.1"
 Setfield = "0.6, 0.7, 0.8, 1"
 Shapefile = "0.10, 0.11"
+SortTileRecursiveTree = "0.1.1"
 Statistics = "1"
 StatsBase = "0.34"
 Test = "1"

src/Rasters.jl

@@ -21,6 +21,7 @@ import Adapt,
        FillArrays,
        Flatten,
        GeoInterface,
+       GeometryOps,
        GeometryOpsCore,
        OffsetArrays,
        ProgressMeter,
@@ -29,6 +30,7 @@ import Adapt,
        RecipesBase,
        Reexport,
        Setfield,
+       SortTileRecursiveTree,
        Statistics
 
 Reexport.@reexport using DimensionalData, GeoFormatTypes
@@ -62,7 +64,7 @@ export Planar, Spherical
 export AbstractRaster, Raster
 export AbstractRasterStack, RasterStack
 export AbstractRasterSeries, RasterSeries
-export Projected, Mapped
+export Projected, Mapped, GeometryLookup
 export Band
 export missingval, boolmask, missingmask, replace_missing, replace_missing!,
        aggregate, aggregate!, disaggregate, disaggregate!, mask, mask!,
@@ -76,6 +78,7 @@ export Extent, extent
 const DD = DimensionalData
 const DA = DiskArrays
 const GI = GeoInterface
+const GO = GeometryOps
 const LA = Lookups
 
 # DimensionalData documentation urls
@@ -119,6 +122,9 @@ const RasterSeriesOrStack = Union{AbstractRasterSeries,AbstractRasterStack}
 include("utils.jl")
 include("skipmissing.jl")
 
+include("geometry_lookup/geometry_lookup.jl")
+include("geometry_lookup/methods.jl")
+
 include("table_ops.jl")
 include("create.jl")
 include("read.jl")