Muchas gracias, ya no me da el error, lo que pasa que no me cerca el área de la 
cuenca, no me hace nada la función "geom_polygon()". Os paso los dos ficheros 
que se necesitan más el script, para ver si podéis probar y ver qué pasa.

Un saludo,

> From: [email protected]
> To: [email protected]
> CC: [email protected]
> Subject: Re: Duda ' proj4string'
> Date: Mon, 10 Aug 2015 07:02:02 +0200
> 
> 
> 
> 
> Hola,
> 
> Intuyo que el error se produce porque "area" no tiene definido el
> proj4string. Creo que lo resolverás cambiando el orden de estas líneas
> así:
> 
> proj4string(area) <- "+proj=utm +zone=30 +north +ellps=WGS84 +datum=WGS84 
> +units=m +no_defs"
> area <- spTransform(area, CRS("+proj=longlat +datum=WGS84"))
> 
> También puedes probar a leer el fichero shp con readOGR (paquete rgdal)
> en lugar de con readShapePoly.
> 
> Saludos.
> 
> Oscar.
> -- 
> Oscar Perpiñán Lamigueiro
> Dpto. Ing. Eléctrica, Electrónica, Automática y Física Aplicada
> Escuela Técnica Superior de Ingeniería y Diseño Industrial
> URL: http://oscarperpinan.github.io
> 
> 
                                          

Attachment: area_cuenca.R
Description: Binary data

IDENTIFICA      PARAJE  COTA_METRO      POINT_X POINT_Y
Pluvi¾metro 1   Arromoro        1831.000        362151.00000000000      
4470922.00000000000
Pluvi¾metro 2   Atalaya 1790.000        361075.00000000000      
4470230.00000000000
Pluvi¾metro 3   Dehesa  746.000 358566.00000000000      4475056.00000000000 
Pluvi¾metro 4   Gasolinera      736.000 358875.00000000000      
4475400.00000000000
Pluvi¾metro 5   Morales 1287.000        361922.00000000000      
4473439.00000000000
Pluvi¾metro 6   Pena Parda      1420.000        360098.00000000000      
4470743.00000000000
Pluvi¾metro 7   Trampalones     920.000 360004.00000000000      
4472896.00000000000

Attachment: venero_Project.shp
Description: Binary data

_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es

Responder a