> 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
Jonathan, If your data are in lon/lat r <- raster() x <- init(r, v='x') y <- init(r, v='y') Or "by hand", for lon/lat do lat <- yFromRow(r, 1:nrow(r)) if the raster is not lon/lat xy <- SpatialPoints(xyFromCell(r, cellFromRowCol(1:nrow(r), 1))) lat <- coordinates( spTransform(xy, CRS( +proj=longlat +datum=WGS84 )))[,2] In cases like this, I stop here. Rather than making a latitude raster with lots of repeated values, I would work with the lat vector. But you could do something like this: x <- t((t(r) * lat)) plot(x) Best, Robert -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Generate-a-coordinate-latitude-image-tp6375697p6378473.html Sent from the R-sig-geo mailing list archive at Nabble.com. _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
