Articulo de referencia

radar de apertura sintética

Esta imagen de radar, obtenida por el radar SIR-C/X-SAR a bordo del transbordador espacial Endeavour, muestra el volcán Teide . La ciudad de Santa Cruz de Tenerife se distingue ...

Esta imagen de radar, obtenida por el radar SIR-C/X-SAR a bordo del transbordador espacial Endeavour, muestra el volcán Teide . La ciudad de Santa Cruz de Tenerife se distingue como la zona púrpura y blanca en el extremo inferior derecho de la isla. Los flujos de lava en el cráter de la cumbre aparecen en tonos verdes y marrones, mientras que las zonas de vegetación se observan como áreas púrpuras, verdes y amarillas en los flancos del volcán.

El radar de apertura sintética ( SAR ) es un tipo de radar que se utiliza para crear imágenes bidimensionales o reconstrucciones tridimensionales de objetos, como paisajes. [ 1 ] El SAR utiliza el movimiento de la antena de radar sobre una región objetivo para proporcionar una resolución espacial más fina que los radares convencionales de exploración de haz estacionarios. El SAR se monta normalmente en una plataforma móvil, como un avión o una nave espacial, y tiene su origen en una forma avanzada de radar aerotransportado de barrido lateral (SLAR). La distancia que recorre el dispositivo SAR sobre un objetivo durante el período en que la escena objetivo está iluminada crea la gran apertura sintética de la antena (el tamaño de la antena). Por lo general, cuanto mayor sea la apertura, mayor será la resolución de la imagen, independientemente de si la apertura es física (una antena grande) o sintética (una antena móvil); esto permite que el SAR cree imágenes de alta resolución con antenas físicas relativamente pequeñas. Para un tamaño y orientación de antena fijos, los objetos más lejanos permanecen iluminados durante más tiempo; por lo tanto, el SAR tiene la propiedad de crear aperturas sintéticas más grandes para objetos más distantes, lo que da como resultado una resolución espacial constante en un rango de distancias de visualización.

Para crear una imagen SAR, se transmiten pulsos sucesivos de ondas de radio para "iluminar" una escena objetivo, y el eco de cada pulso se recibe y registra. Los pulsos se transmiten y los ecos se reciben mediante una única antena formadora de haces , con longitudes de onda desde un metro hasta varios milímetros. A medida que el dispositivo SAR a bordo de la aeronave o nave espacial se mueve, la posición de la antena con respecto al objetivo cambia con el tiempo. El procesamiento de la señal de los ecos de radar registrados sucesivamente permite combinar las grabaciones de estas múltiples posiciones de antena. Este proceso forma la apertura sintética de la antena y permite la creación de imágenes de mayor resolución que la que sería posible con una antena física determinada. [ 2 ] [ 3 ]

Motivación y aplicaciones

La superficie de Venus , tal como fue captada por la sonda Magellan mediante radar de apertura sintética (SAR), coloreada con falso color .

El SAR es capaz de realizar teledetección de alta resolución , independientemente de la altitud de vuelo y de las condiciones meteorológicas, [ 4 ] ya que puede seleccionar frecuencias para evitar la atenuación de la señal causada por el clima. El SAR tiene capacidad de obtención de imágenes tanto de día como de noche, ya que proporciona iluminación. [ 5 ] [ 6 ] [ 7 ]

Las imágenes SAR tienen amplias aplicaciones en la teledetección y el mapeo de superficies de la Tierra y otros planetas. Las aplicaciones de SAR son numerosas. Algunos ejemplos incluyen topografía, oceanografía, glaciología , geología (por ejemplo, discriminación de terreno e imágenes del subsuelo). SAR también se puede utilizar en silvicultura para determinar la altura del bosque, la biomasa y la deforestación. El monitoreo de volcanes y terremotos utiliza interferometría diferencial . SAR también se puede aplicar para monitorear la estabilidad de infraestructura civil como puentes. [ 8 ] SAR es útil en el monitoreo ambiental como derrames de petróleo, inundaciones, [ 9 ] [ 10 ] crecimiento urbano, [ 11 ] y vigilancia militar: incluyendo política estratégica y evaluación táctica. [ 7 ] SAR se puede implementar como SAR inverso al observar un objetivo en movimiento durante un tiempo sustancial con una antena estacionaria.

Principio básico

Principio básico

Un radar de apertura sintética es un radar de imágenes montado sobre una plataforma móvil. [ 12 ] El SAR es una técnica Doppler . Se basa en el hecho de que "las reflexiones de radar de objetos discretos en un campo de haz de radar pasante tienen cada una un pequeño desplazamiento Doppler, o de velocidad, con respecto a la antena". [ 13 ]

Carl Wiley, que trabajaba en Goodyear, Arizona (que más tarde se convirtió en Goodyear Aerospace y, finalmente, en Lockheed Martin Corporation) en 1951, sugirió el principio de que, debido a que cada objeto en el haz del radar tiene una velocidad ligeramente diferente con respecto a la antena, cada objeto tendrá su propio desplazamiento Doppler. Un análisis preciso de la frecuencia de las reflexiones del radar permitirá, por lo tanto, la construcción de una imagen detallada. [ 14 ]

Para materializar este concepto, se transmiten ondas electromagnéticas secuencialmente, se recogen los ecos y la electrónica del sistema digitaliza y almacena los datos para su posterior procesamiento. Como la transmisión y la recepción ocurren en momentos diferentes, se corresponden con distintas posiciones pequeñas. La combinación ordenada de las señales recibidas crea una apertura virtual mucho mayor que el ancho físico de la antena. De ahí proviene el término «apertura sintética», que le confiere la propiedad de un radar de imágenes. [ 7 ] La dirección de alcance es perpendicular a la trayectoria de vuelo y perpendicular a la dirección de acimut, también conocida como dirección longitudinal , ya que coincide con la posición del objeto dentro del campo de visión de la antena.

El procesamiento 3D se realiza en dos etapas. La dirección de acimut y alcance se enfoca para la generación de imágenes 2D (acimut-alcance) de alta resolución, después de lo cual se utiliza un modelo digital de elevación (MDE) [ 15 ] [ 16 ] para medir las diferencias de fase entre imágenes complejas, que se determina a partir de diferentes ángulos de visión para recuperar la información de altura. Esta información de altura, junto con las coordenadas de acimut-alcance proporcionadas por el enfoque SAR 2D, da la tercera dimensión, que es la elevación. [ 5 ] El primer paso solo requiere algoritmos de procesamiento estándar, [ 16 ] para el segundo paso, se utiliza un preprocesamiento adicional como el corregistro de imágenes y la calibración de fase. [ 5 ] [ 17 ]

Además, se pueden utilizar múltiples líneas base para extender la imagen 3D a la dimensión temporal . La imagen SAR 4D y multi-D permite obtener imágenes de escenarios complejos, como áreas urbanas, y ha mejorado el rendimiento con respecto a las técnicas interferométricas clásicas, como la interferometría de dispersores persistentes (PSI). [ 18 ]

Algoritmo

Los algoritmos SAR modelan la escena como un conjunto de objetivos puntuales que no interactúan entre sí (la aproximación de Born ).

Aunque los detalles de los distintos algoritmos SAR difieren, el procesamiento SAR en cada caso consiste en la aplicación de un filtro adaptado a los datos brutos, para cada píxel de la imagen de salida, donde los coeficientes del filtro adaptado son la respuesta de un único objetivo puntual aislado. [ 19 ] En los inicios del procesamiento SAR, los datos brutos se registraban en película y el posprocesamiento mediante filtro adaptado se realizaba ópticamente utilizando lentes de forma cónica, cilíndrica y esférica. El algoritmo Range-Doppler es un ejemplo de un enfoque más reciente.

Enfoques de estimación espectral existentes

El radar de apertura sintética determina la reflectividad 3D a partir de datos SAR medidos. Básicamente, se trata de una estimación espectral, ya que para una celda específica de una imagen, las mediciones SAR de valor complejo de la pila de imágenes SAR son una versión muestreada de la transformada de Fourier de la reflectividad en la dirección de elevación, pero la transformada de Fourier es irregular. [ 20 ] Por lo tanto, las técnicas de estimación espectral se utilizan para mejorar la resolución y reducir el moteado en comparación con los resultados de las técnicas convencionales de imágenes SAR de transformada de Fourier. [ 21 ]

Métodos no paramétricos

FFT

La FFT ( Transformada Rápida de Fourier , es decir, periodograma o filtro adaptado ) es uno de estos métodos, que se utiliza en la mayoría de los algoritmos de estimación espectral, y existen muchos algoritmos rápidos para calcular la transformada discreta de Fourier multidimensional . El álgebra computacional de matriz de núcleo de Kronecker [ 22 ] es un algoritmo popular que se utiliza como una nueva variante de los algoritmos FFT para el procesamiento en sistemas de radar de apertura sintética (SAR) multidimensionales. Este algoritmo utiliza un estudio de las propiedades teóricas de los conjuntos de indexación de datos de entrada/salida y grupos de permutaciones.

Se utiliza una rama del álgebra lineal multidimensional finita para identificar similitudes y diferencias entre diversas variantes del algoritmo FFT y para crear nuevas variantes. Cada cálculo de DFT multidimensional se expresa en forma matricial. La matriz DFT multidimensional, a su vez, se descompone en un conjunto de factores, denominados primitivas funcionales, que se identifican individualmente con un diseño computacional de software/hardware subyacente. [ 7 ]

La implementación de la FFT es esencialmente una realización del mapeo del marco matemático mediante la generación de variantes y la ejecución de operaciones matriciales. El rendimiento de esta implementación puede variar de una máquina a otra, y el objetivo es identificar en qué máquina funciona mejor. [ 23 ]

Ventajas
  • Las propiedades aditivas de la teoría de grupos de los conjuntos de indexación de entrada/salida multidimensionales se utilizan para las formulaciones matemáticas, por lo tanto, es más fácil identificar el mapeo entre las estructuras de computación y las expresiones matemáticas, lo que resulta mejor que los métodos convencionales. [ 24 ]
  • El lenguaje del álgebra CKA ayuda al desarrollador de la aplicación a comprender cuáles son las variantes de FFT más eficientes computacionalmente, reduciendo así el esfuerzo computacional y mejorando su tiempo de implementación. [ 24 ] [ 25 ]
Desventajas
  • La FFT no puede separar sinusoides con frecuencias cercanas. Si la periodicidad de los datos no coincide con la FFT, se observan efectos de borde. [ 23 ]

Método Capon

