RU2274877C2 - Panoramic radar method of determining condition of ocean's layer surface from satellite - Google Patents
Panoramic radar method of determining condition of ocean's layer surface from satellite Download PDFInfo
- Publication number
- RU2274877C2 RU2274877C2 RU2003124710/09A RU2003124710A RU2274877C2 RU 2274877 C2 RU2274877 C2 RU 2274877C2 RU 2003124710/09 A RU2003124710/09 A RU 2003124710/09A RU 2003124710 A RU2003124710 A RU 2003124710A RU 2274877 C2 RU2274877 C2 RU 2274877C2
- Authority
- RU
- Russia
- Prior art keywords
- cell
- elementary
- pulses
- radar
- satellite
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 29
- 239000006185 dispersion Substances 0.000 claims abstract description 21
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 239000002344 surface layer Substances 0.000 claims abstract description 11
- 230000005855 radiation Effects 0.000 claims description 17
- 239000000523 sample Substances 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 239000002245 particle Substances 0.000 abstract 4
- 239000000126 substance Substances 0.000 abstract 1
- 238000004891 communication Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 230000009897 systematic effect Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- OMPJBNCRMGITSC-UHFFFAOYSA-N Benzoylperoxide Chemical compound C=1C=CC=CC=1C(=O)OOC(=O)C1=CC=CC=C1 OMPJBNCRMGITSC-UHFFFAOYSA-N 0.000 description 1
- 206010038743 Restlessness Diseases 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 229940060894 topex Drugs 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- 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
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
Изобретение относится к радиолокации, а именно к радиолокационным методам определения параметров морского волнения, и может быть использовано в метеорологии и океанологии для дистанционного зондирования приповерхностного слоя океанов со спутника.The invention relates to radar, and in particular to radar methods for determining the parameters of sea waves, and can be used in meteorology and oceanology for remote sensing of the surface layer of the oceans from a satellite.
Обеспечение оперативного мониторинга состояния приповерхностного слоя океана со спутника в возможно более широкой полосе обзора является важным для составления достоверных прогнозов погоды, контроля за глобальными изменениями климата, обеспечения безопасности жизнедеятельности в прибрежных районах, изучения Мирового океана и решения множества других задач. При этом для действительно оперативного мониторинга желательно обеспечить построение с достаточным пространственным разрешением двумерного изображения водной поверхности, для чего необходимо определить (восстановить) в каждой "i"-ой элементарной ячейке полосы обзора дисперсию наклонов σ2 i водной поверхности и направление распространения крупномасштабного волнения φwi, а также восстановить в каждой ячейке скорость ветра V.Providing on-line monitoring of the state of the near-surface layer of the ocean from a satellite in the widest possible viewing range is important for making reliable weather forecasts, monitoring global climate changes, ensuring life safety in coastal areas, studying the World Ocean, and solving many other problems. Moreover, for truly operational monitoring, it is desirable to ensure the construction with a sufficient spatial resolution of a two-dimensional image of the water surface, for which it is necessary to determine (restore) the dispersion of the slopes σ 2 i of the water surface and the direction of propagation of large-scale waves φ wi in each “i” elementary cell of the viewing band , and also restore wind speed V in each cell.
Известен способ определения параметров водной поверхности со спутника, оснащенного антенной радиолокатора с двухлучевой диаграммой направленности (D.Hauser et all., "Swimsat: A Real-Aperture Radar to Measure Directional Spectra of Ocean Waves from Space - Main Characteristics and Performance Simulation", J. Atmospheric and Oceanic Technology, 2001, v.18, p.p.421-437). Данный способ состоит в том, что с помощью двухлучевой антенны излучают две независимые последовательности коротких зондирующих импульсов в разных направлениях: первый луч (последовательность импульсов) направляют по нормали (т.е. в надир) к подстилающей водной поверхности, направление второго луча составляет 10° по отношению к первому лучу и этот второй луч вращают при движении спутника по отношению к надиру. На водной поверхности оба луча освещают пятна (оставляют след) с характерными размерами 18×18 км, а расстояние между пятнами (т.е. радиус участка, с которого собирают информацию) при высоте орбиты спутника порядка 500 км, составляет 88 км. Отраженную назад от водной поверхности часть мощности каждого зондирующего импульса обоих лучей принимают соответствующей антенной и записывают форму принимаемых отраженных импульсов. При сеансах связи эту информацию передают на станцию слежения, на которой с помощью компьютера проводят обработку принятых импульсов. Принятую последовательность импульсов первого луча обрабатывают по известному алгоритму, как в известных спутниковых альтиметрах (см., например, Alfred R. Zieger at all., NASA Radar Altimeter for TOPEX/POSEIDON Project, Proceedings of the IEEE, Vol.79, NO.6, June 1991), в результате чего по максимальному значению принятой мощности Рmax каждого отраженного импульса (по форме импульса) определяют сечение обратного рассеяния σо, по которому определяют (восстанавливают) величину скорости V приповерхностного ветра. При этом скорость ветра V восстанавливают с систематической ошибкой, обусловленной неоднозначностью связи между отраженной мощностью Рmax и скоростью ветра V, поскольку известно, что отраженная мощность зависит также и от дисперсии наклонов σ2 i водной поверхности, которая в свою очередь зависит от присутствия волн зыби, не связанных с ветром в точке измерения. По наклону переднего фронта каждого принятого импульса определяют высоту Н значительного волнения в освещенном пятне. Принятую последовательность импульсов второго луча обрабатывают как в самолетном измерителе спектра волн, однако перенос данного радиолокатора на спутник существенно ухудшает его разрешение, поскольку позволяет измерять лишь спектр волн, длина которых более 70 м, а высота - более 2 м, т.е. теряется информация о высоте и наклонах волн короче 70 м. К недостаткам данного способа относится и то, что в основу всех расчетов положено предположение об однородности волнений в пятне (элементе разрешения) с радиусом 88 км, что является некорректным с точки зрения океанологов, поскольку современные стандартные модели открытого океана используют сетку с элементом разрешения (элементарной ячейкой) для открытого океана 50×50 км и для прибрежных районов 28×28 км (см., например, "Satellites, Oceanography and Society" edited by D.Halpern. Elsevier, Amsterdam, p.p.35-56, 2000).A known method for determining water surface parameters from a satellite equipped with a dual-beam radar antenna (D. Hauser et all., "Swimsat: A Real-Aperture Radar to Measure Directional Spectra of Ocean Waves from Space - Main Characteristics and Performance Simulation", J Atmospheric and Oceanic Technology, 2001, v. 18, pp421-437). This method consists in the fact that using a two-beam antenna emit two independent sequences of short probe pulses in different directions: the first beam (a sequence of pulses) is directed along the normal (i.e. to nadir) to the underlying water surface, the direction of the second beam is 10 ° in relation to the first ray and this second ray rotate when the satellite moves in relation to the nadir. On the water surface, both beams illuminate spots (leave a mark) with characteristic dimensions of 18 × 18 km, and the distance between the spots (i.e. the radius of the area from which information is collected) at a satellite orbit altitude of about 500 km is 88 km. The part of the power of each probe pulse of both beams reflected back from the water surface is received by the corresponding antenna and the shape of the received reflected pulses is recorded. During communication sessions, this information is transmitted to a tracking station, where the received pulses are processed using a computer. The received sequence of pulses of the first beam is processed according to a well-known algorithm, as in well-known satellite altimeters (see, for example, Alfred R. Zieger at all., NASA Radar Altimeter for TOPEX / POSEIDON Project, Proceedings of the IEEE, Vol. 79, NO.6 , June 1991), as a result of which the backscattering cross section σ о is determined from the maximum value of the received power P max of each reflected pulse (by the shape of the pulse), from which the velocity V of the surface wind is determined (restored). In this case, the wind speed V is restored with a systematic error due to the ambiguity of the relationship between the reflected power P max and the wind speed V, since it is known that the reflected power also depends on the variance of the slopes σ 2 i of the water surface, which in turn depends on the presence of swell waves not related to the wind at the measurement point. The slope of the leading edge of each received pulse determines the height H of significant excitement in the illuminated spot. The adopted sequence of pulses of the second beam is processed as in an airplane wave spectrum meter, however, the transfer of this radar to the satellite significantly affects its resolution, since it only measures the spectrum of waves with a length of more than 70 m and a height of more than 2 m, i.e. information about wave heights and inclinations shorter than 70 m is lost. The disadvantages of this method include the fact that all calculations are based on the assumption that waves are uniform in a spot (resolution element) with a radius of 88 km, which is incorrect from the point of view of oceanologists, since modern standard open ocean models use a grid with a resolution element (unit cell) for the open ocean 50 × 50 km and for coastal areas 28 × 28 km (see, for example, "Satellites, Oceanography and Society" edited by D. Halpern. Elsevier, Amsterdam , pp35-56, 2000).
Наиболее близким к заявленному способу по технической сущности является панорамный способ определения величины и направления скорости приповерхностного ветра над водной поверхностью со спутника с помощью скаттерометра (радиолокатора для определения скорости ветра), который выбран в качестве прототипа (пат. США №6137437 МПК7 G 01 S 13/60, публ. 24.10.2000). Способ - прототип заключается в том, что с помощью однолучевой вращающейся антенны с ножевой диаграммой направленности, т.е. узкой (1°-3°) в поперечном направлении и широкой протяженной (20°-25°) в продольном направлении, осуществляют наклонное зондирование (при углах падения от 28° до 51°) водной поверхности, что позволяет получить широкую (порядка 1500 км) полосу обзора, т.е. осуществить панорамный режим работы радиолокатора. Этой же антенной радиолокатора принимают последовательность отраженных импульсов, с помощью регистрирующего устройства записывают их форму и передают информацию на станцию слежения. Для обработки полученной информации используют временную селекцию по дальности, с помощью которой формируют элемент разрешения 50×50 км. Для восстановления скорости ветра (и его направления) применяется псевдооднопараметрический алгоритм, который получен методом регрессионного анализа измеряемых со спутника скаттерометрических данных о величине одного параметра - сечения обратного рассеяния σo и буйковых данных о соответствующих скоростях ветра V, при этом второй параметр - дисперсия наклонов σ2 i водной поверхности, также влияющий на зависимость σo от скорости ветра V не учитывается, поскольку в настоящее время его еще не умеют измерять из космоса.Closest to the claimed method in technical essence is a panoramic method for determining the magnitude and direction of the surface wind speed above the water surface from a satellite using a scatterometer (radar for determining wind speed), which is selected as a prototype (US Pat. No. 6,137,437 IPC 7 G 01 S 13/60, publ. 24.10.2000). The prototype method consists in using a single-beam rotating antenna with a knife radiation pattern, i.e. narrow (1 ° -3 °) in the transverse direction and a wide extended (20 ° -25 °) in the longitudinal direction, oblique sounding (at incidence angles from 28 ° to 51 °) of the water surface is carried out, which allows to obtain a wide (about 1500 km ) the span, i.e. to carry out a panoramic mode of operation of the radar. The same radar antenna receives a sequence of reflected pulses, using a recording device to record their shape and transmit information to the tracking station. To process the obtained information, temporary range selection is used, with the help of which a resolution element of 50 × 50 km is formed. To restore the wind speed (and its direction), a pseudo-one-parameter algorithm is used, which is obtained by the method of regression analysis of the scatterometric data measured from the satellite about the value of one parameter - the backscattering cross section σ o and the buoy data on the corresponding wind speeds V, and the second parameter is the slope dispersion 2 i of the water surface, which also affects the dependence of σ o on the wind speed V is not taken into account, because at present they still cannot measure it from space.
Таким образом, недостатком способа-прототипа является то, что он позволяет измерять (восстанавливать) только параметры ветра и не дает информацию о параметрах волнения водной поверхности. При этом для восстановления скорости ветра V используется полученная опытным путем зависимость V=f(σo), которая не является однозначной и поэтому дает систематическую ошибку.Thus, the disadvantage of the prototype method is that it allows you to measure (restore) only the wind parameters and does not provide information about the parameters of the waves of the water surface. In this case, to restore the wind speed V, the experimentally obtained dependence V = f (σ o ) is used, which is not unique and therefore gives a systematic error.
Задачей, на решение которой направлено настоящее изобретение, является разработка панорамного радиолокационного способа определения состояния приповерхностного слоя океана, обеспечивающего измерение в каждой "i"- ой элементарной ячейке широкой полосы обзора следующих параметров: дисперсии наклонов σi 2 водной поверхности, направления распространения φwi крупномасштабного волнения, а также более точное восстановление скорости V приповерхностного ветра за счет исключения вышеуказанной систематической ошибки.The problem to which the present invention is directed is the development of a panoramic radar method for determining the state of the near-surface layer of the ocean, ensuring the measurement of the following parameters in each "i" -th unit cell of a wide viewing band: the dispersion of the slopes σ i 2 of the water surface, the propagation direction φ wi of large-scale unrest, as well as a more accurate restoration of the velocity V of the surface wind due to the elimination of the above systematic error.
Технический результат в разработанном способе достигается тем, что разработанный панорамный радиолокационный способ определения параметров состояния приповерхностного слоя океана со спутника, включает излучение зондирующих импульсов микроволнового диапазона доплеровским радиолокатором, снабженным однолучевой вращающейся антенной с ножевой диаграммой направленности, прием отраженных от водной поверхности импульсов, регистрацию их формы определение сечения обратного рассеяния σO и последующее вычисление по алгоритму скорости приповерхностного ветра V.The technical result in the developed method is achieved by the fact that the developed panoramic radar method for determining the state parameters of the surface layer of the ocean from a satellite includes the emission of probe pulses of the microwave range with a Doppler radar equipped with a single-beam rotating antenna with a knife radiation pattern, reception of pulses reflected from the water surface, registration of their shape determining the backscatter cross section σ O and subsequent calculation using the algorithm speed n ipoverhnostnogo wind V.
Новым в разработанном панорамном радиолокационном способе является то, что ножевую диаграмму направленности антенны вращают вокруг вертикальной оси симметрии антенны, проходящей через ее центр, зондирующие импульсы направляют к поверхности океана в надир и каждым импульсом освещают на водной поверхности пятно с размерами, например, 14×355 км, а при приеме отраженных импульсов используют одновременно временную и доплеровскую селекции по дальности для выделения в упомянутом пятне 14×355 км элементарных рассеивающих ячеек с размерами, например 14×14 км. Затем определяют для каждых двух последовательных "i"-ой и "i+1"-ой ячеек сечения обратного рассеяния σ0 (θi) и σ0 (θi+1), которые корректируют с учетом гауссовой диаграммы направленности антенны в соответствии с выражениемNew in the developed panoramic radar method is that the knife antenna radiation pattern is rotated around the vertical axis of symmetry of the antenna passing through its center, the probe pulses are sent to the nadir to the ocean surface and each spot illuminates a spot with a size of, for example, 14 × 355 km, and when receiving reflected pulses, both time and Doppler distance selection are used simultaneously to select 14 × 355 km of elementary scattering cells with sizes, for example, in the said spot measures 14 × 14 km. Then, for every two consecutive “i” and “i + 1” cells, the backscatter cross sections σ 0 (θ i ) and σ 0 (θ i + 1 ) are determined, which are adjusted taking into account the Gaussian antenna pattern in accordance with expression
δxi - ширина диаграммы направленности на уровне 0,5 по мощности вдоль направления зондирования, θi - угол падения для "i"-ой элементарной ячейки, σ0 (θi) - сечение обратного рассеяния "i"-ой элементарной ячейки. Затем определяют дисперсию наклонов σi 2(φj) вдоль азимутального направления зондирования φj для каждой "i"-ой ячейки:δ xi is the width of the radiation pattern at the power level 0.5 along the sounding direction, θ i is the incidence angle for the “i” th unit cell, σ 0 (θ i ) is the backscattering cross section of the “i” th unit cell. Then determine the variance of the slopes σ i 2 (φ j ) along the azimuthal direction of sounding φ j for each "i" -th cell:
где θi и θi+1 - углы падения для двух последовательных элементарных ячеек, σ0к (θi) и σ0к (θi+1) скорректированные сечения обратного рассеяния этих ячеек соответственно, j - номер азимутального направления, под которым наблюдается каждая элементарная ячейка.where θ i and θ i + 1 are the angles of incidence for two consecutive unit cells, σ 0к (θ i ) and σ 0к (θ i + 1 ) are the adjusted backscattering cross sections of these cells, respectively, j is the number of the azimuthal direction under which each unit cell.
После чего полную дисперсию наклонов для "i"-ой ячейки σi 2 определяют из соотношения σi 2=σi 2(φj)+σ2 i(φj+90°). Далее по графику азимутальной зависимости дисперсии наклонов σi 2(φj) определяют направление φwi распространения в "i"-ой ячейке крупномасштабного волнения, а упомянутую скорость приповерхностного ветра V определяют по сечению обратного рассеяния σ0 и дисперсиям наклонов σ2 i(φj) и σ2 i(φj+90°) с помощью алгоритма V=F[σ0, σ2 i(φj), σ2 i(φj), σ2 i(φj+90°)], полученного стандартным методом регрессии.After that, the total variance of the slopes for the “i” cell σ i 2 is determined from the relation σ i 2 = σ i 2 (φ j ) + σ 2 i (φ j + 90 °). Further, according to the graph of the azimuthal dependence of the slope dispersion σ i 2 (φ j ), the direction of propagation φ wi in the “i” cell of large-scale waves is determined, and the surface wind speed V is determined by the backscattering cross section σ 0 and the slope dispersions σ 2 i (φ j ) and σ 2 i (φ j + 90 °) using the algorithm V = F [σ 0 , σ 2 i (φ j ), σ 2 i (φ j ), σ 2 i (φ j + 90 °)] obtained by the standard regression method.
Целесообразно для более детального изучения отдельных участков водной поверхности сначала с помощью дополнительной временной селекции по дальности сформировать элементарные ячейки с размерами, например, 14×1 км, затем с помощью дополнительной доплеровской селекции по азимутальному углу φj выделить элементарные ячейки с размерами, например, 5×1 км.It is advisable for a more detailed study of individual sections of the water surface, first using additional time selection in range to form elementary cells with dimensions, for example, 14 × 1 km, then using additional Doppler selection in azimuthal angle φ j, select elementary cells with dimensions, for example, 5 × 1 km.
В другом частном случае для определения средней длины волны Lm крупномасштабного волнения в каждой "i"-ой элементарной ячейке полосы обзора целесообразно создать по доплеровским скоростям с помощью дополнительных частотных фильтров, установленных в радиолокаторе, начальную точку отсчета для регистрации отраженных импульсов в каждой "i"-ой ячейке. После этого прием отраженных от водной поверхности импульсов и регистрацию их формы провести для каждой "i"-ой ячейки отдельно и стандартным алгоритмом по тангенсу угла наклона в средней точке переднего фронта импульса, зарегистрированного в "i"-ой ячейке, определить высоту Н значительного волнения в этой ячейке. Уточненное значение высоты Н значительного волнения в "i"-ой ячейке полосы обзора целесообразно получить путем вычисления среднего значения высоты Н по нескольким упомянутым измерениям в данной ячейке. После чего по упомянутому измеренному максимальному значению дисперсии наклонов σw 2 в "i"-ой ячейке на графике фиг.6 и по уточненному значению высоты Н значительного волнения среднюю длину Lm крупномасштабного волнения следует определить из соотношения Lm=H/σw.In another particular case, in order to determine the average wavelength L m of large-scale waves in each "i" -th unit cell of the span, it is advisable to create the starting point for recording reflected pulses in each "i using Doppler speeds using additional frequency filters installed in the radar "th cell. After that, the reception of the pulses reflected from the water surface and the registration of their shape should be carried out separately for each "i" -th cell and using the standard algorithm for determining the slope at the midpoint of the leading edge of the pulse recorded in the "i" -th cell, determine the height H of significant waves in this cell. It is advisable to obtain an updated value of the height H of significant disturbance in the "i" -th cell of the field of view by calculating the average value of the height H from the several mentioned measurements in this cell. Then, according to the aforementioned measured maximum value of the slope dispersion σ w 2 in the "i" -th cell in the graph of Fig. 6 and the specified value of the height H of significant waves, the average length L m of large-scale waves should be determined from the relation L m = H / σ w .
На фиг.1 представлен вариант блок-схемы устройства для реализации разработанного способа.Figure 1 presents a variant of a block diagram of a device for implementing the developed method.
На фиг.2 представлена схема наблюдения поверхности океана в соответствии с разработанным способом. Направление полета выбрано вдоль оси Y и введена вторая система координат X1 и Y1, связанная с антенной и вращающаяся вместе с ней. Ось X1 ориентирована вдоль продольной оси освещенного пятна на водной поверхности. Угол поворота φj оси X1 относительно первой системы координат XY является азимутальным углом.Figure 2 presents a diagram of the observation of the surface of the ocean in accordance with the developed method. The flight direction is chosen along the Y axis and a second coordinate system X 1 and Y 1 is introduced, connected with the antenna and rotating with it. Axis X 1 is oriented along the longitudinal axis of the illuminated spot on the water surface. The rotation angle φ j of the axis X 1 relative to the first coordinate system XY is the azimuthal angle.
Фиг.3 иллюстрирует вращение антенны по азимутальному углу φj во время движения спутника и дискретные измерения в элементарной "i"-ой ячейке полосы обзора под тремя азимутальными углами. Полоса обзора показана пунктирными линиями.Figure 3 illustrates the rotation of the antenna along the azimuthal angle φ j during satellite motion and discrete measurements in the elementary "i" -th cell of the span at three azimuthal angles. The line of sight is indicated by dashed lines.
Фиг.4 иллюстрирует наблюдение двух последовательных "i"-ой и "i+1"-ой элементарных ячеек вдоль взаимно перпендикулярных направлений зондирования под углом φj и φj+90° во время движения спутника для определения дисперсии σ2 i(φj) и σ2 i(φj+90°).Figure 4 illustrates the observation of two consecutive "i" and "i + 1" -th unit cells along mutually perpendicular sensing directions at an angle φ j and φ j + 90 ° during satellite motion to determine the dispersion σ 2 i (φ j ) and σ 2 i (φ j + 90 °).
На фиг.5 представлен пример разбиения полосы обзора на элементарные рассеивающие ячейки.Figure 5 presents an example of dividing the span into elementary scattering cells.
На фиг.6 представлен график азимутальной зависимости дисперсии наклонов σ2 i(φj) в "i"-ой элементарной ячейке.Figure 6 presents a graph of the azimuthal dependence of the slope variance σ 2 i (φ j ) in the "i" -th unit cell.
Один из вариантов устройства для реализации разработанного способа, представленный на фиг.1, содержит антенну 1, выполненную в виде фазированной антенной решетки с ножевой вращающейся диаграммой направленности. Антенна 1 через блок 2 электронного управления вращением диаграммы направленности соединена с импульсным доплеровским радиолокатором 3 сантиметрового диапазона, который в свою очередь соединен с регистрирующим устройством 4, выход которого соединен со входом приемопередатчика 5, обеспечивающего связь со станцией слежения (на чертеже не показана) на Земле. При этом управление всеми элементами устройства для реализации разработанного способа осуществляется с помощью блока управления 6, который имеет соединение с блоком 2 электронного управления, с радиолокатором 3, регистрирующим устройством 4 и приемопередатчиком 5.One of the variants of the device for implementing the developed method, presented in figure 1, contains an
В качестве антенны 1 с вращающейся ножевой диаграммой направленности может быть использована, например, фазированная антенная решетка, разработанная Опытно-конструкторским бюро Московского энергетического института. В другом варианте изготовления устройства для реализации разработанного способа может быть использована вращающаяся однолучевая, щелевая антенна МИУС отечественного производства или однолучевая, щелевая антенна PR-5 чешского производства (фирма TESLA). Ножевая диаграмма направленности антенны обеспечивается за счет выбора размеров щели, например может быть сформирована ножевая диаграмма направленности с угловыми размерами 1°×25°. При высоте полета 800 км такая антенна освещает пятно на поверхности океана с размерами 14×355 км. Вращение диаграммы направленности антенны 1 вокруг вертикальной оси осуществляют с угловой скоростью порядка 6 об/мин.As the
В качестве доплеровского радиолокатора 3 может быть использован, например, доплеровский измеритель скорости и сноса ДИИС (г.Каменск-Уральский) или доплеровский радиолокатор, изготавливаемый в ЦНИИ "Комета" г.Москва. В качестве регистрирующего устройства 4 и приемопередатчика 5 могут быть использованы любые стандартные устройства подобного назначения, эксплуатируемые на спутниках в настоящее время для записи информации и передачи ее на Землю. Импульсный доплеровский радиолокатор 3 обеспечивает длительность зондирующих импульсов от 100 нс до 2000 нс с частотой повторения порядка 5 кГц при длине волны радиолокатора 2,1 см.As a Doppler
Разработанный панорамный радиолокационный способ определения параметров состояния приповерхностного слоя океана осуществляют следующим образом.The developed panoramic radar method for determining the state parameters of the surface layer of the ocean is as follows.
В ходе полета спутника над акваторией океана посредством антенны 1 с вращающейся ножевой диаграммой направленности, управляемой блоком 2 электронного управления, и доплеровского радиолокатора 3, управляемого блоком 6 (см. фиг.1), излучают последовательность зондирующих импульсов по нормали к водной поверхности (см. фиг.2). При вращении ножевой диаграммы направленности с указанными угловыми размерами 1°×25° и скоростью вращения 6 об/мин на водной поверхности освещают полосу обзора шириной 355 км (для указанной высоты полета 800 км), что позволяет получить изображение происходящего на поверхности океана процесса. При этом каждый отдельный импульс излучения радиолокатора 3 освещает при указанной высоте полета пятно с размерами порядка 14×355 км (см. фиг.3). Этой же антенной 1 с радиолокатором 3 принимают последовательность отраженных от водной поверхности импульсов, которые обрабатывают и запоминают в регистрирующем устройстве 4 до очередного сеанса связи со станцией слежения на Земле. С помощью приемопередатчика 5 осуществляют обмен информацией со станцией слежения, при этом принятые команды поступают в блок управления 6, с помощью которого координируют работу всего устройства для реализации разработанного способа.During the flight of the satellite over the ocean by means of an
Обработка с помощью регистрирующего устройства 4 отраженных от водной поверхности зондирующих импульсов заключается в том, что с помощью известных временной и доплеровской селекции по дальности делят упомянутое освещенное пятно с размерами 14×355 км на элементарные "i"-ые (i меняется от 1 до N) рассеивающие ячейки, т.е. деление освещенного пятна проводят вдоль упомянутой (см. фиг.2) продольной оси X1 этого пятна, которая вращается вместе с антенной 1 по азимутальному углу φj (см. фиг.3 и 4). Для каждой "i"-ой ячейки определяют мощность отраженного сигнала и сечение обратного рассеяния σ0 (θi), где θi - угол падения зондирующего излучения на "i"-ую ячейку (см. фиг.2). Эта информация, как было отмечено выше, хранится в регистрирующем устройстве 4 до очередного сеанса связи со станцией слежения. Выбрав шаг по углу падения θ примерно равным 1°, получают деление полосы обзора на элементарные ячейки размером, например, 14×14 км (см. фиг.5). Далее уже на станции слежения проводят коррекцию по мощности принятого сигнала от каждой "i"-ой ячейки с учетом формы диаграммы направленности антенны 1, которая принята гауссовой и задается следующим выражением:The processing by means of a recording device 4 of probe pulses reflected from the water surface consists in the fact that using the known time and Doppler selection by distance, the mentioned illuminated spot with dimensions 14 × 355 km is divided into elementary “i” -ths (i varies from 1 to N ) scattering cells, i.e. division of the illuminated spot is carried out along the aforementioned (see FIG. 2) longitudinal axis X 1 of this spot, which rotates with the
и - ширина диаграммы направленности на уровне 0,5 по осям X1 и Y1 соответственно, R0 - высота полета спутника. Для коррекции умножают мощность отраженного сигнала в "i"-ой элементарной ячейке (сечение обратного рассеяния σ0 (θi)) на коэффициент, связанный с диаграммой направленности антенны 1. Формула перерасчета имеет следующий вид: and - the width of the radiation pattern at the level of 0.5 along the axes X 1 and Y 1, respectively, R 0 - the height of the satellite. For correction, multiply the power of the reflected signal in the "i" -th unit cell (backscattering cross section σ 0 (θ i )) by the coefficient associated with the radiation pattern of
После выполнения коррекции, как установлено авторами, дисперсия наклонов вдоль направления зондирования (оси X1) вычисляется по следующей формуле:After performing the correction, as established by the authors, the variance of the slopes along the sounding direction (axis X 1 ) is calculated by the following formula:
θi и θi+1 - углы падения для двух последовательных "i"-ой и "i+1"-ой элементарных ячеек вдоль направления зондирования (азимутального угла φj) см. фиг.4; σ0к (θi) и σ0к (θi+1) - скорректированные сечения обратного рассеяния этих ячеек соответственно, j - номер азимутального направления, под которым наблюдается каждая элементарная ячейка.θ i and θ i + 1 are the angles of incidence for two consecutive "i" and "i + 1" th unit cells along the sounding direction (azimuthal angle φ j ), see Fig. 4; σ 0к (θ i ) and σ 0к (θ i + 1 ) are the adjusted backscatter cross sections of these cells, respectively, j is the azimuthal direction number under which each unit cell is observed.
Как видно из приведенной формулы (в знаменателе формулы сечение обратного рассеяния σ0к(θi+1) в "i+1"-ой ячейке делится на сечение обратного рассеяния σ0к(θj) в "i"-ой ячейке) при обработке экспериментальных данных со спутника в разработанном способе важны (используются) не абсолютные значения отраженной мощности (сечения обратного рассеяния) в каждой элементарной отражающей ячейке, а относительные изменения отраженной мощности - в соседних ячейках. За счет этого в разработанном способе удается избавиться от таких существенных проблем, как необходимость регулярной калибровки радиолокатора по мощности, необходимость учета ослабления отраженного сигнала дождевыми облаками, необходимость поддержания стабильной мощности радиолокатора в течение всего срока эксплуатации.As can be seen from the above formula (in the denominator of the formula, the backscattering cross section σ 0к (θ i + 1 ) in the "i + 1" cell is divided by the backscattering cross section σ 0к (θ j ) in the "i" th cell) during processing The experimental data from the satellite in the developed method are important (used) not the absolute values of the reflected power (backscattering cross section) in each elementary reflecting cell, but the relative changes in the reflected power in neighboring cells. Due to this, in the developed method, it is possible to get rid of such significant problems as the need for regular radar calibration by power, the need to take into account the attenuation of the reflected signal by rain clouds, the need to maintain stable radar power over the entire life cycle.
Известно, что с достаточной степенью точности полная дисперсия наклонов поверхности представляет собой сумму дисперсий наклонов, измеренных в двух взаимно перпендикулярных направлениях. Поэтому для определения полной дисперсии наклонов σi 2 поверхности в каждой "i"-ой элементарной ячейке измеряют дисперсию наклонов в этой ячейке под разными азимутальными углами φj, используя вращение антенны 1 во время движения спутника (см. фиг.3 и 4). Квадраты на фиг.4 изображают "i"-ую и "i+1"-ую ячейки вдоль направлений зондирования под углами φj и φj+90°. Таким образом, определив дисперсию наклонов для двух взаимно перпендикулярных направлений, определяют полную дисперсию наклонов для "i"-ой ячейки σi 2 из соотношения: σ2 i=σ2 i(φj)+σ2 i(φj+90°).It is known that, with a sufficient degree of accuracy, the total variance of the slopes of the surface is the sum of the variances of the slopes measured in two mutually perpendicular directions. Therefore, to determine the total dispersion of the slopes σ i 2 of the surface in each "i" -th unit cell, the variance of the slopes in this cell is measured at different azimuthal angles φ j using the rotation of
Следующим шагом при обработке данных является определение направления распространения φwi крупномасштабного волнения в "i"-ой ячейке. Поскольку в каждой "i"-ой ячейке дисперсия наклонов σ2 i(φj) определена для нескольких (разных) азимутальных углов φj, то по экспериментально измеренным точкам с помощью алгоритма, используемого в скаттерометрии, восстанавливают азимутальную зависимость дисперсии наклонов (см. фиг.6). Максимальное значение дисперсии σw 2 на данном графике указывает искомое направление распространения крупномасштабного волнения φwi в "i"-ой ячейке.The next step in data processing is to determine the propagation direction φ wi of large-scale waves in the "i" -th cell. Since the slope dispersion σ 2 i (φ j ) in each “i” cell is determined for several (different) azimuthal angles φ j , the azimuthal dependence of the slope dispersion is restored using experimentally measured points using the algorithm used in scatterometry (see 6). The maximum dispersion value σ w 2 in this graph indicates the desired direction of propagation of large-scale waves φ wi in the "i" -th cell.
Скорость приповерхностного ветра V в "i"-ой ячейке определяют по измеренному с помощью регистрирующего устройства 4 в каждой ячейке сечению обратного рассеяния σ0 (θi) и по измеренным во взаимно перпендикулярных направлениях в каждой ячейке дисперсиям наклонов σ2 i(φj) и +σ2 i(φj+90°) (см. выше). Для вычисления искомой скорости ветра V используют новый регрессионный алгоритм V=F[σ0, σ2 i(φj), σ2 i(φj+90°)], который должен быть получен на начальном этапе калибровки радиолокатора спутника по стандартной методике (см., например, Witter and Chelton, A Geosat altimeter wind speed algorithm and a method for altimeter wind speed algorithm development, J. Geophysical Research, v.96, NC5, 1995, pp.8853-8860) с учетом упомянутых измеренных дисперсий наклонов σ2 i(φj) и σ2 i(φj+90°). Дисперсии наклонов σ2 i(φj) и σ2 i,(φj+90°) зависят не только от скорости ветра V, но и от присутствия в "i"-ой ячейке волн зыби, не связанных с ветром в данной ячейке, но влияющих на связь σ0 со скоростью ветра V. Известные и используемые на сегодняшний день алгоритмы вычисления скорости приповерхностного ветра V учитывают только одну переменную - сечение обратного рассеяния σ0 и не учитывают дисперсии наклонов водной поверхности в "i"-ой ячейке во взаимно перпендикулярных направлениях, которые существенно влияют на корреляцию σO со скоростью ветра V и поэтому, если их не учитывать, постоянно вносят ошибку в измерения. Точность определения скорости ветра V по новому алгоритму с учетом упомянутых дисперсий наклонов будет выше, что следует из математической статистики. Методика получения регрессионного алгоритма также известна (см. например, Г.Корн и Т.Корн "Справочник по математике для научных работников и инженеров", Москва, издательство "Наука", 1977 г. стр.553-557).The surface wind speed V in the "i" -th cell is determined by the backscattering cross section σ 0 (θ i ) measured using a recording device 4 in each cell and the slope dispersions σ 2 i (φ j ) measured in mutually perpendicular directions in each cell and + σ 2 i (φ j + 90 °) (see above). To calculate the desired wind speed V use the new regression algorithm V = F [σ 0 , σ 2 i (φ j ), σ 2 i (φ j + 90 °)], which should be obtained at the initial stage of satellite radar calibration according to the standard method (see, for example, Witter and Chelton, A Geosat altimeter wind speed algorithm and a method for altimeter wind speed algorithm development, J. Geophysical Research, v.96, NC5, 1995, pp.8853-8860) taking into account the mentioned measured variances slopes σ 2 i (φ j ) and σ 2 i (φ j + 90 °). The variances of the slopes σ 2 i (φ j ) and σ 2 i , (φ j + 90 °) depend not only on the wind speed V, but also on the presence of swell waves in the “i” cell that are not related to the wind in this cell but affecting the relation between σ 0 and wind speed V. The currently known and currently used algorithms for calculating near-surface wind speed V take into account only one variable — the backscattering cross section σ 0 and do not take into account the variance of the slopes of the water surface in the “i” cell perpendicular directions, which substantially affect the σ O correlation with the speed ve Mr. V, and therefore, if they are not taken into account, constantly making the error in the measurement. The accuracy of determining the wind speed V according to the new algorithm, taking into account the above-mentioned slope dispersions, will be higher, which follows from mathematical statistics. The methodology for obtaining the regression algorithm is also known (see, for example, G. Korn and T. Korn, "A Handbook of Mathematics for Scientists and Engineers", Moscow, Nauka Publishing House, 1977, pp. 553-557).
Таким образом, разработанный панорамный радиолокационный способ определения состояния приповерхностного слоя океана обеспечивает измерение в каждой "i"-ой элементарной ячейке широкой полосы обзора (шириной 355 км) следующих параметров: дисперсии наклонов σi 2, направления распространения φwi крупномасштабного волнения, а также более точное восстановление скорости приповерхностного ветра V, т.е. позволяет решить поставленную задачу и получить двумерное изображение рассеивающей водной поверхности в широкой полосе обзора.Thus, the developed panoramic radar method for determining the state of the near-surface layer of the ocean provides a measurement in each “i” elementary cell of a wide viewing band (355 km wide) of the following parameters: slope dispersion σ i 2 , propagation direction φ wi of large-scale waves, as well as more accurate restoration of the surface wind speed V, i.e. allows us to solve the problem and obtain a two-dimensional image of a scattering water surface in a wide viewing range.
В одном частном случае реализации разработанного способа, когда необходимо повысить пространственное разрешение, т.е. когда необходимо более детальное исследование состояния приповерхностного слоя океана и требуется рассмотреть более мелкие элементарные ячейки чем, 14×14 км, целесообразно с помощью дополнительной временной селекции по дальности сначала сформировать элементарные ячейки с размером, например, 14×1 км. Такая элементарная рассеивающая ячейка будет представлять собой сегмент (отрезок) кольца шириной 1 км и длиной 14 км, благодаря чему оказывается зафиксированным для данной ячейки (сегмента кольца) угол падения θ0, при этом разность доплеровских смещений Vдопл. внутри данной ячейки оказывается зависящей, как это следует из известной формулы для доплеровского смещения Vдопл=2V0sin θi·sin φj, только от изменений азимутального угла φj внутри данной ячейки 14×1 км:In one particular case of the implementation of the developed method, when it is necessary to increase the spatial resolution, i.e. when a more detailed study of the state of the surface layer of the ocean is required and it is necessary to consider smaller unit cells than 14 × 14 km, it is advisable to first form elementary cells with a size of, for example, 14 × 1 km using additional time selection in range. Such an elementary scattering cell will be a segment (segment) of a
Vдопл=2V0sinθ0[sin(φj+1°)-sin(φj)], гдеV add = 2V 0 sinθ 0 [sin (φ j + 1 °) -sin (φ j )], where
V0 - скорость перемещения спутника, оснащенного устройством для реализации разработанного способа (ширина диаграммы направленности антенны 1 в азимутальном направлении выбрана 1°). Как установлено авторами, сильное изменение разности доплеровских смещений Vдопл внутри ячейки 14×1 км происходит даже при изменении азимутального угла φj всего на 0,5-1° в следующих интервалах углов φj при каждом обороте антенны 1: в интервале φj от -45° до 45° и в интервале от 135° до 225°. Поэтому целесообразно использовать эту дополнительную доплеровскую селекцию по азимутальному углу φj для деления ячейки 14×1 км на 2-3 части по азимуту, т.е. для формирования элементарных ячеек с размерами, например, 5×1 км.V 0 - the speed of movement of a satellite equipped with a device for implementing the developed method (the beam width of the
В другом частном случае реализации разработанного способа, когда необходимо измерить среднюю длину волны Lm крупномасштабного волнения в каждой элементарной "i"-ой ячейке, размерами, например, 14×14 км, вначале по доплеровским скоростям с помощью дополнительных частотных фильтров, установленных в радиолокаторе 3, создают начальную точку отсчета для регистрации отраженных импульсов в каждой "i"-ой ячейке. После этого прием отраженных от водной поверхности импульсов и регистрацию их формы с помощью упомянутой временной селекции проводят для каждой "i"-ой ячейки отдельно и стандартным алгоритмом по тангенсу угла наклона в средней точке переднего фронта импульса, зарегистрированного в "i"-ой ячейке, определяют высоту H значительного волнения для нескольких азимутальных направлений в "i"-ой ячейке. Затем определяют уточненное значение высоты Н значительного волнения в "i"-ой ячейке путем определения среднего значения высоты Н по нескольким упомянутым измерениям в данной ячейке. После чего по упомянутому измеренному максимальному значению дисперсии наклонов σw 2 в "i"-ой ячейке на графике фиг.6 и по уточненному значению высоты Н значительного волнения определяют среднюю длину волны Lm крупномасштабного волнения из соотношения Lm=H/σw. Высота Н и длина волны Lm крупномасштабного волнения имеют важное значение для численных моделей приповерхностного слоя океана, кроме того, от длины волны Lm напрямую зависит скорость распространения крупномасштабного волнения, т.е., например, насколько быстро зародившийся шторм достигнет прибрежной зоны.In another particular case of the implementation of the developed method, when it is necessary to measure the average wavelength L m of large-scale waves in each unit "i" -th cell, with dimensions, for example, 14 × 14 km, first by Doppler speeds using additional frequency filters installed in the
Таким образом, разработанный панорамный радиолокационный способ определения параметров состояния приповерхностного слоя океана позволяет построить двумерное изображение водной поверхности, что открывает возможности для анализа волновых процессов на поверхности океана, изучения их структуры и временной динамики при повторных наблюдениях, что необходимо для составления долгосрочных прогнозов погоды, для изучения Мирового океана и решения множества других прикладных задач.Thus, the developed panoramic radar method for determining the state parameters of the surface layer of the ocean allows us to construct a two-dimensional image of the water surface, which opens up the possibility of analyzing wave processes on the surface of the ocean, studying their structure and temporal dynamics during repeated observations, which is necessary for making long-term weather forecasts, for studying the oceans and solving many other applied problems.
Claims (3)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2003124710/09A RU2274877C2 (en) | 2003-08-07 | 2003-08-07 | Panoramic radar method of determining condition of ocean's layer surface from satellite |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2003124710/09A RU2274877C2 (en) | 2003-08-07 | 2003-08-07 | Panoramic radar method of determining condition of ocean's layer surface from satellite |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| RU2003124710A RU2003124710A (en) | 2005-02-10 |
| RU2274877C2 true RU2274877C2 (en) | 2006-04-20 |
Family
ID=35208502
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2003124710/09A RU2274877C2 (en) | 2003-08-07 | 2003-08-07 | Panoramic radar method of determining condition of ocean's layer surface from satellite |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2274877C2 (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2404434C2 (en) * | 2008-07-11 | 2010-11-20 | Морской гидрофизический институт Национальной академии наук Украины (МГИ НАН Украины) | Method for remote determination of wind speed over water surface |
| RU2449312C1 (en) * | 2010-12-27 | 2012-04-27 | Учреждение Российской академии наук Институт прикладной физики РАН | Panoramic radar method of determining parameters of state of ocean surface layer from satellite |
| RU2597148C1 (en) * | 2015-03-20 | 2016-09-10 | федеральное государственное бюджетное образовательное учреждение высшего образования "Российский государственный гидрометеорологический университет" | Method of measuring vector field speed of ocean and river streams in space sar |
| RU2630412C1 (en) * | 2016-11-29 | 2017-09-07 | Федеральное государственное бюджетное научное учреждение "Федеральный исследовательский центр Институт прикладной физики Российской академии наук" (ИПФ РАН) | Method of remote determination of sea current speed |
| RU2794213C2 (en) * | 2018-07-30 | 2023-04-12 | Федеральное государственное бюджетное учреждение науки Институт океанологии им.П.П.Ширшова Российской академии наук | Method for operational search and detection of underwater sound source |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0455842A1 (en) * | 1990-05-07 | 1991-11-13 | Dornier Gmbh | Scatterometer |
| US5234852A (en) * | 1990-10-10 | 1993-08-10 | Sgs-Thomson Microelectronics, Inc. | Sloped spacer for MOS field effect devices comprising reflowable glass layer |
| RU2047874C1 (en) * | 1993-08-18 | 1995-11-10 | Черный Игорь Владимирович | Method of remote diagnosis of condition of system "ocean-atmosphere" and device for its realization |
| US6137437A (en) * | 1999-03-24 | 2000-10-24 | Agence Spatiale Europeenne | Spaceborne scatterometer |
| RU2235344C2 (en) * | 2002-06-03 | 2004-08-27 | Институт прикладной физики РАН | Radiolocation method for determining parameters of condition of surface-adjacent ocean layer from a satellite |
-
2003
- 2003-08-07 RU RU2003124710/09A patent/RU2274877C2/en not_active IP Right Cessation
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0455842A1 (en) * | 1990-05-07 | 1991-11-13 | Dornier Gmbh | Scatterometer |
| US5234852A (en) * | 1990-10-10 | 1993-08-10 | Sgs-Thomson Microelectronics, Inc. | Sloped spacer for MOS field effect devices comprising reflowable glass layer |
| RU2047874C1 (en) * | 1993-08-18 | 1995-11-10 | Черный Игорь Владимирович | Method of remote diagnosis of condition of system "ocean-atmosphere" and device for its realization |
| US6137437A (en) * | 1999-03-24 | 2000-10-24 | Agence Spatiale Europeenne | Spaceborne scatterometer |
| RU2235344C2 (en) * | 2002-06-03 | 2004-08-27 | Институт прикладной физики РАН | Radiolocation method for determining parameters of condition of surface-adjacent ocean layer from a satellite |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2404434C2 (en) * | 2008-07-11 | 2010-11-20 | Морской гидрофизический институт Национальной академии наук Украины (МГИ НАН Украины) | Method for remote determination of wind speed over water surface |
| RU2449312C1 (en) * | 2010-12-27 | 2012-04-27 | Учреждение Российской академии наук Институт прикладной физики РАН | Panoramic radar method of determining parameters of state of ocean surface layer from satellite |
| RU2597148C1 (en) * | 2015-03-20 | 2016-09-10 | федеральное государственное бюджетное образовательное учреждение высшего образования "Российский государственный гидрометеорологический университет" | Method of measuring vector field speed of ocean and river streams in space sar |
| RU2630412C1 (en) * | 2016-11-29 | 2017-09-07 | Федеральное государственное бюджетное научное учреждение "Федеральный исследовательский центр Институт прикладной физики Российской академии наук" (ИПФ РАН) | Method of remote determination of sea current speed |
| RU2794213C2 (en) * | 2018-07-30 | 2023-04-12 | Федеральное государственное бюджетное учреждение науки Институт океанологии им.П.П.Ширшова Российской академии наук | Method for operational search and detection of underwater sound source |
Also Published As
| Publication number | Publication date |
|---|---|
| RU2003124710A (en) | 2005-02-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hauser et al. | Satellite remote sensing of surface winds, waves, and currents: Where are we now? | |
| Hajj et al. | Theoretical description of a bistatic system for ocean altimetry using the GPS signal | |
| CN110824510B (en) | Method for increasing number of sea surface reflection signals received by GNSS-R height measurement satellite | |
| Hauser et al. | SWIMSAT: A real-aperture radar to measure directional spectra of ocean waves from space—Main characteristics and performance simulation | |
| Li et al. | An algorithm for sea-surface wind field retrieval from GNSS-R delay-Doppler map | |
| CN110609287B (en) | Double-frequency radar scatterometer and method for simultaneously measuring sea surface wind field and flow field | |
| CN113050134A (en) | Sea surface wind field inversion observation method based on satellite navigation information | |
| Pozderac et al. | $ X $-Band Beacon-Receiver Array Evaporation Duct Height Estimation | |
| Martín-Neira et al. | Altimetry precision of 1 cm over a pond using the wide-lane carrier phase of GPS reflected signals | |
| Zhang et al. | Wind direction retrieval using spaceborne GNSS-R in nonspecular geometry | |
| Hauser et al. | CFOSAT: A new Chinese-French satellite for joint observations of ocean wind vector and directional spectra of ocean waves | |
| Caudal et al. | KuROS: A new airborne Ku-band Doppler radar for observation of surfaces | |
| RU2466425C1 (en) | Method of measuring characteristics of wavy water surface | |
| Motte et al. | Applications of GNSS-R in continental hydrology | |
| RU2449312C1 (en) | Panoramic radar method of determining parameters of state of ocean surface layer from satellite | |
| RU2274877C2 (en) | Panoramic radar method of determining condition of ocean's layer surface from satellite | |
| RU2235344C2 (en) | Radiolocation method for determining parameters of condition of surface-adjacent ocean layer from a satellite | |
| Wang et al. | Validation and intercomparison of sea state parameter estimation with multisensors for OSMAR-S high-frequency radar | |
| Wang et al. | A passive technique to monitor evaporation duct height using coastal GNSS-R | |
| Rius et al. | Reflectometry | |
| Bernhardt et al. | Bistatic observations of the ocean surface with HF radar, satellite and airborne receivers | |
| Shibo et al. | Inversion method of significant wave height based on time delay window using the ship-borne Beidou reflected signal | |
| Ermoshkin et al. | Determination of the Significant Wave Height from Doppler Radar Images of the Sea Surface | |
| Guo et al. | A new model to estimate significant wave heights with ERS-1/2 scatterometer data | |
| Zhou et al. | Research on the parameter design method and system simulation of multimode microwave remote sensors operating in scatterometer modes |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| MM4A | The patent is invalid due to non-payment of fees |
Effective date: 20110808 |