Hello,

I have the same problem as Costas Douvis, reported in
https://stat.ethz.ch/pipermail/r-sig-geo/2008-June/003713.html

I have a matrix of values, a matrix of latitudes and a matrix of longitudes.
The map is "rectangular" in the Lambert projection. So I did the following:

library(sp)
df <- data.frame(lat=as.numeric(xlat), lon=as.numeric(xlon),
z=as.numeric(hgt))
coordinates(df) <- ~lon+lat         # SpatialPointsDataFrame
proj4string(df) <- CRS("+proj=lcc +a=637000 +lat_1=30 +lat_2=60 +lat_0=40
+lon_0=-7.5")
spplot(df, scales=list(draw=TRUE))

But the map is no projected (it should be a rectangle),

As I want to have lat/lon on the axis of the map, and automatically produce
lots of maps, converting everything to lcc is not straightforward.

I am missing something? Any other packages?

Thank you.
Kind regards,

Rosa Trancoso

-------------
PhD Student
Lisbon Technical University, Portugal
email: arosa@ist.utl.pt
http://meteo.ist.utl.pt

	[[alternative HTML version deleted]]