El método espectral de Capon , también llamado método de mínima varianza, es una técnica de procesamiento de matrices multidimensionales. [ 26 ] Es un método no paramétrico basado en la covarianza, que utiliza un enfoque de banco de filtros adaptado y sigue dos pasos principales:

  1. Pasar los datos a través de un filtro de paso de banda 2D con frecuencias centrales variables (ω1,ω2{\displaystyle \omega _ {1}, \ omega _ {2}}).
  2. Estimando la potencia en (ω1,ω2{\displaystyle \omega _ {1}, \ omega _ {2}}) para todosω1[0,2π),ω2[0,2π){\displaystyle \omega _{1}\in [0,2\pi ),\omega _{2}\in [0,2\pi )}de interés a partir de los datos filtrados.

El filtro de paso de banda adaptativo Capon está diseñado para minimizar la potencia de la salida del filtro, así como para dejar pasar las frecuencias (ω1,ω2{\displaystyle \omega _ {1}, \ omega _ {2}}) sin ninguna atenuación, es decir, para satisfacer, para cada (ω1,ω2{\displaystyle \omega _ {1}, \ omega _ {2}}),

minhhω1,ω2Rhω1,ω2{\displaystyle \min _{h}h_{\omega _{1},\omega _{2}}^{*}Rh_{\omega _{1},\omega _{2}}}sujeto ahω1,ω2aω1,ω2=1,{\displaystyle h_{\omega _{1},\omega _{2}}^{*}a_{\omega _{1},\omega _{2}}=1,}

donde R es la matriz de covarianza ,hω1,ω2{\ Displaystyle h _ {\ omega _ {1}, \ omega _ {2}} ^ {*}}es la transpuesta conjugada compleja de la respuesta impulsional del filtro FIR,aω1,ω2{\ Displaystyle a _ {\ omega _ {1}, \ omega _ {2}}}es el vector de Fourier 2D, definido comoaω1,ω2aω1aω2{\displaystyle a_{\omega _{1},\omega _{2}}\triangleq a_{\omega _{1}}\otimes a_{\omega _{2}}},{\displaystyle \otimes }denota el producto de Kronecker. [ 26 ]

Por lo tanto, transmite una sinusoide 2D a una frecuencia determinada sin distorsión, minimizando la varianza del ruido de la imagen resultante. El objetivo es calcular la estimación espectral de manera eficiente. [ 26 ]

La estimación espectral se da como

Sω1,ω2=1aω1,ω2R1aω1,ω2,{\displaystyle S_{\omega _{1},\omega _{2}}={\frac {1}{a_{\omega _{1},\omega _{2}}^{*}R^{-1}a_{\omega _{1},\omega _{2}}}},}

donde R es la matriz de covarianza yaω1,ω2{\displaystyle a_{\omega _{1},\omega _{2}}^{*}}es la transpuesta conjugada compleja bidimensional del vector de Fourier. El cálculo de esta ecuación sobre todas las frecuencias consume mucho tiempo. Se observa que el estimador de Capon hacia adelante y hacia atrás proporciona una mejor estimación que el enfoque clásico de Capon solo hacia adelante. La razón principal es que, mientras que el Capon hacia adelante y hacia atrás utiliza los vectores de datos hacia adelante y hacia atrás para obtener la estimación de la matriz de covarianza, el Capon solo hacia adelante utiliza únicamente los vectores de datos hacia adelante para estimar la matriz de covarianza. [ 26 ]

Ventajas
  • Capon puede proporcionar estimaciones espectrales más precisas con lóbulos laterales mucho más bajos y picos espectrales más estrechos que el método de transformada rápida de Fourier (FFT). [ 27 ]
  • El método de Capon puede proporcionar una resolución mucho mejor.
Desventajas
  • La implementación requiere el cálculo de dos tareas intensivas: la inversión de la matriz de covarianza R y la multiplicación por laaω1,ω2{\displaystyle a_{\omega _{1},\omega _{2}}}matriz, que debe hacerse para cada punto(ω1,ω2){\displaystyle \left(\omega _{1},\omega _{2}\right)}. [ 5 ]

Método APES

El método APES (estimación de amplitud y fase) es también un método de banco de filtros adaptados, que supone que los datos del historial de fase son una suma de sinusoides 2D con ruido.

El estimador espectral APES tiene una interpretación de filtrado de 2 pasos:

  1. Procesamiento de datos a través de un banco de filtros pasabanda FIR con frecuencia central variable.ω{\displaystyle \omega }.
  2. Obtención de la estimación del espectro paraω[0,2π){\displaystyle \omega \in [0,2\pi )}a partir de los datos filtrados. [ 28 ]

Empíricamente, el método APES produce picos espectrales más amplios que el método Capon, pero estimaciones espectrales más precisas para la amplitud en SAR. [ 29 ] En el método Capon, aunque los picos espectrales son más estrechos que en el APES, los lóbulos laterales son más altos. Como resultado, se espera que la estimación de la amplitud sea menos precisa con el método Capon que con el método APES. El método APES requiere aproximadamente 1,5 veces más cálculos que el método Capon. [ 30 ]

Ventajas
  • El filtrado reduce el número de muestras disponibles, pero cuando se diseña tácticamente, el aumento de la relación señal-ruido (SNR) en los datos filtrados compensará esta reducción, y la amplitud de un componente sinusoidal con frecuenciaω{\displaystyle \omega }se puede estimar con mayor precisión a partir de los datos filtrados que a partir de la señal original. [ 31 ]
Desventajas
  • La matriz de autocovarianza es mucho mayor en 2D que en 1D, por lo tanto, está limitada por la memoria disponible. [ 7 ]

método SAMV

El método SAMV es un algoritmo de reconstrucción de señales dispersas sin parámetros. Logra una superresolución y es robusto ante señales altamente correlacionadas. Su nombre enfatiza su fundamento en el criterio de varianza mínima asintótica (AMV). Es una herramienta poderosa para la recuperación de las características de amplitud y frecuencia de múltiples fuentes altamente correlacionadas en entornos complejos (por ejemplo, número limitado de instantáneas, baja relación señal -ruido ). Entre sus aplicaciones se incluyen la formación de imágenes de radar de apertura sintética y la localización de diversas fuentes.

Ventajas

El método SAMV es capaz de lograr una resolución superior a la de algunos métodos paramétricos establecidos, como MUSIC , especialmente con señales altamente correlacionadas.

Desventajas

La complejidad computacional del método SAMV es mayor debido a su procedimiento iterativo.

Métodos de descomposición de subespacios paramétricos

Método de vectores propios

Este método de descomposición de subespacio separa los autovectores de la matriz de autocovarianza en aquellos que corresponden a señales y a ruido. [ 7 ] La amplitud de la imagen en un punto (ωincógnita,ωy{\displaystyle \omega _{x},\omega _{y}}) viene dado por:

ϕ^miV(ωincógnita,ωy)=1WH(ωincógnita,ωy)(desorden1λivi_vi_H)W(ωincógnita,ωy){\displaystyle {\hat {\phi }}_{EV}\left(\omega _{x},\omega _{y}\right)={\frac {1}{W^{\mathsf {H}}\left(\omega _{x},\omega _{y}\right)\left(\sum _{\text{clutter}}{\frac {1}{\lambda _{i}}}{\underline {v_{i}}}\,{\underline {v_{i}}}^{\mathsf {H}}\right)W\left(\omega _{x},\omega _{y}\right)}}}

dóndeϕ^miV{\displaystyle {\hat {\phi }}_{EV}}es la amplitud de la imagen en un punto(ωincógnita,ωy){\displaystyle \left(\omega _{x},\omega _{y}\right)},vi_{\displaystyle {\underline {v_{i}}}}es la matriz de coherencia yvi_H{\displaystyle {\underline {v_{i}}}^{\mathsf {H}}}es el hermitiano de la matriz de coherencia,1λi{\displaystyle {\frac {1}{\lambda _{i}}}}es el inverso de los valores propios del subespacio de desorden,W(ωincógnita,ωy){\displaystyle W\left(\omega _{x},\omega _{y}\right)}son vectores definidos como [ 7 ]

W(ωincógnita,ωy)=[1exp(jωincógnita)exp(j(METRO1)ωincógnita)][1exp(jωy)exp(j(METRO1)ωy)]{\displaystyle W\left(\omega _{x},\omega _{y}\right)=\left[1\exp \left(-j\omega _{x}\right)\ldots \exp \left(-j(M-1)\omega _{x}\right)\right]\otimes \left[1\exp \left(-j\omega _{y}\right)\ldots \exp \left(-j(M-1)\omega _{y}\right)\right]}

donde ⊗ denota el producto de Kronecker de los dos vectores.

Ventajas
  • Muestra las características de la imagen con mayor precisión. [ 7 ]
Desventajas
  • Alta complejidad computacional. [ 17 ]

Método MUSICAL

MUSIC detecta frecuencias en una señal mediante la descomposición en valores propios de la matriz de covarianza de un vector de datos de las muestras obtenidas de la señal recibida. Cuando todos los vectores propios se incluyen en el subespacio de ruido (orden del modelo = 0), el método EV se vuelve idéntico al método Capon. Por lo tanto, la determinación del orden del modelo es fundamental para el funcionamiento del método EV. El valor propio de la matriz R determina si su vector propio correspondiente pertenece al subespacio de ruido o al de la señal. [ 7 ]

El método MUSIC se considera de bajo rendimiento en aplicaciones SAR. Este método utiliza una constante en lugar del subespacio de ruido. [ 7 ]

En este método, el denominador se iguala a cero cuando una señal sinusoidal correspondiente a un punto en la imagen SAR se alinea con uno de los vectores propios del subespacio de la señal, que es el pico en la estimación de la imagen. Por lo tanto, este método no representa con precisión la intensidad de dispersión en cada punto, sino que muestra los puntos particulares de la imagen. [ 7 ] [ 32 ]

Ventajas
  • La MÚSICA blanquea o ecualiza los autovalores del ruido. [ 21 ]
Desventajas
  • Pérdida de resolución debido a la operación de promediado. [ 12 ]

Algoritmo de retroproyección

El algoritmo de retroproyección tiene dos métodos: retroproyección en el dominio del tiempo y retroproyección en el dominio de la frecuencia . La retroproyección en el dominio del tiempo tiene más ventajas que la retroproyección en el dominio de la frecuencia y, por lo tanto, es la preferida. La retroproyección en el dominio del tiempo forma imágenes o espectros al hacer coincidir los datos adquiridos del radar con lo que espera recibir. Puede considerarse un filtro adaptado ideal para radares de apertura sintética. No es necesario un paso de compensación de movimiento diferente debido a su capacidad para manejar el movimiento/muestreo no ideal. También puede utilizarse para diversas geometrías de imagen. [ 33 ]

