Abrir archivos BIL, BIP o BSQ en QGIS

Cuanto trate con conjuntos de datos de teledetección y científicos, a menudo uno se encuentra con datos en formatos como BIL, BIP o BSQ. La librería GDAL - que es usada por QGIS para leer archivos raster - tiene soporte para estos formatos, pero no puede abrir estos archivos por si mismos. Iremos por el proceso para crear archivos de soporte para que estos formatos puedan ser leídos por QGIS.

Band interleaved by line (BIL), band interleaved by pixel (BIP), and band sequential (BSQ) son métodos comunes de organización de datos de imagen para imágenes multibanda. (Lea más acerca de estos formatos en)

Típicamente, estos archivos están acompañados por un archivo .hdr. Si tu conjunto de datos vino con un archivo .hdr, asegúrese que el nombre raíz del archivo .bil, .bsq o .bip coincidan y que estén en el mismo directorio. Por ejemplo, si el archivo se llama image.bil, el archivo asociado deberá llamarse image.hdr y estar presente en el mismo directorio que el archivo image.bil. Entonces vaya a Capa ‣ Agregar Capa Raster, seleccione el archivo image.bil y abrirá sin problemas.

Muchas veces, los archivos no vienen con un archivo asociado .hdr. En tales casos, debes crear este archivo a mano como se muestra en este tutorial.

Obtener los datos

Usaremos el AVHRR Global Land Cover Classification data del Global Land Cover Facility como ejemplo.

Los conjuntos de datos Cobertura Global son distribuidos como archivos BSQ. Descargue el conjunto de datos 1 Degree pixel resolution.

Para su comodidad, puede descargar directamente una copia del conjunto de datos del siguiente enlace:

gl-latlong-1deg-landcover.bsq.gz

Fuente de Datos [GLCF]

Procedmiento

  1. Descomprima y extraiga el archivo .bsq. En Windows, puede utilizar la excelente utilidad 7-Zip para leer y extraer un archivo .gz. Verá que sólo tiene un archivo .bsq llamado gl-latlong-1deg-landcover.bsq. No existe archivo hdr.

../_images/1171.png
  1. Note que si intenta abrir el archivo gl-latlong-1deg-landcover.bsq en QGIS tal como está, obtendrá un mensaje de error.

../_images/2131.png
  1. Para superar este error, crearemos un archivo encabezado con extensión .hdr. El archivo encabezado contiene información acerca del conjunto de datos y cómo está organizado. Usualmente, esta información es proporcionada como parte de los Metadatos para el conjunto de datos. Si no tienes los metadatos, revisa en el sitio web o la documentación para pistas. Algo de la información puede ser adivinada si no la conoces. En el caso de este conjunto de datos, la página de descarga de datos enlaza a los metadatos. Descargue los metadatos y ábralos.

../_images/370.png
  1. El archivo .hdr necesita ser un archivo de texto plano en el siguiente formato. Alguno de los parámetros nos son dados y algunos necesitan ser investigados. Aprenda más acerca del formato en.

ncols <number of columns or width of the raster>
nrows <number of rows or height of the raster>
cellsize <pixel size or resolution>
xllcorner <X coordinate of lower-left corner of the raster>
yllcorner <Y coordinate of the lower-left corner of the raster>
nodata_value <pixel value to be ignored>
nbits <number of bits per pixel>
pixeltype <type of values stored in a pixel, typically float or integer>
byteorder <byte order in which image pixel values are stored, msb or lsb>
  1. Abra un editor de texto y cree un archivo en el formato especificado en el paso previo. Guarde el archivo como gl-latlong-1deg-landcover.hdr. Asegúrese que el archivo no tiene .txt al final. Alguno de los valores en los archivos texto son fáciles de comprender. Las ncols y nrows vienen de los metadatos como el Número de Líneas y Número de Píxeles por Línea. El cellsize es 1 como la resolución de Píxel de los metadatos. La coordenada X,Y de la esquina inferior izquierda necesita ser identificada por nosotros. Debido a que el archivo cubre todo el mundo y las unidades son lat/long, xllcorner y yllcorner son -180 y -90 respectivamente. No tenemos ninguna información acerca de nodata_value, así que -9999 es un valor seguro. De los metadatos de nuevo, Pixel Format es Byte, así que nbits será igual a 8 y pixeltype será byte_unsigned. No tenemos información acerca de byteorder, así que lo dejamos como msbfirst. Puede descargar el archivo HDR correctamente formateado de aquí.

../_images/550.png
  1. Ahora usted tiene el archivo encabezado, póngalo en el mismo directorio que gl-latlong-1deg-landcover.bsq. Luego en QGIS, vaya a Capa ‣ Agregar Capa Raster. Seleccione gl-latlong-1deg-landcover.bsq como su entrada y clic Abrir.

../_images/648.png
  1. En la siguiente pantalla, puede ser que se te pida la CRS. Debido a que los datos están en Lat/Lon, elige WGS84 EPSG:4326 como tu CRS. Ahora verás el conjunto de datos cargado en QGIS.

../_images/747.png
If you liked tutorials on this site and do check out spatialthoughts.com for more free resources.
comments powered by Disqus

This work is licensed under a Creative Commons Attribution 4.0 International License