cont-limno / lagosusgis Goto Github PK

file <- "list_locality_postcode_meander_valley.tab"
mvfile <- system.file(file.path("extdata/tab", file), package="vapour")

## doesn't need names, they aren't used
## order is assumed
e <- c(xmin = 469750, xmax = 500840, 
       ymin = 5381892, ymax = 5393322)

library(vapour)
## no change
length(vapour_read_geometry(mvfile, extent = e))
#> [1] 58

## needs sql to be run as well
g <- vapour_read_geometry(mvfile, extent = e, sql = "SELECT * FROM list_locality_postcode_meander_valley")
length(g)  ## 9, they intersect the bounding box / extent
#> [1] 9


## compare with sf plot
library(sf)
#> Linking to GEOS 3.6.2, GDAL 2.3.2, PROJ 4.9.3
x <- sf::read_sf(mvfile)
plot(x[1], reset = F)
plot(sf::st_cast(sf::st_as_sfc(g)), add = TRUE, lwd =5)
rect(e[1], e[3], e[2], e[4])

ll_poly <- LAGOSUSgis::query_gis("LAGOS_US_All_Lakes_1ha", "lagoslakeid", 7922)

start_time <- Sys.time()
test <- LAGOSUSgis::query_gis("LAGOS_US_All_Lakes_1ha", "lagoslakeid", 10)
start_time - Sys.time()