[R-sig-Geo] How to translate coordinates from cylindrical to a Lambert Azimuthal Equal Area Projection

Ivailo ubuntero.9161 at gmail.com
Mon Apr 22 12:14:40 CEST 2013

Dear Listmembers,

I have a 40 000 x 40 000 GeoTIFF file that I need to crop to an area
of interest (perhaps about 1000 x 1000 pixels, 25 x 25 m grid). The
file seems to be in the Lambert Azimuthal Equal Area Projection
(coord. ref. : +proj=laea +lat_0=52 +lon_0=10 +x_0=4321000
+y_0=3210000 +ellps=GRS80 +units=m +no_defs) while I use Google Earth
(using a Cylindrical projection with a WGS84 datum) to navigate and
delineate the extent of the area I'd like to crop from the original

As I don't have much experience working with Raster* objects as
defined in the raster package, I'd like to ask you for help in
"translating" the coordinates I get from Google Earth (i.e., decimal
degrees) to the ones needed to define the extent of interest for the
RasterLayer object I have created by reading the GeoTiff file via

Any insight would be much appreciated.

UBUNTU: a person is a person through other persons.

More information about the R-sig-Geo mailing list