Spatialite: Speedup your query with spatial indexing

Spatialite, like any good spatial data management system, can build a spatial index for your layers. Using this index in your spatial queries will dramatically shorten the runtime for that query. The latest version of Spatialite offers a nice compact format for using a spatial index. To demonstrate, I created a point layer of 20,000 theoretical store locations, with sales data for each store, and a polygon layer of over 280 “local councils”. My mission it to sum up the total sales in each local council. So I need to find which stores are located in each local council and aggregate sales for those stores.
Continue reading

Sampling points along a line with R

Creating a random set of points is a standard GIS technique, used often for setting up sampling or monitoring locations. Among the tools in FOSS GIS software which offer this function are QGIS (in Vector -> Reseach Tools -> Random points) or GRASS (using the v.random module). But suppose you need points spread randomly along a line feature? The R-project package of spatial functions called ‘sp‘ can do just that.

Continue reading

My Acacia seedling

Last week’s bike ride was out to the wadi where “my” acacia seedling is struggling to survive. I was very satisfied to see that, with the cool winter weather, the plant has become green and seems to be doing well. On my last visit, during the scorching summer months, the seedling looked totally dried out and dead, and I wasn’t sure it would pull thru. We haven’t had any rainfall yet, just the lower temperature was enough to allow the plant to start growing again. These thorny buggers have eons of evolution behind them, so I guess they know how to deal with the arid season. I’ll visit again later during the winter and take measurements.