Determinando elevaciones máximas y mínimas (estadística zonal) de un ráster utilizando rpy objetos en ambiente de PyQGIS | El Blog de José Guerrero

https://joseguerreroa.wordpress.com/2015/12/06/determinando-elevaciones-maximas-y-minimas-estadistica-zonal-de-un-raster-utilizando-rpy-objetos-en-ambiente-de-pyqgis/[07-12-2015 5:19:19]

← Uniendo rásteres con proyección sinusoidal en QGIS

Determinando elevaciones máximas y mínimas (estadística zonal) de un ráster utilizando rpy objetos en ambiente de PyQGISPublicado el 6 diciembre, 2015

En el GDAL/OGR Cookbook existe código para acometer la estadística zonal de un ráster empleando las librerías GDAL/OGR de Python. Sin embargo, la extensión de este contrasta con lo conciso que puede llegar a ser el mismo procedimiento con lenguaje R.

Para usar R con Python es necesario instalar en Debian los paquetes python-rpy y python-rpy2. Después de ello, utilizando el vectorial y el ráster dem de la imagen siguiente:

desarrollé el código a continuación, en ambiente de PyQGIS, que permite obtener los valores máximo y mínimo (estadística zonal) de un raster con base en un vectorial tipo polígono e introducir en la tabla atributiva de este los valores obtenidos en sus campos respectivos.

Check pagerank

Visitantes:1,346,170 visitas

Mi Perfil en Gis Stack Exchange

Próximo Curso 20 de Octubre de 2015

Suscripción por correo electrónicoEscribe tu dirección de correo electrónico para suscribirte a este blog, y recibir notificaciones de nuevos mensajes por correo.

Únete a otros 550 seguidores

Enlaces RSSRSS - EntradasRSS - Comentarios

Nube de CategoríasAutoCAD Climatología Código C++

El Blog de José Guerrero Especialmente sobre GNU/Linux, código Python y Sistemas de Información Geográfica (SIG)

Home Sobre el autor Curso GIS Enlaces Tips Downloads Inventarios ERDAS Videos ArcGis Comandos Linux GDAL Python

Procedimientos Windows Comandos R GRASS Scripts Ingenio Plugins

Seguir

Seguir “El Blog de José Guerrero”Recibe cada nueva publicación en tu buzón de correo electrónico.

Únete a otros 550 seguidores

Determinando elevaciones máximas y mínimas (estadística zonal) de un ráster utilizando rpy objetos en ambiente de PyQGIS | El Blog de José Guerrero

https://joseguerreroa.wordpress.com/2015/12/06/determinando-elevaciones-maximas-y-minimas-estadistica-zonal-de-un-raster-utilizando-rpy-objetos-en-ambiente-de-pyqgis/[07-12-2015 5:19:19]

from PyQt4.QtCore import *

import rpy2.robjects as robjectsimport os

r = robjects.r

mapcanvas = iface.mapCanvas()layers = mapcanvas.layers()

vp = layers[0].dataProvider()vpath_tmp = vp.dataSourceUri().split('|')vpath = vpath_tmp[0]

rp = layers[1].dataProvider()rpath = rp.dataSourceUri()

vroot, vfile = os.path.split(vpath)rroot, rfile = os.path.split(rpath)

r.setwd(vroot)

r.library('raster')

arg1 = 's <- shapefile(' + "\"" + vfile + "\"" + ')'arg2 = 'd <- raster(' + "\"" + rfile + "\"" + ')'

r(arg1)r(arg2)

print "Wait..."

demmin = r('s$demmin <- extract(d, s, fun = min, na.rm = TRUE)')demmax = r('s$demmax <- extract(d, s, fun = max, na.rm = TRUE)')

fields = [ QgsField('min', QVariant.Double), QgsField('max', QVariant.Double) ]vp.addAttributes( fields )layers[0].updateFields()

idx1 = layers[0].fieldNameIndex('min')idx2 = layers[0].fieldNameIndex('max')

n = len(demmin)

for i in range(n): new_values = { idx1 : demmin[i], idx2 : demmax[i] } vp.changeAttributeValues( {i:new_values} )

print "Done!"

Después de ejecutado el código anterior en la Python Console de QGIS se obtuvo el resultado esperado; tal como se expresa en la tabla de atributos del vectorial tipo polígono de la imagen siguiente:

Código Python Debian Debian testing Documento ArcGis GDAL GDAL/OGR GeoGebra gnuplot

GPS GRASS Lenguaje R Linux OpenStreetMap Plugin

Postgres+postgis PyQGIS QGIS SAGA GIS SAGA SIG

SEXTANTE SIG Software

Freeware Software Libre Teledetección Tips y trucos Transformar Coordenadas WordPress

Entradas recientesDeterminando elevaciones máximas y mínimas (estadística zonal) de un ráster utilizando rpy objetos en ambiente de PyQGISUniendo rásteres con proyección sinusoidal en QGISPlugin de QGIS para producir vectoriales por selección de rasgos con base en consultas tipo SQLDisplay de estadísticas de campos de atributos en un plugin de QGISDiferencia de polígonos que se solapan mediante PyQGIS

