Let's say one wants to find out a box around a given latitude/longitude point. How to do this without GIS software using an artisan UNIX approach? Hell, if it works for coffee! The steps are as follows:

1. Convert the centre point into UTM so we have meters rather ...

This note details how to use parallel processing in python accessing shared memory. The usage case is when you have a process that reads in a lot of data, and where the data can be processed in parallel, such as when reading and processing images/stacks of images in a ...

Quite often, one wants to generate some data at high resolution (say process some image or images) and then calculate some relevant spatial statistics at some other resolution. For example, you might want to process Landsat TM data at 30m resolution, and might want to aggregate it to a resolution ...

Supervised classification of EO data uses a set of samples from known patterns (usually reflectance spectra) to decide whether a given pattern belongs to one class or another. In landcover applications, one goes to the field, and observes that a given location is indeed class $\omega_c$. These ground observations usually ...

Usually, we remote sensing types ignore country boundaries: they don't really make much sense, as they are not aligned with MODIS pixels ;-) However, I was asked what's the easiest way to produce a global mask of countries, so that all th grid cells (say at 0.5 degree ...