[Biomod-commits] Importing georgraphical coordinates in R

Wilfried Thuiller wilfried.thuiller at ujf-grenoble.fr
Sat Sep 11 14:59:47 CEST 2010


You have different way.
Either your data are already in a shapefile with associated coordinate files, or they are simply a txt files with coordinates and you know the projection system

1- Use the rgdal package.
library(rgdal)
poly <- readOGR("path_to_the_right_folder", layer="name_of_the_files_without_extension")
summary(poly)


2- In a text file you import in R with the name 'plots'
Let's say you use the WGS84 projections and the coordinates are in columns X_WGS84 and Y_WGS84
 pointsWGS84 <- plots[,c("X_WGS84", "Y_WGS84")]
 colnames(pointsWGS84) <- c("x","y")
 coordinates(pointsWGS84) <- c("x", "y")
### stipulate to R what is the projection. Look at the help file for proj4string
 proj4string(pointsWGS84) <- CRS("+proj=latlong +datum=WGS84") # WGS84 projection string
 head(pointsWGS84 at coords)

#if you want to transform the coordinates in another projection system. For instance here in Lambert for France II étendu.
# first create the projection file: 

proj4NTFII <- "+proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs"
  #EPSG Projection 27582 - NTF (Paris) / France II


pointsNTFII <- spTransform(pointsWGS84, CRS(proj4NTFII))
head(pointsNTFII at coords) 

You have most of the projection system here:

http://www.remotesensing.org/geotiff/proj_list/




Le 10 sept. 2010 à 22:36, Manuel Spínola a écrit :

> Dear list members,
> 
> How do I import data in geographical coordinates to R?
> 
> ref1                       ref2
> N10 41 08.7     W84 10 46.3
> N10 41 06.8     W84 10 46.4
> N10 41 06.4     W84 10 44.6
> N10 41 05.6     W84 10 46.0
> 
> Best,
> 
> Manuel
> 
> -- 
> Manuel Spínola, Ph.D.
> Instituto Internacional en Conservación y Manejo de Vida Silvestre
> Universidad Nacional
> Apartado 1350-3000
> Heredia
> COSTA RICA
> mspinola at una.ac.cr
> mspinola10 at gmail.com
> Teléfono: (506) 2277-3598
> Fax: (506) 2237-7036
> 
> _______________________________________________
> Biomod-commits mailing list
> Biomod-commits at lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/biomod-commits

--------------------------
Dr. Wilfried Thuiller
Laboratoire d'Ecologie Alpine, UMR CNRS 5553
Université Joseph Fourier
BP53, 38041 Grenoble cedex 9, France
tel: +33 (0)4 76 51 44 97
fax: +33 (0)4 76 51 42 79

Email: wilfried.thuiller at ujf-grenoble.fr
Home page: http://www.will.chez-alice.fr
Website: http://www-leca.ujf-grenoble.fr/equipes/tde.htm

FP6 European MACIS project: http://www.macis-project.net
FP6 European EcoChange project: http://www.ecochange-project.eu





-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/biomod-commits/attachments/20100911/c23eab61/attachment.htm>


More information about the Biomod-commits mailing list