Ventajas

  • Es invariante al modo de imagen : lo que significa que utiliza el mismo algoritmo independientemente del modo de imagen presente, mientras que los métodos del dominio de la frecuencia requieren cambios dependiendo del modo y la geometría. [ 33 ]
  • El aliasing azimutal ambiguo suele producirse cuando las frecuencias superan los requisitos de muestreo espacial de Nyquist. El aliasing inequívoco se produce en geometrías oblicuas donde el ancho de banda de la señal no excede los límites de muestreo, pero ha sufrido un "envolvimiento espectral". El algoritmo de retroproyección no se ve afectado por este tipo de efectos de aliasing. [ 33 ]
  • Coincide con el filtro espacio-temporal: utiliza la información sobre la geometría de la imagen para producir un filtro adaptado que varía píxel a píxel y aproxima la señal de retorno esperada. Esto suele proporcionar una compensación de la ganancia de la antena. [ 33 ]
  • En relación con la ventaja anterior, el algoritmo de retroproyección compensa el movimiento. Esto se convierte en una ventaja en áreas con altitudes bajas. [ 33 ]

Desventajas

  • El algoritmo de retroproyección requiere mayor esfuerzo computacional en comparación con otros métodos en el dominio de la frecuencia.
  • Requiere un conocimiento muy preciso de la geometría de la imagen. [ 33 ]

Aplicación: radar de apertura sintética en órbita geoestacionaria (GEO-SAR)

En GEO-SAR, para centrarse especialmente en la trayectoria relativa en movimiento, el algoritmo de retroproyección funciona muy bien. Utiliza el concepto de procesamiento de acimut en el dominio del tiempo. Para la geometría satélite-tierra, GEO-SAR desempeña un papel importante. [ 34 ]

El procedimiento de este concepto se desarrolla de la siguiente manera. [ 34 ]

  1. Los datos brutos adquiridos se segmentan o se dividen en subaperturas para simplificar la realización rápida del procedimiento.
  2. El rango de los datos se comprime a continuación, utilizando el concepto de "filtrado adaptado" para cada segmento/subepersión creado. Se define mediante:s(t,τ)=exp(j4πλR(t))desde(τ2doR(t)){\textstyle s(t,\tau )=\exp \left(-j\cdot {\frac {4\pi }{\lambda }}\cdot R(t)\right)\cdot \operatorname {sinc} \left(\tau -{\frac {2}{c}}\cdot R(t)\right)}donde τ es el tiempo de alcance, t es el tiempo azimutal, λ es la longitud de onda, c es la velocidad de la luz .
  3. La precisión en la "Curva de Migración de Rango" se logra mediante la interpolación de rango.
  4. La ubicación de los píxeles del terreno en la imagen depende del modelo geométrico satélite-tierra. La división de la cuadrícula se realiza ahora según el tiempo de acimut.
  5. Los cálculos para el " alcance oblicuo " (distancia entre el centro de fase de la antena y el punto en el suelo) se realizan para cada instante de tiempo azimutal utilizando transformaciones de coordenadas.
  6. La compresión de acimut se realiza después del paso anterior.
  7. Los pasos 5 y 6 se repiten para cada píxel, para cubrir cada píxel, y se lleva a cabo el procedimiento en cada subapertura.
  8. Por último, todas las subaperturas de la imagen creadas a lo largo del proceso se superponen entre sí y se genera la imagen HD final.

Comparación entre los algoritmos

Los métodos Capon y APES pueden proporcionar estimaciones espectrales más precisas con lóbulos laterales mucho más bajos y picos espectrales más estrechos que el método de transformada rápida de Fourier (FFT), que también es un caso especial de los enfoques de filtrado FIR. Se observa que, si bien el algoritmo APES produce picos espectrales ligeramente más anchos que el método Capon, el primero proporciona estimaciones espectrales generales más precisas que el segundo y el método FFT. [ 29 ]

El método FFT es rápido y sencillo, pero tiene lóbulos laterales más grandes. Capon tiene alta resolución, pero alta complejidad computacional. EV también tiene alta resolución y alta complejidad computacional. APES tiene mayor resolución, es más rápido que Capon y EV, pero alta complejidad computacional. [ 12 ]

El método MUSIC no es generalmente adecuado para imágenes SAR, ya que el blanqueamiento de los valores propios del ruido destruye las inhomogeneidades espaciales asociadas con el ruido del terreno u otra dispersión difusa en las imágenes SAR. Sin embargo, ofrece una mayor resolución de frecuencia en la densidad espectral de potencia (PSD) resultante que los métodos basados ​​en la transformada rápida de Fourier (FFT). [ 35 ]

El algoritmo de retroproyección es computacionalmente costoso. Resulta especialmente atractivo para sensores de banda ancha, gran angular y/o con aperturas coherentes largas y un movimiento fuera de trayectoria considerable. [ 36 ]

Operación multiestática

El sistema SAR requiere que se tomen capturas de eco en múltiples posiciones de antena. Cuantas más capturas se tomen (en diferentes ubicaciones de antena), más fiable será la caracterización del objetivo.

Se pueden obtener múltiples capturas moviendo una sola antena a diferentes ubicaciones, colocando varias antenas fijas en diferentes ubicaciones o mediante combinaciones de ambas.

La ventaja de una antena móvil única es que se puede colocar fácilmente en múltiples posiciones para generar diversas formas de onda monoestáticas. Por ejemplo, una antena montada en un avión capta múltiples señales por segundo mientras el avión vuela.

Las principales ventajas de utilizar varias antenas estáticas son que se puede caracterizar un objetivo en movimiento (siempre que la electrónica de captura sea lo suficientemente rápida), que no se necesita ningún vehículo ni maquinaria de movimiento, y que las posiciones de las antenas no tienen por qué derivarse de otra información, a veces poco fiable. (Un problema del SAR a bordo de un avión es conocer las posiciones precisas de las antenas durante el vuelo).

Para múltiples antenas estáticas, son posibles todas las combinaciones de capturas de formas de onda de radar monoestáticas y multiestáticas . Sin embargo, tenga en cuenta que no es ventajoso capturar una forma de onda para cada una de las dos direcciones de transmisión para un par de antenas dado, porque esas formas de onda serán idénticas. Cuando se utilizan múltiples antenas estáticas, el número total de formas de onda de eco únicas que se pueden capturar es

norte2+norte2{\displaystyle {\frac {N^{2}+N}{2}}}

donde N es el número de posiciones únicas de la antena.

Modos de escaneo

SAR aerotransportado en modo Stripmap

Ilustración del modo de funcionamiento del mapa de franjas SAR.

La antena permanece en una posición fija. Puede ser ortogonal a la trayectoria de vuelo, o puede estar ligeramente inclinada hacia adelante o hacia atrás. [ 7 ]

Cuando la apertura de la antena se desplaza a lo largo de la trayectoria de vuelo, se transmite una señal a una frecuencia igual a la frecuencia de repetición de pulsos (PRF). El límite inferior de la PRF está determinado por el ancho de banda Doppler del radar. La retrodispersión de cada una de estas señales se suma de forma conmutativa píxel a píxel para obtener la resolución azimutal precisa deseada en las imágenes de radar. [ 37 ]

SAR en modo foco

Representación del modo de imagen Spotlight

La apertura sintética del foco viene dada por [ 38 ].

Lsa=r0Δθa{\displaystyle Lsa=r_{0}\Delta \theta _{a}}[ 32 ]

dóndeΔθa{\displaystyle \Delta \theta _{a}}es el ángulo formado entre el principio y el final de la formación de imágenes, como se muestra en el diagrama de formación de imágenes puntuales yr0{\displaystyle r_{0}}es la distancia del alcance.

El modo de foco proporciona una mejor resolución, aunque para una zona terrestre más pequeña. En este modo, el haz de radar de iluminación se dirige continuamente a medida que la aeronave se mueve, de modo que ilumina la misma zona durante un período de tiempo más prolongado. Este modo no es un modo de imagen de franja continua tradicional; sin embargo, tiene una alta resolución azimutal. [ 32 ] En [ 39 ] se ofrece una explicación técnica del SAR de foco desde sus principios fundamentales.

Modo de escaneo SAR

Representación del modo de imagen ScanSAR

Al operar como un SAR en modo de escaneo, el haz de la antena barre periódicamente y, por lo tanto, cubre un área mucho mayor que los modos de foco y mapa de franjas. Sin embargo, la resolución azimutal se vuelve mucho menor que en el modo de mapa de franjas debido a la disminución del ancho de banda azimutal. Claramente se logra un equilibrio entre la resolución azimutal y el área de escaneo del SAR. [ 40 ] Aquí, la apertura sintética se comparte entre las subfranjas y no está en contacto directo dentro de una subfranja. Se requiere una operación de mosaico en las direcciones de azimut y alcance para unir las ráfagas azimutales y las subfranjas de alcance. [ 32 ]

  • ScanSAR hace que el haz de cobertura sea enorme.
  • La señal de acimut presenta muchas ráfagas.
  • La resolución azimutal es limitada debido a la duración de la ráfaga.
  • Cada objetivo contiene frecuencias variadas que dependen completamente de dónde se encuentre el acimut. [ 32 ]

Técnicas especiales

Polarimetría

Representación cromática de diferentes polarizaciones.

Las ondas de radar tienen polarización . Los distintos materiales reflejan las ondas de radar con diferentes intensidades, pero los materiales anisotrópicos , como el césped, suelen reflejar diferentes polarizaciones con diferentes intensidades. Algunos materiales también convierten una polarización en otra. Al emitir una mezcla de polarizaciones y utilizar antenas receptoras con una polarización específica, se pueden obtener varias imágenes a partir de la misma serie de pulsos. Con frecuencia, se utilizan tres de estas polarizaciones RX-TX (HH-pol, VV-pol, VH-pol) como los tres canales de color en una imagen sintetizada. Esto es lo que se ha hecho en la imagen de la derecha. La interpretación de los colores resultantes requiere pruebas exhaustivas con materiales conocidos.

Los avances recientes en polarimetría incluyen el uso de los cambios en la polarización aleatoria reflejada por ciertas superficies (como hierba o arena) y entre dos imágenes de la misma ubicación tomadas en momentos diferentes para determinar dónde ocurrieron cambios no visibles para los sistemas ópticos. Algunos ejemplos son la construcción de túneles subterráneos o las trayectorias de vehículos que circulan por la zona analizada. Se ha desarrollado una observación mejorada de manchas de petróleo en el mar mediante SAR, gracias a la modelización física adecuada y al uso de mediciones totalmente polarimétricas y de doble polarización.

Polarimetría SAR

Imagen SAR del Valle de la Muerte coloreada mediante polarimetría.

La polarimetría SAR es una técnica utilizada para obtener información física cualitativa y cuantitativa para aplicaciones terrestres, de nieve y hielo, oceánicas y urbanas, basada en la medición y exploración de las propiedades polarimétricas de dispersores artificiales y naturales. La clasificación del terreno y el uso del suelo es una de las aplicaciones más importantes del radar de apertura sintética polarimétrico (PolSAR). [ 41 ]