Archivosdiciembre 2015noviembre 2015octubre 2015septiembre 2015agosto 2015julio 2015junio 2015mayo 2015abril 2015marzo 2015febrero 2015enero 2015diciembre 2014noviembre 2014octubre 2014septiembre 2014agosto 2014julio 2014junio 2014mayo 2014abril 2014marzo 2014febrero 2014enero 2014diciembre 2013noviembre 2013octubre 2013septiembre 2013agosto 2013julio 2013junio 2013mayo 2013abril 2013marzo 2013febrero 2013noviembre 2012octubre 2012septiembre 2012agosto 2012julio 2012junio 2012

Construye un sitio web con WordPress.com

Determinando elevaciones máximas y mínimas (estadística zonal) de un ráster utilizando rpy objetos en ambiente de PyQGIS | El Blog de José Guerrero

https://joseguerreroa.wordpress.com/2015/12/06/determinando-elevaciones-maximas-y-minimas-estadistica-zonal-de-un-raster-utilizando-rpy-objetos-en-ambiente-de-pyqgis/[07-12-2015 5:19:19]

← Uniendo rásteres con proyección sinusoidal en QGIS

Share this:

Twitter Google Facebook Correo electrónico Imprimir

Esta entrada fue publicada en Lenguaje R, PyQGIS, QGIS, SIG, Software Libre. Guarda el enlace permanente.

mayo 2012abril 2012marzo 2012febrero 2012enero 2012diciembre 2011noviembre 2011octubre 2011septiembre 2011agosto 2011julio 2011mayo 2011abril 2011marzo 2011febrero 2011enero 2011diciembre 2010noviembre 2010octubre 2010septiembre 2010agosto 2010julio 2010junio 2010mayo 2010abril 2010febrero 2010enero 2010noviembre 2008marzo 2008febrero 2008junio 2007

CategoríasArcGis10.1AutoCADCódigo C++Código PythonClimatologíaDebianDebian testingDocumento ArcGisDocumento ArcViewFirefoxFlockGarminGDALGDAL/OGRGeoGebragnuplotGoogle mapsGPSGRASSIdiomasImpacto AmbientalJavaKindlekmplotLenguaje RLibreOfficeLinuxMapSourcemodisNewsletterNokia 5800OpenStreetMapPluginPostgres+postgisProcessingPyQGIS

Cargando...

Notificarme los nuevos comentarios por correo electrónico.

Recibir nuevas entradas por email.

Deja un comentario

Correo (requerido) (La dirección no se hará pública)

Nombre (requerido)

Sitio web

Determinando elevaciones máximas y mínimas (estadística zonal) de un ráster utilizando rpy objetos en ambiente de PyQGIS | El Blog de José Guerrero

https://joseguerreroa.wordpress.com/2015/12/06/determinando-elevaciones-maximas-y-minimas-estadistica-zonal-de-un-raster-utilizando-rpy-objetos-en-ambiente-de-pyqgis/[07-12-2015 5:19:19]

QGISSAGA GISSAGA SIGSEXTANTESIGSoftwareSoftware FreewareSoftware LibreSuelos y AguasTeledetecciónTips y trucosTransformar CoordenadasUbuntu 11.04Ubuntu 11.11Ubuntu 12Ubuntu 14.04UncategorizedVirtualizaciónWindowsWordPressWubi

MetaRegistrarseAccederRSS de las entradasRSS de los comentariosWordPress.com

Posts Más VistosManejo de arrays y matrices en pythonDeterminando elevaciones máximas y mínimas (estadística zonal) de un ráster utilizando rpy objetos en ambiente de PyQGISSerie de Fibonacci en C/C++Cómo crear una rejilla regularmente espaciada en AutoCAD 2013Convertir shapefiles en kml y viceversa con QGISConversión masiva de coordenadas UTM a geográficas decimales usando los elipsoides Hayford y WGS 84Determinación del Indice de Area Foliar (IAF) a partir de imágenes de satélite con GRASS en QGISConversión masiva de coordenadas geográficas decimales a UTM usando los elipsoides Hayford y WGS 84Convertir Lineas de Azimut y distancia a CoordenadasIntegración numérica: C++ Linux

diciembre 2015M T W T F S S

« Nov

1 2 3 4 5 6

7 8 9 10 11 12 13

14 15 16 17 18 19 20

21 22 23 24 25 26 27

28 29 30 31

Esta obra está bajo una licencia de Creative Commons.

Determinando elevaciones máximas y mínimas (estadística zonal) de un ráster utilizando rpy objetos en ambiente de PyQGIS | El Blog de José Guerrero

https://joseguerreroa.wordpress.com/2015/12/06/determinando-elevaciones-maximas-y-minimas-estadistica-zonal-de-un-raster-utilizando-rpy-objetos-en-ambiente-de-pyqgis/[07-12-2015 5:19:19]

El Blog de José GuerreroEl tema Twenty Ten. Funciona con WordPress.