Gestión e integración de datos altimétricos y batimétricos en la costa andaluza: uso del “model builder”

Tecnologías de la Información Geográfica: la Información Geográfica al Servicio de los Ciudadanos. Secretariado de Publicaciones de la Universidad de Sevilla. Sevilla 2010 956 Ojeda Zújar, J., Álvarez Francoso, J.I., Fraile Jurado, P., Márquez Pérez, J. y Sánchez Rodríguez, e. (2010): Gestión e integración de datos altimétricos y batimétricos en la costa andaluza: el uso del “model builder”. En: Ojeda, J., Pita, M.F. y Vallejo, I. (eds.), Tecnologías de la información geográfica: la información geográfica al servicio de los ciudadanos. Secretariado de publicaciones de la Universidad de Sevilla. Sevilla. pp. 956-970. ISBN: 978-84-472-1294-1 GESTIÓN E INTEGRACIÓN DE DATOS ALTIMÉTRICOS Y BATIMÉTRICOS EN LA COSTA ANDALUZA: EL USO DEL “MODEL BUILDER” Ojeda Zújar, J1., Álvarez Francoso, J.I1., Fraile Jurado, P1., Márquez Pérez, J1. y Sánchez Rodríguez, E1. (1) Departamento de Geografía Física y AGR. Universidad de Sevilla. C/Maria de Padilla s/n. 41004 Sevilla. Email: [email protected], [email protected], [email protected], [email protected] & [email protected]. RESUMEN En esta comunicación se presenta la estrategia adoptada en la incorporación de la información altimétrica y batimétrica de la costa andaluza al Subsistema del Litoral y Medio Marino (SSLyMM), realizado para la Consejería de Medio Ambiente. En primer lugar se opto por modificar el objetivo final: no se trata de generar un único MDE continuo, sino de mantener los datos altimétricos y batimétricos permanentemente actualizados en un base de datos espacial y diseñar una herramienta interactiva que produzca los MDEs a demanda en función de las aplicaciones en las que se usará. Para ello, se diseño una herramienta utilizando la aplicación ModelBuilder de ArcGis 9.3 que encadena varios procesos de análisis espacial a partir de la información original almacenada en la base de datos (curvas, cotas, sondas o datos raster). Estos procesos se concretan en algoritmos que primero interpolan la batimetría y posteriormente ajustan el datum del MDE resultante al datum altimétrico (mareógrafo Alicante). En segundo lugar se interpola la altimetría y se ajusta a la línea de costa del SSLyMM y, en tercer lugar, se extrae el 0 hidrográfico ya ajustado al 0 de Alicante y se interpola la zona intermareal con los datos existentes desde dicho cero hasta la línea de costa. La posibilidad de configurar variables y parámetros que se deciden interactivamente por el usuario antes de ejecutar la herramienta permite definir el área de análisis y la resolución espacial del modelo resultante, el algoritmo de interpolación a utilizar y la paleta de color a aplicar al MDE resultante entre otras opciones. Por otra parte, esta estrategia permite incorporar de forma inmediata la nueva información disponible (muy variable en el medio submarino) y ajustar las características del MDE (resolución espacial, interpolador, etc.) resultante a las demandas de la aplicación temática para la que va a utilizar (modelado de oleaje, navegación marítima, evaluación de riesgos, generación sombreados para fondos cartográficos, etc.). Palabras Clave: Costa andaluza, altimetría, batimetría, model builder, SIG ABSTRACT This paper presents the methodology used for the integration of the altimetric and bathymetric data of the andalusiam coast into the Littoral and Marine Environment Subsystem (SSLyMM), made for the Environmental Regional Ministry. The aim of this work was not the generation of a only continuous digital elevation model (DEM), as it was firstly required, but also the design of a strategy for keeping the altimetric and bathymetric Tecnologías de la Información Geográfica: la Información Geográfica al Servicio de los Ciudadanos. Secretariado de Publicaciones de la Universidad de Sevilla. Sevilla 2010 957 data always updated in a geodatabase and the development of an interactive tool to produce different DEMs depending on the requirements of the potential user. The tool was designed using the application ModelBuilder in ArcGis 9.3 that uses different spatial analysis processes applied to the original data stored in the spatial database. Those processes reflects a set of algorithms that firstly interpolated the bathymetry data and then adjust the resulting DEM to the national altimetric datum, located in the Mean Sea Level at the Alicante tide gauge in Spain (NMMA). The second step was the interpolation of the altimetry data and the adjustment of the shoreline contained in the SSLyMM. The third step was the extraction of the corrected hydrographic zero contour (once referred to the NMMA) and the interpolation of the intertidal area elevation by means of the available data of this area. The area of the DEM, its spatial resolution, the algorithm that would be used for spatial interpolation and other variables and parameters can be chosen by the user before executing the tool. This strategy let the incorporation of potential new information available (what is quite usual on bathymetry due to the scarcity of data), and the adjustment of the DEM characteristic (spatial resolution, interpolation methods...) to the requirements of the use and applications that the DEM might have. Key Words: Andalucía Coast, altimetry, bathymetry, model builder, GIS INTRODUCCIÓN La integración de información altimétrica y batimétrica se ha convertido en una exigencia tanto para la mayoría de las disciplinas científicas que abordan las zonas costeras como para su planificación y gestión medioambiental y territorial (Batier et alia, 2007). Casi todos los procesos de integración persiguen la obtención de un MDE continuo de las zonas emergidas y sumergidas. Esta necesidad encuentra numerosos problemas en su puesta en práctica: diferentes sistemas de referencia de coordenadas, diferente datum altimétrico/batimétrico, diferente precisión y resolución espacial en los datos de partida, etc. Existen experiencias de integración a escala global (GEBCO, ETOPO 1 Global Relief Model…) y serios intentos en algunos estados (Whitfield, et alia, 2003) como el proyecto ICZMap en Reino Unido a escalas mas detalladas, siendo especialmente complejo la armonización de datum (Turner et alia, 2010). Estos problemas hacen que el MDE final (generalmente en formato raster) sea de difícil obtención y claramente orientado a aplicaciones específicas (modelado de oleaje, riesgo de inundaciones, navegación marítima, etc.) lo que condiciona su producción y los algoritmos de interpolación/integración. Por todo ello, a diferencia de otras experiencias, el objetivo de esta comunicación no es obtener un único MDE que integre las zonas emergidas y sumergidas y que luego pueda ser reutilizado por cualquier aplicación posterior (Fig 1), ya que este sería deudor de los datos de partida y los algoritmos de interpolación empleados para su generación. Conscientes de estas restricciones, en esta comunicación se presenta la estrategia adoptada en la incorporación de la información altimétrica y batimétrica de la costa andaluza al Subsistema del Litoral y Medio Marino (SSLyMM), realizado para la Consejería de Medio Ambiente. Tecnologías de la Información Geográfica: la Información Geográfica al Servicio de los Ciudadanos. Secretariado de Publicaciones de la Universidad de Sevilla. Sevilla 2010 958 Figura 1. Esquema tradicional del proceso de generación de un MDE que integre zonas emergidas y sumergidas en las zonas costeras para su posterior utilización de forma genérica. OBJETIVOS Y AREA DE ESTUDIO Objetivos Esta comunicación se deriva del proyecto de investigación (ver agradecimientos), realizado para la Consejería de Medio Ambiente que incluía la necesidad de incorporar la información altimétrica y batimétrica disponible de la costa andaluza al Subsistema del Litoral y Medio Marino (SSLyMM) para integrarse, posteriormente, en la Red de Información Ambiental de Andalucía (REDIAM). Como resultado de las dificultades ya descritas en el párrafo anterior de esta tarea, se adopta una estrategia diferente a la comúnmente utilizada: no se trata de generar un único MDE continuo, sino que se considera que lo importante es mantener los datos altimétricos y batimétricos permanentemente actualizados en una base de datos espacial y diseñar una herramienta interactiva que produzca los MDEs a demanda en función de las aplicaciones a las que irán destinados (Fig.2). Tecnologías de la Información Geográfica: la Información Geográfica al Servicio de los Ciudadanos. Secretariado de Publicaciones de la Universidad de Sevilla. Sevilla 2010 959 Figura 2. Esquema de la estrategia seguida para generación de MDEs que integren zonas emergidas y sumergidas en las zonas costeras, manteniendo actualizados los datos y utilizando una herramienta para su generación a demanda en función de las aplicaciones. Por ello, el principal objetivo de esta comunicación es:  Diseño una herramienta utilizando la aplicación ModelBuilder de ArcGis 9.3 que encadena varios procesos de análisis espacial a partir de la información original almacenada en la base de datos (curvas, cotas, sondas o datos raster)  Generación de capas de apoyo actualizadas que garanticen la integración, especialmente a partir de la línea de costa y las Unidades fisiográficas del SSLyMM.  Mostrar los resultados a varias escalas en diferentes zonas del litoral andaluz Área de estudio El área de estudio es toda la costa andaluza con más de 800 km de frente costero y con dos fachadas (la atlántica y la mediterránea) con diferentes características de relieve litoral, tanto el emergido como el sumergido (Fig. 3). Tecnologías de la Información Geográfica: la Información Geográfica al Servicio de los Ciudadanos. Secretariado de Publicaciones de la Universidad de Sevilla. Sevilla 2010 960 Figura 3. Área de estudio. DATOS Y METODOLOGÍA Datos  Datos altimétricos: MDT de Andalucía (2005) generado por estereocorrelación y líneas de estructura a partir del vuelo fotogramétrico 1:20.000. El MDT se distribuye en formato raster con una resolución espacial de 10 m.  Datos batimétricos: datos procedentes de las cartas náuticas del Instituto Hidrográfico de la Marina a escala 1:50.000 que proporcionan: el cero hidrográfico, curvas batimétricas y sondas puntuales.  Datos extraídos del repositorio del SSLyMM (Subsistema del Litoral y Medio Marino): - Línea de costa de 2007 a escala 1/2.500, levantada utilizando como criterio de interpretación la marca húmeda de la marea en la playa (límite inferior de la playa seca –backshore-). - Línea de bajamar extraida los vuelos fotogramétricos recientes, utilizando aquel con las condiciones de marea mas extremas (lo mas aproximado al 0 hidrográfico ante la ausencia de éste a escalas de detalle) - Línea de contacto entre marismas mareales con vegetación y las planicies intermareales desnudas. El elemento fotointerpretable mas próximo al nivel medio del mar en las zonas marismeñas. - Líneas direccionadas de los ejes de los principales cañones submarinos extraidas de los Mapas Fisiográficos del Litoral de Andalucía. Tecnologías de la Información Geográfica: la Información Geográfica al Servicio de los Ciudadanos. Secretariado de Publicaciones de la Universidad de Sevilla. Sevilla 2010 961 Metodología El complejo modelo de proceso de la información espacial, diseñado e implementado mediante la aplicación ModelBuilder (Fig. 4) sigue, sintéticamente, los siguientes pasos: Figura 4. Esquema gráfico del encadenamiento de procesos del “model builder” diseñado. 1.- Interpolación espacial de los datos originales procedentes de las cartas náuticas mediante el método Topo to Raster para conseguir un MDE raster (Fig. 5) sobre el datum hidrográfico que necesita los siguientes datos de entrada:  Línea del 0 hidrográfico  Sondas puntuales y/o isóbatas  Líneas direccionadas de los cañones submarinos Tecnologías de la Información Geográfica: la Información Geográfica al Servicio de los Ciudadanos. Secretariado de Publicaciones de la Universidad de Sevilla. Sevilla 2010 962 Figura 5. Superficie batimétrica extraída de la información de las cartas náuticas. 2.- Ajuste del Datum Altimétrico: Para la corrección del Datum hidrográfico del modelo anterior respecto al Datum Altrimético oficial (nivel medio del mar en Alicante) se genera una superficie de las diferencias entre el 0 hidrográfico de los mareógrafos andaluces (tabla 1) y el 0 de Alicante. Para ello:  Se trazan líneas ortogonales a la costa desde el emplazamiento de los mareógrafos y se les asigna el valor de la diferencia altimétrica con Alicante. Tabla 1. Diferencias entre el 0 Hidrográfico local y el nivel medio del mar en Alicante (Puertos del Estado e Instituto Español de Oceanografía) MAREÓGRAFO CERO HIDROG / NMMA MÁLAGA IEO 0,811 MALAGA PE 0,415 BONANZA 1,419 HUELVA 1,576 TARIFA 0,9172 ALGECIRAS 0,6337 Tecnologías de la Información Geográfica: la Información Geográfica al Servicio de los Ciudadanos. Secretariado de Publicaciones de la Universidad de Sevilla. Sevilla 2010 963  Con estas líneas se genera una superficie para la que hemos empleado un interpolador tipo TIN (MDT de diferencias Datum, figura 6). Figura 6. Superficie con los valores de diferencias de Datum interpolados.  Al MDT obtenido de las cartas náuticas se le suma el MDT de diferencias Datum para obtener un MDT raster de la batimetría ajustado al 0 de Alicante, que solo alcanza hasta la línea del 0 hidrográfico de las cartas nauticas: A este modelo lo llamamos Batimetría. 3.- Interpolación de la zona intermareal (0 hidrográfico/linea de costa).  Este paso consiste en la interpolación espacial de la superficie situada entre el cero hidrográfico y la línea de costa. La línea de costa puede ser utilizada como tal, ya que en todo su recorrido conserva el mismo valor. Sin embargo, el cero hidrográfico, al haber sufrido la corrección del Datum, y ser éste variable a lo largo de la costa andaluza, debe ser transformado en puntos que permitan reflejar dicha discontinuidad. Para la obtención de los puntos, se extrae un buffer de 100 metros en torno al 0 hidrográfico del modelo anterior (Batimetría) y las celdillas correspondientes son después transformadas en puntos. Serán estos puntos con batimetría corregida, variable espacialmente, los que servirán de límite inferior para la zona intermareal.  Adicionalmente se integra de la información altimétrica/batimétrica disponible de la zona intermareal corregida al Datum de Alicante si la hubiese. CARTAGENA 0,784 ALICANTE 0,219 VALENCIA 0,98 CADIZ 1,887 Tecnologías de la Información Geográfica: la Información Geográfica al Servicio de los Ciudadanos. Secretariado de Publicaciones de la Universidad de Sevilla. Sevilla 2010 964  Se utiliza el interpolador TopoToRaster para la interpolación la zona intermareal y obtener el MDT raster: MDT intermareal (figura 7). Figura 5. Zona intermareal interpolada. 4.- Finalmente, se realiza un mosaico (Fig.6) de los tres MDT corregidos a Alicante: MDT de Andalucía (altimetria), MDT batimetría (Batimetría) y zona intermareal (MDT intermareal). Figura 7. MDT final resultante del mosaico Km 01234 µ Km 01234 µ