¿Por qué cada cálculo de trama que intento vuelve con valores NAN?

8

Estoy tratando de trabajar con rásteres DEM del USGS usando la calculadora ráster. He tratado de cambiar las unidades de metros a pies usando una expresión como:

"ned10m45111h8@1" * 3.28

Sin embargo, esto devuelve valores NAN para cada celda. También lo he intentado sin las comillas alrededor del nombre de la trama como se muestra en este sitio web: http://spatialgalaxy.net/2012/01/25/using-the-qgis-raster-calculator/

Del mismo modo, probé el código de máscara:

("ned10m45111h8@1" <= 1328.96)*"ned10m45111h8@1"

Con y sin comillas y obtener NAN devuelto para todos los valores. ¿Hay alguna configuración predeterminada que necesito agregar, o algún error fundamental que estoy cometiendo? Soy nuevo en QGIS, pero no en GIS y podría hacer esto fácilmente con ARCGIS y el analista espacial.

Hayden
fuente
1
Creo recordar haber leído en alguna parte que este es un problema con la proyección ... La calculadora ráster en QGIS no se adhiere al "proyecto sobre la marcha". Aunque solo esté utilizando un ráster en la operación, quizás intente exportarlo a la misma proyección que su documento QGIS.
Tangnar
1
¿Te importaría compartir tu versión de QGIS? Si es <2.10, entonces vale la pena probar con la versión 2.10. La calculadora ráster tenía un montón de correcciones para 2.10, incluidas muchas que parecen posibles causas de esto ...
ndawson
¿Ha considerado usar GDAL_Calc gdal.org/gdal_calc.html Gdal_Calc.py -A ned10m45111h8 --outfile = result.tif --calc = "A * 3.28084"
Michael Stimson
simplemente haciendo eco de @ndawson para cualquier otra persona atrapada aquí. Estaba ejecutando 2.4 y tuve este problema - la actualización a 2.14 lo resolvió de inmediato
Owen

Respuestas:

5

Esta es una publicación antigua, pero la encontré mientras tenía el mismo problema. Resulta que mi problema era que las extensiones del cálculo se configuraron de manera predeterminada a las extensiones de capa actuales de la primera banda ráster de la lista. Intenté hacer mi cálculo en otra banda de trama, pero las extensiones no se actualizan cuando se selecciona una diferente.

Tuve que hacer clic en la banda de trama en la que estaba haciendo el cálculo y presionar el botón "extensión de capa actual". Voila

iskionrocks
fuente
Esto funcionó para mí también. Primero haga clic en una de las bandas que desea trabajar, luego en la extensión de la capa actual. Si lo haces, la extensión de la capa actual puede generar NaNs
DarkCygnus el
2

Espero ayudar a cualquier otra persona que busque respuestas para tratar con el ráster QGIS calc que devuelve rásteres NaN. Usé la calculadora ráster GDAL en Processing Toolbox en QGIS 3.4.

Probé todas las soluciones que encontré en el intercambio de pila y nada funcionó para mi ráster, como se muestra a continuación: Screenie de información de trama

Actualicé QGIS, configuré CRS desde la capa seleccionada, forcé a QGIS a proyectar la capa en el CRS correcto, cambié el nombre del archivo de mi ráster y cambié los valores de NaN a 0 en el ráster. Ninguna de estas soluciones me permitió realizar ninguna función de cálculo de trama, ni siquiera una simple "twentysix@1" * 1.

El uso de GDAL Raster Calculator resolvió cualquier problema que QGIS Raster Calculator tuviera con el ráster. ¡Espero que la gente pueda ver esto y lo necesite! Me molestó por días.

Screenie de la ubicación de GDAL calc

Mara
fuente
0

Tuve el mismo problema, probé todo lo anterior sin alegría.

Finalmente lo resolví exportando el .shp como .csv, copiando y pegando los campos relevantes en un nuevo .csv que luego utilicé para realizar la interpolación.

Michael_01
fuente