Измерение влажности и высоты посевов сельскохозяйственных полей с использованием приемника сигналов ГЛОНАСС и GPS

Полный текст

Введение. В настоящее время активно развивается метод дистанционного измерения физических характеристик покровов Земли с использованием сигналов глобальных навигационных спутниковых систем (ГНСС) [1-13]. Сигналы, излучаемые ГНСС в круглосуточном режиме, покрывают поверхность всей Земли и не зависят от погодных условий. Метод измерения основан на регистрации приемником сигналов ГНСС интерференционной картины, возникающей в результате суперпозиции поля прямой волны, распространяющейся от спутника ГНСС к приёмной антенне, расположенной на высоте нескольких метров над поверхностью почвы, с полем отраженной волны от поверхности почвы. Характер интерференционной картины, сформированной многолучевым распространением сигналов ГНСС, зависит от электро-, геофизических характеристик почвенного и надпочвенного покрова (влажность почвы, влажность и высота растительности). В результате эти параметры могут быть восстановлены из измеренной приемником ГНСС интерференционной картины с использованием теоретической модели, описывающей интерференцию амплитуд сигналов в приемнике ГНСС. В работах [1-3; 5; 6] исследовалась принципиальная возможность измерения влажности поверхностного слоя почвенного покрова с использованием сигналов ГНСС. Измерение влажности осуществлялось как с использованием промышленно выпускаемых образцов приемников [1; 4; 6-10], так и с опытными образцами [2; 3]. При этом приемники оснащаются как антеннами с правой круговой [1; 4-10], так и вертикальной поляризацией [2; 3]. В самое последнее время [1] была экспериментально доказана возможность применения для измерения влажности почвенного покрова стандартных приемников сигналов ГНСС, составляющих сеть Earthscope Plate Boundary Observatory, которая создана на территории США для мониторинга смещений тектонических плит в зонах разломов. Предложены методы и проведены эксперименты [12; 13] по одновременному определению влажности почвенного и растительного надпочвенного покрова. Особенностью подхода стало использование специально разработанного приемника и антенн с вертикальной поляризацией. В России впервые исследования по измерению влажности почвенного покрова с применением сигналов ГНСС были начаты в 2011 г. [4]. В работе [10] проведены измерения интерференционных ГНСС сигналов с использованием стандартного приемника МРК-32 производства ОАО « НІ II I «Радиосвязь», г. Красноярск, оснащенного антенной с правой круговой поляризацией; получены оценки влажности с применением алгоритма, предложенного в [2]. В результате этого исследования установлено, что применение стандартной приемной антенны с правой круговой поляризацией позволяет восстанавливать влажность с погрешностью больше чем 0,04 см3/см3, которая принята предельно допустимой для спутников дистанционного зондирования Земли. В дальнейшем в работе [5] был создан оригинальный алгоритм обработки интерференционного сигнала, позволивший при использовании антенны с линейной вертикальной поляризацией восстановить влажность почвы с погрешность меньше чем 0,04 см3/см3. При этом в модели интерференционного сигнала была использована модель диэлектрической проницаемости почвы, обладающая существенно меньшей погрешностью по сравнению с моделями, примененными в [2; 3]. Прием на антенну с вертикальной поляризацией подтвердил свою эффективность при тестировании на большем объеме экспериментальных данных [6]. В данной работе проводится анализ экспериментальных данных [7] и исследуются потенциальные возможности сигналов ГНСС, принятых аппаратурой МРК-32Р, оснащенной антенной с правой круговой поляризацией, для измерения влажности почвы, влажности и высоты растительного покрова сельскохозяйственных полей. Автономный измерительно-регистрирующий комплекс на базе серийного приемника сигналов ГЛОНАСС/GPS МРК-32Р и методика измерений. Для проведения полевых экспериментов по дистанционному зондированию земных покровов с помощью сигналов ГНСС GPS и ГЛОНАСС был собран автономный измерительно-регистрирующий комплекс (АИРК) на базе серийного комплекта МРК-32Р производства НПО «Радиосвязь», г. Красноярск. Блок-схема измерительной установки АИРК представлена на рис. 1. Сигналы ГНСС принимаются с помощью антенн 1, 2, 3, 4 в различных конфигурациях (рис. 1). В данном эксперименте измерения проводились на антенну 1, представляющую собой антенный модуль АМ-7 комплекта МРК-32Р (рис. 2). Антенна размещалась на переносной металлической треноге с удлинительной штангой (рис. 2). Максимальная высота размещения фазового центра антенны в полевых условиях может достигать 3,3 м. Сигналы навигационных спутников, принятые антенной, поступают через смеситель в блок МРК-11 ПрМ-16. Здесь сигналы обрабатываются и через конвертер RS-232/USB2.0 89 Вестник СибГАУ. 2014. N 4(56) направляются в ноутбук. С помощью специального программного обеспечения обеспечивается визуализация процесса приема и обработки сигналов, а также их запись в бинарный файл. Методика проведения эксперимента включала в себя две стадии: полевые измерения и первичная камеральная обработка данных. При полевых работах на месте проведения сеанса наблюдения производился монтаж антенной системы, выставлялась необходимая высота антенны и ее ориентация по азимуту и углу места. После подключения всех разъемов осуществляется проверка работоспособности аппаратуры, прогрев и настройка в течение 10-15 мин. Затем осуществляется включение на запись принимаемого сигнала на жесткий диск ноутбука. Первичная камеральная обработка заключается в извлечении и преобразовании полученных данных из бинарного файла в серию файлов в формате csv (обработка программой Excel). Для удобства дальнейшей обработки все записи сигналов компонуются стандартным образом в один проект в среде программы OriginLab, куда также добавляется информация о существенных обстоятельствах проведенного сеанса (место, время, характеристики антенной конфигурации, объекта и т. п.). Полученные в таком виде экспериментальные данные проходили дальнейшую обработку и анализ. проницаемости будем использовать среднее значение глинистой фракции Wc = 0,312 г/г. Рис. 2. Антенна правой круговой поляризации приемника ГНСС, снабженная поворотным устройством Таблица 1 Содержание глинистой фракции почвы участков полей № Наименование поля Массовое содержание глинистой фракции, г/г 1 Рожь 0,312 2 Пшеница 0,316 3 Рапс 0,309 Рис. 1. Блок-схема приемника сигналов ГНСС В качестве объектов исследования были выбраны поля ФГУП «Минино» Российской академии сельскохозяйственных наук, засеянные рапсом, пшеницей, ячменем и рожью. Во время эксперимента на каждом из полей были проведены измерения высоты растительности, среднее количество стеблей на квадратный метр, отобраны образцы растительности и почвы (с глубины 5 и 10 см) для определения влажности термостатно-весовым методом; проведен гранулометрический анализ образцов почв на содержание глинистой фракции. Результаты гранулометрического анализа почвенных образцов представлены в табл. 1. Как следует из табл. 1, содержание глинистой фракции в почве практически не изменяется для всех тестовых участков. В дальнейших расчетах диэлектрической Поля с различными видами культур находились на достаточно ровной поверхности и имели форму прямоугольника, стороны которого ориентированы по направлениям света. Антенна приемника располагалась на северной границе выбранного поля. Максимум диаграммы направленности антенны ориентировался по горизонту в южном направлении. В результате такой ориентации антенны большая часть зон, существенных для отражения сигналов спутников GPS и ГЛОНАСС, проходили по исследуемому полю. Характеристики исследуемых зерновых культур и условия проведенных сеансов измерений представлены в табл. 2. Диэлектрическая модель растительного покрова и почвы. Модель комплексной диэлектрической проницаемости посевов была построена на основе измерения комплексной диэлектрической проницаемости (КДП) образцов растительности с поля, засеянного рожью. Измерения проводились с использованием методики, разработанной ранее [14], на частоте 1,51 ГГц и в диапазоне влажностей от 0,11 до 0,52 г/г, при комнатной температуре 25 °С. Для измерений был использован векторный анализатор цепей ZVK, изготовленный фирмой Rohde-Schwarz. Исследуемый образец помещается в измерительный контейнер, который представляет собой отрезок коаксиального волновода, подключенный к векторному анализатору цепей. Срезанные на поле колосья ржи вместе со стеблями в лабораторных условиях были высушены 90 Математика, механика, информатика и размолоты до порошкообразного состояния. Далее образец увлажнялся дистиллированной водой и выдерживался в герметичном состоянии 5 суток. По истечении 5 суток проводились диэлектрические измерения, после которых определялась влажность каждого образца. Образец взвешивался и помещался в тепловой шкаф, в котором при температуре 105 °С он выдерживался в течение 24 ч, далее производилось очередное взвешивание. Таким образом, определялась влажность термостатно-весовым способом. На рис. 3 показана зависимость приведенных коэффициентов преломления (ns-1) / pd и поглощения к / pd образца растительности, где pd - плотность сухого образца от его весовой влажности mg. Зависимости, показанные на рис. 3, можно описать рефракционной моделью смеси сухого органического вещества растений и воды в растениях [14]: n = 1+И-аЛ Pd + n - 1)W, к =-Pd +KW, P m (1) где nm, nw, Km, kw - показатели преломления и поглощения соответственно; индексы m и w относятся к органическому веществу сухих растений и воде в растениях соответственно; pd - плотность органического вещества сухих растений, которая измерялась в эксперименте. С помощью регрессионного анализа данных, представленных на рис. 3, найдены следующие параметры модели (1): 1 = 0,13, (nw -1) = 7,69, Km -Pm P, = 1,065, Kw = 0. (2) Комплексный показатель преломления для слоя влажной растительности в полевых условиях оценим исходя из рефракционной формулы смеси (1), учитывая долю воды, содержащейся в объеме биомассы. с Ï. >5 3 т X £ а> ш з а. Р T=25rtÜ 1=1.5-1 ГГц ~т~ Û.1 -і- С 2 I 03 I =1= 05 “Г“ о,з ВыСиЫЫГН ШіЗЖІ-ЮСТЬ, U, I 0£ O.L! L.-1 BtKDSaa ВЛЭИНЭСТЬ. І-'З б Рис. 3. Зависимость приведенного коэффициента преломления образца (n5-1)/pd (a) и приведенного коэффициента затухания ks / Pd (б) от весовой влажности mg образца растительности m Таблица 2 Характеристики зерновых культур и условия сеанса измерения Условия сеансов Ед. изм. Рожь озимая Пшеница яровая Ячмень яровой Рапс яровой Дата проведения 1 сеанса - 20.08.2013 26.08.2013 27.08.2013 28.08.2013 Широта точки измерения с. ш. 56°04'43".2 56°05'02".8 56°04'38".1 56°04'54".4 Долгота точки измерения в. д. 92°40'26".5 92°40'35".5 92°40'42".5 92°40'28".6 Высота фазового центра антенны м 3,08 3,08 3,07 3,06 Ориентация оси антенны град 197 174 185 175 Средняя высота растительного слоя м 1,03 0,87 0,77 0,66 Средняя длина растений м 1,18 0,74 0,81 0,59 Плотность посадки растений шт/м2 380 370 340 150 Содержание сухой биомассы в слое г/см3-10-3 1,01 0,66 0,66 0,50 91 Вестник СибГАУ. 2014. № 4(56) Окончание табл. 2 Условия сеансов Ед. изм. Рожь озимая Пшеница яровая Ячмень яровой Рапс яровой Объемная влажность растительного слоя см3/см3-10-3 0,59 0,79 1,49 2,31 Объемная влажность почвы на глубине 5 см см3/см3 0,27 0,23 0,29 0,19 Объемная влажность почвы на глубине 10 см см3/см3 0,30 0,33 0,31 0,27 Тогда эффективные показатели преломления nl и поглощения к растительного покрова могут быть записаны в следующем виде: п, = 1 + 2 ^ Pl + n - 1)W, Pm (3) Kt = 2 ^ p, +kwW, , Pm где p, - эффективная плотность слоя, заполненного сухими растениями. Эффективная плотность слоя была оценена как содержание сухой биомассы в кубическом метре посевов (см. табл. 2). Далее модель (3) будет использоваться при теоретическом расчете ин-терферограмм сигналов ГНСС. КДП почв сельскохозяйственных земель как функция частоты электромагнитного поля f, объемной влажности Ws и весового содержания глины Wc в почве рассчитывалась на основе рефракционной диэлектрической модели смеси [15]. Частота электромагнитного поля задается передатчиком ГНСС, объемная влажность Ws почвы задавалась значениями из табл. 2. Теоретическая модель интерферограмм сигналов ГЛОНАСС и GPS. Геометрия задачи изображена на рис. 4. Антенна приемника ГНСС расположена на высоте h над поверхностью земли, покрытой растительным слоем толщиной d. Максимум диаграммы направленности антенны ориентирован на горизонт. Излучение ГНСС имеет правую круговую поляризацию. Амплитуда сигнала U(9), регистрируемая приемником ГНСС, оснащенным антенной с правой круговой поляризацией, может быть описана следующей формулой [1]: U(0) = F(0) • |1 + Г(0)e2ik(A cos012, (4) где F(9) - амплитуда на выходе приемника ГНСС в отсутствие волны, рассеянной растительным и почвенным покровом; Г(9) - комплексный коэффициент отражения от почвенного и растительного покрова волны, имеющей правую круговую поляризацию; 9 -угол падения волны на поверхность почвенного и растительного покрова; he - эффективная высота фазового центра антенны над растительным покровом; k0 = а / c - волновое число свободного пространства; ю = 2nf где f- частота электромагнитного поля, излучаемая спутником ГНСС. Будем пренебрегать явлениями дифракции и рассеяния электромагнитных волн на элементах растительного покрова, представленного зерновыми культурами. В этом случае будем рассматривать растительный покров в виде однородного диэлектрического слоя, расположенного на почве. Обозначим через es и єl относительные комплексные диэлектрические проницаемости почвы и растительности соответственно. Почва Рис. 4. Многолучевое распространение волн от спутников ГНСС при расположении приемной антенны (ориентированной на горизонт) приемника ГНСС на высоте h над почвой, покрытой слоем растительности толщиной d: 9 - угол падения, а - угол возвышения, 9 = 90° - а Комплексный коэффициент отражения от слоистой среды «растительность-почва» для плоских волн с горизонтальной Гн(9) и вертикальной ГЦ9) поляризациями можно рассчитать по формуле [2] Г r, h v (0) + RsH V (0)eSem°dewl H ,V ()_1 + R, h v (0)Rs,h v (q)eSe^wl , () где Rl,H,V(9) и RsHiV(9) - коэффициенты отражения Френеля от границы «воздух-растительность» и «растительность-почва»; S = -2[k0acos9]2; ст - среднеквадратичное отклонение высот неровностей поверхности почвы; w, =<Jг, - sin2(0) - волновое число в слое растительности; de - высота слоя растительности. Коэффициенты отражения Френеля RlHV(9), Rs, H>V(9) рассчитываются по формулам г, cos 0 -Jг, - sin2 0 R,v (0) = --V . 2 , г, cos 0 + f г, - sin 0 92 Математика, механика, информатика R V (0.) = £ s cos 8, ~УІ^/ £ s - £, sin2 8, Є s c0S 81 ^V£^V£ s - £l sin2 81 (6) R H (8) =■ RH (81) = cos 8 £, cos 8 V£ і - sin2 8 ■\/Є l - sin2 8 Л/Є s - є, sin2 8 V£ s - є, sin2 8 (7) где 8j = a sin (n 8 ^ ^/£7). Окончательно коэффициент отражения Г(9) в формуле (4) рассчитывается с использованием выражений (5)-(7) по формуле ГR(8)= Гу (8) + Г H (8) 2 ' (8) Угол падения. Частота электромагнитного поля, излучаемого спутниками ГЛОНАСС, рассчитывается по формуле f = 1 602 МГц + n • 0,5625 МГц, где n - это номер частотного канала (n = -7, ..., 6), cпутнники GPS излучают электромагнитные волны на фиксированной частоте: f = 1575,42 МГц. Чувствительность интерферограмм сигналов ГНСС к вариации влажности и высоте растительного слоя. Проведем расчет интерферограмм сигналов ГНСС на основе формул (4)-(8). Амплитуду на выходе приемника ГНСС в отсутствие волны, рассеянной растительным и почвенным покровом, зададим в виде полинома четвертой степени F(9) = -0,0321 + + 0,01499 + 0,000192 - 0,0106-10-493 - 0,0108-10-694 (полином получен в результате регрессионного анализа экспериментальной интерференционной картины одного из сигнала, принятого от спутника GPS, нормированного на средний уровень сигнала при угле падения 80о). При расчете КДП почвы содержание глины и объемная влажность задавалась равными 0,312 г/г и 0,25 см3/см3 соответственно. При расчете КДП растительного слоя плотность сухой биомассы задавалась равной 1,01 • 10-3 г/см3. Высота приемника была задана 3,05 м, а шероховатость поверхности почвы 0 см. На рис. 5 приведены рассчитанные ин-терферограммы сигналов ГНСС в зависимости от влажности растительного слоя в сравнении с интер-ферограммами сигналов ГНСС в отсутствии растительного слоя (влажность слоя равна нулю). Из результатов расчетов, приведенных на рис. 5, видно, что угловые зависимости амплитуд сигналов ГНСС существенно изменяются по величине при увеличении влажности растительного слоя по сравнению с сигналом в отсутствие растительного слоя на поверхности почвы. При фиксированной влажности слоя растительности и увеличении его высоты от 0,1 до 0,9 м (рис. 6) угловые зависимости амплитуд сигналов ГНСС испытывают значительные фазовые вариации. Фазовая вариация интерференционной картины в большей степени чувствительна к изменению высоты растительного слоя, амплитудные вариации - к влажности. Проведенный анализ показывает принципиальную возможность восстановления влажности и высоты растительного слоя на основе измерения угловых амплитудных зависимостей сигналов ГНСС. 40 50 60 70 Угол падения, sa б Рис. 5. Интерферограммы сигналов ГНСС в зависимости от влажности растительного слоя. Высота растительного слоя задавалась равной 0,9 м Определение влажности и высоты растительного покрова с использованием измеренных интерферограмм сигналов ГЛОНАСС и GPS. Метод восстановления, влажности и высоты растительности сельскохозяйственных полей из интерференционной картины сигналов ГНСС основан на минимизации в диапазоне углов зондирования 50о< 9 < 80о функционала, образующего норму невязки между измеренной амплитудой сигнала ГНСС Um(9) и модельным значением U(9) (см. формулу (4)). В общем виде минимизируемый функционал является функцией объемной влажности Ws и весового содержания глины Wc в почве, влажности W, и высоты de растительного покрова, высоты антенны приемника ГНСС he и среднеквадратичной высоты неровностей поверхности почвы ст. Проведенные исследования показали, что коэффициент отражения волны с правой круговой поляризацией при падающей волне от спутника ГНСС правой круговой поляризации слабо меняется (менее 3 %) при вариациях весового содержания глины (0,20-0,49 г/г) и объемного содержания влаги (0,18-0,40 см3/см3) в почве. В связи с этим при минимизации функционала объемное содержание влаги задавалось средним є, cos 93 Вестник СибГАУ. 2014. N 4(56) значением между сухой почвой и влажностью при полной влагоемкости Ws = 0,23 см3/см3, весовое содержание глины задавалось средним значением по содержанию глинистой фракции на исследуемых полях Wc = 0,312 г/г. Среднеквадратичная высота неровностей поверхности почвы задавалась равной ст = 0,0 м. а 40 50 00 Угол падения, 70 ео б Рис. 6. Интерферограммы сигналов ГНСС в зависимости от высоты растительного слоя. Влажность растительного слоя задавалась 1,5 г/см3 Метод восстановления влажности и высоты растительности из интерференционной картины сигналов ГНСС поэтапно реализовывался по следующему алгоритму: 1. В формуле (4) средний уровень амплитуды сигналов на выходе приемника ГНСС U(0) может быть аппроксимирован полиномом четвертой степени U(0) = a0 + aq + a202 + a303 + a404. Воспользуемся этой функцией в регрессионном анализе экспериментального интерференционного сигнала и найдем коэф фициенты ai1 полинома U(0). Полученную экспериментальную функцию U(0) будем использовать в дальнейшем в качестве первого приближения для теоретической модели амплитуды F (0) = U(0) на выходе приемника ГНСС в отсутствие отраженной волны растительным и почвенным покровом. 2. На втором этапе путем минимизации функционала нормы невязки находим первые приближения для высоты антенны приемника ГНСС he,1, высоты de,1 и объемной влажности Wl1 растительного слоя. 3. На заключительном этапе проводится минимизация функционала и определяются следующие восстанавливаемые параметры: aiF, heF, deF, WlF. При этом в качестве начальных величин для этих параметров используются значения, определенные на предыдущих этапах: a,;1, he1, de1, Wl,1. Данная методика была реализована на языке программирования С++ в рабочей среде программы OriginLab, в которой производилась минимизация функционала на основе алгоритма Левенберга-Марк-вардта [16]. Вышеописанный алгоритм был применен для восстановления влажности и высоты посевов следующих культур: рапс, рожь, пшеница и ячмень. На рис. 7 приведены результаты экспериментальной и теоретически рассчитанной интерференционной картины. В табл. 3 сведены восстановленные параметры модели и экспериментально измеренные контактными методами соответствующие величины (обозначены с индексом m). Содержание воды растительного слоя на один квадратный метр Ml рассчитывалось по формуле Ml = = Wl-pw-de, pw = 1 000 кг/м3 - плотность воды. Высота фазового центра антенны hem над растительным покровом рассчитывалась как разница между высотой антенны над почвой h и высотой растительного покрова dm, измеренными на поле во время эксперимента. В экспериментах по измерению влажности и высоты растительности с помощью приемника МРК-32Р показано, что в среднем погрешность измерения высоты и влажности растительного покрова сельскохозяйственных полей с четырьмя различными видами культур составляет величину Ade = 0,17 м и AMl = = 0,21 кг/м2 соответственно (табл. 3). Предложенная теоретическая модель (4)-(8), описывающая амплитуду сигнала, регистрируемого приемником ГНСС в присутствии почвенного покрова, покрытого растительностью, в среднем с коэффициентом корреляции 0,71 коррелирует с измеренными интерференционными сигналами. Таблица 3 Восстановленные и измеренные параметры растительности Поле he,F, м de,F, м Wif10 3, м7м3 MiF, кг/м2 he,m, м dm, м Wl,m10-3, м3 /м3 Ml,m, кг/м2 Рапс 2,56 0,45 3,35 1,51 2,39 0,66 2,30 1,52 Рожь 2,26 0,87 1,10 0,96 2,05 1,03 0,60 0,62 Пшеница 2,23 0,87 1,19 1,04 2,21 0,87 0,79 0,69 Ячмень 2,54 0,49 2,1 1,03 2,30 0,77 1,51 1,16 94 Математика, механика, информатика Угол падения. 245 Угол падения, в 70 - Эксперимент -Теория 57 00 63 65 09 72 75 70 01 Угол падения, в 62 64 66 66 70 72 74 76 76 60 Угол падения, ° Рис. 7. Интерференционная картина, записанная приемником МРК-32 при регистрации сигналов GPS над сельскохозяйственным полем с посевами: а - рапс (спутник № 5, дата 28.08.2013), б - рожь (спутник № 1, дата 20.08.2013), в - пшеница (спутник № 10, дата 26.08.2013), г - ячмень (спутник № 2, дата 27.08.2013) г Заключение. Разработана теоретическая модель интерферограмм сигналов ГЛОНАСС и GPS, в которую интегрирована модель комплексной диэлектрической проницаемости почв с параметрами, характерными для почв тестовых участков сельскохозяйственных земель с посевами пшеницы, ржи, рапса и ячменя в Емельяновском районе Красноярского края. Проведены полевые измерения интерферограмм сигналов ГЛОНАСС и GPS отечественным приемником МРК-32Р на тестовых участках. При сопоставлении теоретической модели интерферограмм с экспериментально измеренными получена высокая корреляция между ними. На основе разработанной теоретической модели интерферограмм и диэлектрической модели почв и растительности создан алгоритм определения влажности и высоты растительного покрова для тестовых участков пшеницы, рапса, ржи и ячменя в Емельяновском районе Красноярского края. Восстановленные значения влажности и высоты растительного покрова тестовых участков в среднем с погреш ностью 0,17 м и 0,21 кг/м2 соответствуют экспериментально измеренным значениям. Для практического применения данных алгоритмов при восстановлении влажности почвы, уточнения погрешностей и условий, ограничивающих применимость данных алгоритмов, необходимо проводить дальнейшую их валидацию на других участках сельскохозяйственных полей при вариации метеорологических условий, шероховатости почвенного покрова. Благодарности. Работа выполнена в рамках проекта РФФИ № 13-05-98089, р_сибирь_а, программы II.12.1. СО РАН «Радиофизические и акустические методы дистанционного зондирования природной среды». Acknowledgements. The work was supported by RFBR, project № 13-05-98089, r_sibir_a, programmy II. 12.1. SO RAN «Radiofizicheskie i akusticheskie metody distancionnogo zondirovaniya prirodnoj sredy». 95 Вестник СибГАУ. 2014. № 4(56) Библиографические ссылки

