<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 12pt;
font-family:Calibri
}
--></style></head>
<body class='hmmessage'><div dir='ltr'><br>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.<br><br>Un saludo,<br><br><div>> From: oscar.perpinan@upm.es<br>> To: public-markbermejo90-PkbjNfxxIARBDgjK7y7TUQ@plane.gmane.org<br>> CC: public-r-help-es-0bNBQ1PAWB4BXFe83j6qeQ@plane.gmane.org<br>> Subject: Re: Duda ' proj4string'<br>> Date: Mon, 10 Aug 2015 07:02:02 +0200<br>> <br>> <br>> <br>> <br>> Hola,<br>> <br>> Intuyo que el error se produce porque "area" no tiene definido el<br>> proj4string. Creo que lo resolverás cambiando el orden de estas líneas<br>> así:<br>> <br>> proj4string(area) <- "+proj=utm +zone=30 +north +ellps=WGS84 +datum=WGS84 +units=m +no_defs"<br>> area <- spTransform(area, CRS("+proj=longlat +datum=WGS84"))<br>> <br>> También puedes probar a leer el fichero shp con readOGR (paquete rgdal)<br>> en lugar de con readShapePoly.<br>> <br>> Saludos.<br>> <br>> Oscar.<br>> -- <br>> Oscar Perpiñán Lamigueiro<br>> Dpto. Ing. Eléctrica, Electrónica, Automática y Física Aplicada<br>> Escuela Técnica Superior de Ingeniería y Diseño Industrial<br>> URL: http://oscarperpinan.github.io<br>> <br>> <br></div> </div></body>
</html>