[R-sig-Geo] alter the lon & lat lines and coastline according to the grid coordinates

[Tomislav Hengl]

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