Об авторах

Валерий Леонидович Миронов

Институт физики им. Л. В. Киренского СО РАН

Email: rsdvm@ksc.krasn.ru. rsdkm@ksc.krasn.ru
доктор физико-математических наук, член-корреспондент РАН, заведующий лабораторией радиофизики дистанционного зондирования, Институт физики СО РАН; заведующий кафедрой КМТ ИКИВТ

Михаил Иванович Михайлов

Институт физики им. Л. В. Киренского СО РАН

Email: rsdvm@ksc.krasn.ru. rsdkm@ksc.krasn.ru
ведущий инженер лаборатории радиофизики дистанционного зондирования

Константин Викторович Музалевский

Институт физики им. Л. В. Киренского СО РАН

Email: rsdvm@ksc.krasn.ru. rsdkm@ksc.krasn.ru
кандидат физико-математических наук, научный сотрудник лаборатории радиофизики дистанционного зондирования, Институт физики СО РАН; старший преподаватель кафедры КМТ ИКИВТ

Анатолий Васильевич Сорокин

Институт физики им. Л. В. Киренского СО РАН

Email: rsdvm@ksc.krasn.ru. rsdkm@ksc.krasn.ru
кандидат физико-математических наук, заместитель заведующего КМТ ИКИВТ, Сибирский государственный аэрокосмический университет имени академика М. Ф. Решетнева; ведущий инженер лаборатории радиофизики дистанционного зондирования

Сергей Викторович Фомин

Институт физики им. Л. В. Киренского СО РАН

Email: rsdvm@ksc.krasn.ru. rsdkm@ksc.krasn.ru
младший научный сотрудник лаборатории радиофизики дистанционного зондирования

Андрей Юрьевич Каравайский

Институт физики им. Л. В. Киренского СО РАН

Email: rsdvm@ksc.krasn.ru. rsdkm@ksc.krasn.ru
ведущий инженер лаборатории радиофизики дистанционного зондирования

Список литературы

Дополнительные файлы