La polarimetría SAR utiliza una matriz de dispersión (S) para identificar el comportamiento de dispersión de los objetos tras interactuar con una onda electromagnética. La matriz se representa mediante una combinación de los estados de polarización horizontal y vertical de las señales transmitidas y recibidas.

S=[SHHSHVSVHSVV]{\displaystyle S={\begin{bmatrix}S_{HH}&S_{HV}\\S_{VH}&S_{VV}\end{bmatrix}}}

donde HH es para transmisión horizontal y recepción horizontal, VV es para transmisión vertical y recepción vertical, HV es para transmisión horizontal y recepción vertical, y VH para transmisión vertical y recepción horizontal.

Las dos primeras combinaciones de polarización se denominan polarización paralela (o copolarización), ya que las polarizaciones de transmisión y recepción son iguales. Las dos últimas combinaciones se denominan polarización cruzada, ya que las polarizaciones de transmisión y recepción son ortogonales entre sí. [ 42 ]

Modelo de potencia de dispersión de tres componentes

El modelo de potencia de dispersión de tres componentes de Freeman y Durden [ 43 ] se utiliza con éxito para la descomposición de una imagen PolSAR, aplicando la condición de simetría de reflexión mediante la matriz de covarianza. El método se basa en mecanismos de dispersión física simples (dispersión superficial, dispersión de doble rebote y dispersión volumétrica). La ventaja de este modelo de dispersión es que es simple y fácil de implementar para el procesamiento de imágenes. Hay 2 enfoques principales para un 3×{\displaystyle \times }3 descomposición de matriz polarimétrica. Una es el enfoque de matriz de covarianza lexicográfica basado en parámetros físicamente medibles, [ 43 ] y la otra es la descomposición de Pauli que es una matriz de descomposición coherente. Representa toda la información polarimétrica en una sola imagen SAR. La información polarimétrica de [S] podría representarse mediante la combinación de las intensidades|SHH|2,|SVV|2,2|SHV|2{\displaystyle |S_{HH}|^{2},|S_{VV}|^{2},2|S_{HV}|^{2}}en una sola imagen RGB donde todas las intensidades anteriores se codificarán como un canal de color. [ 44 ]

Modelo de potencia de dispersión de cuatro componentes

Para el análisis de imágenes PolSAR, puede haber casos en los que no se cumpla la condición de simetría de reflexión. En esos casos, se puede utilizar un modelo de dispersión de cuatro componentes [ 41 ] [ 45 ] para descomponer imágenes de radar de apertura sintética (SAR) polarimétricas. Este enfoque aborda el caso de dispersión no simétrica de reflexión. Incluye y extiende el método de descomposición de tres componentes introducido por Freeman y Durden [ 43 ] a un cuarto componente mediante la adición de la potencia de dispersión helicoidal. Este término de potencia helicoidal generalmente aparece en áreas urbanas complejas, pero desaparece para un dispersor distribuido natural. [ 41 ]

También existe un método mejorado que utiliza el algoritmo de descomposición de cuatro componentes, que se introdujo para los análisis generales de imágenes de datos polSAR. Los datos SAR se filtran primero, lo que se conoce como reducción de moteado, luego cada píxel se descompone mediante un modelo de cuatro componentes para determinar la potencia de dispersión de la superficie (PAGs{\displaystyle P_{s}}), poder de dispersión de doble rebote (PAGd{\displaystyle P_{d}}), poder de dispersión volumétrica (PAGv{\displaystyle P_{v}}), y poder de dispersión helicoidal (PAGdo{\displaystyle P_{c}}). [ 41 ] Los píxeles se dividen en 5 clases (superficie, doble rebote, volumen, hélice y píxeles mixtos) clasificadas según sus potencias máximas. Tras el cálculo, se añade una categoría mixta para los píxeles que presentan dos o tres potencias de dispersión dominantes iguales. El proceso continúa dividiendo los píxeles de todas estas categorías en 20 pequeños grupos de aproximadamente el mismo número de píxeles y fusionándolos según se desee; esto se denomina fusión de clústeres. Se clasifican iterativamente y, a continuación, se asigna automáticamente un color a cada clase. En resumen, este algoritmo permite comprender que los colores marrones representan las clases de dispersión superficial, los rojos las de dispersión de doble rebote, los verdes las de dispersión volumétrica y los azules las de dispersión helicoidal. [ 46 ]

Si bien este método está diseñado para el caso sin reflexión, incluye automáticamente la condición de simetría de reflexión, por lo que puede utilizarse como caso general. Además, conserva las características de dispersión al tener en cuenta la categoría de dispersión mixta, lo que lo convierte en un algoritmo superior.

Interferometría

En lugar de descartar los datos de fase, se puede extraer información de ellos. Si se dispone de dos observaciones del mismo terreno desde posiciones muy similares, se puede realizar una síntesis de apertura para obtener la resolución que proporcionaría un sistema de radar con dimensiones iguales a la separación entre las dos mediciones. Esta técnica se denomina radar de apertura sintética interferométrica (InSAR).

Si las dos muestras se obtienen simultáneamente (quizás colocando dos antenas en la misma aeronave, a cierta distancia una de la otra), cualquier diferencia de fase contendrá información sobre el ángulo desde el que regresó el eco del radar. Combinando esto con la información de distancia, se puede determinar la posición en tres dimensiones del píxel de la imagen. En otras palabras, se puede extraer la altitud del terreno, así como la reflectividad del radar, produciendo un modelo digital de elevación (MDE) con un solo paso del avión. Una aplicación aérea en el Centro Canadiense de Teledetección produjo mapas digitales de elevación con una resolución de 5  m y errores de altitud también de aproximadamente 5  m. La interferometría se utilizó para mapear muchas regiones de la superficie terrestre con una precisión sin precedentes utilizando datos de la Misión Topográfica de Radar del Transbordador .

Si las dos muestras se toman con un intervalo de tiempo diferente, por ejemplo, en dos vuelos sobre el mismo terreno, existen dos posibles causas de desfase. La primera es la altitud del terreno, como se mencionó anteriormente. La segunda es el movimiento del terreno: si este se ha desplazado entre las observaciones, la fase resultante será distinta. El desplazamiento necesario para producir una diferencia de fase significativa es del orden de la longitud de onda utilizada. Esto significa que si el terreno se desplaza unos centímetros, se puede apreciar en la imagen resultante (es necesario disponer de un mapa digital de elevación para diferenciar ambos tipos de desfase; puede ser necesario un tercer paso para obtener una diferencia de fase).

Este segundo método ofrece una herramienta poderosa en geología y geografía . El flujo de los glaciares se puede cartografiar con dos pasadas. Se han publicado mapas que muestran la deformación del terreno tras un terremoto menor o tras una erupción volcánica (mostrando la contracción de todo el volcán en varios centímetros). [ 47 ] [ 48 ] [ 49 ]

Interferometría diferencial

La interferometría diferencial (D-InSAR) requiere la toma de al menos dos imágenes, junto con un modelo digital de elevación (MDE). El MDE puede obtenerse mediante mediciones GPS o generarse por interferometría, siempre que el intervalo entre la adquisición de los pares de imágenes sea corto, lo que garantiza una distorsión mínima de la imagen de la superficie objetivo. En principio, tres imágenes del terreno con una geometría de adquisición similar suelen ser suficientes para D-InSAR. El principio para detectar el movimiento del terreno es bastante simple. Se crea un interferograma a partir de las dos primeras imágenes; este también se denomina interferograma de referencia o interferograma topográfico. Se crea un segundo interferograma que captura la topografía más la distorsión. Al restar este último del interferograma de referencia, se pueden observar franjas diferenciales que indican movimiento. La técnica de generación de D-InSAR de tres imágenes descrita se denomina método de tres pasadas o de doble diferencia.

Las franjas diferenciales que permanecen como franjas en el interferograma diferencial son el resultado de los cambios de alcance SAR de cualquier punto desplazado en el suelo de un interferograma al siguiente. En el interferograma diferencial, cada franja es directamente proporcional a la longitud de onda SAR, que es de aproximadamente 5,6  cm para el ciclo de fase única de ERS y RADARSAT. El desplazamiento de la superficie lejos de la dirección de observación del satélite provoca un aumento en la diferencia de trayectoria (que se traduce en diferencia de fase). Dado que la señal viaja desde la antena SAR hasta el objetivo y de vuelta, el desplazamiento medido es el doble de la unidad de longitud de onda. Esto significa que en la interferometría diferencial un ciclo de franja de −π a + π o una longitud de onda corresponde a un desplazamiento relativo a la antena SAR de solo media longitud de onda (2,8  cm). Existen varias publicaciones sobre la medición del movimiento de subsidencia, el análisis de estabilidad de taludes, deslizamientos de tierra, movimiento de glaciares, etc., utilizando herramientas D-InSAR. Un avance adicional de esta técnica permite utilizar la interferometría diferencial de paso ascendente y descendente del SAR satelital para estimar el movimiento del suelo en 3D. Las investigaciones en este campo han demostrado que se pueden lograr mediciones precisas del movimiento del terreno en 3D con una exactitud comparable a la de las mediciones basadas en GPS.

Tomo-SAR

La tomografía SAR es un subcampo del concepto de interferometría multilínea base. Se ha desarrollado para proporcionar una exposición tridimensional a la imagen, utilizando el concepto de formación de haces. Se puede utilizar cuando se requiere una atención de fase específica entre los componentes de magnitud y fase de los datos SAR durante la recuperación de información. Una de las principales ventajas de Tomo-SAR es que puede separar los parámetros que se dispersan, independientemente de cuán diferentes sean sus movimientos. [ 50 ] Al utilizar Tomo-SAR con interferometría diferencial, se desarrolla una nueva combinación denominada "tomografía diferencial" (Diff-Tomo). [ 50 ]

Tomo-SAR tiene una aplicación basada en imágenes de radar, que es la representación del volumen de hielo y la coherencia temporal del bosque ( la coherencia temporal describe la correlación entre ondas observadas en diferentes momentos en el tiempo). [ 50 ]

SAR de banda ultraancha

Los sistemas de radar convencionales emiten pulsos de energía radioeléctrica con un rango de frecuencias bastante estrecho. Un canal de banda estrecha, por definición, no permite cambios rápidos en la modulación. Dado que es el cambio en la señal recibida lo que revela el momento de su llegada (obviamente, una señal constante no revelaría nada sobre "cuándo" se reflejó en el objetivo), una señal con un cambio lento en la modulación no puede revelar la distancia al objetivo con la misma precisión que una señal con un cambio rápido en la modulación.

