[R-sig-Geo] Generate a coordinate (latitude) image?

[Manuel Schneider]

Jonathan Greenberg <greenberg <at> ucdavis.edu> writes:

> 
> Folks:
> 
> Is there any way to take a raster and generate two outputs where
> each pixel is the center coordinate of that pixel (e.g. you will have
> one x- image and one y-image)?  I'm trying to produce a "latitude map"
> for solar calculations.
> 
> --j
> 

You could use coordinates(as(obj, "SpatialPixels")) but only for relatively
small rasters. And why not calculate by hand?