3 votos

Recepción de valores nulos de la herramienta de muestreo de puntos en QGIS

Estoy utilizando QGIS 2.8.2. Mi proyecto incluye una capa raster de profundidad de flujo que representa la salida de una simulación de inundación (archivo .DAT). También incluye una capa vectorial de puntos que representan los centroides de los edificios dentro de la zona inundada. Estoy tratando de consultar los valores de los píxeles del raster (profundidad del flujo) en las ubicaciones de los centros de los edificios utilizando la herramienta de muestreo de puntos, pero estoy recibiendo valores NULL para cada entrada.

No puedo adjuntar los archivos del proyecto, ya que son datos privados, pero aquí hay un pantallazo:

enter image description here

  • El CRS es el mismo para la capa rasterizada y la capa de puntos (EPSG:2452. Esto también se establece como el CRS para el proyecto, visto en la parte inferior derecha de la capa de trama). Los CRS se han establecido y comprobado mediante el formulario "Selector de sistema de referencia de coordenadas", al que se accede haciendo clic con el botón derecho del ratón en las capas y seleccionando "Establecer CRS de la capa".

  • Las dos capas se superponen definitivamente (véase la captura de pantalla, los puntos negros muestran los centros de los edificios).

  • La capa ráster está correctamente seleccionada mediante la forma de la herramienta de selección de puntos (véase la captura de pantalla. Sólo hay una capa rasterizada y una capa de puntos, por lo que no hay posibilidad de error).

  • El archivo .DAT rasterizado se importa bien (usando Capa -> Añadir capa -> Añadir capa rasterizada), y se visualiza correctamente (los valores de los píxeles son correctos, el renderizado de la banda muestra los valores de los píxeles correctamente -ver captura de pantalla-, el histograma de propiedades de la capa muestra la distribución correcta de los valores de los píxeles).

  • La capa de salida de la herramienta de muestreo de puntos tiene el número correcto de entradas (19.814), por lo que la capa que contiene los puntos de muestreo (centros de construcción) parece utilizarse correctamente.

Todas las consultas anteriores sobre este problema se han resuelto comprobando que el SIR es el mismo para ambas capas. En mi caso son iguales, por lo que mi pregunta es:

Si el CRS es el mismo para la capa de puntos de muestreo y la capa de trama, y las capas se muestran superpuestas, ¿por qué la herramienta de muestreo de puntos devuelve valores NULL para cada punto?

1voto

GriffinHeart Puntos 187

En tu (excelente) descripción del problema, hay una cosa que me preocupa: tu uso del "Set Layer CRS". Esto casi nunca debería usarse, y a menudo causa problemas. Esta opción sólo tiene un propósito (raramente necesario): si se carga una capa que:

  1. Carece de una definición de SRI Y
  2. Sabe de antemano cuál es el SRI correcto

entonces puede adjuntar la definición correcta de CRS a la capa.

Pero, si ambas capas de datos fueran desde el principio en el mismo CRS NO debería ser necesario "Establecer el CRS de la capa".

¿Puede verificar si la reproyección "Sobre la marcha" está activada? Si no es así, actívela y compruebe si las capas siguen solapándose.

\=== Continuación ===

Lo siguiente que yo probaría: hacer el mismo procedimiento usando GRASS. Aquí están los pasos:

  1. Inicie GRASS y cree un nuevo LOCATION basado en el EPSG:2452 CRS
  2. Importe tanto el ráster de profundidad de la inundación como los centroides del edificio a esta UBICACIÓN utilizando los módulos r.in.gdal y v.in.ogr
  3. Obtenga los valores del raster en cada ubicación de punto. En primer lugar, añada una columna que contenga los valores ráster: v.db.addcol map=<centroids> column="flood_depth DOUBLE PRECISION" o bien, si aún no existe una base de datos para los centoides v.db.addtable map=<centroids> column="flood_depth DOUBLE PRECISION"
  4. Ahora corre v.what.rast vector=<centroids> raster=<flood_depth_raster> column=flood_depth

Ahora debería tener una columna de atributos en el vector de centoides con las profundidades de inundación en cada lugar.

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X