[Gvsig_usuarios] Eliminar polígonos dentro de un raster (enclavados)

manuel peixoto manuel.8777b en hotmail.com
Mie Nov 24 19:24:07 CET 2010


Hola Nacho. 
primero darte las gracias por responderme a mi pregunta y segundo volver a darte las gracias por anticipado, aunque sea por leer el coñazo que viene a continuación.

Si. Es esto lo que te comento.... pero te explico.
Yo lo que hago es cortar un mde.asc general con el, por ejemplo, Enclavado1=> Enclavado1_mde.asc. 
Después voy a Sextante/Herramientas de cálculo para capas raster y la divido por si misma "/" (o bien con la calculadora de mapas). Ahora tengo las zonas de los enclavados con valor 1 y el resto con valor -99999 (celdas sin datos).Nota=> Cuando hago esta operación no le puedo dar un formato "raster.asc", si no que lo tengo que guardar como "archivo temporal". Si lo intento guardar en una carpeta dentro de MisDocumentos, con extension *.asc o *.tiff me genera un raster sin ningún tipo de datos. ¿que pasa aquí?. No lo entiendo Esto me lo tiene hecho bien. Aunque no sé si era con Sextante 0.33. Como siga persistiendo este error voy a tener que probar esta versión de sextante en gvSIG 1.9.
Bueno... una vez tengo la máscara.... voy a "Herramientas básicas para capas raster/Invertir máscara" y ahí ya no me genera nada de nada. Ni guardándolo en archivos temporales.
¿Que estoy haciendo mal? 
x------------------------------------------------x
Ahora un tema que en principio está aparte... pero que después se relaciona con la primera pregunta.
Cuando tengo un mde de una zona determinada (generalmente un monte, soy ingeniero técnico forestal) y hago un mde a partir de las curvas de nivel lo hago de la siguiente manera:
1.- Consigo una curvas de nivel que abarquen algo más (bastante más, por si acaso) que la superficie a estudiar.
2.- Rasterizo y Relleno Celdas sin datos (Lo hago en un paso con un "modelo" que me hice yo... aunque lo borré sin querer... bueno... a veces pasa) y obtengo mi "MDE.asc"
3.- Ahora tengo un mde que abarca más que la zona de estudio (que por supuesto tengo en un *shp).
4.- Corto mi mde con mi polígono y obtengo mi "Polígono_mde.asc". Como es un polígono irregular el resto de la superficie regular con celdas me la muestra como celdas sin datos.
5.- Una vez tengo este raster "Polígono_mde.asc" y después de aplicarle una tabla de color y hacer que las celdas sin datos sean transparentes... voy a "Sextante/Herramientas categóricas para capas raster/Estadísticas de clases y obtengo una tabla *dbf la cual abro en Excel (o openOffice).
6.- Si a mi tabla "Poligono_mde.xlsx" le elimino el rango -99999 y le sumo el Area obtengo, con un error ridículo, el área planimétrica de mi polígono.
7.- Si a mi "Polígono_mde.asc" le aplico el algorítmo Sextante/Geomorfometría y análisis del relieve/Ãrea real obtengo un nuevo raster.
8.- si a mi "Polígono_mdeAreaReal.asc" le calculo las estadísticas de clases me da el mismo área que en el raster "Polígono_mde.asc"
la pregunta es ¿Como puedo calcular el área real de mi polígono? Valdría realizar sobre "Polígono_mdeAreaReal.asc" el algoritmo "Sextante/Herramientas básicas para capas raster/Estadísticas básicas?
x--------------------------------------------x
Ahora, y volviendo a la primera pregunta
Podría restarle a los datos de  mi archivo "Poligono_mde.xlsx" los datos de "Enclavado1_mde.xls (y Enclavado 2, 3,4 y 5)  pertenecientes a los mismos valores del campo IDENTIFICA en los raster de mde, pendientes y orientaciones? 
obtendría así un resultado válido sobre la superficie de mi polígono sin los enclavados?
x-------------------------------------------x
Muchas gracias y recibe un cordial saludo. 
Y también gracias a todo el foro... y espero que la información que aporto con mis dudas le pueda valer a alguien.
Date: Wed, 24 Nov 2010 10:33:07 +0100 
From: nachouve en gmail.com
To: gvsig_usuarios en listserv.gva.es
Subject: Re: [Gvsig_usuarios]	Eliminar polígonos dentro de un raster (enclavados)

Hola Manuel,

No he entendido bien el caso que comentas. Lo que quieres es asignar a 5 zonas de tu raster valores "nulos" y esas zonas están definidas en una capa vectorial de polígonos... ¿es así?

Si es esto, podrías rasterizar la capa de polígonos y darle valor 1 a todo menos a las zonas a "eliminar". Luego con la "calculadora de capas" conseguirías el resultado usando como máscara esta capa.


Un saludo,
Nacho



El 23 de noviembre de 2010 20:29, manuel peixoto <manuel.8777b en hotmail.com> escribió:






Hola a tod en s.
Estoy trabajando con gvSIG 1.9 + SEXTANTE 0.5 sobre Windows Xp

Tengo un polígono al cual le he realizado un estudio del terreno que consiste en un raster de modelo digital de elevaciones, un raster de pendientes y un raster de orientaciones.

Ahora necesito quitarle... eliminar.... 5 polígonos (enclavados) dentro de mi superficie... porque no entran dentro del estudio.
como puedo hacer?
Hay algún proceso en sextante que sea inverso a cortar capa raster con polígonos de herramientas básicas para capas raster?

Gracias
Un saludo 		 	   		  

_______________________________________________

gvSIG_usuarios mailing list

gvSIG_usuarios en listserv.gva.es



Para ver histórico de mensajes, editar sus preferencias de usuario o darse de baja en esta lista, acuda a la siguiente dirección:



http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios



_______________________________________________
gvSIG_usuarios mailing list
gvSIG_usuarios en listserv.gva.es

Para ver hist�rico de mensajes, editar sus preferencias de usuario o darse de baja en esta lista, acuda a la siguiente direcci�n:

http://listserv.gva.es/cgi-bin/mailman/listinfo/gvsig_usuarios 		 	   		  
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: http://listserv.gva.es/pipermail/gvsig_usuarios/attachments/20101124/457bf390/attachment.htm 


Más información sobre la lista de distribución gvSIG_usuarios