La banda ultraancha (UWB) se refiere a cualquier transmisión de radio que utiliza un ancho de banda muy amplio, lo que equivale a decir que utiliza cambios de modulación muy rápidos. Si bien no existe un valor de ancho de banda fijo que califique una señal como "UWB", los sistemas que utilizan anchos de banda superiores a una porción considerable de la frecuencia central (normalmente alrededor del diez por ciento) suelen denominarse sistemas "UWB". Un sistema UWB típico podría utilizar un ancho de banda de entre un tercio y la mitad de su frecuencia central. Por ejemplo, algunos sistemas utilizan un ancho de banda de aproximadamente 1  GHz centrado en torno a los 3  GHz.

Los dos métodos más comunes para aumentar el ancho de banda de la señal en radares UWB, incluyendo SAR, son los pulsos muy cortos y la modulación de frecuencia de alta amplitud (chirp). En este artículo se ofrece una descripción general de la modulación de frecuencia. El ancho de banda de un sistema con modulación de frecuencia puede ser tan estrecho o tan amplio como lo deseen los diseñadores. Aquí se describen los sistemas UWB basados ​​en pulsos, el método más común asociado al término "radar UWB".

Un sistema de radar basado en pulsos transmite pulsos muy cortos de energía electromagnética, generalmente de tan solo unas pocas ondas o incluso menos. Un pulso tan corto es, por supuesto, una señal que cambia muy rápidamente y, por lo tanto, ocupa un ancho de banda muy amplio. Esto permite una medición de distancia mucho más precisa y, por consiguiente, una mayor resolución.

La principal desventaja del SAR UWB basado en pulsos es que la electrónica de transmisión y recepción es difícil de diseñar para aplicaciones de alta potencia. Específicamente, el ciclo de trabajo de transmisión es excepcionalmente bajo y el tiempo de pulso excepcionalmente corto, por lo que la electrónica debe ser capaz de proporcionar una potencia instantánea extremadamente alta para rivalizar con la potencia promedio de los radares convencionales. (Si bien es cierto que UWB proporciona una ganancia notable en la capacidad del canal sobre una señal de banda estrecha debido a la relación del ancho de banda en el teorema de Shannon-Hartley y porque el bajo ciclo de trabajo de recepción recibe menos ruido, aumentando la relación señal-ruido , todavía existe una disparidad notable en el presupuesto de enlace porque el radar convencional puede ser varios órdenes de magnitud más potente que un radar típico basado en pulsos). Por lo tanto, el SAR UWB basado en pulsos se usa típicamente en aplicaciones que requieren niveles de potencia promedio en el rango de microvatios o milivatios, y por lo tanto se usa para escanear áreas objetivo más pequeñas y cercanas (varias decenas de metros), o en casos donde es posible una integración prolongada (durante un lapso de minutos) de la señal recibida. Sin embargo, esta limitación se resuelve en los sistemas de radar UWB modulados en frecuencia.

Las principales ventajas del radar UWB son una mejor resolución (unos pocos milímetros utilizando componentes electrónicos comerciales estándar ) y una mayor información espectral sobre la reflectividad del objetivo.

Agudeza del haz Doppler

El realce del haz Doppler se refiere comúnmente al método de procesar el historial de fase del haz real desenfocado para lograr una mejor resolución que la que se obtendría procesando el haz real sin él. Debido a que la apertura real de la antena de radar es muy pequeña (en comparación con la longitud de onda utilizada), la energía del radar se dispersa en un área amplia (generalmente de varios grados de ancho en una dirección ortogonal (en ángulo recto) a la dirección de la plataforma (aeronave)). El realce del haz Doppler aprovecha el movimiento de la plataforma, ya que los objetivos que se encuentran delante de la plataforma devuelven una señal Doppler desplazada hacia arriba (con una frecuencia ligeramente mayor) y los objetivos que se encuentran detrás de la plataforma devuelven una señal Doppler desplazada hacia abajo (con una frecuencia ligeramente menor).

La magnitud del desplazamiento varía con el ángulo hacia adelante o hacia atrás respecto a la dirección ortonormal. Conociendo la velocidad de la plataforma, la señal de retorno del objetivo se ubica en un intervalo angular específico que cambia con el tiempo. Las señales se integran a lo largo del tiempo, reduciendo sintéticamente el haz del radar a una apertura mucho menor; o, más precisamente (y basándose en la capacidad de distinguir desplazamientos Doppler más pequeños), el sistema puede tener cientos de haces muy estrechos simultáneamente. Esta técnica mejora drásticamente la resolución angular; sin embargo, es mucho más difícil aprovecharla para la resolución de alcance. (Véase radar Doppler pulsado ).

Radares de pulsos comprimidos (chirp)

Una técnica común en muchos sistemas de radar (generalmente también en sistemas SAR) consiste en modular la señal mediante una señal chirriante . En un radar con modulación chirriante, el pulso puede ser mucho más largo. Un pulso más largo permite emitir y, por lo tanto, recibir más energía, pero suele dificultar la resolución de alcance. Sin embargo, en un radar con modulación chirriante, este pulso más largo también presenta un desplazamiento de frecuencia durante su duración (de ahí la modulación chirriante o desplazamiento de frecuencia). Cuando se recibe la señal modulada, debe correlacionarse con el pulso emitido. Clásicamente, en los sistemas analógicos, se envía a una línea de retardo dispersiva (a menudo un dispositivo de ondas acústicas superficiales ) cuya velocidad de propagación varía según la frecuencia. Esta técnica comprime el pulso en el tiempo, lo que produce el efecto de un pulso mucho más corto (mejora la resolución de alcance) a la vez que ofrece la ventaja de una mayor duración del pulso (mayor energía de señal recibida). Los sistemas más modernos utilizan la correlación digital de pulsos para detectar el pulso recibido en la señal.

Operación típica

El instrumento AirSAR de la NASA está acoplado al lateral de un DC-8.

Recopilación de datos

En una aplicación típica de radar de apertura sintética (SAR), se instala una única antena de radar en una aeronave o nave espacial de manera que una parte sustancial del haz radiado por la antena tenga una dirección de propagación perpendicular a la trayectoria de vuelo. El haz debe ser amplio en la dirección vertical para que ilumine el terreno desde casi debajo de la aeronave hasta el horizonte.

Resolución de imagen y ancho de banda

La resolución en la dimensión de rango de la imagen se logra mediante la creación de pulsos que definen intervalos de tiempo muy cortos, ya sea emitiendo pulsos cortos que consisten en una frecuencia portadora y las bandas laterales necesarias, todo dentro de un ancho de banda determinado, o utilizando pulsos de modulación de frecuencia más largos en los que la frecuencia varía (a menudo linealmente) con el tiempo dentro de ese ancho de banda. Los diferentes tiempos en los que regresan los ecos permiten distinguir puntos a diferentes distancias.

La resolución de imagen del SAR en su coordenada de alcance (expresada en píxeles de imagen por unidad de distancia) es principalmente proporcional al ancho de banda de radio del tipo de pulso utilizado. En la coordenada transversal, la resolución es proporcional al ancho de banda del desplazamiento Doppler de las señales reflejadas dentro del ancho del haz. Dado que la frecuencia Doppler depende del ángulo de la dirección del punto de dispersión con respecto a la dirección perpendicular, el ancho de banda Doppler disponible dentro del ancho del haz es el mismo en todos los alcances. Por lo tanto, los límites teóricos de resolución espacial en ambas dimensiones de la imagen permanecen constantes con la variación del alcance. Sin embargo, en la práctica, tanto los errores que se acumulan con el tiempo de adquisición de datos como las técnicas específicas utilizadas en el posprocesamiento limitan aún más la resolución transversal a grandes distancias.

Resolución de imagen y ancho del haz

Antena SAR de los satélites SAOCOM .

La señal total es la de una porción del suelo del tamaño del ancho del haz. Para producir un haz estrecho en la dirección transversal , los efectos de difracción requieren que la antena sea ancha en esa dimensión. Por lo tanto, la distinción entre puntos co-alcance simplemente por la intensidad de las señales de retorno que persisten mientras se encuentran dentro del ancho del haz es difícil con antenas transportables por aeronaves, porque sus haces pueden tener anchos lineales solo unas dos órdenes de magnitud (cientos de veces) menores que el alcance. (Las transportables por naves espaciales pueden hacerlo 10 o más veces mejor). Sin embargo, si se registran tanto la amplitud como la fase de las señales de retorno, entonces la porción de esa señal de retorno de múltiples objetivos que fue dispersada radialmente por cualquier elemento más pequeño de la escena puede extraerse mediante la correlación del vector de fase de la señal de retorno total con la forma de la señal de retorno esperada de cada uno de dichos elementos.

El proceso puede concebirse como la combinación de una serie de observaciones distribuidas espacialmente, como si todas se hubieran realizado simultáneamente con una antena del tamaño del ancho del haz y enfocada en ese punto específico. La "apertura sintética" simulada al alcance máximo del sistema mediante este proceso no solo es más larga que la antena real, sino que, en aplicaciones prácticas, es mucho más larga que la aeronave de radar y muchísimo más larga que la nave espacial de radar.

Aunque algunas referencias a los radares de apertura sintética (SAR) los han caracterizado como "telescopios de radar", su analogía óptica real es el microscopio, ya que el detalle en sus imágenes es menor que la longitud de la apertura sintética. En términos de ingeniería de radar, mientras que el área objetivo se encuentra en el " campo lejano " de la antena iluminadora, se encuentra en el "campo cercano" de la simulada. Un diseño y operación cuidadosos pueden lograr una resolución de objetos menores a una millonésima parte del alcance, por ejemplo, 30 centímetros (12 pulgadas) a 300 kilómetros (190 millas) .  

Transmisión y recepción de pulsos

La conversión del tiempo de retardo de retorno a distancia geométrica puede ser muy precisa debido a la constancia natural de la velocidad y la dirección de propagación de las ondas electromagnéticas. Sin embargo, para una aeronave que vuela a través de la atmósfera, que nunca es uniforme ni está en calma, la relación entre los tiempos de transmisión y recepción de pulsos y las posiciones geométricas sucesivas de la antena debe ir acompañada de un ajuste constante de las fases de retorno para compensar las irregularidades detectadas en la trayectoria de vuelo. Los radares de apertura sintética (SAR) en naves espaciales evitan este problema atmosférico, pero aun así deben realizar correcciones por los movimientos conocidos de la antena debidos a las rotaciones de la nave, incluso aquellos que son reacciones a los movimientos de la maquinaria a bordo. Localizar un SAR en una nave espacial tripulada puede requerir que los humanos permanezcan inmóviles con respecto a la nave durante los períodos de recopilación de datos.

