La colocación es un procedimiento utilizado en teledetección para hacer coincidir mediciones de dos o más instrumentos diferentes. Esto se hace por dos razones principales: con fines de validación al comparar mediciones de la misma variable, y para relacionar mediciones de dos variables diferentes, ya sea para realizar recuperaciones o para predicción. En el segundo caso, los datos se introducen posteriormente en algún tipo de método estadístico inverso , como una red neuronal artificial , un algoritmo de clasificación estadística , un estimador de núcleo o un mínimo cuadrado lineal . En principio, la mayoría de los problemas de colocación se pueden resolver mediante una búsqueda del vecino más cercano , pero en la práctica hay muchas otras consideraciones involucradas y el mejor método es altamente específico para la comparación particular de instrumentos. Aquí abordamos algunas de las consideraciones más importantes junto con ejemplos específicos.
Hay al menos dos consideraciones principales al realizar colocaciones. El primero es el patrón de muestreo del instrumento. Las mediciones pueden ser densas y regulares, como las de un instrumento satelital de escaneo transversal. En este caso, puede ser apropiada alguna forma de interpolación . Por otro lado, las mediciones pueden ser escasas, como una campaña de campo única diseñada para algún ejercicio de validación particular. La segunda consideración es la huella del instrumento, que puede variar desde algo cercano a una medición puntual, como la de una radiosonda , o puede tener varios kilómetros de diámetro, como la de un radiómetro de microondas montado en un satélite. En el último caso, es apropiado tener en cuenta el patrón de antena del instrumento al hacer comparaciones con otro instrumento que tenga una huella más pequeña y un muestreo más denso, es decir, varias mediciones de un instrumento encajarán en la huella del otro.
Así como el instrumento tiene una huella espacial, también tendrá una huella temporal, a menudo llamada tiempo de integración. Si bien el tiempo de integración suele ser inferior a un segundo, lo que para aplicaciones meteorológicas es esencialmente instantáneo, hay muchos casos en los que algún tipo de promedio de tiempo puede facilitar considerablemente el proceso de colocación.
Las colocaciones deberán examinarse basándose tanto en las escalas de tiempo como de duración del fenómeno de interés. Esto facilitará aún más el proceso de colocación, ya que los datos de teledetección y otras mediciones casi siempre se agrupan de alguna manera. Ciertos fenómenos atmosféricos, como las nubes o la convección, son bastante transitorios, por lo que no es necesario considerar colocaciones con un error temporal de más de una hora aproximadamente. El hielo marino, por otra parte, se mueve y evoluciona con bastante lentitud, por lo que las mediciones separadas por hasta un día o más podrían seguir siendo útiles.
Los satélites que más nos preocupan son los que tienen una órbita polar cercana a la Tierra , ya que los satélites geoestacionarios ven el mismo punto durante toda su vida. El diagrama muestra mediciones de instrumentos AMSU-B montados en tres satélites durante un período de 12 horas. Esto ilustra tanto la trayectoria orbital como el patrón de escaneo que se ejecuta transversalmente. Dado que la órbita de un satélite es determinista , salvo maniobras orbitales , podemos predecir la ubicación del satélite en un momento determinado y, por extensión, la ubicación de los píxeles de medición. En teoría, las colocaciones se pueden realizar invirtiendo las ecuaciones determinantes a partir del período de tiempo deseado. En la práctica, los datos parcialmente procesados (generalmente denominados nivel 1b, 1c o nivel 2) contienen las coordenadas de cada uno de los píxeles de medición y es común simplemente enviar estas coordenadas a una búsqueda del vecino más cercano. Como se mencionó anteriormente, los datos satelitales siempre se agrupan de alguna manera. Como mínimo, los datos se organizarán en franjas que se extenderán de un polo a otro. Las franjas se etiquetarán por período de tiempo y se conocerá la ubicación aproximada.
Las radiosondas son particularmente importantes para los estudios de colocación porque miden variables atmosféricas con mayor precisión y más directamente que los satélites u otros instrumentos de detección remota. Además, las muestras de radiosonda son efectivamente mediciones puntuales instantáneas. Un problema con las radiosondas transportadas por globos meteorológicos es la deriva del globo. En [1] esto se maneja promediando todos los píxeles del satélite dentro de un radio de 50 km desde el lanzamiento del globo.
Si se utilizan datos de sonda de alta resolución, que normalmente tienen una frecuencia de muestreo constante o incluyen el tiempo de medición, entonces se puede rastrear el movimiento lateral a partir de los datos del viento. Incluso con datos de baja resolución, el movimiento aún se puede aproximar suponiendo una velocidad de ascenso constante. Excepto un poco hacia el final, el ascenso lineal se puede ver claramente en la figura de arriba. Podemos demostrar que la velocidad de ascenso de un globo viene dada por la siguiente ecuación
donde g es la aceleración gravitacional, k relaciona la altura, h , y el área de superficie, A , del globo con su volumen: V = khA ; R s es la "constante de gas" equivalente del globo, R a es la constante de gas del aire y c D es el coeficiente de resistencia del globo. Sustituyendo algunos valores sensibles para cada una de las constantes, k =1. (el globo es un cilindro perfecto), h =2. m, c D = 1. y Ra es la constante del gas del helio, arroja una velocidad de ascenso de 4,1 m/s. Compárese esto con los valores mostrados en el histograma que recopila todos los lanzamientos de radiosondas desde el buque de investigación Polarstern durante un período de once años entre 1992 y 2003.
Para datos cuadriculados, como datos de asimilación o reanálisis , la interpolación es probablemente el método más apropiado para realizar cualquier tipo de comparación. Un punto específico tanto en posición física como en tiempo es fácil de ubicar dentro de la cuadrícula y se realiza la interpolación entre los vecinos más cercanos. La interpolación lineal ( bilineal , trilineal , etc.) es la más común, aunque también se utiliza la cúbica, pero probablemente no valga la pena la sobrecarga computacional adicional. Si la variable de interés tiene una tasa de cambio relativamente suave (la temperatura es un buen ejemplo de esto porque tiene un mecanismo de difusión, transferencia radiativa , no disponible para otras variables atmosféricas), entonces la interpolación puede eliminar gran parte del error asociado con la colocación.
La interpolación también puede ser apropiada para muchos tipos de instrumentos satelitales, por ejemplo, un instrumento de escaneo transversal como Landsat . En [3] los datos derivados de la Unidad Avanzada de Sondeo por Microondas (AMSU) se interpolan (aunque no con fines de colocación) utilizando una ligera variación de la interpolación trilineal. Dado que las mediciones dentro de una única pista de escaneo se presentan en una cuadrícula aproximadamente rectangular, se puede realizar una interpolación bilineal. Al buscar la pista de escaneo superpuesta más cercana tanto hacia adelante como hacia atrás en el tiempo, las interpolaciones espaciales se pueden interpolar en el tiempo. Esta técnica funciona mejor con cantidades derivadas que con temperaturas de brillo sin procesar, ya que el ángulo de escaneo ya se habrá tenido en cuenta.
Para instrumentos con un patrón de muestreo más irregular, como el instrumento Advanced Microwave Scanning Radiometer-EOS (AMSR-E) que tiene un patrón de escaneo circular, necesitamos una forma más general de interpolación, como la estimación del núcleo . Un método comúnmente utilizado para este instrumento en particular, así como para SSM/I , es un promedio diario simple dentro de contenedores espaciales cuadriculados regularmente. [4]
Para colocar mediciones de un trazador atmosférico de vida media a larga con un segundo instrumento, las trayectorias pueden mejorar considerablemente la precisión. También simplifica un poco el análisis: se recorre una trayectoria tanto hacia adelante como hacia atrás desde el lugar de medición y entre la ventana de tiempo deseada. Tenga en cuenta que la ventana de tiempo aceptable ahora se ha alargado porque se elimina el error de los cambios inducidos por el transporte en el trazador: la vida útil del trazador sería una buena ventana para usar. Dado que las trayectorias proporcionan una ubicación para cada punto en el tiempo dentro de la ventana de tiempo, no es necesario verificar múltiples mediciones desde el segundo instrumento. Cada vez que se realiza dentro de la trayectoria se verifica el criterio de distancia, pero dentro de una ventana muy estrecha. Alternativamente, los tiempos exactos de las mediciones del segundo instrumento se interpolan dentro de la trayectoria. Sólo se utiliza el error de distancia más pequeño por debajo del umbral y, como consecuencia, el criterio de distancia puede reducirse.
Las colocaciones del espesor del hielo marino y las temperaturas de brillo tomadas durante la Campaña Pol-Ice son un excelente ejemplo, ya que ilustran muchos de los principios más importantes y demuestran la necesidad de tener en cuenta el caso individual. La campaña Pol-Ice se llevó a cabo en el Báltico Norte en marzo de 2007 como parte del proyecto SMOS-Ice en preparación para el lanzamiento del satélite de humedad del suelo y salinidad del océano . Debido a la baja frecuencia del instrumento SMOS, se espera que proporcione información sobre el espesor del hielo marino, por lo que la campaña comprendió mediciones tanto del espesor del hielo marino como de la temperatura de brillo emitida. Las temperaturas de brillo se midieron con el radiómetro de microondas de banda L EMIRAD [5] que se lleva a bordo de un avión. El espesor del hielo se midió con el medidor de espesor de hielo EM Bird, transportado por un helicóptero. El EM Bird mide el espesor del hielo con una combinación de mediciones de inductancia para determinar la ubicación de la interfaz hielo-agua y un altímetro láser para medir la altura de la superficie del hielo. [6] El mapa de arriba muestra las trayectorias de vuelo de ambos instrumentos, que fueron aproximadamente coincidentes pero obviamente sujetas a errores del piloto.
Dado que las trayectorias de vuelo de ambos aviones eran aproximadamente lineales, el primer paso en el proceso de colocación fue convertir todos los vuelos coincidentes a coordenadas cartesianas, siendo el eje x la distancia lateral y el eje y la distancia transversal. De esta manera, las colocaciones se pueden realizar de dos maneras: de manera tosca, haciendo coincidir solo las distancias x y, más precisamente, haciendo coincidir ambas coordenadas.
Más importante aún, el tamaño de la huella del radiómetro es muchas veces mayor que el del medidor EM Bird. La figura de la izquierda muestra la función de respuesta de la antena del radiómetro. El ancho total a la mitad del máximo es de 31 grados. [5] Dado que el avión volaba a aproximadamente 500 m, esto se traduce en una huella de 200 mo más. Mientras tanto, el tamaño de la huella del EM Bird fue de aproximadamente 40 m con un espacio entre muestras de sólo 2 a 4 m. [6] En lugar de mirar a los vecinos más cercanos, lo que habría producido malos resultados, se realizó un promedio ponderado de las mediciones de espesor para cada medición del radiómetro. Los pesos se calcularon basándose en la función de respuesta del radiómetro, que es casi una gaussiana perfecta hasta aproximadamente 45 grados. Se podrían excluir puntos según la distancia a lo largo de la ruta de vuelo. Para la validación de los cálculos del modelo directo de emisividad del hielo marino, esto se perfeccionó aún más realizando un cálculo de emisividad para cada medición de espesor y promediando la huella del radiómetro. [7]
La siguiente figura ilustra las ubicaciones de medición relativas de cada uno de los instrumentos utilizados en la campaña Pol-Ice. Se muestran dos pasos elevados: uno desde el avión que lleva el radiómetro EMIRAD y otro desde el helicóptero que lleva el instrumento EM Bird. El eje x está a lo largo de la línea de la trayectoria de vuelo. Las huellas de EMIRAD se dibujan con líneas, las mediciones de inductancia de EM Bird se representan con círculos y las mediciones de LIDAR con puntos.
