Page 8 - MAPPING_182
P. 8

Abderrahim Nemmaoui, Manuel Ángel Aguilar , Antonio Novelli, Mª Carmen Vicente , Fernando José Aguilar , Malgorzata Betlej, Piotr Cichón




          al., 2015) la segmentación se realizó de forma manual, evi-  (MS) respectivamente. Este satélite fue el primer sensor co-
          tando tener que encontrar los parámetros óptimos en la fase   mercial en captar 8 bandas en su imagen MS. Estas bandas
          de segmentación.                                     son: coastal (C, 400–450 nm), blue (B, 450–510 nm), green (G,
             Este trabajo se centra en la fase de segmentación de in-  510–580 nm), yellow (Y, 585–625 nm), red (R, 630–690 nm),
          vernaderos de plástico en entornos OBIA, intentando estimar   red edge (RE, 705–745 nm), near infrared-1 (NIR1, 760–895
          los parámetros óptimos (i.e., Escala, Forma y Compacidad)   nm) y near infrared-2 (NIR2, 860–1040 nm).
          del algoritmo de segmentación multi-resolución incluido   Una imagen de WV2, tomada el 30 de septiembre de
          en eCognition a partir de una ortoimagen multiespectral del   2013 sobre la zona de estudio, fue usada para este trabajo.
          satélite WorldView-2. Se intentará también relacionar la bon-  Fue recibida en formato  Ortho  Ready  Standard  Level-2A
          dad de la segmentación con la precisión en la fase posterior   (ORS2A), conteniendo las imágenes PAN y MS. Esta imagen
          de clasificación de objetos.                         de satélite presentó un ángulo off-nadir de 11.8° y un 0% de
                                                               nubes. El tamaño de píxel de las imágenes PAN y MS fue de
                                                               0.4 m y 1.6 m respectivamente. El formato ORS2A de WV2
          2. ZONA DE ESTUDIO Y DATOS DE                        presenta correcciones radiométricas y geométricas, y ya está
                                                               georreferenciada respecto a una proyección cartográfica.
          PARTIDA                                              Incluye los correspondientes datos de coeficientes raciona-
                                                               les para poder reproducir el modelo de cámara. El producto
             2.1. Zona de estudio                              empleado tenía un rango dinámico de 11 bit.
             El trabajo fue realizado en el poniente almeriense, en el   A partir de las imágenes PAN y MS se obtuvo una imagen
          lugar del mundo con mayor concentración de invernade-  pansharpened o fusionada con 0.4 m GSD y conteniendo
          ros conocido como «Mar de Plástico». La zona de estudio   todas las bandas RGB de la imagen MS. Para tal fin se utili-
          comprende un área rectangular de 8 000 ha (8 km x 10 km)   zó el módulo PANSHARP del  software Geomática v. 2014
          centrado en coordenadas geográficas WGS84 36.7824°N y   (PCI Geomatics, Richmond Hill, Canadá). Para obtener las
          2.6867°W (Figura 1). Comprende parte de los Términos Mu-  ortoimágenes finales se emplearon las coordenadas de 7
          nicipales de Vícar, El Ejido, Roquetas de Mar y La Mojonera,   puntos de control de campo (ground control points, GCPs) y
          siendo una zona principalmente de uso agrícola, presen-  32 puntos de comprobación independientes (independent
          tando parcelas con y sin invernaderos, además de algunos   check points, ICPs) que fueron obtenidas mediante GPS dife-
          núcleos de población.                                rencial (GPS Topcon HiPer PRO station) trabajando en modo
                                                               RTK. Estos puntos de campo fueron medidos con referencia
             2.2.  Imagen de satélite VHR                      al datum European Terrestrial Reference System 1989 (ETRS89)
             WorldView-2 (WV2) es un satélite VHR lanzado en Oc-  y con proyección UTM Zona 30. Una ortoimagen fusionada
          tubre de 2009. Este sensor es capaz de adquirir imágenes   con tres bandas (RGB) y con 0.4 m GSD fue generada usando
          ópticas de la superficie terrestre con un tamaño de píxel o   funciones racionales refinadas de grado cero y un modelo di-
          Ground Sample Distance (GSD) de 0.46 m y 1.84 m en su posi-  gital de elevaciones (MDE) con 10 m de resolución espacial y
          ción más nadiral en modo pancromático (PAN) y multiesptral   1.34 m de error vertical de la Junta de Andalucía. La precisión
                                                                             planimétrica medida como error cuadrático
                                                                             medio (root mean square error, RMSE) en los
                                                                             32 ICPs de esta ortoimagen resultó ser de
                                                                             0.59 m.
                                                                               Por otra parte, se generó también una
                                                                             ortoimagen MS con 1.6 m GSD con las 8
                                                                             bandas  originales  de  WV2.  Se  usaron  los
                                                                             mismos 7 GCPs, el mismo modelo de sen-
                                                                             sor y MDE empleados para la obtención de
                                                                             la  ortoimagen  pansharpened o fusionada.
                                                                             Esta ortoimagen fue la empelada para los
                                                                             ensayos OBIA realizados en este trabajo.
                                                                             Por ello, la imagen MS de WV2 original fue
                                                                             corregida atmosféricamente aplicándole
                                                                             el módulo ATCOR (atmospheric correction)
          Figura 1. Situación de la zona de estudio sobre una imagen Landsat 8 de agosto de 2013. Sistema de   incluido  en  Geomatica  v.  2014.  Mediante
          coordenadas ETRS89 UTM Zona 30N                                    esta corrección atmosférica absoluta trans-





          Pág. 06                                    REVISTA MAPPING  VOL.26   Nº182   MARZO-ABRIL 2017  ISSN: 1131-9100
   3   4   5   6   7   8   9   10   11   12   13