ES2455225T3 - Procedimiento y sistema de formación de imágenes radiométrica de síntesis de apertura espacio-temporal - Google Patents
Procedimiento y sistema de formación de imágenes radiométrica de síntesis de apertura espacio-temporal Download PDFInfo
- Publication number
- ES2455225T3 ES2455225T3 ES10723692.9T ES10723692T ES2455225T3 ES 2455225 T3 ES2455225 T3 ES 2455225T3 ES 10723692 T ES10723692 T ES 10723692T ES 2455225 T3 ES2455225 T3 ES 2455225T3
- Authority
- ES
- Spain
- Prior art keywords
- antenna
- elementary
- sampling
- satellite
- equal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
- G01S13/955—Radar or analogous systems specially adapted for specific applications for meteorological use mounted on satellite
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S3/00—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
- G01S3/02—Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
- G01S3/04—Details
- G01S3/06—Means for increasing effective directivity, e.g. by combining signals having differently oriented directivity characteristics or by sharpening the envelope waveform of the signal derived from a rotating or oscillating beam antenna
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
- G01S7/4026—Antenna boresight
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
- G01S7/4026—Antenna boresight
- G01S7/403—Antenna boresight in azimuth, i.e. in the horizontal plane
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
- G01S7/4026—Antenna boresight
- G01S7/4034—Antenna boresight in elevation, i.e. in the vertical plane
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
Procedimiento de formación de imágenes radiométrica destinado a cartografiar por lo menos un parámetro característico de objetos distantes dispuestos sobre una superficie y detectable por la radiación electromagnética multiangular emitida por los objetos, utilizando el procedimiento una antena de red (8; 432; 532) formada por unas antenas elementales (A(k)) dispuestas a bordo de por lo menos un satélite (4; 404, 406; 504, 506) móvil con respecto a la superficie, estando el procedimiento caracterizado porque comprende las etapas que consisten en: identificar (202) cada antena elemental (A(k)) por un índice k, aplicar (206) una indicación de tiempo según una base de tiempos (60), a posiciones sucesivas (P(t(i))) de una antena elemental (A(1)) predeterminada de la antena de red (8; 432, 532) que sirve como referencia espaciotemporal, muestreándose las posiciones sucesivas de la antena elemental de referencia (A(1)) según un período de muestreo (Dtmue) en instantes sucesivos (t(i)) ordenados según un índice (i) y referenciándose dentro de un marco de referencia espacial fijo con respecto a la superficie sobrevolada, estando un reloj local (H(k)) asociado a cada antena elemental (A(k)), para cada posición sucesiva muestreada (P(t(i))) de la antena elemental de referencia A(1), sincronizar (208) cada reloj local asociado (H(k)) con la base de tiempos (60), para cada posición sucesiva muestreada (P(t(i))) de la antena elemental de referencia (A(1)), determinar (210) las posiciones espaciales respectivas de las antenas elementales A(k) con respecto a la antena elemental de referencia (A(1)), en cada instante de muestreo (t(i)) correspondiente a una posición muestreada (P(t(i))) de la antena elemental (A(1)), medir (212) de manera analítica un campo electromagnético recibido por cada antena elemental (A(k)) y suministrar, para cada una de ellas, una componente en fase y una componente en cuadratura del campo, asignar (214) a cada medición analítica del campo electromagnético (E(k,t(i))) medido el índice (k) de la antena elemental correspondiente y la fecha (t(i)) de muestreo correspondiente, para un índice de imagen (im) predeterminado, definir (220) un dominio de imagen temporal (DIM(im)) constituido por los instantes de muestreo (t(j)) tales que la diferencia t(im)-t(j) es inferior o igual a una duración temporal de dominio de imagen DtIM predeterminada, para cada par (k, j) de índice (k) de antena elemental (A(k)) y de índice de instante de muestreo (j) correspondiente al dominio de imagen temporal DIM(im), suministrar (222) un vector de medición M(k, j) que contiene ordenadas temporalmente las muestras de las mediciones del campo electromagnético recibido por la antena elemental (A(k)) en los instantes de muestreo comprendidos entre el instante de muestreo t(j) y el instante de muestreo )**Fórmula** , designando Dtc una duración de integración de la correlación entre dos vectores de medición M(k,j), M(k',j'), construir (224) una imagen (IM(im)) a partir de una función de visibilidad espacio-temporal que comprende un conjunto de coeficientes de visibilidades, siendo igual por lo menos un coeficiente de visibilidad a la correlación de dos vectores de medición M(k,j), M(k',j') cuyos índices de muestreo j, j' son diferentes.
Description
Procedimiento y sistema de formación de imágenes radiométrica de síntesis de apertura espacio-temporal.
La presente invención se refiere a un procedimiento de formación de imágenes de síntesis de apertura que utiliza un radiómetro de formación de imágenes llevado a bordo de, por lo menos, un satélite y se refiere también al sistema de formación de imágenes correspondiente de realización del procedimiento.
Se han desarrollado satélites de observación de la Tierra, como, por ejemplo, el satélite SMOS (denominación en inglés de Soil Moisture and Ocean Salinity) para medir las temperaturas de brillo en forma de radiaciones electromagnéticas polarizadas, multiangulares, del suelo y de los océanos, con una alta posibilidad de revisita, típicamente un período de tres días, y de realizar así un seguimiento de la evolución de la humedad de los suelos y la salinidad de los océanos.
El satélite SMOS recorre una órbita heliosincrónica a la altitud de 755 km y a la velocidad de 7 km/s según un período de 100 minutos.
El radiómetro de formación de imágenes llevado a bordo del satélite SMOS funciona en la banda L dentro de una banda de 17 MHz en torno a la frecuencia correspondiente a una longitud de onda de 21 cm.
El radiómetro de formación de imágenes está compuesto por 54 antenas elementales dispuestas por grupos de dieciocho en tres palas de un rotor en forma de “Y”.
Las antenas elementales están asociadas, cada una de ellas, a un receptor diferente y están conectadas, a través de los receptores y de un arnés compuesto por fibras ópticas, a un correlador de salida codificada en un bit, que permite correlacionar los campos eléctricos recibidos por cada par de antenas elementales.
Las longitudes de las fibras ópticas que forman el arnés son prácticamente idénticas y la precisión de posicionamiento de las antenas en las palas es de 3 mm. Cuando el satélite está en servicio, el mismo está inclinado un ángulo de 31,2 grados con respecto a la dirección del nadir.
Al ser de aproximadamente 1.000 km la instantánea, denominada snap-shot, el pixel de la imagen construida por un procedimiento de síntesis de apertura espacial es igual a 40 km y el tiempo de integración para calcular la correlación para un par cualquiera de antenas elementales es igual a 1,2 segundos.
El tamaño del píxel es igual a la resolución geométrica del dispositivo de formación de imágenes y el tiempo de integración constitutivo de una instantánea.
Así, cien instantáneas consecutivas permiten la adquisición del mismo píxel para diferentes ángulos de incidencia comprendidos entre 0 y 55 grados.
Además, la sensibilidad radiométrica, denominada también resolución radiométrica, está comprendida entre 3,5 y 5,8º K.
El tiempo de muestreo de cada medición elemental del campo eléctrico recibido por una antena elemental es igual a 10-6 segundos.
Debido a la estructura de tres palas del radiómetro, la superficie de apertura obtenida es, por una parte, más pequeña que la correspondiente del disco en el cual están circunscritas las palas, y, por otra parte, anisótropa.
El aumento del número de palas permite incrementar la superficie de apertura y mejorar la isotropía de su contorno, aunque esto da como resultado un aumento del número de antenas elementales y consecuentemente un aumento de la masa del satélite.
El problema técnico consiste en incrementar la superficie de apertura sintetizada del radiómetro sin aumentar el número de antenas elementales, la resolución geométrica y sin deteriorar la resolución radiométrica.
A este efecto, la invención tiene como objetivo un procedimiento de formación de imágenes radiométrica destinado a cartografiar por lo menos un parámetro característico de objetos distantes dispuestos sobre una superficie y detectable por la radiación electromagnética multiangular emitida por los objetos, de manera que el procedimiento utiliza una antena de red formada por antenas elementales dispuestas a bordo de, por lo menos, un satélite móvil con respecto a la superficie, estando el procedimiento caracterizado porque comprende las etapas que consisten en:
identificar cada antena elemental A(k) por un índice k,
aplicar una indicación de tiempo según una base de tiempos, a posiciones sucesivas P(t(i)) de una antena
elemental A(1) predeterminada de la antena de red que sirve como referencia espacio-temporal, muestreándose las posiciones sucesivas de la antena elemental de referencia A(1) según un período de muestreo Δtmue en instantes sucesivos t(i) ordenados según un índice i y referenciándose dentro de un marco de referencia espacial fijo con respecto a la superficie sobrevolada,
estando un reloj local H(k) asociado a cada antena elemental A(k), para cada posición sucesiva muestreada P(t(i)) de la antena elemental de referencia A(1), sincronizar cada reloj local asociado H(k) con la base de tiempos,
para cada posición sucesiva muestreada P(t(i)) de la antena elemental de referencia A(1), determinar las posiciones espaciales respectivas de las antenas elementales A(k) con respecto a la antena elemental de referencia A(1),
en cada instante de muestreo t(i) correspondiente a una posición muestreada P(t(i)) de la antena elemental A(1), medir de manera analítica un campo electromagnético recibido por cada antena elemental A(k) y suministrar, para cada una de ellas, una componente en fase y una componente en cuadratura del campo,
asignar a cada medición analítica del campo electromagnético (E(k,t)) medido el índice k de la antena elemental correspondiente y la fecha t(i) de muestreo correspondiente,
para un índice de imagen im predeterminado, definir un dominio de imagen temporal DIM(im) constituido por los instantes de muestreo t(j) tales que la diferencia t(im)-t(j) es inferior o igual a una duración temporal de dominio de imagen ΔtIM predeterminada,
para cada par (k, j) de índice k de antena elemental A(k) y de índice de instante de muestreo j correspondiente al dominio de imagen temporal DIM(im), suministrar un vector de medición M(k, j) que contiene ordenadas temporalmente las muestras de las mediciones del campo electromagnético recibido por la antena elemental A(k) en los instantes de muestreo comprendidos entre el instante de muestreo t(j) y el instante de muestreo
Δt
c
t( j − ) , designando Δtc una duración de integración de la correlación entre dos vectores de medición
Δ
tmue
M(k,j), M(k’,j’),
construir una imagen IM(im) a partir de una función de visibilidad espacio-temporal que comprende un conjunto de coeficientes de visibilidades, siendo por lo menos un coeficiente de visibilidad igual a la correlación de dos vectores de medición M(k,j), M(k’,j’) cuyos índices de muestreo j, j’ son diferentes.
Según unos modos particulares de realización, el procedimiento de formación de imágenes radiométrica comprende una o varias de las siguientes características:
- -
- estando un tamaño de píxel Tpíxel predeterminado y presentando una velocidad de desplazamiento con respecto a la superficie observada Vsat por lo menos un satélite portador de la antena de referencia elemental A(1), la duración de integración de correlación Δtc es igual al producto de un coeficiente fraccionario α por el tamaño Tpíxel de un píxel dividido por la velocidad Vsat del por lo menos un satelite;
- -
- el coeficiente fraccionario α es inferior o igual a 1/3 y preferentemente igual a 1/5;
- -
- la duración temporal de dominio de imagen ΔtIM es igual a la relación de la más grande de las distancias D que separan un par cualquiera de antenas elementales A(k) con respecto a la velocidad de desplazamiento Vsat del por lo menos un satélite, y
- -
- una banda pasante B de filtrado de las frecuencias en torno a una frecuencia central v0 del campo electromagnético medido por cada una de las antenas elementales A(k) es inferior o igual a la relación de la velocidad de desplazamiento Vsat del por lo menos un satélite con respecto a la más grande de las distancias D que separan un par cualquiera de antenas elementales A(k);
- -
- una banda de frecuencias de medición del campo electromagnético recibido por las antenas elementales A(k) está comprendida dentro de la banda L y preferentemente dentro de la banda de las frecuencias comprendidas entre 1.400 y 1.427 MHz;
- -
- la precisión de localización de una antena elemental A(k) es igual al producto de un coeficiente fraccionario β por la longitud de onda λ0 correspondiente a la frecuencia central v0, siendo el coeficiente fraccionario β inferior a 1/5 y preferentemente inferior a 3/20;
- -
- la precisión de la base de tiempos es inferior a la relación de la precisión de localización con respecto a la
velocidad de la luz;
- -
- el período de muestreo Δtmue es inferior al mínimo de una primera duración mínima de muestreo ΔtShannon correspondiente a un criterio de Shannon y de una segunda duración ΔTB correspondiente al cumplimiento de la exigencia de una resolución radiométrica predeterminada ΔTB,
5 siendo la primera duración ΔtShannon igual a la mitad del período correspondiente a una longitud de onda λ0 de una frecuencia de medición,
siendo la segunda duración ΔtB función de la resolución radiométrica predeterminada ΔTB, de una temperatura de
10 antena TA, de una temperatura de receptor Trec, de un ancho de banda pasante de filtrado B, de un número total (N+1) de antenas elementales, de una longitud D de la antena de red, de una superficie eficaz Se de una antena elemental, del tamaño de un píxel Tpíxel, de la velocidad Vsat del por lo menos un satelite;
- -
- la segunda duración (ΔtB) cumple la ecuación: 15
3 ΔT1
1)
TS
píxel V
B TT
Δ
e
2 (Nt
+
=
B 2
+
D B
A sat rec
La invención tiene también por objetivo un sistema de formación de imágenes radiométrica por satélite destinado a cartografiar por lo menos un parámetro característico de objetos distantes dispuestos sobre una superficie y 20 detectable por la radiación electromagnética multiangular emitida por los objetos, caracterizado porque comprende
una antena de red formada por antenas elementales A(k) dispuestas a bordo de por lo menos un satélite móvil con respecto a la superficie,
25 por lo menos un conjunto de recepción conectado a la antena de red apto para
aplicar una indicación de tiempo, según una base de tiempos de alta precisión, a posiciones sucesivas P(t(i)) de una antena elemental A(1) predeterminada de la antena de red que sirve como referencia espacio-temporal, muestreándose las posiciones sucesivas de la antena elemental de referencia A(1) según un período de
30 muestreo ΔTmue en instantes sucesivos t(i) ordenados según un índice i y referenciándose dentro de un marco de referencia espacial fijo con respecto a la superficie sobrevolada,
estando un reloj local H(k) asociado a cada antena elemental A(k), para cada posición sucesiva muestreada P(t(i)) de la antena elemental de referencia A(1), sincronizar cada reloj local asociado H(k) con la base de 35 tiempos,
para cada posición sucesiva muestreada P(t(i)) de la antena elemental de referencia A(1), determinar las posiciones espaciales respectivas de las antenas elementales A(k) con respecto a la antena elemental de referencia A(1),
40 en cada instante de muestreo t(i) correspondiente a una posición muestreada P(t(i)) de la antena elemental A(1), medir de manera analítica un campo electromagnético recibido por cada antena elemental A(k) y suministrar, para cada una de ellas, una componente en fase y una componente en cuadratura,
45 asignar a cada medición analítica del campo electromagnético E(k,t(i)) medido el índice k de la antena elemental correspondiente y la fecha t(i) de muestreo correspondiente,
y una unidad de tratamiento apta para
50 para un índice de imagen im predeterminado, definir un dominio de imagen temporal DIM(im) constituido por los instantes de muestreo t(j) tales que la diferencia t(im)-t(j) es inferior o igual a una duración temporal de dominio de imagen ΔtIM predeterminada,
para cada par (k, j) de índice k de antena elemental A(k) y de índice de instante de muestreo j correspondiente al
55 dominio de imagen temporal DIM(im), suministrar un vector de medición M(k, j) que contiene ordenadas temporalmente las muestras de las mediciones del campo electromagnético recibido por la antena elemental A(k) en los instantes de muestreo comprendidos entre el instante de muestreo t(j) y el instante de muestreo
Δt
c
t( j − ) , designando Δtc una duración de integración de la correlación entre dos vectores de medición
Δ
tmue
M(k,j), M(k’,j’), y 60 construir una imagen IM(im) a partir de una función de visibilidad espacio-temporal que comprende un conjunto
de coeficientes de visibilidades, siendo por lo menos un coeficiente de visibilidad igual a la correlación de dos vectores de medición M(k,j), M(k’,j’) cuyos índices de muestreo j, j’ son diferentes.
Según unos modos particulares de realización, el sistema de formación de imágenes radiométrica comprende una o 5 varias de las siguientes características:
- -
- la antena de red comprende un primer segmento de antenas elementales alineadas, separado de un segundo segmento de antenas elementales o de una antena aislada única por una distancia de separación inferior o igual a un tamaño de la antena de red,
- -
- presentando cada antena elemental A(k) un centro respectivo C(k), los centros C(k), C(k+1) de dos antenas adyacentes de un mismo segmento están separados por una semilongitud de onda, y
- -
- el primer segmento y el segundo segmento, respectivamente el primer segmento y la antena aislada, están 15 orientados segun una misma direccion;
- -
- el sistema comprende un solo satélite y una sola antena de red de extensión unidimensional de longitud D,
- -
- la antena de red comprende un solo segmento de antenas elementales alineadas, presentando el segmento una longitud D1 que es la mitad de la longitud D de la antena de red, y una antena elemental aislada, separada por un espacio vacío de la misma longitud D1 que la del segmento y alineada con las antenas
elementales del segmento;
- -
- las antenas elementales de la antena de red están dispuestas en dos satélites que describen, cada uno de
25 ellos, una órbita polar heliosincrónica y que forman un vuelo en formación, estando el primer segmento dispuesto en el primer satélite y estando el segundo segmento dispuesto en el segundo satelite;
- -
- las altitudes de las órbitas polares son diferentes de una altura inferior al doble de la longitud de un segmento; y
- -
- el primer satélite y el segundo satélite están situados sobre dos órbitas polares heliosincrónicas de la misma altitud, estando los satélites separados en la dirección longitudinal por una distancia comprendida entre una y dos longitudes de un segmento.
35 La invención se comprenderá mejor al leer la descripción de una forma de realización que se ofrece a continuación, ofrecida únicamente a título de ejemplo y realizada en referencia a los dibujos, en los cuales:
- -
- la figura 1 es una vista general de una primera forma de realización de un sistema de formación de imágenes radiométrica según la invención,
- -
- la figura 2 es una vista detallada de la antena de red radiométrica utilizada en el sistema de formación de imágenes de la figura 1;
- -
- la figura 3 es una vista detallada por bloques funcionales de un conjunto de recepción de señales
45 radiométricas captadas por la antena de red de la figura 2 y llevado a bordo del satélite del sistema de la figura 1;
- -
- la figura 4 es una vista detallada por bloques de un módulo de medición del conjunto de recepción de la figura
3;
- -
- la figura 5 es un diagrama de flujo de un procedimiento de formación de imágenes radiométrica llevado a la práctica por el sistema de formación de imágenes descrito en las figuras 1 a 4;
- -
- la figura 6 es una representación espacio-temporal de la red interferométrica correspondiente al sistema 55 descrito en la figura 1;
- -
- la figura 7 es una vista general de una segunda forma de realización de un sistema de formación de imágenes radiométrica según la invención, basada en una formación en vuelo de dos satélites que se reparten el soporte de las antenas elementales que constituyen la antena de red,
- -
- la figura 8 es una vista detallada de una tercera forma de realización de un sistema de formación de imágenes radiométrica según la invención.
Según la figura 1, un sistema de formación de imágenes radiométrica 2 comprende un satélite de observación 4 en 65 órbita en torno a la Tierra, y una unidad de tratamiento 6 de los datos de observación recogidos y retransmitidos al suelo por el satélite 4.
El satélite de observación 4 describe una órbita polar heliosincrónica de 755 kilómetros de altitud, de periodo igual a
100 minutos y de velocidad lineal designada por Vsat cuyo módulo es igual a 7 km/s.
El satélite 4 está inclinado con respecto a la dirección del Nadir según un ángulo de cabeceo designado por θ en la
figura 1, e igual a 31,2 grados.
Un marco de referencia triaxial vinculado a la plataforma del satélite 4, designado por (XS, YS, ZS) comprende un eje
YS orientado según el sentido de la velocidad Vsat y un eje XS resultante de la rotación del ángulo θ en torno al eje
10 YS del eje de nadir.
El satélite 4 comprende una antena de red 8 de observación radiométrica formada por una antena de red elementales 9, un conjunto 10 de recepción y de transformación de señales radioeléctricas, transmitidas a través de un enlace radioeléctrico 11 y captadas por la antena de red 8, conectado a la salida de la antena de red 8.
15 El satélite 4 comprende también una antena de telemedición 12 conectada a la salida del conjunto 10 de recepción y apta para retransmitir hacia la unidad de tratamiento 6 datos de mediciones radiométricas a través de un enlace radioeléctrico de telemedición 14.
20 La antena de red 8 funciona dentro de una banda de frecuencias dedicada únicamente a la observación, designada por “banda L”, cuyas frecuencias están comprendidas entre 1.400 y 1.427 MHz.
Esta banda de observación está centrada en torno a una frecuencia v0 correspondiente a una longitud de onda λ0 igual a 21 cm y que presenta un ancho de banda igual a 27 MHz.
25 Esta banda de observación permite medir temperaturas de brillo de objetos en forma de radiaciones electromagnéticas polarizadas y multiangulares emitidas por los objetos, siendo las temperaturas de brillo representativas particularmente de la humedad de los suelos y de la salinidad de los océanos.
30 La antena de red 8 está compuesta en este caso por un conjunto de antenas elementales 9 de estructura idéntica, aunque dispuestas en emplazamientos diferentes, designadas por A(k) para diferenciarlas, siendo k un índice entero diferente asignado a cada antena A(k).
Las antenas elementales A(k) están dispuestas en alineación según un eje de extensión longitudinal de la antena de
35 red 8 designado por YA y de la misma dirección que el eje YS del satélite 4. Un marco de referencia triaxial vinculado a la antena de red 8, designado por (XA, YA, ZA), comprende ejes XA, YA, ZA cuyas orientaciones respectivas son idénticas a las de los ejes XS, YS, ZS del marco de referencia correspondiente vinculado al satélite 4.
De forma más precisa según la figura 2, la antena de red 8 forma una red interferométrica de antenas elementales y 40 de un tamaño designado por D.
La antena de red 8 comprende un segmento 20, que tiene una longitud total D1, igual en este caso a la mitad de D, es decir, 25 metros, y compuesto por un número entero N de antenas elementales 20 sucesivas, alineadas.
45 La antena de red comprende también una única antena aislada A(N+1) separada del segmento 20 por un espacio vacío de longitud D2 sustancialmente igual a D1, es decir, 25 metros, estando la antena aislada A(N+1) alineada según la alineación de las antenas elementales que componen el segmento 20.
El tamaño D de la antena de red 8 se selecciona en este caso de manera que la resolución geométrica, indicada 50 como Rg, de una imagen obtenida por las mediciones radiométricas, es decir, el tamaño Tpíxel de un píxel de imagen, sea igual a 4 km.
Al cumplir la resolución geométrica la ecuación:
λ
55 R = 0H (rel. 1)
g
D
designando H la altitud del satélite, se deduce que el valor D es igual a aproximadamente 50 metros, y la longitud D1 es igual a 25 metros.
Las antenas elementales 9 del segmento 20 son, por ejemplo, antenas de tipo “patch” clásicas con cavidad 60 resonante dieléctrica que tienen una superficie radiante o de recepción en forma de disco.
Cada antena elemental A(k) de índice k comprende un centro asociado C(k).
Dos antenas elementales A(k), A(k+1) del segmento 20 dispuestas adyacentes están separadas por sus centros respectivos C(k), C(k+1) una distancia dA inferior o igual a una semilongitud de onda λ0/2 para cumplir la condición de muestreo de Shannon.
En este caso, la distancia dA se selecciona igual a 10 cm y el número total N de antenas elementales del segmento 20 es por lo tanto igual a 250.
El tamaño D de la antena de red 8, el número y la disposición de las antenas elementales permiten de este modo obtener una resolución geométrica según la dirección ZA perpendicular a la velocidad del satélite 4, igual a 4 km, y cumplir los criterios de reconstrucción de una imagen por una transformada de Fourier según la dirección ZA después de la utilización de todas las mediciones de pares de antenas, definiendo cada par una base espacial de la red interferométrica formada por la antena de red elementales A(k).
La precisión de localización de cada antena elemental A(k) es en este caso igual a 3 mm.
En una variante y de forma general, la precisión de localización es igual al producto de un coeficiente fraccionario β de la longitud de onda λ0 por la longitud de onda λ0, siendo el coeficiente fraccionario β inferior a 1/3 y preferentemente inferior o igual a 3/20.
En una variante, la antena de red está provista de receptores de posicionamiento de tipo GPS (Global Positioning System) que permiten estimar la deformación de la antena de red 8 en función de fluctuaciones de temperatura a velocidad lenta con respecto a la velocidad del satélite 4.
Según la figura 3, el conjunto de recepción 10 comprende terminales de entradas I(k), variando el índice de antena elemental k de 1 a N+1, estando cada terminal de entrada I(k) conectado a una antena elemental A(k) asociada distinta, y comprende un terminal de salida (50) conectado a la antena de telemedición 12.
El conjunto de recepción 10 comprende una batería 52 de módulos de mediciones designados por Mod(k), variando k entre 1 y N+1, una memoria masiva 54 conectada por la entrada a cada módulo de medición Mod(k) por medio de un par asociado distinto L(k) de enlaces, y un único emisor de datos de telemedición 56 conectado por la entrada a la memoria masiva 54 y por la salida a la antena de telemedición 12.
El conjunto de recepción 10 comprende también una base de tiempos 60 ó reloj maestro de alta precisión que está conectado por entradas distintas S(k), variando k entre 1 y N+1, correspondiéndose cada entrada S(k) con un módulo de medición Mod(k) distinto asociado. Un enlace diferente de sincronización LS(k) conecta la base de tiempos 60 a cada entrada de sincronización S(k).
La precisión de la base de tiempos 60 es en este caso igual a 3/20 de la longitud de onda λ0 dividido por la velocidad de la luz c, es decir igual a 10-11 segundos. Esta precisión de la base de tiempos 60 se corresponde con la precisión de localización de cada antena elemental A(k) igual a 3 mm.
En una variante y de manera general, la precisión absoluta de la base de tiempos 60 es igual al producto de un coeficiente fraccionario β de la longitud de onda λ0 por la longitud de onda λ0 dividido por la velocidad de la luz, siendo el coeficiente fraccionario β inferior a 1/3 y preferentemente inferior o igual a 3/20.
Según la figura 4, un módulo de medición Mod(k) cualquiera comprende de manera genérica, conectados en serie, un filtro de encabezamiento 102, un amplificador de bajo ruido 104, un filtro de canal 106 ó de ruido de banda pasante designada por B y un mezclador reductor 108 de frecuencias.
El mezclador reductor 108 comprende una primera entrada 110 conectada al filtro de canal 106, una segunda entrada 112 apta para recibir una señal de frecuencia de desplazamiento o de traslación, un primer terminal de salida 114, respectivamente un segundo terminal de salida 116 apto para suministrar una señal de salida en fase de frecuencia intermedia (FI), respectivamente una señal de salida en cuadratura de la misma frecuencia intermedia (FI) correspondiente a la parte real, respectivamente la parte imaginaria de una señal de salida denominada analítica
o incluso sintetizada.
El módulo de medición Mod(k) comprende también, conectada al primer terminal de salida 114, respectivamente el segundo terminal de salida 116 del mezclador reductor 108, una primera cadena de tratamiento 118 de frecuencia intermedia (FI), respectivamente una segunda cadena de tratamiento 120 de frecuencia intermedia (FI).
El módulo de medición Mod(k) comprende también un reloj propio H(k) designado por 122, controlado por la base de tiempos 60 de alta precisión, conectado por la entrada S(k) a la base de tiempos 60, y también conectado a cada cadena de tratamiento de frecuencia intermedia 118, 120 así como al mezclador reductor 108.
La primera cadena de tratamiento 118 de frecuencia intermedia, respectivamente la segunda cadena de tratamiento
120 de frecuencia intermedia, comprenden, conectados en serie, un primer filtro de FI 124, un primer amplificador de FI 126, una primera unidad de corrección y de regulación 128, una primera unidad de muestreo 130, un primer digitalizador 132, respectivamente un segundo filtro de FI 134, un segundo amplificador de FI 136, una segunda unidad de corrección y de regulación 138, una segunda unidad de muestreo 140, un segundo digitalizador 142.
El primer digitalizador 132, respectivamente el segundo digitalizador 134, están conectados a la primera salida Oi(k) en fase, respectivamente la segunda salida Oq(k) en cuadratura, y son aptos para suministrar la parte real, respectivamente la parte imaginaria, cuantificada según un número entero q predeterminado de bits de la señal de medición analítica del campo electromagnético recibido por la antena correspondiente A(k). El número entero q es un entero que puede adoptar un valor cualquiera incluyendo el valor 1.
El reloj 122 propio del módulo de medición Mod(k) comprende una entrada 146 conectada a la entrada de sincronización S(k), una primera salida 148 conectada al segundo terminal de entrada 112 del mezclador reductor 108 y una segunda salida 150 conectada a una entrada de control 152 del primer muestreador 130 y a una entrada de control 154 del segundo muestreador 142.
El reloj 122 propio del módulo de medición Mod(k) comprende un generador de tiempo local 156, por ejemplo un oscilador de cuarzo, y una unidad de elaboración de señales de reloj 158 apta para controlar la señal de reloj local con la base de tiempos 60 de alta precisión y entregar señales de reloj síncrono de la base de tiempos 60.
Así, el módulo de medición Mod(k) es apto para aplicar una indicación de tiempo, con una alta precisión, a las señales medidas localmente por la antena elemental A(k) al mismo tiempo que la suministrada por la base de tiempos 60.
La banda pasante B del filtro de canal 106 se selecciona de tal manera que se verifica la siguiente relación:
Vsat
B ≤ (rel. 2)
D
designando Vsat el módulo de la velocidad del satélite 4 y designando D el tamaño de la antena de red elementales, dicho de otra manera, la distancia más grande que separa dos antenas elementales de la antena de red 8.
En este caso, se selecciona un valor de banda pasante B igual a 100 Hz.
Ya se han desarrollado filtros integrados de banda estrecha tales que funcionan en las longitudes de onda próximas a λ0.
La duración de muestreo Δtmue se selecciona de tal manera que sea inferior a la vez a una primera duración de muestreo ΔTShannon que cumple un criterio de Shannon y a una segunda duración ΔtB correspondiente al cumplimiento de la exigencia de una resolución radiométrica deseada, indicada como ΔTB, es decir cumple la inecuación:
Δtmue ≤ Min(Δtshannon ,ΔtB ) (rel. 3)
La primera duración de muestreo ΔTShannon que cumple un criterio de Shannon es igual a la mitad del periodo correspondiente a la longitud de onda λ0, es decir, en este caso aproximadamente 1,25.10-5.
La segunda duración ΔtB correspondiente al cumplimiento de exigencia de una resolución radiométrica deseada ΔTB cumple la ecuación siguiente:
ΔTB Se 1 Tpíxel 2 B
Δt =
2
(N + 1) (rel. 4)
T + T V
A recD
B sat
Designando TA la temperatura de antena, designando Trec la temperatura del receptor, B el ancho de banda pasante de filtrado, N+1 el número total de antenas elementales, D la longitud de la antena de red 8, Se la superficie eficaz de una antena elemental, Tpíxel el tamaño de un píxel, Vsat la velocidad del satélite 4.
B
Para un valor de la relación ΔT idéntico al del sistema SMOS, se determina el valor de 2*10-5 segundos, lo
T + T
A rec
cual se correspondería con un desplazamiento del satélite igual a 14 cm.
Aquí parece, en este caso concreto, que la primera duración de muestreo ΔTShannon que cumple un criterio de
Shannon es la duración de dimensionamiento del periodo de muestreo Δtmue.
En una variante, la frecuencia del filtro de canal 106 es programable, lo cual permite seleccionar una temperatura de brillo de observación.
En una variante, el conjunto de recepción es apto para medir varias temperaturas de brillo utilizando por lo menos dos baterías de filtros de canal centrados en frecuencias centrales diferentes.
Según la figura 5, un procedimiento de formación de imágenes radiométrica 200, destinado a cartografiar por lo menos un parámetro característico de objetos distantes dispuestos sobre una superficie y detectable por la radiación electromagnética multiangular emitida por los objetos, y que se lleva a la práctica por medio del sistema descrito en las figuras 1 a 4, comprende un conjunto de etapas 202, 204, 206, 208, 210, 212, 214, 216, 218, 220, 222, 224.
En una primera etapa 202, a cada antena elemental A(k) se le asocia un índice entero distinto k que permite identificar la antena elemental A(k).
En una etapa sucesiva de inicialización 204, se selecciona una antena elemental A(k) de la antena de red 8 como referencia espacial con respecto a las otras antenas. Por ejemplo, se selecciona la antena elemental A(1).
En la misma etapa de inicialización 204, se fija un instante inicial t(0) de puesta en marcha del muestreo de la posición inicial P(t(0)) de la antena elemental de referencia A(1) a través de la puesta a cero de un contador de índice de muestreo i.
En la misma etapa de inicialización 204, para cada antena elemental A(k), el reloj asociado (H(k)) se sincroniza con la base de tiempos 60 y se determinan con precisión las posiciones espaciales respectivas de las antenas elementales A(k) con respecto a la antena elemental de referencia A(1).
En una etapa sucesiva 206, para un valor actual i del índice del contador de muestreo, la posición muestreada P(t(i)) de la antena elemental A(1) seleccionada de la antena de red elementales se referencia en un marco de referencia espacial fijo con respecto a la superficie sobrevolada por el satélite y se le aplica una indicación temporal por medio de la base de tiempos 60 de alta precisión según una coordenada temporal t(i).
En la siguiente etapa 208, para cada antena elemental A(k), se mantiene la sincronización del reloj asociado (H(k)) con la base de tiempos 60 por un bucle de seguimiento de sincronización.
Esta etapa de mantenimiento de sincronización que viene a continuación de una fase de adquisición dentro de la etapa de inicialización 204 se efectúa, por ejemplo, con la ayuda de intercambio de señales de señalización y de un bucle de enganche de fase.
Cabe destacar que, en el caso en el que el reloj local sea muy estable, las correcciones de sincronización se efectúan según períodos de tiempo que pueden ser prolongados.
Cabe destacar también que, si las posiciones relativas entre las antenas elementales A(k) no varían en el transcurso del tiempo, pueden ponerse en práctica estrategias de sincronización de grupos de antenas.
En una etapa sucesiva 210, se actualiza, llegado el caso, la determinación de las posiciones espaciales respectivas de las antenas elementales A(k) con respecto a la antena elemental de referencia A(1).
En el caso de la antena descrita en la figura 2, las posiciones relativas de las antenas elementales A(k) se consideran en este caso invariables en el tiempo. La etapa de determinación propiamente dicha de las posiciones relativas de las antenas elementales no se pone en práctica más que una sola vez, reduciéndose las etapas de determinación siguientes a una recuperación de valores de posición relativos, guardados.
Sin embargo, si las posiciones relativas de las antenas elementales A(k) varían en el transcurso del tiempo con una velocidad reducida con respecto a la velocidad del satélite debido, por ejemplo, a fluctuaciones de temperatura experimentadas por la antena, puede ponerse en práctica un mecanismo de determinación de las posiciones relativas, basado en receptores de posicionamiento global y un modelo de deformación.
Cabe destacar que también se pone en práctica un mecanismo de determinación de las posiciones relativas de las antenas elementales cuando la antena de red está repartida en satélites en vuelo de formación, determinándose las posiciones y actitudes relativas de los satélites entre sí, que varían débilmente con respecto a la velocidad del satélite.
En una etapa sucesiva 212, para el instante de muestreo actual (t(i)) correspondiente a la posición muestreada P(t(i)) de la antena elemental A(1), se mide de manera analítica el campo electromagnético, por ejemplo eléctrico, recibido por cada antena elemental A(k), con el fin de suministrar, para cada antena elemental A(k) una componente en fase y una componente en cuadratura del campo medido E(k, t(i)).
En una etapa sucesiva 214, el índice k de la antena elemental y la fecha t(i) de muestreo correspondiente a la medición de campo electromagnético E(t, k) se asigna al valor complejo del campo medido.
5 A continuación, en una etapa sucesiva 216, después del transcurso de un período de muestreo ΔTmue a partir de t(i), se incrementa en una unidad el contador de los instantes de muestreo, y se pone en práctica la etapa 206 con el instante de muestreo incrementado. En este caso, el valor del periodo de muestreo Δtmue se toma igual a 1,25 10-5 segundos.
10 De forma paralela a la ejecución de las etapas de bucle 206, 208, 210, 212, 214, 216, se lleva a cabo una etapa de prueba 218 para verificar si el índice de muestreo actual i supera un valor de umbral iumbral predeterminado. El valor de umbral iumbral predeterminado se define por la relación
T
píxel
15 iumbral = (rel. 5)
V ΔT
sat mue
En este caso, con la velocidad del satélite igual a 7 km/s, el tamaño del píxel Tpíxel igual a 4 km, y el período de muestreo Δtmue igual a 1,25.10-5, el valor de umbral es igual a aproximadamente 4715.
20 Si el índice actual i no supera el valor de umbral predeterminado, entonces en la etapa de temporización 220, la ejecución ulterior de la etapa 218 se retarda en un período fijado.
Si el índice de muestreo i actual supera el valor de umbral predeterminado iumbral, en una etapa 220, para valores de im separados, por ejemplo, por Tpíxel/αVsat/Δtmue, se define un dominio de imagen temporal DIM(im) que está 25 constituido por los instantes de muestreo (t(j)) tales que la diferencia t(im)-t(j) es inferior o igual a una duración temporal de dominio de imagen ΔtIM predeterminada.
En este caso, la duración temporal de dominio de imagen ΔtIM es igual a la relación de la más grande de las distancias D que separan un par cualquiera de antenas elementales (A(k)) con respecto a la velocidad de 30 desplazamiento Vsat del satélite 4.
En la etapa siguiente 222, se suministra a cada par (k, j) de índice (k) de antena elemental A(k) y de índice de instante de muestreo j del dominio de imagen temporal DIM(im), un vector de medición M(k, j) que contiene, ordenadas temporalmente, las muestras de las mediciones del campo electromagnético recibido por la antena 35 elemental A(k) en los instantes de muestreo comprendidos entre el instante de muestreo t(j) y el instante de
Δt
c
muestreo t( j − ) , designando Δtc una duración de integración de la correlación entre dos vectores de medición
Δt
mue
M(k,j), M(k’,j’), y caracterizando el grado de aceptabilidad del efecto de desenfoque de movimiento debido al movimiento del satélite.
40 La duración de integración de la correlación entre dos mediciones es igual al producto de un coeficiente fraccionario α del tamaño de un píxel dividido por la velocidad del satélite Vsat, siendo el valor del coeficiente fraccionario inferior
o igual a 1/3 y preferentemente igual a 1/5.
Considerando en este caso el valor del coeficiente fraccionario α igual a 1/5, la duración de integración de la 45 correlación Δtc es igual a aproximadamente 0,12 segundos.
En la siguiente etapa 224, se construye una imagen a partir de una función de visibilidad espaciotemporal calculada a partir:
50 por una parte, de una matriz de bases espacio-temporales que se corresponden con un interferómetro espaciotemporal, definiéndose una base espacio-temporal por el dato de un doblete de pared (k, j), (k’, j’) con, en este caso, k el índice de una primera antena elemental A(k) y j un índice de muestreo de la primera antena elemental A(k), k’ el índice de una segunda antena elemental A(k’) y j’ un índice de muestreo de la segunda antena elemental A(k’), variando los índices k, k’ entre 1 y 251, y variando los índices de muestreo j entre i e i-ΔtIM/Δtmue,
55 por otra parte, de los productos de correlación con respecto a la duración Δtc entre los vectores de medición M(k,j), M(k’,j’) correspondientes a cada base ((k,j), (k’,j’)).
Un producto de correlación define un coeficiente de visibilidad indicado como Vim((k,j), (k’,j’)), y el conjunto de los 60 coeficientes de visibilidad para k, k’ que varían de 1 a N+1 y j, j’ que varían sobre DIM(im) define una función de visibilidad Vi.
En la etapa 224, se reconstruye así el dominio de Fourier de la imagen IM(im) a través de la función de visibilidad Vim, y la imagen IM(im) se obtiene finalmente mediante un tratamiento clásico de inversión.
Entre estos pares distintos, existen pares para los cuales los instantes de muestreo son diferentes. Así, la imagen 5 reconstruida a partir de la función de visibilidad es una verdadera imagen obtenida por síntesis de apertura espaciotemporal utilizando las relaciones de coherencia espacio-temporales entre las diferentes mediciones M(k, j).
La función de visibilidad se expresa por la siguiente expresión genérica de un coeficiente de visibilidad:
10 V((k, j),(k', j' )) =
l=0
de manera que el operador * representa el operador de conjugación y el índice de sumatorio l varía de 0 al número de periodos de muestreo Δtmue contenidos dentro de una duración de correlación Δtc.
15 El tratamiento de inversión se basa en una relación de constitución que relaciona los coeficientes de visibilidad Vi((k,j), (k’,j’)) y la fuente de temperatura de brillo de la superficie observada designada por S.
El procedimiento permite, con la misma resolución geométrica y la misma resolución radiométrica, incrementar la superficie de apertura sintetizada del radiómetro.
20 El procedimiento y el sistema permiten también mejorar la resolución geométrica sin aumentar el número de antenas elementales, sin disminuir la superficie de apertura y sin deteriorar la resolución radiométrica, utilizando una antena unidimensional más sencilla de realizar y de montar en un satélite y sin presentar ningún efecto de aumento drástico de masa en el satélite.
25 El procedimiento permite mejorar la resolución radiométrica cuando aumenta la frecuencia de muestreo 1/Δtmue.
En una variante, el procedimiento utiliza una red interferométrica de antenas elementales A(k) dispuestas a bordo de, por lo menos, dos satélites móviles con respecto a la superficie observada, describiendo los satélites un vuelo en 30 formación, es decir, sus posiciones relativas están coordinadas de forma precisa.
Según la figura 6, se describe una representación espacio-temporal de la red de interferometría espacio-temporal 300 obtenida con ayuda de la antena de red 8 espacio-temporal, en relación con un primer eje 302 orientado según el desplazamiento del satélite, representativo de la dimensión temporal, y un segundo eje 304, orientado según la
35 dirección de extensión unidimensional de la antena de red 8 física, constituida por el conjunto de las antenas elementales, y representativo de una dimensión espacial.
El primer eje 302 es representativo de la extensión de la imagen IM(im) según el sentido de desplazamiento del satélite, mientras que el segundo eje 304 es representativo de la extensión de la imagen IM(im) en una dirección 40 perpendicular al desplazamiento del satélite.
Se representan diferentes vistas 310, 312, 314, 316 de la antena 8 en diferentes instantes de muestreo.
La primera vista 310 de la antena 8 se corresponde con la posición de la antena 8 en el instante de muestreo t(im) 45 referenciado por el índice im.
La segunda vista 312 de la antena 8 se corresponde con la posición de la antena 8 en el instante de muestreo t(im-1) que precede inmediatamente al instante de muestreo de la primera vista 310.
50 La vista 314 de la antena 8 se corresponde con la posición de la antena 8 en el instante de muestreo t(im) -ΔtIM, es decir, el último instante de muestreo que se tiene en cuenta en la provisión de un vector de medición M(k,j), siendo k un índice cualquiera de una antena elemental A(k) y siendo j tal que t(j) pertenece al dominio de imagen temporal DIM(im).
55 La vista 318 se corresponde con la posición de la antena más alejada temporalmente de la vista 310 en el horizonte del tiempo de integración de correlación Δtc.
A título de ejemplo, se representa un vector de base 320 que se corresponde con el doblete de par ((k,j), (k’,j’)) con k igual a 1, j igual a i, k’ igual a 256 y j’ igual a i-ΔtIM/Δtmue.
60 Según la figura 7, una variante del sistema de formación de imágenes radiométrica 402 comprende un primer satélite 404 y un segundo satélite 406 en órbita en torno a la Tierra, aptos para sobrevolar la superficie observada de manera coordinada, y una unidad de tratamiento 408 de los datos de observación recogidos y retransmitidos al suelo
por los dos satélites.
El primer satélite de observación 404 describe una primera órbita polar 410 heliosincrónica de 755 kilómetros de
altitud, de período igual a 100 minutos y de velocidad lineal designada por Vsat cuyo módulo es igual a 7 km/s. 5 El segundo satélite de observación 406 describe una segunda órbita polar 412 desplazada en longitud aproximadamente 16,60 metros en el nivel del ecuador con respecto a la primera órbita, estando la altitud modulada según una amplitud de ±25 metros en torno a una altitud central idéntica a la de la primera órbita, para evitar una colisión en los polos y garantizar un mantenimiento de los satélites uno junto a otro en una misma latitud. 10 Se supone que en este caso los satélites 404 y 406 no están inclinados con respecto a su dirección respectiva de Nadir.
El primer satélite 404 comprende un primer segmento 424 de antenas elementales 426 y, no representados, un
15 primer conjunto de recepción, y una primera antena de telemedición apta para retransmitir hacia la unidad de tratamiento 408 los datos de mediciones radiométricas de los campos electromagnéticos recibidos por las antenas elementales 426 del primer segmento 424.
El segundo satélite 406 comprende un segundo segmento 428 de antenas elementales 430 y, no representados en
20 la figura 7, un segundo conjunto de recepción, y una segunda antena de telemedición apta para retransmitir hacia la unidad de tratamiento 408 los datos de mediciones radiométricas de los campos electromagnéticos recibidos por las antenas elementales 430 del segundo segmento 428.
El primer satélite 404 y el segundo satélite comprenden también medios de determinación, de alta precisión, de sus 25 posiciones y actitudes relativas, utilizando, por ejemplo, láseres.
Las antenas elementales 426 y 430 son idénticas a las descritas en la figura 2.
Las longitudes del primer segmento 424 y del segundo segmento 428 son idénticas e iguales a 16,6 metros.
30 La arquitectura del primer conjunto de recepción difiere de la arquitectura del conjunto de recepción descrito en la figura 3 en que una primera interfaz de enlace entre satélites está conectada a la base de tiempos de alta precisión del primer conjunto de recepción.
35 La arquitectura del segundo conjunto de recepción difiere de la arquitectura del conjunto de recepción descrito en la figura 3 en que la base de tiempos se sustituye por una segunda interfaz de enlace entre satélites compatible de la primera interfaz de enlace entre satélites.
Los dos segmentos 424 y 428 forman así una antena de red 432 análoga a la antena de red 8 descrita en la figura 2.
40 Sin embargo, la antena de red 432 difiere de la antena de red 8 de la figura 2 en las posiciones relativas de las antenas elementales del segundo segmento con respecto a las antenas elementales del primer segmento debido a un desplazamiento relativo entre el primer satélite.
45 Incluso si este desplazamiento relativo permanece confinado dentro de un volumen de 25 m3 aproximadamente, varía lentamente con respecto a la velocidad del satélite y permanece insignificante mientras dura el tiempo de integración de correlación, es necesario corregir regularmente las posiciones de las antenas elementales A(k) con respecto a una antena elemental de referencia A(1).
50 Esta corrección se pone en práctica a partir del conocimiento de las posiciones precisas de las antenas elementales con respecto a su segmento de pertenencia y con el conocimiento de la posición relativa de los segmentos 424 y
428.
La provisión de una indicación precisa de tiempo para cada antena elemental se lleva a cabo mediante la 55 sincronización del reloj local asociado a la antena elemental con la base de tiempos única, de alta precisión, que se lleva a bordo del primer satélite 404.
La antena de red 432 funciona de manera análoga a la antena de red 8 y se puede aplicar el procedimiento de formación de imágenes radiométrica descrito en la figura 5.
60 La ventaja proporcionada por la segmentación de la antena de red en dos es la realización de satélites más pequeños en los cuales se montan antenas más pequeñas, incluso si estas antenas más pequeñas no constituyen más que segmentos de la antena de red 432.
65 Según la figura 8, una variante 502 del sistema de formación de imágenes radiométrica comprende un primer satélite 504, respectivamente un segundo satélite 506, idéntico al primer satélite 404, respectivamente al segundo satélite 406 de la figura 7, siendo los dos satélites 504 y 506 en órbita en torno a la Tierra aptos para sobrevolar la superficie observada, de manera coordinada, y una unidad de tratamiento 508 idéntica a la unidad 408 de la figura 7.
5 La órbita 510 del primer satélite 504 es idéntica a la del primer satélite 404 de la figura 7.
El segundo satélite 506 describe una órbita polar 512 heliosincrónica, de la misma altitud que la altitud de la órbita 510 del primer satélite 504 y desplazada 25 metros permanentemente en latitud.
10 La órbita polar 512 está desplazada en longitud aproximadamente 16,60 metros en el nivel del ecuador con respecto a la primera órbita.
Así se evita una colisión de los satélites en los polos.
15 La antena de red 532 correspondiente funciona de forma análoga a la antena de red 432 y se puede aplicar el procedimiento de formación de imágenes radiométrica descrito en la figura 5.
La configuración de la figura 8 presenta la ventaja, con respecto a la de la figura 7, de no tener que ajustar la altura del segundo satélite con respecto al primer satélite. 20
Claims (13)
- REIVINDICACIONES1. Procedimiento de formación de imágenes radiométrica destinado a cartografiar por lo menos un parámetro característico de objetos distantes dispuestos sobre una superficie y detectable por la radiación electromagnética multiangular emitida por los objetos, utilizando el procedimiento una antena de red (8; 432; 532) formada por unas antenas elementales (A(k)) dispuestas a bordo de por lo menos un satélite (4; 404, 406; 504, 506) movil con respecto a la superficie, estando el procedimiento caracterizado porque comprende las etapas que consisten en:identificar (202) cada antena elemental (A(k)) por un índice k,aplicar (206) una indicación de tiempo según una base de tiempos (60), a posiciones sucesivas (P(t(i))) de una antena elemental (A(1)) predeterminada de la antena de red (8; 432, 532) que sirve como referencia espaciotemporal, muestreándose las posiciones sucesivas de la antena elemental de referencia (A(1)) según un período de muestreo (Δtmue) en instantes sucesivos (t(i)) ordenados según un índice (i) y referenciándose dentro de un marco de referencia espacial fijo con respecto a la superficie sobrevolada,estando un reloj local (H(k)) asociado a cada antena elemental (A(k)), para cada posición sucesiva muestreada (P(t(i))) de la antena elemental de referencia A(1), sincronizar (208) cada reloj local asociado (H(k)) con la base de tiempos (60),para cada posición sucesiva muestreada (P(t(i))) de la antena elemental de referencia (A(1)), determinar (210) las posiciones espaciales respectivas de las antenas elementales A(k) con respecto a la antena elemental de referencia (A(1)),en cada instante de muestreo (t(i)) correspondiente a una posición muestreada (P(t(i))) de la antena elemental (A(1)), medir (212) de manera analítica un campo electromagnético recibido por cada antena elemental (A(k)) y suministrar, para cada una de ellas, una componente en fase y una componente en cuadratura del campo,asignar (214) a cada medición analítica del campo electromagnético (E(k,t(i))) medido el índice (k) de la antena elemental correspondiente y la fecha (t(i)) de muestreo correspondiente,para un índice de imagen (im) predeterminado, definir (220) un dominio de imagen temporal (DIM(im)) constituido por los instantes de muestreo (t(j)) tales que la diferencia t(im)-t(j) es inferior o igual a una duración temporal de dominio de imagen ΔtIM predeterminada,para cada par (k, j) de índice (k) de antena elemental (A(k)) y de índice de instante de muestreo (j) correspondiente al dominio de imagen temporal DIM(im), suministrar (222) un vector de medición M(k, j) que contiene ordenadas temporalmente las muestras de las mediciones del campo electromagnético recibido por la antena elemental (A(k)) en los instantes de muestreo comprendidos entre el instante de muestreo t(j) y el instanteΔtcde muestreo t( j − ) , designando Δtc una duración de integración de la correlación entre dos vectores deΔtmuemedición M(k,j), M(k’,j’),construir (224) una imagen (IM(im)) a partir de una función de visibilidad espacio-temporal que comprende un conjunto de coeficientes de visibilidades, siendo igual por lo menos un coeficiente de visibilidad a la correlación de dos vectores de medición M(k,j), M(k’,j’) cuyos índices de muestreo j, j’ son diferentes.
-
- 2.
- Procedimiento de formación de imágenes radiométrica según la reivindicación 1, caracterizado porque un tamaño de píxel (Tpíxel) está predeterminado y por lo menos un satelite (4; 404; 504) que lleva la antena de referencia elemental (A(1)) presenta una velocidad de desplazamiento con respecto a la superficie observada (Vsat), la duración de integración de correlación (Δtc) es igual al producto de un coeficiente fraccionario (α) por el tamaño (Tpíxel) de un píxel dividido por la velocidad (Vsat) de dicho por lo menos un satélite (4; 404; 504).
-
- 3.
- Procedimiento de formación de imágenes radiométrica según la reivindicación 2, caracterizado porque el coeficiente fraccionario (α) es inferior o igual a 1/3 y preferentemente igual a 1/5.
-
- 4.
- Procedimiento de formación de imágenes radiométrica según cualquiera de las reivindicaciones 2 a 3, caracterizado porque
la duración temporal de dominio de imagen (ΔtIM) es igual a la relación de la más grande de las distancias (D) que separan un par cualquiera de antenas elementales (A(k)) con respecto a la velocidad de desplazamiento (Vsat) de dicho por lo menos un satélite (4; 404; 504), yporque una banda pasante (B) de filtrado de las frecuencias en torno a una frecuencia central (v0) del campo electromagnético medido por cada una de las antenas elementales (A(k)) es inferior o igual a la relación de lavelocidad de desplazamiento (Vsat) de dicho por lo menos un satélite (4; 404; 504) con respecto a la más grande de las distancias (D) que separan un par cualquiera de antenas elementales (A(k)). -
- 5.
- Procedimiento de formación de imágenes radiométrica según cualquiera de las reivindicaciones 1 a 4, caracterizado porque una banda de frecuencias de medición del campo electromagnético recibido por las antenas elementales (A(k)) está comprendida dentro de la banda L y preferentemente dentro de la banda de las frecuencias comprendidas entre 1.400 y 1.427 MHz.
-
- 6.
- Procedimiento según cualquiera de las reivindicaciones 4 a 5, caracterizado porque la precisión de localización de una antena elemental (A(k)) es igual al producto de un coeficiente fraccionario (β) por la longitud de onda (λ0) correspondiente a la frecuencia central (v0), siendo el coeficiente fraccionario (β) inferior a 1/5 y preferentemente inferior a 3/20.
-
- 7.
- Procedimiento según la reivindicación 6, caracterizado porque la precisión de la base de tiempos (60) es inferior a la relación de la precisión de localización con respecto a la velocidad de la luz.
-
- 8.
- Procedimiento según cualquiera de las reivindicaciones 1 a 7, caracterizado porque el período de muestreo (Δtmue) es inferior o igual al mínimo de una primera duración mínima de muestreo (ΔtShannon) correspondiente a un criterio de Shannon y de una segunda duración (ΔtB) correspondiente al cumplimiento de la exigencia de una resolución radiométrica predeterminada (ΔTB),
siendo la primera duración (ΔtShannon) igual a la mitad del período correspondiente a una longitud de onda (λ0) de una frecuencia de medición,siendo la segunda duración (ΔtB) función de la resolución radiométrica predeterminada (ΔTB), de una temperatura de antena (TA), de una temperatura de receptor T(rec), de un ancho de banda pasante de filtrado (B), de un número total (N+1) de antenas elementales, de una longitud (D) de la antena de red (8), de una superficie eficaz (Se) de una antena elemental, del tamaño de un píxel (Tpíxel), de la velocidad (Vsat) de dicho por lo menos un satélite (4; 404; 504). - 9. Procedimiento según la reivindicación 8, caracterizado porque la segunda duración (ΔtB) cumple la ecuación:ΔTB Se 1 Tpíxel 2 BΔt = 2 (N + 1)TA +TrecD2 B Vsat
- 10. Sistema de formación de imágenes radiométrica destinado a cartografiar por lo menos un parámetro característico de objetos distantes dispuestos sobre una superficie y detectable por la radiación electromagnética multiangular emitida por los objetos, caracterizado porque comprendeuna antena de red (8, 432, 532) formada por unas antenas elementales (A(k)) dispuestas a bordo de por lo menos un satélite (4; 404, 406; 504, 506) movil con respecto a la superficie,por lo menos un conjunto de recepción conectado a la antena de red (8; 432; 532) apto paraaplicar una indicación de tiempo, según una base de tiempos (60) de alta precisión, a posiciones sucesivas (P(t(i))) de una antena elemental (A(1)) predeterminada de la antena de red (8, 432, 532) que sirve como referencia espacio-temporal, muestreándose las posiciones sucesivas de la antena elemental de referencia (A(1)) según un período de muestreo ΔTmue en instantes sucesivos (t(i)) ordenados según un índice (i) y referenciándose dentro de un marco de referencia espacial fijo con respecto a la superficie sobrevolada,estando un reloj local (H(k)) asociado a cada antena elemental (A(k)), para cada posición sucesiva muestreada (P(t(i))) de la antena elemental de referencia A(1), sincronizar cada reloj local asociado (H(k)) con la base de tiempos (60),para cada posición sucesiva muestreada (P(t(i))) de la antena elemental de referencia (A(1)), determinar las posiciones espaciales respectivas de las antenas elementales A(k) con respecto a la antena elemental de referencia (A(1)),en cada instante de muestreo (t(i)) correspondiente a una posición muestreada (P(t(i))) de la antena elemental (A(1)), medir de manera analítica un campo electromagnético recibido por cada antena elemental (A(k)) y suministrar, para cada una de ellas, una componente en fase y una componente en cuadratura,asignar a cada medición analítica del campo electromagnético (E(k,t(i))) medido el índice (k) de la antena elemental correspondiente y la fecha (t(i)) de muestreo correspondiente,y una unidad de tratamiento (6) apta parapara un índice de imagen (im) predeterminado, definir (220) un dominio de imagen temporal (DIM(im)) constituido 5 por los instantes de muestreo (t(j)) tales que la diferencia t(im)-t(j) es inferior o igual a una duración temporal de dominio de imagen ΔtIM predeterminada,para cada par (k, j) de índice (k) de antena elemental (A(k)) y de índice de instante de muestreo (j) correspondiente al dominio de imagen temporal (DIM(im)), suministrar (222) un vector de medición M(k, j) que10 contiene ordenadas temporalmente las muestras de las mediciones del campo electromagnético recibido por la antena elemental (A(k)) en los instantes de muestreo comprendidos entre el instante de muestreo t(j) y el instante Δtcde muestreo t( j − ) , designando Δtc una duración de integración de la correlación entre dos vectores deΔtmuemedición M(k,j), M(k’,j’),15 construir (224) una imagen (IM(im)) a partir de una función de visibilidad espacio-temporal que comprende un conjunto de coeficientes de visibilidades, siendo por lo menos un coeficiente de visibilidad igual a la correlación de dos vectores de medición M(k,j), M(k’,j’) cuyos índices de muestreo j, j’ son diferentes.
- 11. Sistema de formación de imágenes radiométrica según la reivindicación 10, caracterizado porque20 la antena de red (8; 432; 532) comprende un primer segmento (20; 424), de antenas elementales alineadas, separado de un segundo segmento (428) de antenas elementales o de una antena aislada única (A(N+1)) por una distancia de separación inferior o igual a la mitad de un tamaño de la antena de red (8; 432; 532),25 porque cada antena elemental (A(k)) tiene un centro respectivo (C(k)), los centros (C(k), C(k+1)) de dos antenas adyacentes de un mismo segmento (20, 424, 428) están separados por una semilongitud de onda (λ0/2), y porqueel primer segmento (424) y el segundo segmento (428), respectivamente el primer segmento (20) y la antena 30 aislada (A(N+1)), están orientados según una misma dirección.
- 12. Sistema de formación de imágenes radiométrica según la reivindicación 11, caracterizado porque comprende un solo satélite (4) y una sola antena de red (8) de extensión unidimensional de longitud (D),35 porque la antena de red (8) comprende un solo segmento (20) de antenas elementales alineadas, presentando el segmento (20) una longitud (D1) que es la mitad de la longitud (D) de la antena de red (8), y una antena elemental aislada (A(N+1)), separada por un espacio vacío (30) de la misma longitud (D1) que la del segmento(20) y alineada con las antenas elementales del segmento (20).40 13. Sistema según la reivindicación 11, caracterizado porque las antenas elementales (426; 436) de la antena de red (432; 532) están dispuestas en dos satelites (404, 406; 504, 506) que describen, cada uno de ellos, una orbita polar heliosincrónica y que forman un vuelo en formación, estando el primer segmento (424) dispuesto en el primer satélite (404; 504) y estando el segundo segmento (428) dispuesto en el segundo satelite (406; 506).45 14. Sistema según la reivindicación 13, caracterizado porque las altitudes de las órbitas polares difieren en una altura inferior al doble de la longitud de un segmento.
- 15. Sistema según la reivindicación 13, caracterizado porque el primer satélite (504) y el segundo satélite (506) estánsituados sobre dos órbitas polares heliosincrónicas de la misma altitud, estando los satélites (504, 506) separados 50 en la dirección longitudinal por una distancia comprendida entre una y dos longitudes de un segmento.
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0953787 | 2009-06-08 | ||
FR0953787A FR2946434B1 (fr) | 2009-06-08 | 2009-06-08 | Procede et systeme d'imagerie radiometrique a synthese d'ouverture spatio-temporelle |
PCT/FR2010/050739 WO2010142876A1 (fr) | 2009-06-08 | 2010-04-16 | Procédé et système d'imagerie radiométrique à synthèse d'ouverture spatio-temporelle |
Publications (1)
Publication Number | Publication Date |
---|---|
ES2455225T3 true ES2455225T3 (es) | 2014-04-15 |
Family
ID=41682692
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
ES10723692.9T Active ES2455225T3 (es) | 2009-06-08 | 2010-04-16 | Procedimiento y sistema de formación de imágenes radiométrica de síntesis de apertura espacio-temporal |
Country Status (4)
Country | Link |
---|---|
EP (1) | EP2440946B1 (es) |
ES (1) | ES2455225T3 (es) |
FR (1) | FR2946434B1 (es) |
WO (1) | WO2010142876A1 (es) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102914774B (zh) * | 2012-09-18 | 2014-06-04 | 华中科技大学 | 一种综合孔径辐射计图像反演方法 |
FR3010798B1 (fr) * | 2013-09-19 | 2017-05-05 | Centre Nat D'etudes Spatiales (Cnes) | Procede d'imagerie radiometrique destine a mesurer au moins un parametre physique d'au moins une cible distante disposee sur une surface |
CN111811660B (zh) * | 2020-07-14 | 2021-06-04 | 中国科学院空天信息创新研究院 | 二维综合孔径微波辐射计亮温数据的优化处理方法 |
CN114994088B (zh) * | 2022-08-03 | 2022-10-25 | 北京国信航宇科技有限公司 | 星载大气探测仪及其大气参数测量方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4724439A (en) * | 1984-05-07 | 1988-02-09 | Hughes Aircraft Company | Microwave radiometer using fanbeam inversion |
FR2902526B1 (fr) * | 2006-06-16 | 2008-09-12 | Agence Spatiale Europeenne | Radiometre interferometrique |
-
2009
- 2009-06-08 FR FR0953787A patent/FR2946434B1/fr not_active Expired - Fee Related
-
2010
- 2010-04-16 WO PCT/FR2010/050739 patent/WO2010142876A1/fr active Application Filing
- 2010-04-16 ES ES10723692.9T patent/ES2455225T3/es active Active
- 2010-04-16 EP EP10723692.9A patent/EP2440946B1/fr active Active
Also Published As
Publication number | Publication date |
---|---|
WO2010142876A1 (fr) | 2010-12-16 |
EP2440946B1 (fr) | 2014-01-08 |
FR2946434A1 (fr) | 2010-12-10 |
FR2946434B1 (fr) | 2011-07-08 |
EP2440946A1 (fr) | 2012-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fu et al. | Remote sensing satellites for digital earth | |
JP7391429B2 (ja) | 衛星編隊ベースの遠隔検知システム及びコンステレーションシステム | |
Taylor et al. | The Canadian galactic plane survey | |
US20180364393A1 (en) | System and method for high-resolution radio occultation measurement through the atmosphere | |
ES2320599T3 (es) | Radiometro interferometrico. | |
Fallourd et al. | Monitoring temperate glacier displacement by multi-temporal TerraSAR-X images and continuous GPS measurements | |
EP3469395B1 (en) | Software-defined radio earth atmosphere imager | |
ES2455225T3 (es) | Procedimiento y sistema de formación de imágenes radiométrica de síntesis de apertura espacio-temporal | |
Li et al. | Observation of hurricane-generated ocean swell refraction at the Gulf Stream north wall with the RADARSAT-1 synthetic aperture radar | |
US9816863B1 (en) | Wind field vertical profile estimation using spectral radiance of multiband imagery and temperature profiles | |
CN113253233B (zh) | 基于全天空流星雷达信号的分析处理方法、系统 | |
ES2927925T3 (es) | Correlación de datos térmicos de imágenes satelitales para generar mapas térmicos de alta resolución espacial | |
Tadono et al. | Calibration and validation of PRISM onboard ALOS | |
Li et al. | Comprehensive calibration and validation site for information remote sensing | |
Kurum et al. | Integration of smartphones into small unmanned aircraft systems to sense water in soil by using reflected GPS signals | |
Cabot et al. | ULID: An unconnected L-band interferometer demonstrator | |
Chambon et al. | On-ground evaluation of MTG image navigation and registration (INR) performances | |
Fomalont | Astrometry and Geodesy | |
Cramer et al. | Data capture | |
Fallourd et al. | Alpine glacier 3D displacement derived from ascending and descending TerraSAR-X images on Mont-Blanc test site | |
Tampubolon et al. | Quality Assessment of an Extended Interferometric Radar Data Processing Approach | |
Li et al. | FPASMR: A new instrument for future sea surface salinity measurement | |
Aggarwal | Earth resource satellites | |
Sahawneh | Evaluation of the Hurricane Imaging Radiometer (HIRAD) Brightness Temperatures | |
Tian et al. | Comparison of Advanced Technology Microwave Sounder Biases Estimated Using Radio Occultation and Hurricane Florence (2018) Captured by NOAA-20 and S-NPP |