ES2392512T3 - Método para la reducción del desenfoque angular en imágenes de radar - Google Patents
Método para la reducción del desenfoque angular en imágenes de radar Download PDFInfo
- Publication number
- ES2392512T3 ES2392512T3 ES04445004T ES04445004T ES2392512T3 ES 2392512 T3 ES2392512 T3 ES 2392512T3 ES 04445004 T ES04445004 T ES 04445004T ES 04445004 T ES04445004 T ES 04445004T ES 2392512 T3 ES2392512 T3 ES 2392512T3
- Authority
- ES
- Spain
- Prior art keywords
- radar
- angular
- angular position
- shape
- antenna
- 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.)
- Expired - Lifetime
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000001514 detection method Methods 0.000 claims abstract description 4
- 101100234408 Danio rerio kif7 gene Proteins 0.000 claims abstract description 3
- 101100221620 Drosophila melanogaster cos gene Proteins 0.000 claims abstract description 3
- 101100398237 Xenopus tropicalis kif11 gene Proteins 0.000 claims abstract description 3
- 238000013459 approach Methods 0.000 claims abstract description 3
- 238000012937 correction Methods 0.000 description 12
- 238000002592 echocardiography Methods 0.000 description 6
- JTJMJGYZQZDUJJ-UHFFFAOYSA-N phencyclidine Chemical compound C1CCCCN1C1(C=2C=CC=CC=2)CCCCC1 JTJMJGYZQZDUJJ-UHFFFAOYSA-N 0.000 description 5
- 238000005259 measurement Methods 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 239000012141 concentrate Substances 0.000 description 2
- 101150118300 cos gene Proteins 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
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
- 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/4034—Antenna boresight in elevation, i.e. in the vertical 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
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
- G01S13/426—Scanning radar, e.g. 3D radar
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Método para la reducción del desenfoque angular debido a la anchura del haz del radar en imágenes de radar logradas en sistemas de detección basados en bins de distancias en los que la amplitud de la potencia S (x) es medida como una función de una posición angular x a una cierta distancia, comprendiendo: la determinación en cada bin de distancias, en el caso de una serie de posiciones angulares x, de una expresión del tipo x0 >= x + q/c / (rc (1 + q/c) + 1) ó x0 >= x + (l/d) ·arctan (q/ (2·d) ), en las que c es una constante que determina la forma de un polinomio de segundo orden utilizado para realizar la aproximación de la forma del haz, d es una constante que determina la forma de una aproximación a la forma del haz según S (x) >= k·cos2 (d· (x -x0) ), siendo k una amplitud, q es la relación q >= S’ (x) / S (x), siendo S’ (x) la primera derivada dS/dx de S (x), y x0 es una posición angular corregida; y se desplaza S (x) desde la posición angular x a la posición angular corregida x0 mediante la adición a S (x0) .
Description
Método para la reducción del desenfoque angular en imágenes de radar
La presente invención se refiere a un método para la reducción del desenfoque angular en imágenes de radar.
Un radar funciona mediante la radiación de ondas electromagnéticas procedentes de una antena, y detectando las ondas reflejadas (los ecos) por los objetos próximos al radar. Muy a menudo, el radar puede distinguir ecos de objetos a distancias diferentes, y de este modo obtiene una resolución de la distancia. A menudo, la medición de la distancia es discontinua, en el sentido de que los valores medidos de las distancias no varían de forma continua, sino que adoptan más bien un valor dentro de un conjunto predeterminado de valores, muy a menudo equidistantes. La explicación siguiente está limitada a dichos radares, con una medición con valores discontinuos de las distancias (“bins” de distancias).
Para vigilar una zona amplia, muchos radares de búsqueda tienen una antena de exploración, muy a menudo giratoria, para abarcar todas las direcciones, o que se desplaza hacia atrás y hacia adelante para cubrir un sector, por ejemplo, en sentido acimutal. En este caso, el radar recibe ecos en cada posición angular desde todos los bins de distancias estudiados (en algunos bins de distancias puede no existir eco). Si el radar es pulsante, la detección se realiza solamente en posiciones angulares discontinuas, una para cada impulso.
En este caso, la amplitud recibida (o potencia) de los ecos del radar puede ser mostrada como una imagen de radar con un conjunto finito de elementos de la imagen del radar (píxeles). En la mayor parte de casos, la imagen de radar de un sector cartesiano se presenta como un rectángulo en coordenadas polares con un número finito de posiciones angulares a lo largo del eje horizontal y un número finito de bins de distancias a lo largo del eje vertical (esta presentación se denomina presentación en “ámbito B”). Para mayor comodidad, muy a menudo se presenta el logaritmo de las amplitudes (potencias) (escala en decibelios, dB).
La explicación siguiente está limitada a un radar con datos presentados como dichas imágenes rectangulares de ámbito B con un número finito de píxeles. La coordenada angular se considera que en general es un acimut al realizar la antena de radar la exploración alrededor de un eje vertical, pero asimismo es posible en elevación.
La sensibilidad del radar es máxima en la dirección determinada por la antena. En la mayor parte de casos, esta dirección es frente a la antena. Dado que la anchura física de la antena es un número limitado de longitudes de onda que corresponden a la frecuencia utilizada, el haz principal de la antena, en el que la sensibilidad es máxima, no es infinitamente estrecho, sino que tiene una cierta anchura angular. El radar recibe un eco de un objeto, no solo cuando la antena está orientada directamente hacia el objeto, sino también cuando la antena está orientada ligeramente a la derecha del objeto, cuando el objeto está iluminado por la parte izquierda del haz principal de la antena y, cuando la antena está orientada ligeramente a la izquierda del objeto, cuando el objeto está iluminado por la parte derecha del haz principal de la antena. Un pequeño reflector del radar (en forma de punta), proporciona de este modo ecos en todo el alcance angular del haz principal del radar.
Dado que existen (en caso de duplicidad, en caso de análisis Doppler, o por otros motivos) diversas posiciones angulares utilizadas (diversos impulsos transmitidos) para cada anchura del haz principal de la antena del radar, la consecuencia es que en la coordenada angular se consiguen imágenes desenfocadas del radar. Este desenfoque hace que el análisis de la imagen del radar sea más difícil. Resulta que no es trivial la determinación del ángulo con respecto a un blanco del radar, incluso si el blanco tiene forma puntual, y resulta difícil la comparación entre la imagen del radar y una imagen de la misma zona por medio de otro sensor con características diferentes (posiblemente otro radar).
Existe una serie de métodos para reducir el desenfoque angular, y a continuación se comentarán algunos rasgos importantes de dichos métodos.
Según una primera técnica principal los métodos están basados en la desconvolución. Dado que el desenfoque puede ser considerado el resultado de la convolución de una imagen sin desenfoque con el haz de la antena, en teoría, puede ser eliminado realizando una transformación de Fourier de la imagen en dirección angular, dividiendo el esquema de la antena con la transformación de Fourier y realizando, a continuación, una transformación de Fourier inversa. Lo racional es que la transformación de Fourier convierta la convolución en multiplicación.
Según una segunda técnica principal, se comentan los métodos dentro de un radar general. Frecuentemente, datos tales como los de la figura 1, que muestran una imagen de radar medida de un sector cartesiano en escala de dB con posiciones angulares a lo largo del eje horizontal y “bins” de distancias a lo largo del eje vertical, pasan a través de un “extractor de blancos” que a menudo trabaja mediante la definición de un “blanco” como un máximo local con una amplitud que sobrepasa un umbral definido de una manera adecuada. A continuación, en un “bin” de distancias dado, se define la posición angular del blanco como la posición en la que se consigue el máximo. Puede calcularse una posición angular más exacta ajustando una curva adecuada tal como una parábola, con la posición angular como la variable independiente, a los datos próximos al máximo local, y calcular la posición angular del máximo de la parábola ajustada. En particular, no es difícil ajustar una parábola a los tres puntos más próximos al máximo local, con las mayores amplitudes, en particular si las posiciones angulares son equidistantes. Para ajustar una parábola a más de tres puntos, puede utilizarse un cálculo con el procedimiento de los mínimos cuadrados.
Este procedimiento, en general, funciona bien para el seguimiento de un número limitado de blancos, pero no siempre es practicable para eliminar el desenfoque de una imagen entera, para su presentación o para otros propósitos.
Según una tercera técnica principal, se utiliza un radar móvil. En el caso de radares aerotransportados, los métodos avanzados tales como el radar de apertura sintética (Synthetic Aperture Radar), (SAR, [LHMN cap. 15]) en su forma más simple “Adelgazamiento del haz mediante Doppler” (Doppler Beam Sharpening) [LHMN p. 260], utiliza un radar móvil con medición de la velocidad radial del radar (Doppler) y un tratamiento avanzado de la señal para obtener una “imagen” muy delgada de la parte controlada del terreno (no siempre descrita como un rectángulo en coordenadas polares). En esta descripción, se utiliza LHMN como una referencia abreviada a “Una introducción a sistemas de radar aéreos y aerotransportados Lacomme/Hardange/Marchais/Normant”, (Lacomme/Hardange/Marchais/Normant, Air and Spaceborne Radar Systems. An Introduction). ISBN 1 -891121 -13 -8, SciTech Publishing 2001.
Según una cuarta técnica principal, los métodos para la reducción del desenfoque angular presuponen un tipo particular de antena de radar, tal como una antena monoimpulso o una agrupación de antenas.
Los radares de seguimiento son radares que pueden medir en cuanto difiere la dirección de la antena del radar (la alineación) de una dirección deseada (la diferencia se denomina generalmente ángulo de error). Los radares de seguimiento se utilizan para apuntar cañones, guiar misiles, mantener una antena siguiendo un satélite de comunicaciones, etc. En general, un radar de seguimiento utiliza una antena monoimpulso [Sherman]. En esta descripción se utiliza Sherman como una referencia abreviada a “Sherman, Principios y técnicas de los monoimpulsos” (Monopulse Principles and Techniques), ISBN 0-89006-137-8, Artech House 1984. Una antena monoimpulso puede utilizar, por lo menos, dos disposiciones de antena, una disposición de suma con su direccionalidad máxima en la dirección de la antena y una disposición de diferencia que es sensible a blancos ligeramente desviados de la dirección de la antena. Existe una disposición de diferencia para cada variable angular utilizada (acimut, elevación). La figura 3 muestra disposiciones típicas de suma y de antena ([Sherman p. 138]). En la figura 3, las cantidades (fuerzas de campo) tienen, ambas, signos positivos y negativos. La figura 4 muestra el valor absoluto de la disposición de diferencia, comparada con la disposición de suma.
Algunos radares, a menudo radares aerotransportados multifuncionales, utilizan una antena de seguimiento monoimpulso en inspecciones y búsquedas. En este caso, puede utilizarse el canal de diferencia de diversas formas para reducir el desenfoque de la imagen debido a la anchura del haz principal de la antena (suma). A continuación, se indican dos tipos de métodos comprobados. Ambos tipos pueden estar incluidos bajo la denominación de “Haz monoimpulso adelgazado” (Monopulse Beam Sharpening).
1) Supongamos que S y D indican el canal suma y diferencia, respectivamente. Por tanto S -k·abs (D), con una constante pequeña positiva adecuada k, es un modelo que es = S en la dirección de la antena (alineación), en la que D = 0, pero que es menor que S fuera de la alineación, y por lo tanto más estrecha que S. Dicho modelo, puede ser utilizado para reducir el desenfoque angular en las imágenes de radar. Para evitar valores negativos de la amplitud (potencia), puede utilizarse una combinación max (0, S -k·abs(D)).
2) Un método mucho mejor y más sofisticado es el de calcular el ángulo de error del monoimpulso en cada punto en la imagen y determinar una posición angular corregida como “posición angular + ángulo de error”. La potencia detectada en un bin de distancias en una posición angular determinada no se almacena en la posición angular original, sino más bien en la posición angular corregida [LHMN páginas 205 -206]:
posición angular corregida = posición angular + ángulo de error
Esto proporciona un “adelgazamiento del haz” muy efectivo para la visualización o para otros propósitos.
En general, en la imagen en la que se debe eliminar el desenfoque, solamente existe un número limitado de posiciones angulares permitidas (denominadas “sub-haces” (“sub beams”) en [LHMN, página 206]) y, a continuación, el lado derecho de la fórmula deberá ser redondeado a la posición permitida más próxima. El “almacenamiento” puede ser realizado mediante la adición a valores ya almacenados (especialmente cuando se añade potencia) o mediante la utilización de una operación del tipo max(nuevo_valor, valores_ya_almacenados). Una alternativa, que en algunos casos proporciona resultados más exactos, es distribuir la magnitud (potencia) en las dos posiciones angulares permitidas, más próximas, posiblemente mediante interpolación lineal.
Si el radar utiliza una disposición de antena con varios elementos independientes de transmisión o de recepción, se dispone asimismo de métodos especiales.
A partir del documento U.S.A.-A-5 831 570 era conocido previamente un sistema de detección que comprende técnicas monoimpulso en acimut basadas en señales procedentes de dos canales, un canal de suma y un canal de diferencia.
El objetivo de la invención es obtener un método general para reducir de forma eficiente el desenfoque angular, y que puede ser utilizado, asimismo, para grandes cantidades de datos. El objetivo se alcanza mediante un método como el reivindicado en la reivindicación 1. Por medio de características tales como las determinadas en las reivindicaciones independientes, se obtienen ventajas adicionales. El método de la invención no utiliza la desconvolución ni presupone mediciones Doppler, ni ningún tipo particular de antena.
La invención se describirá a continuación con mayor detalle haciendo referencia a los dibujos adjuntos, en los cuales:
La figura 1 muestra una imagen de radar medida de un sector cartesiano en escala de dB.
La figura 2 muestra un corte horizontal a través de la figura 1 para el bin de distancias #69 en una escala de dB.
La figura 3 muestra una disposición monoimpulso típica de suma (simétrica) y una disposición de diferencia.
La figura 4 muestra una disposición monoimpulso de suma y el valor absoluto de la disposición de diferencia.
La figura 5 muestra una imagen de radar, según la figura 1, después de la reducción del desenfoque, según el método de la invención.
La figura 6 muestra un corte horizontal, según la figura 2, después de la reducción del desenfoque, según el método de la invención.
A continuación, se presenta un procedimiento para reducir el desenfoque angular en imágenes de radar, del tipo definido anteriormente. En la explicación, se supone que el método es puesto en práctica mediante cálculos en un procesador digital.
En la figura 1, se presenta una típica imagen medida de radar de un sector cartesiano. La imagen muestra la magnitud de la señal de radar recibida (video) en 230 bins de distancias (valores de distancias) para cada una de las 140 posiciones angulares. Si la magnitud está indicada mediante S, los datos pueden ser descritos de este modo como una disposición S(i,j), 1 : i : 140, 1 : j : 230. Los puntos oscuros significan una magnitud grande. Tanto los bins de distancias como las posiciones angulares se considerarán en este caso equidistantes.
Es evidente que la figura 1 contiene muchos objetos distintos, pero el desenfoque angular debido a la anchura del haz del radar (acimut) es muy importante. En él, parece que no existan blancos en forma de punto; en vez de ello, todos los blancos están prolongados en la dirección angular. El radar utilizado tiene una anchura del haz acimutal 3 dB de unas 10 posiciones angulares, lo que significa que un blanco puntual con un contraste con respecto al fondo de 20 dB o mayor puede parecer que tenga una anchura de unas 20 posiciones angulares.
La figura 2 muestra un corte horizontal (distancia constante) en la figura 1, a través del bin de distancias #69, que contiene un reflector del radar muy potente (blanco). Los datos de la figura 2 se pueden describir de este modo como S(i,69), 1: i : 140. El haz principal (la forma del lobulo pri ncipal de la disposición de direccionalidad de la antena del radar) es claramente visible. Cerca de un máximo local, las amplitudes medidas siguen la disposición del haz.
En la figura 2, el tamaño del haz principal del acimut de la antena del radar proporciona retornos apreciables del blanco en posiciones angulares, por lo menos entre #20 y #40, aunque el blanco parece estar situado aproximadamente en la posición angular #28 ó #29.
La forma de la figura 2 es realmente la imagen simétrica de izquierda a derecha de la forma del haz de la antena (a la izquierda de la posición correcta, la parte derecha del haz de la antena ilumina el blanco, etc.). Sin embargo, la mayor parte de los haces reales de la antena son aproximadamente simétricos y, en consecuencia, la diferencia entre izquierda y derecha es apenas perceptible.
A continuación, se mostrará la reducción del desenfoque angular en los datos de las figuras 1, 2.
Antes de una presentación más detallada del método sugerido para la reducción del desenfoque angular en imágenes de radar siguen, a continuación, como antecedentes, algunas explicaciones sobre ecuaciones diferenciales. Supongamos que una variable medida S(x) como la de la figura 2 tiene un máximo local en el ángulo x = x0. Cerca del máximo, se puede aproximar la curva mediante la disposición del haz de la antena, asimismo indicado con S(x), con un ligero abuso de la notación. En general, S es una función uniforme semejante a una parábola, simétrica alrededor de x = x0 y con una tangente horizontal en su máximo en x = x0.
El objetivo deseado es utilizar los datos medidos para determinar el máximo local x0 más próximo y concentrar la energía desenfocada del haz en este punto. Se puede determinar x0 mediante algún tipo de ajuste de curva local, ajustando localmente una expresión para la función S de la forma del haz con relación a los datos medidos, por ejemplo, mediante el método de los mínimos cuadrados (como en el “extractor de blancos” mencionado anteriormente).
Sin embargo, se puede determinar asimismo x0 mediante la utilización de una ecuación diferencial, utilizando una o varias derivadas (parciales) (en la dirección angular) de S(x), estimadas a partir de los datos. Dado que la primera derivada proporciona la pendiente de S(x) y la segunda derivada proporciona la curvatura, esto es geométricamente razonable. Dicho procedimiento es ventajoso en muchos casos, y se demostrará a continuación mediante algunos ejemplos.
Ejemplo 1. Una aproximación muy sencilla de la forma del haz próximo a x = x0 es un polinomio de segundo orden: S(x) = k·(1 -c·(x -x0)2) lo que da una curva parabólica con un máximo = k en x = x0. La constante k da la amplitud en el máximo, y la
constante c determina la forma: una c positiva grande proporciona un haz delgado, en punta. La aproximación solamente debe ser utilizada para valores pequeños de (x -x0) que dan valores positivos de S(x). La derivada (parcial) de S, como una función de la variable angular x, viene dada por
S’(x) = -2·k·c·(x -x0) La variable de perturbación k se elimina realizando el cociente entre S’ y S: q(x) = S’(x) / S(x) = -2·c·(x -x0) / (1 -c·(x -x0)2) que da la ecuación diferencial
(1 -c·(x -x0)2) · S’(x) + 2·c(x -x0) · S(x) = 0 En el caso de x = x0, 2·c·(x -x0) · S(x) = 0, pero asimismo S’(x0) = 0 (tangente horizontal). Si para un cierto ángulo x, S’(x) y S(x) pueden ser calculadas o estimadas a partir de los datos del radar, q(x) es un
valor conocido. La relación en el caso de q(x) da una ecuación de segundo orden para la incógnita (x0 -x):
(x0-x)2 + 2 / q·(x0 -x) -1/c = 0 habiendo anotado q(x) simplemente como q. Solamente interesa una de las raíces cuadradas (rc): x0-x =-1/q + rc(1/q2+ 1/c) = q/c / (rc(1 + q/c) + 1) (*) (utilizando la regla de la conjugada) y de este modo x0 = x + q/c / (rc(1 + q/c) + 1)
que, para cualquier valor del ángulo x cerca del máximo x = x0, da una estimación de x0. El algoritmo es (observar que c es un parámetro conocido en la forma supuesta del haz): En cada bin de distancias, para un valor angular x, se obtiene S(x) a partir de los datos medidos y de la estimación
de S’(x). Calcular la relación q = S’(x) / S(x) y sustituir este valor en el lado de la derecha de la última ecuación anterior. A continuación el lado de la derecha, evaluado numéricamente, da una estimación del ángulo x0 más próximo a x, donde el haz tiene un máximo.
El algoritmo es razonable. Si S es una magnitud o potencia, S(x) es siempre positiva. Si q(x), y de este modo la pendiente S’(x), es > 0, entonces el máximo x0 más próximo debería ser > x. Si q(x) < 0, x0 debería ser < x. Dado que el denominador rc(1 + q/c) + 1 anterior es siempre > 0, el algoritmo es cualitativamente tal como se espera.
El algoritmo es de interés principalmente para pequeños valores de q(x), (cerca del máximo, S(x) es grande y S’(x)
es numéricamente pequeño).
Con x = posición angular y x0 interpretado como “posición angular corregida”, la fórmula utilizada es del tipo
posición angular corregida = posición angular + corrección del ángulo
La fórmula proporciona asimismo una expresión para la “corrección del ángulo” (“ángulo de error”)
Final del Ejemplo 1.
Dependiendo de los supuestos sobre la forma del haz S(x), las fórmulas pueden ser ligeramente diferentes. Esto se
mostrará a continuación. Ejemplo 2. Otra aproximación común a la forma del haz en la literatura sobre radares es S(x) = k·cos2(d·(x -x0)) = k/2·(1+ cos(2·d·(x -x0))) En realidad, la aproximación utilizada [Sherman, página 138] para el lóbulo suma de la figura 3, es de este tipo. Al
igual que en Ejemplo 1, la constante k da la amplitud del máximo y la constante d determina la forma del haz (máximo estrecho o amplio). La derivada S’(x) viene dada por
S’(x) = -2·k·d·cos(d·(x -x0))·sen(d·(x -x0)) Eliminando la variable k de perturbación al tomar el cociente q(x) = S’(x) / S(x), se obtiene
q(x) = S’(x) / S(x) = -2·d·tan(d·(x -x0)) = 2·d·tan(d·(x0 -x)) (que da la ecuación diferencial S’(x) + 2·d·tan(d·(x -x0))·S(x) = 0) Resolviendo la incógnita x0 -x se obtiene
d·(x0 -x) = arctan(q(x) / 2/d) y finalmente x0 = x + 1/d·arctan(q/(2·d)) que es la relación deseada, una fórmula del tipo
posición angular corregida = posición angular + corrección del ángulo Final del Ejemplo 2. Las características comunes de los ejemplos son: En cada bin de distancias, para cada valor del ángulo x, la derivada (parcial) S’(x) (y las derivadas de orden superior
S”(x), etc.) en la dirección angular, se determinan (se estiman) a partir de los datos. Se utiliza una fórmula del tipo posición angular corregida = posición angular + corrección del ángulo en la que la expresión “corrección del ángulo” se calcula mediante una fórmula que utiliza S(x), S’(x), … y
parámetros conocidos de la forma del haz. Muy a menudo, en la fórmula solamente aparecen cocientes del tipo S’(x)
/ S(x).
Estas características se utilizan en el método para reducir el desenfoque angular, las cuales serán descritas a
continuación.
El método sugerido no presupone una antena monoimpulso, sino solamente una disposición de una antena ordinaria
(de búsqueda). La amplitud medida (potencia) en el ángulo x y la distancia y se indican mediante S(x, y). Cuando el
valor (fijado) de la distancia no es importante, simplemente se anota S(x). Se supone que el ángulo x y la distancia y
solamente se producen como valores diferenciados (discontinuos). Para incrementar esto, se pueden anotar los
datos medidos como S(i,j), i1 : i : i2, j1 : j : j2.
En cada bin de distancias y para la posición angular x, se utiliza la fórmula de la sección anterior:
posición angular corregida = posición angular + corrección del ángulo
en la que la “corrección del ángulo” viene dada por una fórmula que implica primeras derivadas o derivadas de orden superior (derivadas parciales en la dirección angular) de S(x):
corrección del ángulo = función de S(x), S’(x), … y parámetros conocidos de la forma del haz
Las fórmulas para la “corrección del ángulo” pueden ser ligeramente diferentes, dependiendo de la aproximación de la disposición del haz utilizada. El método sugerido comprende la totalidad de esta clase de fórmulas. En particular, el método comprende la clase de fórmulas que utilizan solamente relaciones del tipo S’(x) / S(x).
El método propuesto puede ser realizado de diversas formas. Una realización ventajosa utiliza el método siguiente:
En cada bin de distancias las derivadas (parciales) dS/dx, S”(x), etc. están calculadas de forma aproximada, utilizando generalmente diferencias finitas en la dirección angular. Si la amplitud en la imagen de radar en la posición angular i y en el “bin” de distancias j se denomina S(i, j), esto significa en el caso de la primera derivada dS/dx (en realidad una derivada parcial a lo largo de la dirección angular “i”) en el “bin” de distancias j:
dS/dx (i,j) = (aproximadamente) (S(i +1,j) -S(i -1,j) / (posición angular (i+1) -posición angular (i -1))
En este caso se utiliza una diferencia simétrica, pero existen otras posibilidades que pueden ser utilizadas en diferentes realizaciones del método. Si las posiciones angulares son aproximadamente equidistantes, el denominador es constante y los cálculos se pueden simplificar.
De manera simbólica, esta realización significa que dS/dx se estima mediante diferencias tales como fS/fu.
Es ventajoso que el número de posiciones angulares por ancho de haz no sea demasiado pequeño, ya que entonces es posible una cierta uniformización numérica (realizando una media, con un filtro contra el paso de bajas frecuencias) de los cocientes diferenciales fS/fu, etc., para obtener una estimación numérica mejor.
Asimismo, es ventajoso evitar utilizar la fórmula cuando el denominador S es demasiado pequeño. Una posibilidad es utilizar la fórmula de corrección solamente en los puntos en los que S es mayor que un umbral adecuado. Otra posibilidad es utilizar la fórmula en todos los puntos de la imagen, pero sustituir el denominador S por S + t, en que t es un umbral positivo adecuado (pequeño).
A continuación, la figura 5 muestra la figura 1 después de una corrección utilizando una de las fórmulas contempladas por el método.
Es evidente, por la figura 5, que el método de corrección sugerido “adelgaza el haz de la antena” y concentra los ecos del radar alrededor de las posiciones angulares correctas.
La figura 6 muestra la puerta de distancia #69 de la figura 5. Comparar con la figura 2. En la figura 6, la mayor parte de la potencia del haz ancho de la figura 2 es elevada y añadida a la posición angular correcta del blanco, posiciones angulares 28 -29.
Otras características del método sugerido:
El método no presupone que el radar sea móvil (no se utiliza una apertura sintetizada).
El método no utiliza un proceso Doppler.
El método sugerido funciona particularmente bien cuando el número de posiciones angulares por ancho de haz es grande. El método es menos adecuado en el caso de números pequeños (<<10) de posiciones angulares por ancho de haz.
El método proporciona una cierta mejora en la resolución angular, interpretada en el sentido de la posibilidad de resolver blancos muy próximos en ángulo y en el mismo bin de distancias. No se reivindica que esta mejora sea grande. La mejora es más bien limitada y tiene sus características particulares.
El método funciona muy bien cuando concentra blancos aislados (puntos) (figuras 5, 6).
La invención no está limitada a los ejemplos facilitados anteriormente, pero puede estar sometida a modificaciones dentro del ámbito de las reivindicaciones adjuntas.
Claims (5)
- REIVINDICACIONES1. Método para la reducción del desenfoque angular debido a la anchura del haz del radar en imágenes de radar logradas en sistemas de detección basados en bins de distancias en los que la amplitud de la potencia S(x) es medida como una función de una posición angular x a una cierta distancia, comprendiendo:la determinación en cada bin de distancias, en el caso de una serie de posiciones angulares x, de una expresión del tipox0 = x + q/c / (rc(1 + q/c) + 1)óx0 = x + (l/d)·arctan (q/(2·d)),en las quec es una constante que determina la forma de un polinomio de segundo orden utilizado para realizar la aproximación de la forma del haz,d es una constante que determina la forma de una aproximación a la forma del haz según S(x) = k·cos2 (d·(x -x0)), siendo k una amplitud,q es la relación q = S’(x) / S(x), siendo S’(x) la primera derivada dS/dx de S(x), y x0 es una posición angular corregida; yse desplaza S(x) desde la posición angular x a la posición angular corregida x0 mediante la adición a S(x0).
-
- 2.
- Método, según la reivindicación anterior, caracterizado porque las derivadas son estimaciones numéricas a partir de los valores medidos de S(x).
-
- 3.
- Método, según la reivindicación 2, caracterizado porque para la estimación de las derivadas S’(x) se utilizan diferencias finitas en la dirección angular fS/fu.
-
- 4.
- Método, según la reivindicación 3, caracterizado porque se utilizan diferencias finitas promediadas o filtradas contra el paso de bajas frecuencias.
-
- 5.
- Método, según cualquiera de las reivindicaciones anteriores, caracterizado porque la información de la variable S(x) solamente se utiliza para unos valores tales de x en los que la variable S(x) sobrepasa un cierto umbral t dado.
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP04445004A EP1560038B1 (en) | 2004-01-27 | 2004-01-27 | Method for reducing angular blur in radar pictures |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| ES2392512T3 true ES2392512T3 (es) | 2012-12-11 |
Family
ID=34639520
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| ES04445004T Expired - Lifetime ES2392512T3 (es) | 2004-01-27 | 2004-01-27 | Método para la reducción del desenfoque angular en imágenes de radar |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US7236125B2 (es) |
| EP (1) | EP1560038B1 (es) |
| ES (1) | ES2392512T3 (es) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8044846B1 (en) * | 2007-11-29 | 2011-10-25 | Lockheed Martin Corporation | Method for deblurring radar range-doppler images |
| JP5186033B1 (ja) * | 2011-10-12 | 2013-04-17 | ダンロップスポーツ株式会社 | ゴルフクラブシャフト及びゴルフクラブ |
| US9182482B2 (en) * | 2011-10-25 | 2015-11-10 | Navico Holding As | Radar beam sharpening system and method |
| US10705201B1 (en) * | 2015-08-31 | 2020-07-07 | Rockwell Collins, Inc. | Radar beam sharpening system and method |
| US11391831B2 (en) * | 2018-04-25 | 2022-07-19 | Qualcomm Incorporated | Association aware radar beamforming |
| US11047955B2 (en) * | 2018-07-06 | 2021-06-29 | Lyft, Inc. | Calibrating a radar antenna |
Family Cites Families (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3710384A (en) * | 1969-11-17 | 1973-01-09 | North American Rockwell | Signal-correlating apparatus for improving the angular resolution of a directionally ranging system |
| US6201496B1 (en) * | 1981-01-12 | 2001-03-13 | Northrop Grumman Corporation | Target angle estimator for search radars |
| USH1181H (en) * | 1985-05-16 | 1993-05-04 | The United States Of America As Represented By The Secretary Of The Air Force | Method for high resolution radar imagery and accurate dimensional measurements |
| US4717916A (en) * | 1986-05-16 | 1988-01-05 | Holodyne Ltd., 1986 | High resolution imaging doppler interferometer |
| US5610610A (en) * | 1988-05-18 | 1997-03-11 | Hughes Electronics | Inverse synthetic array radar system and method |
| US5184133A (en) * | 1991-11-26 | 1993-02-02 | Texas Instruments Incorporated | ISAR imaging radar system |
| US5243351A (en) * | 1992-06-25 | 1993-09-07 | Hughes Aircraft Company | Full aperture image synthesis using rotating strip aperture image measurements |
| US5644508A (en) * | 1995-06-05 | 1997-07-01 | Hughes Electronics | Detection of moving objects |
| US6155704A (en) * | 1996-04-19 | 2000-12-05 | Hughes Electronics | Super-resolved full aperture scene synthesis using rotating strip aperture image measurements |
| US6192322B1 (en) * | 1996-04-19 | 2001-02-20 | Raytheon Company | Moving object and transient event detection using rotation strip aperture image measurements |
| ATE237811T1 (de) * | 1996-05-14 | 2003-05-15 | Honeywell Int Inc | Autonomes landeführungssystem |
| USH1720H (en) * | 1997-03-31 | 1998-04-07 | Chen; Victor C. | Time frequency processor for radar imaging of moving targets |
| US6002983A (en) * | 1997-08-27 | 1999-12-14 | Delphi Technologies, Inc. | Angle extent estimation method for a motor vehicle object detection system |
| US6356227B1 (en) * | 1999-09-16 | 2002-03-12 | Honeywell International, Inc. | Smearing compensation apparatus for a radar system |
| US6661369B1 (en) * | 2002-05-31 | 2003-12-09 | Raytheon Company | Focusing SAR images formed by RMA with arbitrary orientation |
-
2004
- 2004-01-27 ES ES04445004T patent/ES2392512T3/es not_active Expired - Lifetime
- 2004-01-27 EP EP04445004A patent/EP1560038B1/en not_active Expired - Lifetime
-
2005
- 2005-01-25 US US11/041,244 patent/US7236125B2/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| EP1560038A1 (en) | 2005-08-03 |
| EP1560038B1 (en) | 2012-08-15 |
| US7236125B2 (en) | 2007-06-26 |
| US20060066475A1 (en) | 2006-03-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fienup | Detecting moving targets in SAR imagery by focusing | |
| Davidson et al. | Signal properties of spaceborne squint-mode SAR | |
| US8698668B2 (en) | SAR radar system | |
| Borden | Radar Imaging of Airborne Targets: A primer for applied mathematicians and physicists | |
| Smith et al. | Robust through-the-wall radar image classification using a target-model alignment procedure | |
| ES2392512T3 (es) | Método para la reducción del desenfoque angular en imágenes de radar | |
| CN116299227B (zh) | Sar雷达二维天线方向图测量方法、系统、介质及设备 | |
| US5270718A (en) | Method and apparatus for tracking targets from direct and multipath reflected radar signals | |
| McDonald et al. | Elevation angle-of-arrival determination for a standard and a modified SuperDARN HF radar layout | |
| US20180074180A1 (en) | Ultrafast target detection based on microwave metamaterials | |
| Monti Guarnieri et al. | Reduction to monostatic focusing of bistatic or motion uncompensated SAR surveys | |
| CN115166749B (zh) | 一种基于单个相控阵天气雷达获取三维风场的方法及设备 | |
| Smith et al. | Dual carrier probing for spaceborne SAR imaging | |
| Lagovsky et al. | Regression methods of obtaining angular superresolution | |
| Garnier et al. | Coherent interferometric imaging for synthetic aperture radar in the presence of noise | |
| Marques et al. | Optimal detection and imaging of moving objects with unknown velocity | |
| Slattery | Use of mills cross receiving arrays in radar systems | |
| Artyushenko et al. | Measurement error estimation of motion variables for extended objects under changing range conditions | |
| Klochko | Algorithms of 3D radio-wave imaging in airborne Doppler radar | |
| Karaev et al. | Waveform of the Reflected Impulse at the Oblique Sounding of the Sea Surface | |
| Gressani et al. | Wave spectrometer tilt modulation transfert function using near-nadir Ku-and Ka-band GPM radar measurements | |
| Vorobev et al. | System for adjustment of angle coordinates for sea surface surveillance radar | |
| Verdone et al. | Processing algorithms for COSMO-SkyMed SAR sensor | |
| Li et al. | SAR image correction of moving ships in marine scene based on wakes | |
| Kulpa et al. | Forward scattering effect exploitation in passive radars |