Las señales de retorno de los dispersores dentro del alcance de cualquier imagen se distribuyen en un intervalo de tiempo equivalente. El intervalo entre pulsos debe ser lo suficientemente largo para permitir que las señales de retorno de mayor alcance de cualquier pulso lleguen antes de que comiencen a aparecer las de menor alcance del siguiente pulso, de modo que no se superpongan en el tiempo. Por otro lado, la frecuencia entre pulsos debe ser lo suficientemente rápida para proporcionar suficientes muestras para la resolución transversal deseada (o transversal al haz). Cuando el radar se transporta en un vehículo de alta velocidad y debe obtener imágenes de una gran área con alta resolución, estas condiciones pueden entrar en conflicto, lo que da lugar a lo que se conoce como el problema de ambigüedad del SAR. Las mismas consideraciones se aplican también a los radares "convencionales", pero este problema se presenta de forma significativa solo cuando la resolución es tan fina que solo se puede obtener mediante procesos SAR. Dado que la base del problema es la capacidad de transmisión de información del único canal de entrada de señal proporcionado por una antena, la única solución es utilizar canales adicionales alimentados por antenas adicionales. El sistema se convierte entonces en un híbrido entre un radar de apertura sintética (SAR) y una antena de fase, a veces denominado antena Vernier.

Proceso de datos

La combinación de las series de observaciones requiere importantes recursos computacionales, generalmente mediante técnicas de transformada de Fourier . La alta velocidad de procesamiento digital disponible actualmente permite realizar dicho procesamiento prácticamente en tiempo real a bordo de un avión SAR. (Existe necesariamente un retardo mínimo hasta que se hayan recibido todas las partes de la señal). El resultado es un mapa de reflectividad de radar, que incluye tanto la amplitud como la fase.

Datos de amplitud

La información de amplitud, al visualizarse en un mapa, proporciona datos sobre la cobertura del suelo de forma muy similar a como lo hace una fotografía en blanco y negro. También se pueden realizar variaciones en el procesamiento, tanto en estaciones instaladas en vehículos como en estaciones terrestres, con diversos fines, para resaltar ciertas características de la imagen y realizar un análisis detallado del área objetivo.

Datos de fase

Aunque la información de fase de una imagen generalmente no está disponible para un observador humano en un dispositivo de visualización de imágenes, se puede conservar numéricamente y, en ocasiones, permite reconocer ciertas características adicionales de los objetos.

Moteado de coherencia

Desafortunadamente, las diferencias de fase entre los elementos de imagen adyacentes ("píxeles") también producen efectos de interferencia aleatorios denominados " moteado de coherencia ", que consisten en una especie de granularidad con dimensiones del orden de la resolución, lo que hace que el concepto de resolución adquiera un significado sutilmente diferente. Este efecto es el mismo que se observa tanto visual como fotográficamente en escenas ópticas iluminadas con láser. La escala de esta estructura de moteado aleatorio está determinada por el tamaño de la apertura sintética en longitudes de onda y no puede ser menor que la resolución del sistema. La estructura de moteado puede atenuarse a costa de la resolución.

Holografía óptica

Antes de que existieran las computadoras digitales rápidas, el procesamiento de datos se realizaba mediante una técnica de holografía óptica . Los datos analógicos del radar se registraban como un patrón de interferencia holográfica en una película fotográfica a una escala que permitía que la película conservara el ancho de banda de la señal (por ejemplo, 1:1.000.000 para un radar que utilizaba una longitud de onda de 0,6 metros). Luego, la luz, utilizando, por ejemplo, ondas de 0,6 micrómetros (como las de un láser de helio-neón ) que atravesaba el holograma, podía proyectar una imagen del terreno a una escala registrable en otra película a distancias focales razonables del procesador, de alrededor de un metro. Esto funcionaba porque tanto el SAR como los arreglos de fase son fundamentalmente similares a la holografía óptica, pero utilizan microondas en lugar de ondas de luz. Los "procesadores de datos ópticos" desarrollados para este propósito de radar [ 51 ] [ 52 ] [ 53 ] fueron los primeros sistemas informáticos ópticos analógicos eficaces y, de hecho, se idearon antes de que la técnica holográfica se adaptara completamente a la formación de imágenes ópticas. Debido a las diferentes fuentes de alcance y a las estructuras de señal transversales presentes en las señales de radar, los procesadores de datos ópticos para SAR incluían no solo lentes esféricas y cilíndricas, sino también, en ocasiones, lentes cónicas.

Apariencia de la imagen

Las siguientes consideraciones también se aplican a los radares de imágenes de terreno de apertura real, pero son más importantes cuando la resolución en alcance se ajusta a una resolución transversal que solo está disponible en un radar de apertura sintética (SAR).

Imagen SAR de 25 cm de resolución del centro de Cleveland, Ohio, tomada por Umbra.

Alcance, alcance transversal y ángulos

Las dos dimensiones de una imagen de radar son el alcance y la distancia transversal. Las imágenes de radar de zonas limitadas de terreno pueden parecerse a fotografías oblicuas, pero no a las tomadas desde la ubicación del radar. Esto se debe a que la coordenada de alcance en una imagen de radar es perpendicular a la coordenada de ángulo vertical de una fotografía oblicua. Por lo tanto, la posición aparente de la pupila de entrada (o centro de la cámara ) para ver dicha imagen no es como si estuviera en el radar, sino como si estuviera en un punto desde el cual la línea de visión del observador es perpendicular a la dirección de alcance oblicuo que conecta el radar con el objetivo, con el alcance oblicuo aumentando de arriba hacia abajo de la imagen.

Debido a que las pendientes en terrenos llanos varían en ángulo vertical, cada elevación de dicho terreno aparece como una superficie curva, específicamente una curva de coseno hiperbólico . Las verticales a diferentes distancias son perpendiculares a esas curvas. La dirección aparente de la mirada del observador es paralela al eje hipotético de la curva. Los objetos directamente debajo del radar aparecen como si se vieran ópticamente de forma horizontal (es decir, de lado) y los que se encuentran a distancias lejanas como si se vieran ópticamente desde arriba. Estas curvaturas no son evidentes a menos que se observen grandes extensiones de terreno cercano, incluyendo pendientes pronunciadas.

Visibilidad

Cuando se observan como se especificó anteriormente, las imágenes de radar de alta resolución de áreas pequeñas pueden parecerse mucho a las imágenes ópticas habituales, por dos razones. La primera se comprende fácilmente imaginando un mástil en la escena. El alcance oblicuo hasta su extremo superior es menor que hasta su base. Por lo tanto, el mástil solo puede aparecer correctamente con la parte superior hacia arriba cuando se observa en la orientación indicada. En segundo lugar, dado que la iluminación del radar incide hacia abajo, las sombras se ven en su dirección más familiar de "iluminación cenital".

La imagen de la parte superior del poste se superpondrá a la de algún punto del terreno situado en el mismo arco de alcance inclinado, pero a un alcance horizontal menor ("alcance terrestre"). Las imágenes de las superficies de la escena que miran tanto hacia la iluminación como hacia el punto de vista aparente tendrán geometrías similares a las de una escena óptica vista desde ese punto de vista. Sin embargo, las pendientes que miran hacia el radar aparecerán acortadas y las que miran en dirección opuesta se alargarán con respecto a sus dimensiones horizontales (mapa). Por lo tanto, las primeras aparecerán más brillantes y las segundas más oscuras.

Las señales provenientes de pendientes más pronunciadas que la perpendicular al rango de inclinación se superpondrán a las de terrenos de menor elevación a menor altura, siendo ambas visibles pero entremezcladas. Esto ocurre especialmente en superficies verticales como las paredes de los edificios. Otro inconveniente visual que surge cuando una superficie es más pronunciada que la perpendicular al rango de inclinación es que se ilumina por una cara, pero se observa desde la cara opuesta. Entonces, por ejemplo, se ve la pared de un edificio que da al radar como si fuera desde el interior, mientras que el interior del edificio y la pared trasera (la más cercana al observador, y por lo tanto, la que se espera que sea visible ópticamente) desaparecen, ya que carecen de iluminación al estar en la sombra de la pared frontal y el techo. Parte de la señal del techo puede superponerse a la de la pared frontal, y ambas pueden superponerse a la señal del terreno frente al edificio. La sombra visible del edificio incluirá las de todos los elementos iluminados. Las sombras largas pueden presentar bordes borrosos debido al movimiento de la antena de iluminación durante el tiempo de exposición necesario para crear la imagen.

Artefactos y sombras que reflejan

Las superficies que solemos considerar rugosas, si su relieve es menor que la longitud de onda del radar, se comportan como espejos lisos, mostrando, más allá de dicha superficie, imágenes adicionales de los objetos que se encuentran frente a ella. Estas imágenes reflejadas aparecerán dentro de la sombra de la superficie reflectante, llegando incluso a llenarla por completo, impidiendo así su reconocimiento.

La dirección de superposición de cualquier punto de la escena no apunta directamente hacia el radar, sino hacia el punto de la trayectoria actual del radar de apertura sintética (SAR) más cercano al objetivo. Si el SAR se desvía hacia adelante o hacia atrás con respecto a la dirección perpendicular, la dirección de iluminación, y por lo tanto la de la sombra, no será opuesta a la dirección de superposición, sino que estará inclinada hacia la derecha o hacia la izquierda. La imagen aparecerá con la geometría de proyección correcta cuando se visualice de manera que la dirección de superposición sea vertical, la trayectoria de vuelo del SAR esté por encima de la imagen y el alcance aumente ligeramente hacia abajo.

Objetos en movimiento

Los objetos en movimiento dentro de una escena SAR alteran las frecuencias Doppler de los retornos. Por lo tanto, dichos objetos aparecen en la imagen en ubicaciones desplazadas en la dirección transversal en cantidades proporcionales al componente de la dirección transversal de su velocidad. Los vehículos de carretera pueden aparecer fuera de la calzada y, por lo tanto, no ser reconocidos como elementos del tráfico rodado. Los trenes que aparecen lejos de sus vías se reconocen más fácilmente por su longitud paralela a la vía conocida, así como por la ausencia de una longitud igual de firma de la plataforma ferroviaria y de algún terreno adyacente, ambos sombreados por el tren. Si bien las imágenes de embarcaciones en movimiento pueden estar desplazadas de la línea de las partes iniciales de sus estelas, las partes más recientes de la estela, que aún participan de parte del movimiento de la embarcación, aparecen como curvas que conectan la imagen de la embarcación con la estela relativamente tranquila de la parte posterior. En tales casos identificables, la velocidad y la dirección de los elementos en movimiento pueden determinarse a partir de las cantidades de sus desplazamientos. El componente longitudinal del movimiento de un objetivo causa cierta falta de enfoque. Los movimientos aleatorios, como el del follaje de los árboles movido por el viento, el de los vehículos que circulan por terrenos irregulares o el de personas u otros animales que caminan o corren, generalmente hacen que esos objetos no se puedan enfocar, lo que provoca que se vean borrosos o incluso que se vuelvan prácticamente invisibles.

Estas consideraciones, junto con la estructura de moteado debida a la coherencia, requieren un periodo de adaptación para interpretar correctamente las imágenes SAR. Para facilitar este proceso, se han recopilado grandes conjuntos de firmas de objetivos significativas mediante numerosos vuelos de prueba sobre terrenos conocidos y objetos de interés cultural.

