Working with @pudo on a general lonlat to region-code resolver based on point-in-polygon testing in shapefiles. Quite useful for #ddj