Más

¿Qué sistema de coordenadas es este ráster netCDF?

¿Qué sistema de coordenadas es este ráster netCDF?


Obtuve un archivo ráster netCDF, pero no pude obtener ningún metadato para obtener el nombre del sistema de coordenadas en el que se ha construido el ráster. El ráster en sí no tiene ningún sistema de coordenadas incrustado. Pensé que solo sería un WGS84, y se veía así a primera vista, pero con una mayor investigación en ArcMap, vi que es un sistema bastante poco común. Así es como se muestra:

La trama naranja es una trama normal en WGS84 que he insertado aquí para fines de comparación. El violeta es el ráster con el sistema de coordenadas desconocido. ¿Tienes alguna idea de lo que podría ser esto?

Algunas actualizaciones: Aquí está el ráster netCDF: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 También pude obtener algunos metadatos del proveedor de imágenes:

netcdf sic_average_nclimate {dimensiones: nlon = 361; nlat = 90; nseas = 4; variables: float SIC_Change (nlat, nlon); SIC_Change: Title = "Cambio medio anual medio del conjunto de modelos múltiples cuadriculados en la concentración de hielo 21C-20C"; flotar SIC_Season_Change (nseas, nlat, nlon); SIC_Season_Change: Title = "Cambio medio estacional medio del conjunto multimodelo cuadriculado en la concentración de hielo 21C-20C"; float SIC_Change_STD (nlat, nlon); SIC_Change_STD: Título = "Desviación estándar cuadriculada de múltiples modelos del cambio medio anual en la concentración de hielo 21C-20C"; flotar SIC_Season_Change_STD (nseas, nlat, nlon); SIC_Season_Change_STD: Título = "Desviación estándar cuadriculada de múltiples modelos del cambio medio estacional en la concentración de hielo 21C-20C"; flotar LAT (nlat); LAT: Título = "Latitud"; flotar LON (nlon); LON: Título = "Longitud"; // atributos globales:: Título = "Métricas de concentración de hielo para el subconjunto del modelo como en la Figura 1 del documento NCLIMATE";

Muestran el límite de latitud, pero aparentemente no hay información sobre el sistema de coordenadas.


Puede cambiar los datos a la posición correcta utilizando estos comandos GDAL:

gdal_translate -a_srs EPSG: 4326 -a_ullr 0 0360-90 NETCDF: "sic_average_nclimate.nc": SIC_Change change360.tif gdalwarp -t_srs WGS84 change360.tif change180.tif -wo SOURCE_EXTRA = 1000 --config

(con un poco de ayuda de Frank Warmerdam: Cómo reproyectar el ráster de 0360 a -180 180 con el meridiano de corte 180)

y el resultado se ve así en una proyección estereográfica:


(Esta no es una respuesta, sino más bien una sugerencia de proceso que es demasiado larga para ser un comentario)

Puede georreferenciarlo en QGIS utilizando la herramienta GeoReferencer, ya que parece que puede identificar puntos comunes en ambos rásteres. Esto generaría un script GDAL con un conjunto de puntos GCP y parámetros de transformación. Puede agregar este guión con una muestra de estos puntos a su pregunta, ya que podría explicar cuál es la diferencia entre las dos proyecciones.

En el peor de los casos, tendría una estimación aproximada de los valores basada en el ráster reposicionado que podría ayudarlo con su objetivo.


La esquina inferior izquierda está en las coordenadas (0,0), pero para mí (0,0) es la coordenada de la esquina superior derecha en su caso y no conozco ningún sistema de coordenadas geográficas con su origen en -90.

Mientras trabaja en ArcGIS, puede usar la barra de herramientas de georreferenciación para mover sus datos a la ubicación correcta. Mueva la parte superior derecha a (0,0) y la inferior izquierda a (-360, -90) y esto debería producir una transformación exacta, pero aún tendría un problema con la línea de fecha.

Por lo tanto, puede probar este enlace que le ayudará a crear un sistema de coordenadas geográficas personalizado. Para resumir, en su sistema de referencia, cree un WGS 84 personalizado con el primer meridiano en -180, luego una transformación personalizada usando la rotación de longitud.

La georreferenciación en este nuevo sistema consistiría en establecer el origen en (-180, -90). Un punto de control es suficiente porque solo necesita una traducción. Seleccione con cuidado el píxel inferior izquierdo, haga clic con el botón derecho y agregue los valores XY exactos). Entonces debería verse así

Como observación final, no hay certeza de que el sistema de coordenadas esté basado en WGS 84. Pero es seguro que es un sistema de coordenadas geográficas en grados y la resolución de su ráster es probablemente más tosca que el error potencial debido a una selección de datum incorrecta.


NetCDF

NetCDF (Formulario de datos comunes de red) es un conjunto de bibliotecas de software y formatos de datos autodescriptivos e independientes de la máquina que admiten la creación, el acceso y el intercambio de datos científicos orientados a matrices. La página de inicio del proyecto [1] está alojada en el programa Unidata de la Corporación Universitaria para la Investigación Atmosférica (UCAR). También son la fuente principal de software netCDF, desarrollo de estándares, actualizaciones, etc. El formato es un estándar abierto. NetCDF Classic y el formato de compensación de 64 bits son un estándar internacional del Open Geospatial Consortium. [2]

El proyecto se inició en 1988 y todavía cuenta con el apoyo activo de UCAR. El formato binario netCDF original (lanzado en 1990, ahora conocido como "formato clásico netCDF") todavía se usa ampliamente en todo el mundo y continúa siendo totalmente compatible con todas las versiones de netCDF. La versión 4.0 (lanzada en 2008) permitió el uso del formato de archivo de datos HDF5. La versión 4.1 (2010) agregó soporte para el acceso de clientes de C y Fortran a subconjuntos específicos de datos remotos a través de OPeNDAP. La versión 4.3.0 (2012) agregó un sistema de compilación CMake para compilaciones de Windows. La versión 4.7.0 (2019) agregó compatibilidad para leer objetos de Amazon S3. Se planean lanzamientos adicionales para mejorar el rendimiento, agregar funciones y corregir errores.


