Bike-commuter & East Phillips resident focused on environmental (in)justice in so-called "Minneapolis."
I do spatial analysis, community outreach, & mapping with open-source tech
| Website | https://rwhendrickson.github.io/Portfolio |
Just a thought, but I might be totally off - Do you have a spatial index?
I'm not familiar with SQLite, but I remember my friend having a similar issue when they had a lot of points in a postgres database (hour+ long waits for queries). After adding a spatial index, the spatial queries took seconds
Cyclone #Chido strikes #Mayotte with winds of 215 km/h (130 mph).
Now the cyclone is heading for #Mozambique, where it is expected to make landfall south of Pemba early on Sunday.
Global tree percent cover and impervious surfaces (urban) cover for years 1985 and 2022 at 30 m spatial resolution based on the GLC_FCS30D data set produced by the State Key Laboratory of Remote Sensing, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing. The tree cover has been estimated using the 37 class legend (GLC_FCS30D_legend.csv) and is only an approximation of might have been the actual tree cover. Note: the GLC_FCS30D data set is probably over-estimating tree-based vegetation types in arid and semi-arid areas so this would need to be corrected; in addition, several artifacts due to the tiling system / locally fitted classification are visible in the GLC_FCS30D maps. The GLC_FCS30D data set is available for viewing from https://OpenLandMap.org and is explained in detail in: Zhang, X., Liu, L., Chen, X., Gao, Y., Xie, S., and Mi, J. (2021): GLC_FCS30: global land-cover product with fine classification system at 30 m using time-series Landsat imagery, Earth Syst. Sci. Data, 13, 2753–2776, https://doi.org/10.5194/essd-13-2753-2021 It is estimated that world has lost 1/3rd of original pre-industrial forest areas, however, many current forest areas also have a thinner tree cover than before. Similar data sets (not used here) are the JRC's Tropical Moist Forests product (TMFv2023) which is also available for download although it only covers tropical forest areas, and the MapBiomas data set: Souza Jr, C. M., Z. Shimbo, J., Rosa, M. R., Parente, L. L., A. Alencar, A., Rudorff, B. F., ... & Azevedo, T. (2020). Reconstructing three decades of land use and land cover changes in brazilian biomes with landsat archive and earth engine. Remote Sensing, 12(17), 2735. https://doi.org/10.3390/rs12172735 Vancutsem, C., Achard, F., Pekel, J. F., Vieilledent, G., Carboni, S., Simonetti, D., ... & Nasi, R. (2021). Long-term (1990–2019) monitoring of forest cover changes in the humid tropics. Science advances, 7(10), eabe1603. https://doi.org/10.1126/sciadv.abe1603 To convert GLC_FCS30D to tree cover and urban cover I've used the following code: library(terra); library(parallel) g1 = terra::vect("glc_100km_tiles.gpkg")gcl = read.csv("GLC_FCS30D_legend.csv")rcl = gcl[,c("GLC_FCS30D","Tree_cover")] fraction.cover <- function(tif, tile, rcl){ bb = as.vector(ext(tile)) year = substr(strsplit(tif, "_")[[1]][6], 1, 4) out.tif = paste0("./t", year, "/tree.cover_30m_", year, "_", paste(bb, collapse = "."), ".tif") if(!file.exists(out.tif)|!file.exists(gsub("30m", "1km", out.tif))){ r = terra::rast(tif) r.t = terra::crop(r, ext(tile)) m = terra::classify(r.t, rcl) writeRaster(m, filename=out.tif, wopt= list(gdal=c("COMPRESS=DEFLATE")), datatype='INT1U', NAflag=255) system(paste0('gdal_translate ', out.tif, ' ', gsub("30m", "1km", out.tif), ' -co COMPRESS=DEFLATE -tr 0.008333333 0.008333333 -r "average" -q')) m0 = terra::ifel(r.t==190, 100, 0) writeRaster(m0, filename=gsub("tree", "urban", out.tif), wopt= list(gdal=c("COMPRESS=DEFLATE")), datatype='INT1U', NAflag=255) rm(m); rm(m0); gc() tmpFiles(remove=TRUE) } } ## run in parallel: tif = "lc_glc.fcs30d_c_30m_s_20220101_20221231_go_epsg.4326_v20231026.tif" #tif = "lc_glc.fcs30d_c_30m_s_19850101_19851231_go_epsg.4326_v20231026.tif" x = parallel::mclapply(1:length(g1), function(i){try( fraction.cover(tif, tile=g1[i], rcl=rcl) )}, mc.cores=80)
“There’s a certain amount of hubris involved in saying we’re going to stick a 40-mile pipeline through the earth, and if we mess up these wetlands that nature has created over thousands of years, we’ll just put them back together again." #ShutDownLine5
https://insideclimatenews.org/news/29112024/wisconsin-enbridge-line-5-reroute-impacts/
It’s fresh snow season, which means that we can see bicycle tracks and bike lane usage in #Toronto again. A very well-used bike lane, even in snowy weather!
I present a map type @dschep and I are calling the Fog Of War On Cars. It will mask of your area that isn't within 50m of bicycle parking as determined by OpenStreetMap data.
Here's what it looks like for Redmond WA, the supposed "Bicycle Capital of the Northwest".
Pan the map to your town an hit run. Let me know if you find anything interesting!
https://overpass-ultra.us/#query=gist:8aa0f652498160c3a51bde830b680bab
Zoom Earth
NetBlocks
Tomislav Hengl
Streets.mn
Ian R Buck
Three Waters PRT
Majd Al-Shihabi 🏴مجد الشهابي
watmildon
Brian Mitchell