Industria comercial

En los últimos años, la comercialización de la tecnología de radar de apertura sintética (SAR) se ha expandido, y empresas privadas lanzan satélites SAR para proporcionar capacidades de imágenes de alta resolución para diversas aplicaciones. Estas incluyen monitoreo ambiental, respuesta a desastres, defensa e inteligencia, monitoreo de infraestructura y vigilancia marítima.

Capella Space , una empresa estadounidense de observación de la Tierra, opera una constelación de satélites SAR diseñados para proporcionar imágenes de alta resolución en cualquier condición meteorológica. Capella fue la primera empresa estadounidense en desplegar un satélite SAR comercial y continúa ampliando su flota con capacidades SAR avanzadas, que incluyen mayor resolución, tiempos de revisita más rápidos y automatización de tareas. Los datos SAR son utilizados frecuentemente por agencias gubernamentales, organizaciones de defensa y clientes comerciales para monitorear cambios en la Tierra prácticamente en tiempo real.

Otros proveedores comerciales de SAR incluyen a ICEYE , una empresa finlandesa especializada en pequeños satélites SAR, y a Airbus , que opera las misiones TerraSAR-X y PAZ. La industria comercial de SAR ha crecido significativamente gracias a los avances en la miniaturización de satélites, el procesamiento de datos en la nube y la inteligencia artificial, que mejoran la accesibilidad y la utilidad para una amplia gama de usuarios.

Historia

La historia del radar de apertura sintética comienza en 1951 con la invención de la tecnología por el matemático Carl A. Wiley y su desarrollo en la década siguiente. Inicialmente desarrollada para uso militar, la tecnología se ha aplicado desde entonces en el campo de la ciencia planetaria .

Relación con los conjuntos de antenas en fase

Una técnica estrechamente relacionada con el SAR utiliza una matriz (conocida como " matriz de fase ") de elementos de antena reales distribuidos espacialmente en una o dos dimensiones perpendiculares a la dimensión de alcance del radar. Estas matrices físicas son verdaderamente sintéticas, ya que se crean mediante la síntesis de un conjunto de antenas físicas subsidiarias. Su funcionamiento no requiere movimiento con respecto a los objetivos. Todos los elementos de estas matrices reciben simultáneamente en tiempo real, y las señales que las atraviesan pueden someterse individualmente a cambios controlados de fase. Un resultado posible es responder con mayor intensidad a la radiación recibida de una pequeña área específica de la escena, centrándose en dicha área para determinar su contribución a la señal total recibida. El conjunto de señales detectadas coherentemente y recibidas en toda la apertura de la matriz puede replicarse en varios canales de procesamiento de datos y procesarse de forma diferente en cada uno. El conjunto de respuestas así rastreadas a diferentes áreas pequeñas de la escena puede visualizarse conjuntamente como una imagen de la misma.

En comparación, el elemento de antena física (generalmente único) de un radar de apertura sintética (SAR) capta señales en diferentes posiciones y en diferentes momentos. Cuando el radar es transportado por una aeronave o un vehículo en órbita, esas posiciones dependen de una sola variable: la distancia a lo largo de la trayectoria del vehículo, que es una dimensión matemática (no necesariamente la misma que una dimensión geométrica lineal). Las señales se almacenan, convirtiéndose así en funciones, ya no del tiempo, sino de las ubicaciones de registro a lo largo de esa dimensión. Cuando las señales almacenadas se leen posteriormente y se combinan con desfases específicos, el resultado es el mismo que si los datos registrados hubieran sido captados por una matriz de fase de igual longitud y forma. Lo que se sintetiza, por lo tanto, es un conjunto de señales equivalente a lo que podría haber recibido simultáneamente una matriz de fase real de gran apertura (en una dimensión). El SAR simula (en lugar de sintetizar) esa larga matriz de fase unidimensional. Si bien el término del título de este artículo se ha derivado incorrectamente, ahora está firmemente establecido tras medio siglo de uso.

Si bien el funcionamiento de una antena de fase se comprende fácilmente como una técnica completamente geométrica, el hecho de que un sistema de apertura sintética recopile sus datos mientras se mueve (o su objetivo) a cierta velocidad implica que las fases que variaban con la distancia recorrida originalmente variaban con el tiempo, constituyendo así frecuencias temporales. Dado que las frecuencias temporales son las variables que suelen utilizar los ingenieros de radar, sus análisis de los sistemas SAR se formulan habitualmente (y de manera muy productiva) en estos términos. En particular, la variación de fase durante el vuelo a lo largo de la apertura sintética se considera una secuencia de desplazamientos Doppler de la frecuencia recibida con respecto a la frecuencia transmitida. Una vez que los datos recibidos se han registrado y, por lo tanto, se han vuelto atemporales, el procesamiento de datos SAR también se comprende como un tipo especial de antena de fase, que puede tratarse como un proceso completamente geométrico.

La base de las técnicas SAR y de antenas de fase reside en que las distancias que recorren las ondas de radar hacia y desde cada elemento de la escena consisten en un número entero de longitudes de onda más una fracción de una longitud de onda "final". Estas fracciones provocan diferencias entre las fases de la re-irradiación recibida en distintas posiciones del SAR o de la antena. Se requiere detección coherente para capturar la información de fase de la señal, además de la información de amplitud. Este tipo de detección exige encontrar las diferencias entre las fases de las señales recibidas y la fase simultánea de una muestra bien conservada de la iluminación transmitida.

Véase también

Referencias

  1. Kirscht, Martin y Carsten Rinke. "Reconstrucción 3D de edificios y vegetación a partir de imágenes de radar de apertura sintética (SAR)". MVA. 1998.
  2. "Introducción al RADAR aerotransportado", GW Stimson, Capítulo 1 (13 págs.).
  3. WILLIAMS, PDL (15 de noviembre de 1992). "Reseña de " Synthetic Aperture Radar—Systems & Signal Processing ". J. CURLANDER y RN McDONOUGH. (Chichester: John Wiley & Sons Ltd.) [ Pp. 647. ] Precio £46.95" . International Journal of Remote Sensing . 13 (17): 3364–3364 . doi : 10.1080/01431169208928371 . ISSN 0143-1161 . 
  4. "Viendo lo invisible: Cómo el radar de apertura sintética está revolucionando las operaciones espaciales y militares" . www.mrcy.com . Consultado el 28 de agosto de 2025 .
  5. 1 2 3 4 SAR tomográfico. Gianfranco Fornaro. Consejo Nacional de Investigación (CNR). Instituto de Detección Electromagnética del Medio Ambiente (IREA). Via Diocleziano, 328, I-80124 Nápoles, ITALIA.
  6. Oliver, C. y Quegan, S. Comprensión de las imágenes de radar de apertura sintética. Artech House, Boston, 1998.
  7. 1 2 3 4 5 6 7 8 9 10 11 12 Imágenes de radar de apertura sintética mediante técnicas de estimación espectral. Shivakumar Ramakrishnan, Vincent Demarcus, Jerome Le Ny, Neal Patwari, Joel Gussy. Universidad de Michigan.
  8. "Ciencia, ingeniería y sostenibilidad: Monitoreo de puentes con datos satelitales SAR" .
  9. Wu, Xuan; Zhang, Zhijie; Xiong, Shengqing; Zhang, Wanchang; Tang, Jiakui; Li, Zhenghao; An, Bangsheng; Li, Rui (12 de abril de 2023). "Un método de detección de inundaciones casi en tiempo real basado en aprendizaje profundo e imágenes SAR" . Teledetección . 15 (8): 2046. Código Bib : 2023RemS...15.2046W . doi : 10.3390/rs15082046 . hdl : 10150/674478 . ISSN 2072-4292 . 
  10. Garg, Shubhika; Feinstein, Ben; Timnat, Shahar; Batchu, Vishal; Dror, Gideon; Rosenthal, Adi Gerzi; Gulshan, Varun (7 de noviembre de 2023). "Destilación multimodal para el mapeo de la extensión de inundaciones" . Environmental Data Science . 2 e37. arXiv : 2302.08180 . Bibcode : 2023EnvDS...2E..37G . doi : 10.1017/eds.2023.34 . ISSN 2634-4602 . 
  11. A. Maity (2016). "Clasificación supervisada de datos polarimétricos de RADARSAT-2 para diferentes características del terreno". arXiv : 1608.00501 [ cs.CV ].
  12. 1 2 3 Moreira, Alberto; Prats-Iraola, Pau; Younis, Marwan; Krieger, Gerhard; Hajnsek, Irena; P. Papathanassiou, Konstantinos (2013). "Un tutorial sobre radar de apertura sintética" (PDF) . IEEE Geoscience and Remote Sensing Magazine . 1 (1): 6– 43. Bibcode : 2013IGRSM...1a...6M . doi : 10.1109/MGRS.2013.2248301 . S2CID 7487291 . 
  13. "Lockheed Martin: Radar de apertura sintética: "Reconocimiento las 24 horas"" .
  14. Colburn, Robert (10 de marzo de 2009). "Su legado de ingeniería: radar de apertura sintética" . IEEE USA InSight . IEEE.
  15. R. Bamler; P. Hartl (agosto de 1998). "Interferometría de radar de apertura sintética". Problemas inversos . 14 (4): R1– R54. Bibcode : 1998InvPr..14R...1B . doi : 10.1088/0266-5611/14/4/001 . S2CID 250827866 . 
  16. 1 2 G. Fornaro, G. Franceschetti, "SAR Interferometry", Capítulo IV en G. Franceschetti, R. Lanari, Synthetic Aperture Radar Processing, CRC-PRESS, Boca Raton, marzo de 1999.
  17. 1 2 Fornaro, Gianfranco; Pascazio, Vito (2014). "Interferometría y tomografía SAR: teoría y aplicaciones". Academic Press Library in Signal Processing: Volumen 2 – Comunicaciones y procesamiento de señales de radar . Vol. 2. pp. 1043–1117 . doi : 10.1016/B978-0-12-396500-4.00020-X . ISBN   9780123965004.
  18. Reigber, Andreas; Lombardini, Fabrizio; Viviani, Federico; Nannini, Matteo; Martínez Del Hoyo, Antonio (2015). "Imágenes tridimensionales y de orden superior con SAR tomográfico: técnicas, aplicaciones y problemas". Simposio Internacional de Geociencias y Teledetección IEEE 2015 (IGARSS) . págs. 2915–2918 . doi : 10.1109/IGARSS.2015.7326425 . hdl : 11568/843638 . ISBN  978-1-4799-7929-5. S2CID 9589219 . 
  19. Allen, Chris. "Implementación de radar de apertura sintética (SAR)" . Departamento de EECS, Universidad de Kansas . Consultado el 14 de enero de 2023 .
  20. Xiaoxiang Zhu, "Estimación espectral para tomografía de radar de apertura sintética", Ciencia y tecnología espacial orientada a la Tierra – ESPACE, 19 de septiembre de 2008.
  21. 1 2 DeGraaf, SR (mayo de 1998). "Imágenes SAR mediante métodos modernos de estimación espectral 2D". IEEE Transactions on Image Processing . 7 (5): 729– 761. Bibcode : 1998ITIP....7..729D . doi : 10.1109/83.668029 . PMID 18276288 . 
  22. D. Rodríguez. "Un sistema de modelado de generación de datos brutos SAR mediante álgebra de matriz de núcleo de Kronecker computacional". Señales, sistemas y computadoras, 2001. Actas de la trigésimo quinta conferencia de Asilomar del año: 2001. 1 .
  23. 1 2 T. Gough, Peter (junio de 1994). "Un algoritmo rápido de estimación espectral basado en la FFT". IEEE Transactions on Signal Processing . 42 (6): 1317– 1322. Bibcode : 1994ITSP...42.1317G . doi : 10.1109/78.286949 .
  24. 1 2 Datcu, Mihai; Popescu, Anca; Gavat, Inge (2008). "Caracterización de imágenes SAR complejas mediante análisis espectral de variantes espaciales". Conferencia de radar IEEE 2008 .
  25. J. Capo4 (agosto de 1969). "Análisis del espectro de número de onda de frecuencia de alta resolución". Actas del IEEE . 57 (8): 1408– 1418. doi : 10.1109/PROC.1969.7278 .{{cite journal}}: CS1 maint: nombres numéricos: lista de autores ( enlace )
  26. 1 2 3 4 A. Jakobsson; SL Marple; P. Stoica (2000). "Análisis del espectro de Capon bidimensional computacionalmente eficiente". IEEE Transactions on Signal Processing . 48 (9): 2651– 2661. Bibcode : 2000ITSP...48.2651J . CiteSeerX 10.1.1.41.7 . doi : 10.1109/78.863072 . 
  27. I. Yildirim; NS Tezel; I. Erer; B. Yazgan. "Una comparación de estimadores espectrales no paramétricos para imágenes SAR". Avances recientes en tecnologías espaciales, 2003. RAST '03. Conferencia internacional sobre. Actas del año: 2003 .
  28. "Realización iterativa del método Capon 2-D aplicado al procesamiento de imágenes SAR", Conferencia Internacional de Radar IET 2015.
  29. 1 2 R. Alty, Stephen; Jakobsson, Andreas; G. Larsson, Erik. "Implementación eficiente de los estimadores espectrales recursivos en el tiempo de Capon y APES". Conferencia de Procesamiento de Señales, 12.ª Europea, 2004 .
  30. Li, Jian; P. Stoica (1996). "Un enfoque de filtrado adaptativo para la estimación espectral y la obtención de imágenes SAR" . IEEE Transactions on Signal Processing . 44 (6): 1469– 1484. Bibcode : 1996ITSP...44.1469L . doi : 10.1109/78.506612 .
  31. Li, Jian; EG Larsson; P. Stoica (2002). "Estimación del espectro de amplitud para datos bidimensionales con huecos" . IEEE Transactions on Signal Processing . 50 (6): 1343– 1354. Bibcode : 2002ITSP...50.1343L . doi : 10.1109/tsp.2002.1003059 .
  32. 1 2 3 4 5 Moreira, Alberto. "Radar de apertura sintética: principios y aplicaciones" (PDF) .
  33. 1 2 3 4 5 6 Duersch, Michael. "Retroproyección para radar de apertura sintética". Archivo de becarios de BYU .
  34. 1 2 Zhuo, Li; Chungsheng, Li (2011). "Algoritmo de retroproyección para la formación de imágenes GEO-SAR de alta resolución". Simposio Internacional de Geociencias y Teledetección IEEE 2011. IEEE. págs. 336–339 . doi : 10.1109/IGARSS.2011.6048967 . ISBN  978-1-4577-1003-2. S2CID 37054346 . 
  35. Xiaoling, Zhang; Chen, Cheng (2011). "Un nuevo método de imágenes SAR 3D de superresolución basado en el algoritmo MUSIC". 2011 IEEE RadarCon (RADAR) . págs. 525–529 . doi : 10.1109/RADAR.2011.5960592 . ISBN  978-1-4244-8901-5.
  36. Yegulalp, AF (1999). «Algoritmo de retroproyección rápida para radar de apertura sintética». Actas de la Conferencia de Radar IEEE de 1999. Radar en el próximo milenio (Cat. No. 99CH36249) . págs. 60–65 . doi : 10.1109/NRC.1999.767270 . ISBN  0-7803-4977-6.
  37. Mark T. Crockett, "Introducción al radar de apertura sintética: una alternativa de alta resolución a la imagen óptica"
  38. Stimson, George (1998). Introducción al radar aerotransportado . Raleigh: SciTech Publishing, Inc. pág. 13. ISBN  9781891121012.
  39. Bauck, Jerald (19 de octubre de 2019). Una justificación para la retroproyección en la formación de imágenes de radar de apertura sintética Spotlight (Informe). doi : 10.31224/osf.io/5wv2d . S2CID 243085005 . 
  40. C. Romero, Simulación de alta resolución de imágenes de radar de apertura sintética. 2010. [En línea]. Disponible en: http://digitalcommons.calpoly.edu/cgi/viewcontent.cgi?article=1364&context=theses . Consultado el 14 de noviembre de 2016.
  41. 1 2 3 4 Y. Yamaguchi; T. Moriyama; M. Ishido; H. Yamada (2005). "Modelo de dispersión de cuatro componentes para la descomposición de imágenes SAR polarimétricas". IEEE Transactions on Geoscience and Remote Sensing . 43 (8): 1699. Bibcode : 2005ITGRS..43.1699Y . doi : 10.1109/TGRS.2005.852084 . S2CID 10094317 . 
  42. Woodhouse, HI 2009. Introducción a la teledetección por microondas. CRC Press, Taylor & Francis Group, Edición especial para la India.
  43. 1 2 3 A. Freeman; SL Durden (mayo de 1998). "Un modelo de dispersión de tres componentes para datos SAR polarimétricos". IEEE Transactions on Geoscience and Remote Sensing . 36 (3): 963– 973. Bibcode : 1998ITGRS..36..963F . doi : 10.1109/36.673687 .
  44. "PolSARpro v6.0 (Biomass Edition) Toolbox" (PDF) . ESA . Consultado el 20 de noviembre de 2022 .
  45. "Gianfranco Fornaro; Diego Reale; Francesco Serafino,"Imágenes SAR de cuatro dimensiones para la estimación de altura y el monitoreo de dispersores simples y dobles"". IEEE Transactions on Geoscience and Remote Sensing . 47 (1). 2009.
  46. "Haijian Zhang; Wen Yang; Jiayu Chen; Hong Sun," Clasificación mejorada de datos SAR polarimétricos basada en un modelo de dispersión de cuatro componentes"". Conferencia Internacional CIE sobre Radar 2006 .
  47. Bathke, H.; Shirzaei, M.; Walter, TR (2011). "Inflación y deflación en el estratovolcán de laderas empinadas Llaima (Chile) detectadas mediante InSAR" . Geophys. Res. Lett . 38 (10): L10304. Bibcode : 2011GeoRL..3810304B . doi : 10.1029/2011GL047168 .
  48. Dawson, J.; Cummins, P.; Tregoning, P.; Leonard, M. (2008). "Terremotos intraplaca superficiales en Australia Occidental observados mediante radar de apertura sintética interferométrica". J. Geophys. Res . 113 (B11): B11408. Bibcode : 2008JGRB..11311408D . doi : 10.1029/2008JB005807 .
  49. "Monitoreo de volcanes: Uso de InSAR para observar cambios en la forma de los volcanes - Instituciones de Investigación Incorporadas para la Sismología" .
  50. 1 2 3 Lombardini, Fabrizio; Viviani, Federico (2014). "Tomografía SAR multidimensional: avances para aplicaciones urbanas y perspectivas para aplicaciones forestales/de hielo". 11.ª Conferencia Europea de Radar de 2014. pp. 225–228 . doi : 10.1109/EuRAD.2014.6991248 . ISBN  978-2-8748-7037-8. S2CID 37114379 . 
  51. "Radar de apertura sintética", L. J. Cutrona, Capítulo 23 (25 págs.) del "Manual de radar" de McGraw Hill, 1970. (Escrito cuando el procesamiento óptico de datos era todavía el único método viable, por la persona que lideró ese desarrollo).
  52. "Una breve historia del Grupo de Óptica de los Laboratorios Willow Run", Emmett N. Leith, en Tendencias en Óptica: Investigación, Desarrollo y Aplicaciones (libro), Anna Consortini, Academic Press, San Diego: 1996.
  53. "Automatización visual e imágenes de alta resolución", WM Brown, JL Walker y WR Boario, IEEE Transactions on Aerospace and Electronic Systems, vol. 40, n.º 4, octubre de 2004, págs. 1426–1445.

Bibliografía

  • Curlander, John C.; McDonough, Robert N. (1992). Radar de apertura sintética: sistemas y procesamiento de señales . Teledetección y procesamiento de imágenes. Wiley. ISBN 978-0-471-85770-9.
  • Gart, Jason H (2006). Industria electrónica y aeroespacial en Arizona durante la Guerra Fría, 1945–1968: Motorola, Hughes Aircraft, Goodyear Aircraft (Tesis). Universidad Estatal de Arizona.
  • Moreira, A.; Prats-Iraola, P.; Younis, M.; Krieger, G.; Hajnsek, I.; Papathanassiou, KP (2013). "Un tutorial sobre radar de apertura sintética" (PDF) . IEEE Geoscience and Remote Sensing Magazine . 1 (1): 6– 43. Bibcode : 2013IGRSM...1a...6M . doi : 10.1109/MGRS.2013.2248301 . S2CID 7487291 . 
  • Woodhouse, Iain H (2006). Introducción a la teledetección por microondas . CRC Press. ISBN 9780415271233.
  • Mediciones InSAR del transbordador espacial
  • Imágenes del instrumento SAR del transbordador espacial en Wayback Machine (archivadas el 26 de febrero de 2013).
  • El Centro de Satélites de Alaska cuenta con numerosos documentos técnicos, incluyendo un texto introductorio sobre la teoría del radar de apertura sintética (SAR) y sus aplicaciones científicas.
  • SAR Journal realiza un seguimiento de la industria del radar de apertura sintética (SAR).
  • El radar de la NASA revela restos ocultos en la antigua Angkor en la Wayback Machine (archivado el 24 de diciembre de 2014)Laboratorio de Propulsión a Chorro