¿Abrir un archivo vrt?

Question

Me gustaría convertir un vrt en un conjunto de datos ráster mediante gdal_translate . Estoy haciendo lo siguiente:

gdal_translate in.vrt out.tif

Sin embargo, obtengo el siguiente error:

`ERROR 4: 'in_vrt' no reconocido como formato de archivo soportado.

He comprobado y vuelto a comprobar el nombre y la ruta de mi archivo y no consigo averiguar qué estoy haciendo mal. Traté de abrir en Arcmap y recibí este mensaje de error:

Conjunto de datos ráster no válido. Error al crear capa raster

y en QGIS tengo

Capa no válida: Proveedor GDAL No se puede abrir el conjunto de datos GDAL in.vrt: `in.vrt' no se reconoce como formato de archivo compatible. [ ] no es válido (proveedor: gdal, URI: in.vrt

¿Alguien sabe qué estoy haciendo mal? ¿Cómo puedo abrir este archivo vrt sin errores?

EDITAR He aquí algunas líneas de un vrt típico que estoy utilizando

<OGRVRTDataSource>
        <OGRVRTLayer name="gouda_rws_env_dsc_v3_ds_qh">
            <SrcDataSource>gouda_rws_env_dsc_v3_ds_qh.csv</SrcDataSource>
            <GeometryType>wkbPoint</GeometryType>
            <LayerSRS>WGS84</LayerSRS>
            <GeometryField encoding="PointFromColumns" x="pnt_lon" y="pnt_lat"/>
            <Field name="pnt_id" src="pnt_id" type="String" />
            <Field name="pnt_lat" src="pnt_lat" type="Real" />
            <Field name="pnt_lon" src="pnt_lon" type="Real" />
            <Field name="pnt_rdx" src="pnt_rdx" type="Real" />
            <Field name="pnt_rdy" src="pnt_rdy" type="Real" />
            <Field name="pnt_demheight" src="pnt_demheight" type="Real" />
            <Field name="pnt_height" src="pnt_height" type="Real" />

Answer 1

Su VRT es un conjunto de datos vectoriales. gdal_translate se utiliza para convertir entre formatos de trama, no de vector a trama.

En su lugar, puede utilizar:

• ogr2ogr para generar un shapefile u otro formato vectorial,
  • QGIS Processing toolbox | GDAL | Conversión de vectores | Convertir formato
• gdal_rasterize para grabar los puntos en un tiff, o
  • QGIS Processing toolbox | GDAL | Conversión vectorial | Rasterizar
• gdal_grid para interpolar los puntos a un raster de superficie (he observado algunos campos de elevación en tu VRT).
  • QGIS Processing toolbox | GDAL | Raster Analysis | Grid (*)

Si el archivo CSV subyacente está destinado a ser un conjunto de datos rasterizados en Formato XYZ gdal/QGIS debería ser capaz de leerlo directamente como un raster si eliminas los campos adicionales y sólo mantienes "pnt_lon", "pnt_lat" y "pnt_height" (O "pnt_demheight" pero no ambos), ordenados por lat y luego por lon. A continuación, podría utilizar gdal_translate xyz.csv out.tif

GDAL admite la lectura y escritura de conjuntos de datos ráster ASCII XYZ reticulados (es decir, los XYZ no reticulados, LIDAR XYZ, etc. deben abrirse por otros medios. Consulte la documentación de la utilidad gdal_grid).

Estos conjuntos de datos son archivos ASCII con (al menos) 3 contiene las coordenadas X e Y del centro de la celda y el valor de la celda.

El espacio entre cada celda debe ser constante y admitido. Las celdas con las mismas coordenadas Y deben colocarse en líneas consecutivas. Para un mismo valor de coordenada Y, las líneas del conjunto de datos deben estar organizarse por valores X crecientes. El valor de la coordenada Y puede aumentar o disminuir. Los separadores de columna admitidos son espacio, coma, punto y coma y tabulaciones.

El controlador intenta autodetectar una línea de cabecera y buscará 'x', 'lon' o 'este' para detectar el índice de la columna X, 'y', 'lat' o 'norte' para la columna Y y 'z', 'alt' o 'altura' para la columna Z. columna. Si no identificada en la cabecera, las columnas X, Y y Z (en ese orden) se las 3 primeras columnas de cada línea.