alter the lon & lat lines and coastline according to the grid coordinates
For European projects, we commonly use the official European Terrestrial Reference System (www.euref.eu). You can convert your data from the longlat system to the ETRS using:
library(maptools)
proj4string(pointmap) <- CRS("+proj=longlat +datum=WGS84")
pointmap.etrs <- spTransform(pointmap, CRS("+init=epsg:3035"))
spplot(pointmap.etrs[1])
To simply (without any intervention) generate a gridded surface showing the change of values (the pattern), see: https://stat.ethz.ch/pipermail/r-sig-geo/2008-June/003703.html Many European GI data you can obtain from: EEA (http://dataservice.eea.europa.eu); but also take a look at: 2. Land cover maps (http://www-gem.jrc.it/glc2000/) 3. Landsat images (http://image2000.jrc.it/) 4. General type maps (http://www.inspire-geoportal.eu/) Tom Hengl http://spatial-analyst.net -----Original Message----- From: r-sig-geo-bounces at stat.math.ethz.ch [mailto:r-sig-geo-bounces at stat.math.ethz.ch] On Behalf Of Costas Douvis Sent: dinsdag 10 juni 2008 17:43 To: r-sig-geo at stat.math.ethz.ch Subject: [R-sig-Geo] alter the lon & lat lines and coastline according to the grid coordinates Hi everyone My problem might be easy to deal with but I have already spent quite some time and effort with no result. Any help will be appreciated My data come from an RCM that I ran, Regcm3. My domain spans over the wider area of Europe. The model grid is Lambert conformal. This means that the data values that come from the same line or column of a data matrix do not correspond to the exact same latitude or longitude. I have the values of latitude and longitude of each grid point in 2 matrices with the same dimensions as my data So what I need to do is to plot those data on a map. I guess that the best way to do this is to plot them on a rectangle map, preferably using filled.contour (forgetting for a moment that the grid is actually bended), and afterwards overlay the latitude and longitude lines and the coastline as they really are on that map (i.e. the lat and lon lines should not be straight) Is that possible? Or do you have a better suggestion?
Kostas Douvis PhD Student University of Athens - Department of Geography and Climatology Academy of Athens - Research Centre for Atmospheric Physics and Climatology email: cdouvis at geol.uoa.gr tel: +30-210-8832048, +30-210-8847280 fax: +30-210-8842098 _______________________________________________ R-sig-Geo mailing list R-sig-Geo at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo