Diferencias

Muestra las diferencias entre dos versiones de la página.

Enlace a la vista de comparación

manuales:sig:ortofotos_mdt [2016/09/22 00:28] (actual)
jan creado
Línea 1: Línea 1:
  
 +====== Procesar Ortofotos y MDT con GDAL ======
 +
 +===== Ortofotografías =====
 +
 +Fotografías aéreas ortorectificadas o ortofotos, se estima que son del año 2010-2011. ​
 +
 +  * En un formato de datos propietario que requiere liciencia para usarlo (ECW).
 +  * Se desconoce el sistema de referencia espacial de los archivos. Se supone que es SIRGAS 2000 / Zona 18S, tal como los MDT entregados.
 +
 +Las ortofotografías tienen las siguientes resoluciones espaciales:
 +  * Loreto: 0,3m
 +  * Orellana: 0,5m
 +  * Sacha: 0,4m 
 +  * Aguarico: 0,5m
 +
 +Se realizan los siguientes procesamientos:​
 +  - Transformación de coordenadas de SIRGAS 2000, zona 18S a WGS84, Zona 18S.
 +  - Conversión a un formato libre (GeoTIFF) con compresión JPEG con colores YCBCR
 +  - Cortar con máscara y fondo en transparente
 +  - Crear pirámides
 +
 +
 +El siguiente apartado resume los comandos que se han realizado para transformar los imágenes en tres pasos usando los comandos ''​gdal_translate'',​ ''​gdalwarp''​ y ''​gdaladdo''​ de la librería GDAL.
 +
 +<code bash>
 +# Primer paso es convertir el formato del propietario jpeg2000 al libre GeoTIFF, siendo una ortofoto con compresión jpeg que es óptima en el espacio de colores YCBCR.
 +gdal_translate -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR -co TILED=YES --config GDAL_TIFF_INTERNAL_MASK YES JOYA_DE_LOS_SACHAS.jp2 /​media/​jan/​geoportal/​georellana/​geodatos/​sacha_jpeg_ycbcr.tif
 +# En el caso de que sea necesario quitar la capa alfa
 +gdal_translate -b 1 -b 2 -b 3 -co ALPHA=NO -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR -co BIGTIFF=YES -co TILED=YES ORELLANA.jp2 orellana_jpeg_ycbcr.tif
 +gdal_translate -b 1 -b 2 -b 3 -co ALPHA=NO -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR -co BIGTIFF=YES -co TILED=YES LORETO.jp2 loreto_jpeg_ycbcr.tif
 +gdal_translate -b 1 -b 2 -b 3 -co ALPHA=NO -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR -co BIGTIFF=YES -co TILED=YES LORETO.jp2 aguarico_oriental_jpeg_ycbcr.tif
 +gdal_translate -b 1 -b 2 -b 3 -co ALPHA=NO -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR -co BIGTIFF=YES -co TILED=YES LORETO.jp2 aguarico_occidental_jpeg_ycbcr.tif
 +
 +# Segundo paso es cortar el imagen con una masquara que se ha preparado a partir de la capa de los límites cantonales y el distribuidor de la ortofoto 1:5000
 +gdalwarp -s_srs EPSG:31978 -t_srs EPSG:32718 -q -cutline sacha.shp -crop_to_cutline -dstalpha -of GTiff -co COMPRESS=JPEG -co TILED=YES -co BIGTIFF=YES sacha_jpeg_ycbcr.tif sacha_clip.tif
 +gdalwarp -q -cutline aguarico_occidental.shp -crop_to_cutline -dstalpha -of GTiff -co COMPRESS=JPEG -co TILED=YES -co BIGTIFF=YES aguarico_occidental_jpeg_ycbcr.tif aguarico_occidental_clip.tif
 +
 +# Tercer paso: Creando los pirámides
 +gdaladdo --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config INTERLEAVE_OVERVIEW PIXEL --config GDAL_TIFF_INTERNAL_MASK YES -r average sacha.tif 2 4 8 16 128 254
 +gdaladdo --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config INTERLEAVE_OVERVIEW PIXEL --config GDAL_TIFF_INTERNAL_MASK YES -r average orellana.tif 2 4 8 16 32 64 128 254
 +gdaladdo --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config INTERLEAVE_OVERVIEW PIXEL --config GDAL_TIFF_INTERNAL_MASK YES -r average loreto.tif 2 4 8 16 32 64 128 254
 +gdaladdo --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config INTERLEAVE_OVERVIEW PIXEL --config GDAL_TIFF_INTERNAL_MASK YES -r average aguarico_oriental.tif 2 4 8 16 32 64 128 254
 +gdaladdo --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config INTERLEAVE_OVERVIEW PIXEL --config GDAL_TIFF_INTERNAL_MASK YES -r average aguarico_occidental.tif 2 4 8 16 32 64 128 254
 +</​code>​
 +
 +===== Modelo de elevación del terreno (MDT) =====
 +
 +
 +Se desconocen los años de elaboración de los Modelos Digitales de Terreno (MDT)
 +
 +Los MDT entregados no son homogéneas con respecto a su resolución espacial:
 +  * Orellana: 5m
 +  * Loreto: 3m
 +  * Sacha: 4m
 +  * Aguarico: 5m
 +
 +
 +Se realizan los siguientes procesamientos:​
 +  * Los MDT se reproyectan al Sistema de Referencia Especial WGS84 / UTM Zona 18S (EPSG: 32718).
 +  * Se comprimen DEFLATE (sin pérdida) con predictor=1 y zlevel=6
 +  * Se generan pirámides internos
 +
 +<code bash>
 +gdalwarp -s_srs EPSG:31978 -t_srs EPSG:32718 -co COMPRESS=DEFLATE -co PREDICTOR=1 -co ZLEVEL=6 -co BIGTIFF=YES -co TILED=YES -dstnodata -32767 orellana.tif epsg32718/​orellana.tif
 +gdalwarp -s_srs EPSG:31978 -t_srs EPSG:32718 -co COMPRESS=DEFLATE -co PREDICTOR=1 -co ZLEVEL=6 -co BIGTIFF=YES -co TILED=YES -dstnodata -32767 aguarico.tif epsg32718/​aguarico.tif
 +gdalwarp -s_srs EPSG:31978 -t_srs EPSG:32718 -co COMPRESS=DEFLATE -co PREDICTOR=1 -co ZLEVEL=6 -co BIGTIFF=YES -co TILED=YES -dstnodata -32767 loreto.tif epsg32718/​loreto.tif
 +# Es necesario cortar el MDT del Sacha con los límites de la provincia, para quitar los valores erroneos en zonas fuera de la provincia. La máscara han sido editado manualmente para cerrar vacíos entre orellana y sacha.
 +gdalwarp -q -cutline ../​clip/​sacha_mask.shp -crop_to_cutline -s_srs EPSG:31978 -t_srs EPSG:32718 -co COMPRESS=DEFLATE -co PREDICTOR=1 -co ZLEVEL=6 -co BIGTIFF=YES -co TILED=YES -dstnodata -32767 sacha.tif epsg32718/​sacha.tif
 +
 +gdaladdo ​ --config COMPRESS_OVERVIEW DEFLATE --config PREDICTOR_OVERVIEW 1 -r AVERAGE epsg32718/​orellana.tif 2 4 8 16 32 64 128 256
 +gdaladdo --config COMPRESS_OVERVIEW DEFLATE --config PREDICTOR_OVERVIEW 1 -r AVERAGE epsg32718/​aguarico.tif 2 4 8 16 32 64 128 256
 +gdaladdo --config COMPRESS_OVERVIEW DEFLATE --config PREDICTOR_OVERVIEW 1 -r AVERAGE epsg32718/​loreto.tif 2 4 8 16 32 64 128 256
 +gdaladdo --config COMPRESS_OVERVIEW DEFLATE --config PREDICTOR_OVERVIEW 1 -r AVERAGE epsg32718/​sacha.tif 2 4 8 16 32 64 128 256
 +
 +# El MDT se importa a GRASS GIS y se exporta un MDT de la provincia creado con r.patch
 +# Comprimir y crear pyrámides del raster exportado
 +gdal_translate ​ -co COMPRESS=DEFLATE -co PREDICTOR=1 -co ZLEVEL=6 -co BIGTIFF=YES -co TILED=YES mdt_provincia_output.tif mdt_provincia.tif
 +gdaladdo --config COMPRESS_OVERVIEW DEFLATE --config PREDICTOR_OVERVIEW 1 -r AVERAGE mdt_provincia.tif 2 4 8 16 32 64 128 256
 +</​code>​
 +
 +===== Enlaces de interés =====
 +
 +  * Acerca de compresión sin pérdidas [[http://​linfiniti.com/​2011/​05/​gdal-efficiency-of-various-compression-algorithms/​|GDAL. Efficiency of various compression algorithms]]
 +  * Compresión con pérdidas [[http://​blog.cleverelephant.ca/​2015/​02/​geotiff-compression-for-dummies.html|Paul Ramsey: GeoTIFF compression for Dummies]]