Sintaxis de secuencias de comandos

NetCDFConvert2DVariableToArcGISRaster_GeoEco (inputFile, outputRaster, variableName, xLowerLeftCorner, yLowerLeftCorner, cellSize, nodataValue, transpose, mirror, flip, swapHemispheres, coordenadaSystem, proyectadoCoordinateSystem, resumir, transposición, espejo, voltear, intercambiar

En el momento en que se desarrolló esta herramienta, se estaban desarrollando dos versiones principales de netCDF: NetCDF-3, que estaba en la versión 3.6.2, y NetCDF-4, que estaba en la versión 4.0 beta 1. La implementación actual de esta herramienta puede procesar Archivos NetCDF-3 pero no archivos NetCDF-4. La compatibilidad con NetCDF-4 se agregará en un futuro próximo.

Si proporciona un archivo comprimido en un formato de compresión compatible, se descomprimirá automáticamente. Si es un archivo (por ejemplo, .zip o .tar), debe contener exactamente un archivo, que no debe estar en un subdirectorio.

Nombre de una variable en el archivo netCDF. La variable debe tener dos dimensiones y un tipo de datos entero o de punto flotante.

Si no conoce el nombre de la variable, proporcione su mejor estimación. Si una variable no existe con ese nombre, se generará un ValueError. El mensaje de error enumerará los nombres de todas las variables 2D presentes en el archivo de entrada. Seleccione uno de estos nombres y vuelva a intentarlo.

También puede utilizar la herramienta Extraer encabezado NetCDF para volcar el encabezado netCDF en un archivo de texto. Puede encontrar los nombres de las variables buscando en el archivo la sección que comienza con "variables:". Cada variable se enumera como:

dónde tipo de datos es el tipo de datos de la variable (p. ej., "flotante"), nombre es el nombre de la variable, y dimensiones es la lista de dimensiones de la variable. Por ejemplo, el siguiente extracto de un Aviso DT-MADT sobre corrientes geostróficas netCDF enumera dos variables que tienen dos dimensiones, Grid_0001 y Grid_0002:

Esta herramienta extrae la variable netCDF a una cuadrícula ArcInfo ASCII y la convierte en un ráster mediante la herramienta de geoprocesamiento ArcGIS ASCII to Raster. Esa herramienta y el formato ráster de ArcGIS tienen varias limitaciones que restringen los tipos de variables netCDF que se pueden convertir correctamente al formato ráster de ArcGIS.

Las variables NetCDF que utilizan el tipo de datos flotante o doble no deben contener valores "infinito" (INF) o "no un número" (NAN). Se generará un ValueError si se descubren estos valores.

El formato ráster de ArcGIS admite el tipo de datos flotantes de 32 bits, pero no el tipo de datos dobles de 64 bits. Si proporciona una variable netCDF con el tipo de datos doble, se convertirá en un ráster flotante de 32 bits mediante la herramienta de geoprocesamiento ArcGIS ASCII to Raster. El comportamiento de la herramienta en esta situación no está documentado. En ArcGIS 9.1 parece ser:

Los valores en los que el exponente oscila entre -38 y +38 se representan correctamente en el ráster flotante de 32 bits resultante, aunque se pierde algo de precisión debido a la mantisa más pequeña del tipo de datos flotantes de 32 bits.

Los valores donde el exponente es menor que -38 (por ejemplo, -39, -40, etc.) se convierten a 0.

Los valores donde el exponente es mayor que +38 se convierten a -INF o + INF, según el signo del valor (por ejemplo, -5.3083635279597874e-212 aparece como -1. # INF en la GUI de ArcCatalog, mientras que 2.5502286890301497e + 084 aparece como 1. # INF).

La herramienta ArcGIS 9.1 ASCII a ráster también presenta algunas peculiaridades al convertir rásteres enteros:

Para datos enteros de 8 bits, la herramienta creará un ráster de 16 bits si aparece el valor -128, a menos que -128 se designe como valor NODATA. Especificar un valor NODATA diferente, como 0, aún produce un ráster de 16 bits si aparece -128.

De manera similar, para datos enteros de 16 bits, la herramienta creará un ráster de 32 bits si el valor -32768 aparece en el archivo ASCII, a menos que se designe como valor NODATA.

Peor aún, para datos enteros de 32 bits, la herramienta informará un error si aparece el valor -2147483648 a menos que se designe como valor NODATA. Aún más extraño, el valor -2147483647 siempre se traduce a NODATA, pase lo que pase.

Para todos los tipos de datos enteros, la herramienta produce un comportamiento extraño cuando especifica un valor NODATA que no es el valor más pequeño posible para el tipo de datos. Por ejemplo, si el rango de datos va de 0 a 255 y 0 se designa como valor NODATA, la herramienta produce un ráster de salida de 8 bits sin signo. Pero si se designa 1 como valor NODATA, produce un ráster de salida firmado de 16 bits, y ArcCatalog muestra en Propiedades del dataset ráster que el valor NoData es -32768, aunque la herramienta Identificar muestra que las celdas que tenían el valor 1 son en realidad NODATA. Se pueden obtener resultados extraños similares para rásteres enteros de otros tipos de datos, cuando se designa un NODATA que no es el valor más pequeño posible.

Coordenada X de la esquina inferior izquierda del ráster.

La coordenada es para la esquina de la celda inferior izquierda, no el centro de esa celda. Por ejemplo, si el ráster es una proyección geográfica de toda la Tierra, la coordenada de la esquina inferior izquierda sería -180,0, correspondiente a una longitud de 180 grados Oeste.

Coordenada Y de la esquina inferior izquierda del ráster.

La coordenada es para la esquina de la celda inferior izquierda, no el centro de esa celda. Por ejemplo, si el ráster es una proyección geográfica de toda la Tierra, la coordenada de la esquina inferior izquierda sería -90,0, correspondiente a una latitud de 90 grados Sur.

Por ejemplo, si el ráster es una proyección geográfica de toda la Tierra, con 720 columnas y 360 filas, tendría un tamaño de celda de 0.5, correspondiente a 1/2 de un grado geográfico.

El formato de datos subyacente requiere que las celdas sean cuadradas. No es posible especificar un tamaño de celda para cada dimensión.

Valor que indica que una celda no tiene datos.

Si es Verdadero, la imagen se transpondrá (se volteará sobre el eje diagonal) antes de la conversión. Utilice esta opción para corregir una imagen que tenga el eje este / oeste hacia arriba y hacia abajo en lugar de hacia la izquierda y hacia la derecha.

Si es Verdadero, la imagen se volteará sobre el eje vertical antes de la conversión. Utilice esta opción para corregir una imagen que sea la "imagen reflejada" de lo que se supone que es.

Si es Verdadero, la imagen se volteará sobre el eje horizontal antes de la conversión. Utilice esta opción para arreglar una imagen que está al revés.

Si es Verdadero, se intercambiarán los hemisferios este y oeste de la imagen. Utilice esta opción para cambiar la orientación de una imagen global de una orientación de 0 a 360 centrada en el océano Pacífico a una orientación de -180 a +180 centrada en el océano Atlántico, o viceversa.

Sistema de coordenadas a definir para el ráster de salida. Si no se proporciona un valor, el sistema de coordenadas del ráster de salida permanecerá sin definir.

Nuevo sistema de coordenadas para proyectar el ráster de salida.

El ráster solo se puede proyectar en un nuevo sistema de coordenadas si se define la proyección original. Se generará un error si especifica un nuevo sistema de coordenadas sin definir el sistema de coordenadas original.

La herramienta ArcGIS Project Raster se utiliza para realizar la proyección. La documentación de esa herramienta recomienda que también especifique un tamaño de celda para el nuevo sistema de coordenadas.

He notado que para ciertos sistemas de coordenadas, la herramienta ArcGIS 9.2 Project Raster parece recortar el ráster proyectado en una extensión arbitraria que es demasiado pequeña. Por ejemplo, al proyectar una imagen global de clorofila MODIS Aqua de 4 km en coordenadas geográficas a Lambert_Azimuthal_Equal_Area con meridiano central de -60 y latitud de origen de -63, la imagen resultante se recorta para mostrar solo una cuarta parte del planeta. Este problema no ocurre cuando Project Raster se invoca de forma interactiva desde la interfaz de usuario de ArcGIS; solo ocurre cuando la herramienta se invoca mediante programación (el método ProjectRaster_management del geoprocesador). Por lo tanto, es posible que no lo vea cuando use Project Raster usted mismo, pero puede suceder cuando use herramientas MGET que invocan Project Raster como parte de sus operaciones de geoprocesamiento.

Si encuentra este problema, puede solucionarlo de la siguiente manera:

Primero, ejecute esta herramienta sin especificar un nuevo sistema de coordenadas, para obtener el ráster de salida en el sistema de coordenadas original.

En ArcCatalog, use la herramienta Proyectar ráster para proyectar el ráster al nuevo sistema de coordenadas. Verifique que esté presente todo el ráster, que no se haya recortado en una extensión demasiado pequeña.

En ArcCatalog, busque la extensión del ráster proyectado haciendo clic derecho sobre él en el árbol del catálogo, seleccionando Propiedades y desplazándose hacia abajo hasta Extensión.

Ahora, antes de ejecutar la herramienta MGET que proyecta el ráster, establezca la configuración del entorno de Extensión en los valores que buscó. Si está invocando la herramienta MGET de forma interactiva desde ArcCatalog o ArcMap, haga clic en el botón Entornos en el cuadro de diálogo de la herramienta, abra Configuración general, cambie el menú desplegable Extensión a "Como se especifica a continuación" y escriba los valores que buscó. Si lo está invocando desde un modelo de geoprocesamiento, haga clic con el botón derecho en la herramienta en el modelo, seleccione Hacer variable, Desde el entorno, Configuración general, Extensión. Esto colocará Extent como una variable en su modelo, adjunta a la herramienta MGET. Abra la variable Extensión, cámbiela a "Como se especifica a continuación" y escriba los valores que buscó. Si está invocando la herramienta MGET mediante programación, debe establecer la propiedad Extensión del geoprocesador en los valores que buscó. Consulte la documentación de ArcGIS para obtener más información sobre esto y la configuración del entorno en general.

Ejecute la herramienta MGET. La extensión del ráster de salida debería tener ahora el tamaño adecuado.

Un método de transformación utilizado para convertir entre el sistema de coordenadas original y el nuevo sistema de coordenadas.

Este parámetro es una nueva opción introducida por ArcGIS 9.2. Debe tener ArcGIS 9.2 para usar este parámetro.

Este parámetro solo es necesario cuando especifica que el ráster debe proyectarse a un nuevo sistema de coordenadas y ese nuevo sistema usa un datum diferente al sistema de coordenadas original, o existe alguna otra diferencia entre los dos sistemas de coordenadas que requiere una transformación. Para determinar si se necesita una transformación, recomiendo el siguiente procedimiento:

Primero, ejecute esta herramienta sin especificar un nuevo sistema de coordenadas, para obtener el ráster de salida en el sistema de coordenadas original.

A continuación, use la herramienta ArcGIS 9.2 Project Raster en el ráster de salida para proyectarlo en el sistema de coordenadas deseado. Si se necesita una transformación geográfica, esa herramienta le pedirá una. Anote el nombre exacto de la transformación que utilizó.

Finalmente, si se necesitaba una transformación, escriba el nombre exacto en esta herramienta, vuelva a ejecutarla y verifique que el ráster de salida se proyectó como deseaba.

El algoritmo de remuestreo que se utilizará para proyectar el ráster original en un nuevo sistema de coordenadas. La herramienta ArcGIS Project Raster se utiliza para realizar la proyección y acepta los siguientes valores:

MÁS CERCANO: interpolación del vecino más cercano

BILINEAR - interpolación bilineal

Debe especificar uno de estos algoritmos para proyectar a un nuevo sistema de coordenadas. Se generará un error si especifica un nuevo sistema de coordenadas sin seleccionar un algoritmo.

El tamaño de celda del sistema de coordenadas proyectadas. Aunque este parámetro es opcional, para recibir los mejores resultados, la documentación de ArcGIS recomienda que siempre lo especifique al proyectar a un nuevo sistema de coordenadas.

Las coordenadas xey (en el espacio de salida) utilizadas para la alineación de píxeles.

Este parámetro es una nueva opción introducida por ArcGIS 9.2. Debe tener ArcGIS 9.2 para usar este parámetro. Se ignora si no especifica que el ráster debe proyectarse en un nuevo sistema de coordenadas.

Rectángulo al que se debe recortar el ráster.

Si se especificó un sistema de coordenadas proyectadas, el recorte se realiza después de la proyección y las coordenadas del rectángulo deben especificarse en el nuevo sistema de coordenadas. Si no se especificó un sistema de coordenadas proyectadas, las coordenadas deben especificarse en el sistema de coordenadas original.

La herramienta ArcGIS Clip se utiliza para realizar el clip. El rectángulo de recorte debe pasarse a esta herramienta como una cadena de cuatro números separados por espacios. La interfaz de usuario de ArcGIS formatea automáticamente la cadena correctamente al invocar esta herramienta desde la interfaz de usuario de ArcGIS, no necesita preocuparse por el formato. Pero al invocarlo mediante programación, tenga cuidado de proporcionar una cadena con el formato adecuado. Los números están ordenados IZQUIERDA, INFERIOR, DERECHA, SUPERIOR. Por ejemplo, si el ráster está en un sistema de coordenadas geográficas, se puede recortar a 10 W, 15 S, 20 E y 25 N con la cadena:

Pueden proporcionarse números enteros o decimales.

Expresión de álgebra de mapas para ejecutar en el ráster de salida.

ADVERTENCIA: El Generador de modelos de geoprocesamiento de ArcGIS puede eliminar de forma aleatoria y silenciosa el valor de este parámetro. Este es un error en ArcGIS. Antes de ejecutar un modelo que ha guardado, abra esta herramienta y valide que el valor del parámetro todavía existe.

La expresión se ejecuta después de proyectar y recortar el ráster convertido (si se especifican esas opciones). Utilice la cadena inputRaster que distingue entre mayúsculas y minúsculas para representar el ráster sobre el que ahora desea realizar el álgebra de mapas. Por ejemplo, para convertir el ráster en un ráster entero y agregar 1 a todas las celdas, use esta expresión:

La cadena inputRaster distingue entre mayúsculas y minúsculas. Antes de ejecutar la expresión de álgebra de mapas, la cadena se reemplaza con la ruta a un ráster temporal que representa el ráster de salida que se está generando. La expresión final debe tener menos de 4000 caracteres o ArcGIS informará un error.

La herramienta Álgebra de mapas de salida única de ArcGIS se utiliza para ejecutar la expresión de álgebra de mapas. Debe tener una licencia para la extensión ArcGIS Spatial Analyst para realizar álgebra de mapas.

La sintaxis del álgebra de mapas puede ser muy delicada. Aquí hay algunos consejos que le ayudarán a tener éxito con esta herramienta:

Antes de usar esta herramienta, construya y pruebe su expresión de álgebra de mapas con la herramienta Álgebra de mapas de salida única de ArcGIS. Luego pegue la expresión en esta herramienta y edítela para usar la variable inputRaster en lugar del valor de prueba que usó con Single Output Map Algebra.

Si desarrolla su expresión directamente en esta herramienta, comience con una expresión muy simple. Verifique que funcione correctamente, agregue un poco y verifique nuevamente. Repita este proceso hasta que haya construido la expresión completa.

Separe siempre los operadores matemáticos de las rutas ráster utilizando espacios. En el ejemplo anterior, el operador / contiene un espacio a cada lado. Siga este patrón. En algunas circunstancias, ArcGIS fallará al procesar expresiones de álgebra ráster que no separan las rutas ráster de los operadores que utilizan espacios. El mensaje de error informado generalmente no indica que este sea el problema, y ​​rastrearlo puede ser muy frustrante.

Si es True, se construirán pirámides para el ráster de salida, lo que mejorará su velocidad de visualización en la interfaz de usuario de ArcGIS. Este es el último paso realizado en el procesamiento posterior a la conversión.


Parámetros

Las capas de estadísticas geográficas 3D que se exportarán al archivo Output netCDF. Si se proporciona más de una capa, la salida será un archivo netCDF multivariante.

El archivo netCDF de salida que contiene los valores exportados de las capas de estadísticas geográficas 3D de entrada.

Especifica las ubicaciones para exportar desde las capas de estadísticas geográficas 3D de entrada. Puede exportar a puntos de cuadrícula 3D o proporcionar características de puntos 3D personalizadas para representar las ubicaciones de exportación. Si elige puntos cuadriculados 3D, debe proporcionar valores para los parámetros Espaciado X, Espaciado Y y Espaciado de elevación que representan la distancia entre cada punto cuadriculado en todas las dimensiones. Si elige Puntos 3D personalizados, debe proporcionar entidades de puntos 3D en el parámetro de ubicaciones de puntos 3D que representan las ubicaciones para exportar.

  • Puntos cuadriculados 3D: las ubicaciones de predicción son puntos cuadriculados 3D. Este es el predeterminado.
  • Puntos 3D personalizados: las ubicaciones de predicción se definen mediante características de puntos 3D personalizados.

El espacio entre cada punto de la cuadrícula en la dimensión x. El valor predeterminado crea 40 puntos a lo largo de la extensión x de salida.

El espacio entre cada punto de la cuadrícula en la dimensión y. El valor predeterminado crea 40 puntos a lo largo de la extensión y de salida.

El espaciado entre cada punto de la cuadrícula en la dimensión de elevación (z). El valor predeterminado crea 40 puntos a lo largo de la extensión z de salida.

Las entidades de puntos 3D que representan ubicaciones para exportar. Las entidades de puntos deben tener sus elevaciones almacenadas en el atributo de geometría Shape.Z.

Especifica los tipos de salida para las capas de estadísticas geográficas 3D de entrada. Puede especificar uno o más tipos de salida para cada una de las capas o puede aplicar un tipo de salida a todas las capas de estadísticas geográficas de entrada. De forma predeterminada, se exportarán las predicciones para todas las capas.

Para exportar otros tipos de salida, especifique la capa a exportar (o elija Todo para especificar todas las capas) en la primera entrada de la tabla de valores. Especifique el tipo de salida en la segunda entrada de la tabla de valores. Si elige Probabilidad o Cuantil como tipo de salida, especifique el valor de umbral (para probabilidad) o el valor de cuantil (para cuantil) en la tercera entrada de la tabla de valores. Si elige Predicción o Error estándar de predicción como tipo de salida, puede dejar la tercera entrada en la tabla de valores vacía.

Las entidades poligonales que representan el área de estudio. Solo los puntos que están dentro del área de estudio se guardan en el archivo netCDF de salida. Cuando se visualiza como una capa de vóxeles, solo los vóxeles dentro del área de estudio se mostrarán en la escena. Se determina que los puntos están dentro o fuera del área de estudio utilizando solo sus coordenadas x e y.

Las capas de estadísticas geográficas 3D que se exportarán al archivo Output netCDF. Si se proporciona más de una capa, la salida será un archivo netCDF multivariante.

El archivo netCDF de salida que contiene los valores exportados de las capas de estadísticas geográficas 3D de entrada.

Especifica las ubicaciones para exportar desde las capas de estadísticas geográficas 3D de entrada. Puede exportar a puntos de cuadrícula 3D o proporcionar características de puntos 3D personalizadas para representar las ubicaciones de exportación. Si elige puntos cuadriculados 3D, debe proporcionar valores para los parámetros Espaciado X, Espaciado Y y Espaciado de elevación que representan la distancia entre cada punto cuadriculado en todas las dimensiones. Si elige Puntos 3D personalizados, debe proporcionar entidades de puntos 3D en el parámetro de ubicaciones de puntos 3D que representan las ubicaciones para exportar.

  • 3D_GRIDDED_POINTS: las ubicaciones de predicción son puntos cuadriculados en 3D. Este es el predeterminado.
  • CUSTOM_3D_POINTS: las ubicaciones de predicción se definen mediante características de puntos 3D personalizadas.

El espacio entre cada punto de la cuadrícula en la dimensión x. El valor predeterminado crea 40 puntos a lo largo de la extensión x de salida.

El espacio entre cada punto de la cuadrícula en la dimensión y. El valor predeterminado crea 40 puntos a lo largo de la extensión y de salida.

El espacio entre cada punto de la cuadrícula en la dimensión de elevación (z). El valor predeterminado crea 40 puntos a lo largo de la extensión z de salida.

Las entidades de puntos 3D que representan ubicaciones para exportar. Las entidades de puntos deben tener sus elevaciones almacenadas en el atributo de geometría Shape.Z.

Especifica los tipos de salida para las capas de estadísticas geográficas 3D de entrada. Puede especificar uno o más tipos de salida para cada una de las capas o puede aplicar un tipo de salida a todas las capas de estadísticas geográficas de entrada. De forma predeterminada, se exportarán las predicciones para todas las capas.

Para exportar otros tipos de salida, especifique la capa a exportar (o elija Todo para especificar todas las capas) en la primera entrada de la tabla de valores. Especifique el tipo de salida en la segunda entrada de la tabla de valores. Si elige Probabilidad o Cuantil como tipo de salida, especifique el valor de umbral (para probabilidad) o el valor de cuantil (para cuantil) en la tercera entrada de la tabla de valores. Si elige Predicción o Error estándar de predicción como tipo de salida, puede dejar la tercera entrada en la tabla de valores vacía.

Las entidades poligonales que representan el área de estudio. Solo los puntos que están dentro del área de estudio se guardan en el archivo netCDF de salida. Cuando se visualiza como una capa de vóxeles, solo los vóxeles dentro del área de estudio se mostrarán en la escena. Se determina que los puntos están dentro o fuera del área de estudio utilizando solo sus coordenadas x e y.

Muestra de código

Interpole puntos 3D dos veces y convierta las salidas en un archivo netCDF multivariado.

Interpole puntos 3D dos veces y convierta las salidas en un archivo netCDF multivariado. Exportar a puntos 3D personalizados y cuadriculados.


¿Cómo convierto un archivo netcdf proyectado a ráster para trazar en R?

Tengo un archivo netcdf de rayos que contiene variables para x, y y Lat Lon. El sistema de referencia de coordenadas para los datos xy es albers_conical_equal_area y hay una serie de variables adicionales falso este, falso norte, etc. la proyección. Me gustaría convertir el archivo a ráster y debo completar un parámetro crs. Me gustaría que el ráster se muestre en el formato proyectado, pero no sé cómo formatear correctamente la entrada crs.

¿Es esto siquiera posible (el ráster solo funciona con datos de latitud?)

Si es así, ¿cómo convierto la información en el archivo netcdf para que funcione con el parámetro crs en el comando ráster?

Intenté obtener los atributos del archivo netcdf pero no está claro cómo usarlos.

He usado el siguiente código para obtener los atributos de la variable crs pero.

No estoy seguro de cómo aplicar esta información al comando ráster

El archivo que estoy usando es un solo día de descargas eléctricas en los EE. UU. Y preveo importar muchos de estos archivos para agregar ataques en un área de interés durante todo un año, así que me gustaría resolver este problema del proyecto. en este archivo primero.


Visualización de datos temporales de NetCDF con QGIS (escaparate global de temperatura mensual)

Ya agregamos los datos NetCDF en QGIS y vemos la información de atributos de los datos. Ahora visualicemos o animemos el cambio temporal.

Desafortunadamente, después de jugar un poco con esta herramienta, llegué a la conclusión de que la herramienta de controlador temporal no se puede usar directamente para representar una banda con una fecha respectiva. Es diferente de los datos vectoriales que tienen un atributo de tiempo para cada característica. Después de buscar una solución, finalmente encontré el siguiente código de Python que se puede usar para superar este problema. Cambié un poco el código agregando entrada para la fecha / hora de inicio y finalización. El código Python original se puede encontrar aquí.

Volver al escritorio de QGIS. Abierto Consola de Python. Copie el código y péguelo en un editor como en la figura 6. Cambie la fecha / hora de inicio y finalización con los datos NetCDF correspondientes. Puede verificar la fecha / hora de inicio y finalización en la información de propiedades de atributo / capa. La información de tiempo en la información de propiedades está en una marca de tiempo Unix como esta: 1297320 (horas desde 1800: 01: 01 00:00:00), por lo que debemos convertirlo en formato de fecha y hora (aaaa: mm: dd hh: mm: ss) si queremos usarlo. Pero para simplificar la vida, en este caso, miré la página de descarga del archivo y descubrí que el rango de datos es desde 1948 hasta febrero de 2021.


Creación de presentaciones de mapas con datos en un sistema de referencia de coordenadas proyectadas

Este ejemplo ilustra cómo importar y mostrar datos geográficos que contienen coordenadas en un sistema de referencia de coordenadas proyectadas.

En particular, este ejemplo ilustra cómo

Importar conjuntos de datos vectoriales y ráster específicos

Cree visualizaciones de mapas para visualizar los datos

Mostrar múltiples conjuntos de datos en una visualización de mapa

Muestre múltiples conjuntos de datos con coordenadas en sistemas de referencia de coordenadas geográficas y proyectadas en una sola pantalla de mapa

Ejemplo 1: Importar datos ráster en un sistema de referencia de coordenadas proyectadas

Los datos ráster geográficos que contienen coordenadas en un sistema de referencia de coordenadas proyectadas se pueden almacenar en una variedad de formatos diferentes, incluidos formatos de archivo estándar como GeoTIFF, Spatial Data Transfer Standard (SDTS), NetCDF, HDF4 o HDF5. Este ejemplo ilustra la importación de datos desde un archivo GeoTIFF. Los datos del archivo contienen coordenadas en el sistema de referencia de coordenadas del mapa proyectado Sistema de coordenadas de la zona continental del plano del estado de Massachusetts .

Las coordenadas de la imagen en el archivo GeoTIFF, boston.tif, están en un sistema de referencia de coordenadas proyectadas. Puede determinarlo utilizando la función geotiffinfo y examinar los valores de los campos PCS y Proyección.

La unidad de longitud de las coordenadas se define mediante el campo UOMLength en la estructura de información.

Para importar la imagen y el objeto de referencia espacial, use readgeoraster.

Ejemplo 2: Mostrar datos ráster en un sistema de referencia de coordenadas proyectadas

Puede mostrar la imagen en ejes normales de MATLAB usando mapshow, que muestra la imagen y establece los límites de los ejes a los límites definidos por el objeto de referencia, R. Las coordenadas, como se mencionó anteriormente, están al pie de la encuesta de EE. UU. Y son relativas a un origen al suroeste del mapa, por lo que los números son grandes. Las coordenadas siempre son positivas dentro de la zona.

Ejemplo 3: Importar datos vectoriales en un sistema de referencia de coordenadas proyectadas

Los datos vectoriales geográficos que contienen coordenadas en un sistema de referencia de coordenadas proyectadas se pueden almacenar en shapefiles. Este ejemplo ilustra cómo importar datos vectoriales en un sistema de referencia de coordenadas proyectadas desde el shapefile, boston_roads.shp.

Importe datos de líneas vectoriales desde el archivo boston_roads.shp.

Para obtener información sobre el sistema de referencia de coordenadas proyectadas, primero devuelva información sobre el shapefile como una estructura. Luego, consulte el campo CoordinateReferenceSystem.

Ejemplo 4: Mostrar datos vectoriales y ráster en un sistema de referencia de coordenadas proyectadas

Los datos vectoriales y ráster de este ejemplo están en el mismo sistema de referencia de coordenadas proyectadas. Sin embargo, los datos vectoriales están en unidades de longitud en metros, mientras que los datos ráster están en unidades de longitud del pie topográfico. Convierta los datos ráster a unidades de longitud en metros y muestre los datos en los mismos ejes.

Convierta las coordenadas de la imagen ráster de unidades de pies topográficos estadounidenses a metros.

Display the raster image and vector data using mapshow .

Example 5: Display Data in both Geographic and Projected Coordinate Reference Systems

You may have geographic data whose coordinates are in latitude and longitude and other data whose coordinates are in a projected coordinate reference system. You can display these data sets in the same map display. This example illustrates how to display data in a geographic coordinate reference system (latitude and longitude) with data in a projected map coordinate reference system (Massachusetts State Plane Mainland Zone coordinate system).

Read a raster image with a worldfile whose coordinates are in latitude and longitude. Use imread to read the image and worldfileread to read the worldfile and construct a spatial referencing object.

To display the overview image and the GeoTIFF image in the same map display, you need to create a map display using a Mapping Toolbox™ projection structure containing the projection information for the data in the projected coordinate reference system, Massachusetts State Plane Mainland Zone coordinate system . To make a map display in this system, you can use the projection information contained in the GeoTIFF file. Use the geotiff2mstruct function to construct a Mapping Toolbox™ projection structure, from the contents of the GeoTIFF information structure. The geotiff2mstruct function returns a projection in units of meters. Use the projection structure to define the projection parameters for the map display.

Use the latitude and longitude limits of the Boston overview image.

Create a map display by using the projection information stored in the map projection structure and set the map latitude and longitude limits. Display the geographic data in the map axes. geoshow projects the latitude and longitude coordinates.

Since the coordinates of the GeoTIFF image are in a projected coordinate reference system, use mapshow to overlay the more detailed Boston image onto the display. Plot the boundaries of the Boston image in red.

Zoom to the geographic region of the GeoTIFF image by setting the axes limits to the limits of the Boston image and add a small buffer. Note that the buffer size ( delta ) is expressed in meters.

You can overlay the road vectors onto the map display. Use a symbol specification to give each class of road its own color.

You can also overlay data from a GPS stored in a GPX file. Import point geographic vector data from the boston_placenames.gpx file included with the Mapping Toolbox™ software. The file contains latitude and longitude coordinates of geographic point features in part of Boston, Massachusetts, USA. Use gpxread to read the GPX file and return a geopoint vector .

Overlay the placenames onto the map and increase the marker size, change the markers to circles and set their edge and face colors to yellow.

Data Set Information

The files boston.tif and boston_ovr.jpg include materials copyrighted by GeoEye, all rights reserved. GeoEye was merged into the DigitalGlobe corporation on January 29th, 2013. For more information about the data sets, use the commands type boston.txt and type boston_ovr.txt .

The files boston_roads.shp and boston_placenames.gpx are from the Bureau of Geographic Information (MassGIS), Commonwealth of Massachusetts, Executive Office of Technology and Security Services. For more information about the data sets, use the commands type boston_roads.txt and type boston_placenames_gpx.txt .


GeoTIFF file (*.tif), Arc/Info ASCII file (*.asc) (Raster data)¶

When you import GeoTIFF file (*.tif), Arc/Info ASCII file (*.asc) to import raster data, if information about coordinate system is not included in the file, [Warning] dialog ( Figure 62 ) is shown, and [Select Coordinate System] dialog ( Figure 63 ) is shown.

When importing aGeoTIFF file, the coordinate system information included in the file is imported. When importing an Arc/Info ASCII file, the coordinate system information included in *.prj file is imported, if it exists.

Figure 63 [Select Coordinate System] dialog


Creating Map Displays with Data in Projected Coordinate Reference System

This example illustrates how to import and display geographic data that contain coordinates in a projected coordinate reference system.

In particular, this example illustrates how to

Import specific raster and vector data sets

Create map displays for visualizing the data

Display multiple data sets in a map display

Display multiple data sets with coordinates in geographic and projected coordinate reference systems in a single map display

Example 1: Import Raster Data in Projected Coordinate Reference System

Geographic raster data that contains coordinates in a projected coordinate reference system can be stored in a variety of different formats, including standard file formats such as GeoTIFF, Spatial Data Transfer Standard (SDTS), NetCDF, HDF4, or HDF5. This example illustrates importing data from a GeoTIFF file. The data in the file contains coordinates in the projected map coordinate reference system Massachusetts State Plane Mainland Zone coordinate system .

The coordinates of the image in the GeoTIFF file, boston.tif , are in a projected coordinate reference system. You can determine that by using the geotiffinfo function and examine the PCS and Projection field values.

The length unit of the coordinates are defined by the UOMLength field in the info structure.

To import the image and the spatial referencing object, use readgeoraster .

Example 2: Display Raster Data in Projected Coordinate Reference System

You can display the image on a regular MATLAB axes using mapshow , which displays the image and sets the axes limits to the limits defined by the referencing object, R . The coordinates, as mentioned above, are in US survey foot and are relative to an origin to the southwest of the map, which is why the numbers are large. The coordinates are always positive within the zone.

Example 3: Import Vector Data in Projected Coordinate Reference System

Geographic vector data that contains coordinates in a projected coordinate reference system can be stored in shapefiles. This example illustrates how to import vector data in a projected coordinate reference system from the shapefile, boston_roads.shp .

Import vector line data from the boston_roads.shp file.

To get information about the projected coordinate reference system, first return information about the shapefile as a structure. Then, query the CoordinateReferenceSystem field.

Example 4: Display Vector and Raster Data in Projected Coordinate Reference System

The vector and raster data in this example are in the same projected coordinate reference system. However, the vector data is in length units of meter, while the raster data is in length unit of survey foot. Convert the raster data to length units of meter and display the data on the same axes.

Convert the coordinates of the raster image from units of US survey foot to meter.

Display the raster image and vector data using mapshow .

Example 5: Display Data in both Geographic and Projected Coordinate Reference Systems

You may have geographic data whose coordinates are in latitude and longitude and other data whose coordinates are in a projected coordinate reference system. You can display these data sets in the same map display. This example illustrates how to display data in a geographic coordinate reference system (latitude and longitude) with data in a projected map coordinate reference system (Massachusetts State Plane Mainland Zone coordinate system).

Read a raster image with a worldfile whose coordinates are in latitude and longitude. Use imread to read the image and worldfileread to read the worldfile and construct a spatial referencing object.

To display the overview image and the GeoTIFF image in the same map display, you need to create a map display using a Mapping Toolbox™ projection structure containing the projection information for the data in the projected coordinate reference system, Massachusetts State Plane Mainland Zone coordinate system . To make a map display in this system, you can use the projection information contained in the GeoTIFF file. Use the geotiff2mstruct function to construct a Mapping Toolbox™ projection structure, from the contents of the GeoTIFF information structure. The geotiff2mstruct function returns a projection in units of meters. Use the projection structure to define the projection parameters for the map display.

Use the latitude and longitude limits of the Boston overview image.

Create a map display by using the projection information stored in the map projection structure and set the map latitude and longitude limits. Display the geographic data in the map axes. geoshow projects the latitude and longitude coordinates.

Since the coordinates of the GeoTIFF image are in a projected coordinate reference system, use mapshow to overlay the more detailed Boston image onto the display. Plot the boundaries of the Boston image in red.

Zoom to the geographic region of the GeoTIFF image by setting the axes limits to the limits of the Boston image and add a small buffer. Note that the buffer size ( delta ) is expressed in meters.

You can overlay the road vectors onto the map display. Use a symbol specification to give each class of road its own color.

You can also overlay data from a GPS stored in a GPX file. Import point geographic vector data from the boston_placenames.gpx file included with the Mapping Toolbox™ software. The file contains latitude and longitude coordinates of geographic point features in part of Boston, Massachusetts, USA. Use gpxread to read the GPX file and return a geopoint vector .

Overlay the placenames onto the map and increase the marker size, change the markers to circles and set their edge and face colors to yellow.

Data Set Information

The files boston.tif and boston_ovr.jpg include materials copyrighted by GeoEye, all rights reserved. GeoEye was merged into the DigitalGlobe corporation on January 29th, 2013. For more information about the data sets, use the commands type boston.txt and type boston_ovr.txt .

The files boston_roads.shp and boston_placenames.gpx are from the Bureau of Geographic Information (MassGIS), Commonwealth of Massachusetts, Executive Office of Technology and Security Services. For more information about the data sets, use the commands type boston_roads.txt and type boston_placenames_gpx.txt .


FME considers a dataset to be a folder name. The feature type of each dataset is the filename.

FME will write a netCDF file with a single data variable for each feature. Coordinate variables and dimensions can be specified through attributes.

By default, FME will write coordinate variables named y and x to capture the georeferencing of the raster. For geographic coordinate systems, these variables will instead be named lat and lon . Multi-band rasters will get an additional band coordinate variable which stores a band index.

Instead of these default coordinate variables, user-defined coordinate variables may be defined using the netcdf_coordinate_variable<> list attribute. For example, the following attributes could be specified to create a data variable that depends on four coordinate variables:

When user-defined coordinate variables are specified, values must also be specified for each band for each "extra" (non-x/y) coordinate variable. For example, if the year variable defined above has two possible values <2000, 2001>and the month variable has 12 possible values <1, 2, . 12>, then the raster must have 24 bands, and the following attributes should be defined:

Note: The last two coordinate variables are always assumed to be y and x, and their values are based on the georeferencing of the raster. As a result, values do not need to be specified for the last two coordinate variables using the netcdf_band<>.coordinate_variable<> list attribute.


Frequently asked questions

The Add Multidimensional Voxel Layer dialog box does not recognize my netCDF file. Why?

Unsupported netCDF files do not display variable information in the Add Multidimensional Voxel Layer dialog box and cannot be added to a local scene. A warning button appears on the dialog box. Click it to see a report on any issues with the netCDF file. Ensure that the netCDF file is volumetric and contains the necessary x,y,z dimensions. See supported voxel formats for more information. The Python netcdf4 library or ArcPy NetCDFFileProperties class can validate dimension and variable information for a netCDF file.

Which coordinate systems do voxel layers support?

Voxel layers support both geographic and projected coordinate systems in a local scene. Ellipsoidal and gravity-based vertical coordinate systems are also supported.

How are coordinate systems read from a netCDF file?

Coordinate variables in a netCDF file determine which variables represent the x,y,z information. The coordinate system of the data created from a netCDF variable is determined by the units of the coordinate variables and whether the grid_mapping attribute exists.

If the units of the coordinate variables specified as the x- and y-dimensions are the units of longitude and latitude, respectively, the data is in a geographic coordinate system (GCS) and is automatically set as WGS 1984.

If the standard_name of the coordinate variables specified as the x- and y-dimensions are projection_x_coordinate and projection_y_coordinate, respectively, and the variable has a grid_mapping attribute, the data is in a projected coordinate system (PCS). The projection name, type, and other parameters are defined by a grid_mapping variable.

Esri software that creates netCDF files from space-time cubes or geostatistical analysis layers write out a global attribute esri_pe_string . The voxel layer reads the well known text (WKT) from this global attribute to set the coordinate system.

My NetCDF file has the wrong coordinate system. What should I do?

Place a projection file ( .prj ) with the same name as the netCDF file in the same directory to override the coordinate system of the netCDF file. See Save a coordinate system as a projection file to learn how to create a projection file.

I've received a notification that my voxel layer will not draw.

Both the horizontal and vertical coordinate systems of the voxel layer must match the scene's coordinate system. Open the scene properties to confirm both the horizontal and vertical coordinate systems are set properly.

How do I repair a broken voxel layer?

To repair a voxel layer, click the red exclamation point in the Contents pane. The netCDF file used to repair the broken layer must be the same netCDF file that was used to create the original voxel layer. If the netCDF files do not match, an error message appears.

How do I share a voxel layer?

Voxel layers are supported in ArcGIS Pro . Voxel layers can be shared as a layer package, map package, or project package.

Can I add multiple voxel layers to my local scene?

Si. Ensure that the voxel layers do not overlap in the same area, as only one voxel layer displays at a time. Also ensure that your graphics card can handle displaying multiple voxel layers. It is recommended that you have at least 4 GB of dedicated graphics memory. See ArcGIS Pro system requirements for more information.

How do voxel layers support time?

Voxel layers support netCDF files with X,Y,T or X,Y,Z,T dimensions. The time variable must have a units attribute and the time values must be sorted in ascending order. The following formats are supported: YYYY:MM:DD , hh:mm:ss , YYYY:MM:DD , and YYYY:MM:DD hh:mm:ss hh:mm , where the last hh:mm is a time zone. Use the Python netcdf4 library if the time values need to be sorted in ascending order.

What is the .vxc1 that is created next to my netCDF file?

The .vxc1 is created anytime a user adds a voxel layer to a local scene with the Optimize for performance option enabled. This file can be used by any ArcGIS Pro application that is displaying the netCDF file as a voxel layer. The .vxc1 file increases in size as a user displays more variables and creates more isosurfaces and sections. Disable this functionality by unchecking Optimize for performance on the Display page of the voxel layer properties. The .vxc1 file can be manually deleted in the file system.