gis code snippets in r