Hi Edzer, no problem at all - Robert's approach is much less verbose,
but I think this explicit method is easier to understand and share as
it exposes the underlying indexing.
I agree there's a sense that it is image processing, but in the "GIS
raster" sense (image, DEM, data, etc.) it's a useful
Michael, to be hounest, it even surprised me, being the author, that
data(meuse.grid)
gridded(meuse.grid) = ~x+y
fullgrid(meuse.grid) = TRUE
image(meuse.grid[104:1,])
would not flip the image. It just works hard to not destroy the spatial
location of all individual pixels -- it will select all ro
Ergh, sorry for the update - a pox on Gmail for sabotaging my
attempts at plain text! I've attached the functions in a text file to
avoid [ampersand] to " at " conversion ...
-- Forwarded message --
From: Michael Sumner
Date: Tue, Nov 3, 2009 at 8:45 AM
Subject: Re: [R-sig-Ge
