Este es el comando r.in.gdalgrass que se puede ejecutar en el proveedor de alojamiento gratuito de OnWorks utilizando una de nuestras múltiples estaciones de trabajo en línea gratuitas, como Ubuntu Online, Fedora Online, emulador en línea de Windows o emulador en línea de MAC OS.
PROGRAMA:
NOMBRE
r.in.gdal - Importa datos ráster en un mapa ráster GRASS usando la biblioteca GDAL.
PALABRAS CLAVE
ráster, importación
SINOPSIS
r.in.gdal
r.in.gdal --ayuda
r.in.gdal [-ojeflkcp] Las opciones de entrada=nombre salida=nombre [Número=entero[,entero, ...]]
[memoria=entero] [dirigidos=nombre ] [título=redactar] [compensar=entero] [località=nombre ]
[mesa=presentar] [-exagerar] [-ayuda] [-verboso] [-tranquilo] [-ui]
Banderas
-o
Anular la proyección (usar la proyección de la ubicación)
-j
Realice solo la comprobación de proyección y salga
-e
Extienda las extensiones de la región en función de un nuevo conjunto de datos
También actualiza la región predeterminada si está en el conjunto de mapas PERMANENTE
-f
Enumere los formatos admitidos y salga
-l
Forzar los mapas Lat / Lon para que se ajusten a las coordenadas geográficas (90N, S; 180E, W)
-k
Mantenga los números de banda en lugar de usar nombres de colores de banda
-c
Cree la ubicación especificada por el parámetro "ubicación" y salga. No importe el
archivo de trama.
-p
Imprime el número de bandas y sal
--Sobrescribir
Permitir que los archivos de salida sobrescriban los archivos existentes
--ayuda
Resumen de uso de impresión
--verboso
Salida del módulo detallado
--tranquilo
Salida de módulo silencioso
--ui
Forzar el inicio del cuadro de diálogo GUI
parámetros:
Las opciones de entrada=nombre [requerido]
Nombre del archivo ráster que se va a importar
salida=nombre [requerido]
Nombre del mapa ráster de salida
Número=entero [, entero, ...]
Banda (s) para seleccionar (por defecto son todas las bandas)
memoria=entero
Memoria máxima a utilizar (en MB)
Tamaño de caché para filas ráster
Opciones: 0 - 2047
Por defecto: 300
dirigidos=nombre
Nombre de la ubicación de destino de GCP
Nombre de la ubicación para crear o leer la proyección para la transformación de GCP
título=redactar
Título del mapa ráster resultante
compensar=entero
El desplazamiento se agregará al número de banda mientras se crea el nombre del mapa ráster de salida
Por defecto: 0
località=nombre
Nombre de la nueva ubicación para crear
mesa=presentar
Prefijo de archivo para tablas de atributos ráster
El número de banda y ".csv" se agregarán al prefijo del archivo.
DESCRIPCIÓN
r.in.gdal permite a un usuario crear una capa de mapa ráster (binario) GRASS, o un grupo de imágenes,
desde cualquier formato de mapa ráster compatible con GDAL, con un título opcional. El archivo importado puede
también se puede utilizar opcionalmente para crear una nueva ubicación.
Gdal apoyadas raster formatos
Los detalles completos sobre los formatos compatibles con GDAL están disponibles en:
http://www.gdal.org/formats_list.html
Formatos seleccionados de los más de 80 formatos compatibles:
Nombre de formato largo Creación de código Georreferenciación Tamaño máximo de archivo
Arc / Info ASCII Grid AAIGrid Sí Sí Sin límites
Arco / Info Cuadrícula binaria AIG No Sí -
AIRSAR Polarimétrico AIRSAR No No -
Mapa de bits independiente del dispositivo de Microsoft Windows (.bmp) BMP Sí Sí 4GiB
Formato de carta náutica BSB (.kap) BSB No Sí -
Formato de terreno binario VTP (.bt) BT Sí Sí -
CEOS (Spot por ejemplo) CEOS No No -
DOQ (.doq) DOQ1 del USGS de primera generación No Sí -
Nueva etiqueta USGS DOQ (.doq) DOQ2 No Sí -
Datos (.dt0, .dt1) DTED No Sí -
ERMapper Comprimidos Wavelets (.ecw) ECW Sí Sí
ESRI .hdr Etiquetado EHdr No Sí -
ENVI .hdr Raster etiquetado ENVI Sí Sí Sin límites
Producto de imagen Envisat (.n1) Envisat No No -
EOSAT FAST Formato FAST No Sí -
AJUSTA (.ajusta) AJUSTA Sí No
Formato de intercambio de gráficos (.gif) GIF Sí No 2 GB
Cuadrícula binaria Arc / Info (.adf) GIO Sí Sí
GRASS Rásteres GRASS No Sí -
TIFF / GeoTIFF (.tif) GTiff Sí Sí 4GiB
Formato de datos jerárquico versión 4 (HDF4) HDF4 Sí Sí 2GiB
Erdas Imagine (.img) HFA Sí Sí Sin límites
Atlantis MFF2e HKV Sí Sí Sin límites
Visualización y análisis de imágenes (WinDisp) IDA Sí Sí 2 GB
Mapa ráster de ILWIS (.mpr, .mpl) ILWIS Sí Sí -
Japonés DEM (.mem) JDEM No Sí -
JPEG JFIF (.jpg) JPEG Sí Sí 4GiB (dimensiones máximas 65500x65500)
JPEG2000 (.jp2, .j2k) JPEG2000 Sí Sí 2GiB
JPEG2000 (.jp2, .j2k) JP2KAK Sí Sí Sin límites
Conjunto de datos de nivel 1b del orbitador polar NOAA (AVHRR) L1B No Sí -
Erdas 7.x .LAN y .GIS LAN No Sí 2GB
En memoria Raster MEM Sí Sí 2GiB
Atlantis MFF MFF Sí Sí Sin límites
Base de datos de imágenes sin fisuras de resolución múltiple MrSID No Sí -
Formato de datos NDF NLAPS No Sí Sin límites
NITF NITF Sí Sí
NetCDF netCDF Sí Sí 2GB
Puente OGDI OGDI No Sí -
PCI .aux etiquetado PAux Sí No Sin límites
Archivo de base de datos de PCI Geomatics PCIDSK Sí Sí Sin límites
Gráficos de red portátiles (.png) PNG Sí No
PCRaster (.map) PCRaster Sí No
Netpbm (.ppm, .pgm) PNM Sí No Sin límites
RadarSat2 XML (product.xml) RS2 No Sí 4GB
USGS SDTS DEM (* CATD.DDF) SDTS No Sí -
SAR CEOS SAR_CEOS No Sí -
USGS ASCII DEM (.dem) USGSDEM No Sí -
Mapa de píxeles X11 (.xpm) XPM Sí No
Ubicación contenido SEO
r.in.gdal intenta preservar la información de proyección al importar conjuntos de datos si el
El formato de origen incluye información de proyección y si el controlador GDAL lo admite. Si el
la proyección del conjunto de datos de origen no coincide con la proyección de la ubicación actual
r.in.gdal informará un mensaje de error (la proyección del conjunto de datos no parece coincidir
ubicación actual) y luego reportar los parámetros PROJ_INFO del conjunto de datos de origen.
Si el usuario desea ignorar la diferencia entre el aparente sistema de coordenadas del
datos de origen y la ubicación actual, pueden pasar el -o bandera para anular la proyección
comprobar.
Si el usuario desea importar los datos con la definición de proyección completa, es posible
para que r.in.gdal cree automáticamente una nueva ubicación basada en la proyección y las extensiones
del archivo que se está leyendo. Esto se logra pasando el nombre que se utilizará para el nuevo
ubicación a través del località parámetro. Al completar el comando, se abrirá una nueva ubicación.
se han creado (solo con un conjunto de mapas PERMANENTE), y el ráster se habrá importado
con lo indicado salida nombre en el conjunto de mapas PERMANENTE.
Soporte for GCP
En caso de que la imagen contenga GCP, se escriben en un archivo POINTS dentro de un grupo de imágenes.
Se pueden utilizar directamente para i.rectify.
El sistema dirigidos La opción le permite volver a proyectar automáticamente los GCP desde sus propios
proyección en otra proyección leída del archivo PROJ_INFO del nombre de la ubicación
dirigidos.
Si dirigidos la ubicación no existe, se creará una nueva ubicación que coincida con el
definición de proyección de los GCP. El objetivo del grupo de salida se establecerá en el nuevo
location, y i.rectify ahora se puede utilizar sin ninguna preparación adicional.
Algunas imágenes de satélite (por ejemplo, NOAA / AVHRR, ENVISAT) pueden contener cientos o miles de
GCP. En estos casos, se recomienda la transformación de coordenadas de spline de placa delgada, ya sea
antes de importar con gdalwarp -tps o después de importar con rectificar -t.
NOTAS
La importación de archivos grandes puede ser significativamente más rápida al configurar memoria al tamaño de la
fichero de entrada.
El sistema r.in.gdal El comando admite las siguientes funciones, siempre que el
el controlador de formato lo admite:
Tabla de colores
Las bandas con tablas de colores asociadas tendrán las tablas de colores transferidas. Tenga en cuenta que
si la fuente no tiene mapa de colores, r.in.gdal en GRASS 5.0 no emitirá ningún mapa de colores. Usar
r.colors map = ... color = gray para asignar un mapa de colores en escala de grises. En una versión futura de
GRASS r.in.gdal probablemente se actualizará para emitir automáticamente mapas de color en escala de grises.
Tipos de datos
Se admiten la mayoría de los tipos de datos GDAL. Las bandas de tipo Float32 y Float64 se traducen como
Células de punto flotante de GRASS (pero no de doble precisión ... esto podría agregarse si
necesario), y la mayoría de los otros tipos se traducen como celdas enteras de GRASS. Esto incluye
Fuentes de datos enteros de 16 bits. Las bandas de datos complejas (algunos formatos de datos de señal SAR) son
traducido a dos capas de celdas de punto flotante (* .real e * .imaginary).
Georreferenciación
Si el conjunto de datos tiene información de georreferenciación afín, esto se utilizará para establecer la
Bordes norte, sur, este y oeste. Los coeficientes rotacionales se ignorarán, lo que resultará
en un posicionamiento incorrecto para conjuntos de datos rotados.
Proyección
La proyección de los conjuntos de datos se utilizará para comparar con la ubicación actual o para definir una
nueva ubicacion. Internamente, GDAL representa proyecciones en OpenGIS Well Known Text
formato. Se admite un gran subconjunto del conjunto total de proyecciones de GRASS.
Valores nulos
Las bandas ráster para las que GDAL reconoce un valor nulo tendrán píxeles nulos
transformado en nulos de estilo GRASS durante la importación. Muchos formatos genéricos (y formatos
pobremente compatible con GDAL) no tienen una forma de reconocer píxeles nulos, en cuyo caso
r.null debe usarse después de la importación.
GCP
Los conjuntos de datos que tienen puntos de control terrestre se importarán como un archivo de PUNTOS
asociado con el grupo de imágenes. Conjuntos de datos con una sola banda que de otra manera
se han traducido como un mapa ráster simple también tendrá un grupo de imágenes asociado
si hay puntos de control en tierra. El sistema de coordenadas del control de tierra
r.in.gdal informa los puntos, pero no se conservan. Depende del usuario asegurarse
que la ubicación establecida con i.target tiene un sistema de coordenadas compatible antes
usando los puntos con i.rectify.
Mejoras planificadas a r.in.gdal en el futuro incluir soporte para informar de todo
conocido sobre un conjunto de datos si el salida El parámetro no está configurado.
Error Mensajes
"ERROR: Entrada mapa is girada - no puede importar."
En este caso, la imagen debe primero rotarse externamente, aplicando la información de rotación almacenada
en el campo de metadatos del archivo de imagen ráster. Por ejemplo, el software gdalwarp puede ser
utilizado para transformar el mapa a North-up (nota, hay varios parámetros gdalwarp para
seleccione el algoritmo de remuestreo):
gdalwarp rotado.tif northup.tif
"ERROR: Proyección of datos sí no Aparecer a partido de la forma más corriente localización."
Debe crear una ubicación cuya proyección coincida con los datos que desea importar. Tratar
usando località parámetro para crear una nueva ubicación basada en la información de proyección en
el archivo. Si lo desea, puede volver a proyectarlo a otra ubicación con r.proj.
Alternativamente, puede anular este error utilizando el -o bandera.
"ADVERTENCIA: G_set_window (): ilegal latitud for Norte"
Las ubicaciones de latitud / longitud en GRASS no pueden tener regiones que excedan los 90 ° norte o
Sur. Las imágenes no georreferenciadas tendrán coordenadas basadas en el número de imágenes de
píxeles: 0,0 en la parte inferior izquierda; cols, filas en la parte superior derecha. Normalmente, las imágenes serán mucho
más de 90 píxeles de alto, por lo que el SIG se niega a importarlo. Si está seguro de que
Los datos son apropiados para su ubicación Lat / Lon y tienen la intención de restablecer los límites del mapa con
de la forma más r.región módulo directamente después de la importación, puede utilizar el -l bandera para restringir el mapa
coordenadas a valores legales. Si bien los límites y la resolución resultantes probablemente serán
incorrecto para su mapa, los datos del mapa se mantendrán inalterados y seguros. Después de restablecer a conocido
límites con r.región deberías comprobarlos dos veces con r.info, prestando especial atención a
la resolución del mapa. En la mayoría de los casos, querrá importar al archivo nativo del archivo de datos.
proyección, o en una ubicación XY simple y use las herramientas de Georrectificación (rectificar et
al.) para proyectar correctamente en la ubicación de destino. los -l bandera debe only ser usado si tu
sabe que la proyección es correcta pero la georreferenciación interna se ha perdido, y usted
Conozca de antemano cuáles deben ser los límites y la resolución del mapa.
EJEMPLOS
CAD Data
El proyecto European Climate Assessment and Dataset (ECAD) proporciona datos climáticos para
Europa desde 1950 hasta 2010. Para importar los diferentes fragmentos de datos proporcionados por el
proyecto como archivos netCDF, el parámetro de compensación se puede utilizar para obtener ráster numerado diariamente
mapas desde el 1 de enero de 1950 en adelante. Asegúrese de estar en una ubicación LatLong.
# Importar datos de precipitación
r.in.gdal -o entrada = rr_0.25deg_reg_1950-1964_v4.0.nc salida = compensación de precipitación = 0
r.in.gdal -o entrada = rr_0.25deg_reg_1965-1979_v4.0.nc salida = compensación de precipitación = 5479
r.in.gdal -o entrada = rr_0.25deg_reg_1980-1994_v4.0.nc salida = compensación de precipitación = 10957
r.in.gdal -o entrada = rr_0.25deg_reg_1995-2010_v4.0.nc salida = compensación de precipitación = 16436
# Importar datos de presión de aire
r.in.gdal -o input = pp_0.25deg_reg_1950-1964_v4.0.nc salida = air_pressure offset = 0
r.in.gdal -o input = pp_0.25deg_reg_1965-1979_v4.0.nc salida = air_pressure offset = 5479
r.in.gdal -o input = pp_0.25deg_reg_1980-1994_v4.0.nc salida = air_pressure offset = 10957
r.in.gdal -o input = pp_0.25deg_reg_1995-2010_v4.0.nc salida = air_pressure offset = 16436
# Importar datos de temperatura mínima
r.in.gdal -o entrada = tn_0.25deg_reg_1950-1964_v4.0.nc salida = temperatura_min offset = 0
r.in.gdal -o entrada = tn_0.25deg_reg_1965-1979_v4.0.nc salida = temperatura_min offset = 5479
r.in.gdal -o entrada = tn_0.25deg_reg_1980-1994_v4.0.nc salida = temperatura_min offset = 10957
r.in.gdal -o entrada = tn_0.25deg_reg_1995-2010_v4.0.nc salida = temperatura_min offset = 16436
# Importar datos de temperatura máxima
r.in.gdal -o entrada = tx_0.25deg_reg_1950-1964_v4.0.nc salida = temperatura_máx offset = 0
r.in.gdal -o entrada = tx_0.25deg_reg_1965-1979_v4.0.nc salida = temperatura_máx offset = 5479
r.in.gdal -o entrada = tx_0.25deg_reg_1980-1994_v4.0.nc salida = temperatura_máx offset = 10957
r.in.gdal -o entrada = tx_0.25deg_reg_1995-2010_v4.0.nc salida = temperatura_máx offset = 16436
# Importar datos de temperatura media
r.in.gdal -o input = tg_0.25deg_reg_1950-1964_v4.0.nc salida = temperature_mean offset = 0
r.in.gdal -o input = tg_0.25deg_reg_1965-1979_v4.0.nc salida = temperature_mean offset = 5479
r.in.gdal -o input = tg_0.25deg_reg_1980-1994_v4.0.nc salida = temperature_mean offset = 10957
r.in.gdal -o input = tg_0.25deg_reg_1995-2010_v4.0.nc salida = temperature_mean offset = 16436
GTOPO30 DEM
Para evitar que los datos de GTOPO30 se lean incorrectamente, puede agregar una nueva línea "PIXELTYPE
SIGNEDINT "en el .HDR para forzar la interpretación del archivo como firmado en lugar de sin firmar
enteros. Luego, se puede importar el archivo .DEM. Por último, por ejemplo, la tabla de colores de 'terreno' puede
ser asignado al mapa importado con r.colores.
GLOBE DEM
Para importar mosaicos GLOBE DEM (resolución de aproximadamente 1 km, mejor que los datos de GTOPO30 DEM), el usuario
tiene que descargar adicionalmente los archivos HDR relacionados. Finalmente, por ejemplo, el color del 'terreno'
La tabla se puede asignar al mapa importado con r.colores. Consulte también su portal DEM.
Mundoclim.org
Para importar datos de Worldclim, se debe agregar la siguiente línea a cada archivo .hdr:
PIXELTYPE FIRMADO
HDF
La importación de bandas HDF requiere la especificación de las bandas individuales como las ve
GDAL:
# Ejemplo MODIS FPAR
gdalinfo MOD15A2.A2003153.h18v04.004.2003171141042.hdf
...
Subconjuntos de datos:
SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:"MOD15A2.A2003153.h18v04.004.2003171141042.hdf":MOD_Grid_MOD15A2:Fpar_1km
SUBDATASET_1_DESC = [1200x1200] Fpar_1km MOD_Grid_MOD15A2 (entero sin signo de 8 bits)
SUBDATASET_2_NAME=HDF4_EOS:EOS_GRID:"MOD15A2.A2003153.h18v04.004.2003171141042.hdf":MOD_Grid_MOD15A2:Lai_1km
SUBDATASET_2_DESC = [1200x1200] Lai_1km MOD_Grid_MOD15A2 (entero sin signo de 8 bits)
...
# importación de la primera banda, aquí FPAR 1km:
r.in.gdal HDF4_EOS:EOS_GRID:"MOD15A2.A2003153.h18v04.004.2003171141042.hdf":MOD_Grid_MOD15A2:Fpar_1km
out = fpar_1km_2003_06_02
# ... igualmente para otras bandas HDF en el archivo.
Use r.in.gdalgrass en línea usando los servicios de onworks.net