close

Вход

Забыли?

вход по аккаунту

?

Патент BY15950

код для вставкиСкачать
ОПИСАНИЕ
ИЗОБРЕТЕНИЯ
К ПАТЕНТУ
РЕСПУБЛИКА БЕЛАРУСЬ
(46) 2012.06.30
(12)
(51) МПК
НАЦИОНАЛЬНЫЙ ЦЕНТР
ИНТЕЛЛЕКТУАЛЬНОЙ
СОБСТВЕННОСТИ
(54)
(2006.01)
G 01J 1/10
G 01D 18/00 (2006.01)
СПОСОБ ОПРЕДЕЛЕНИЯ КАЛИБРОВОЧНОГО КОЭФФИЦИЕНТА
ДЛЯ КАЛИБРОВКИ ОПТИЧЕСКОГО СПУТНИКОВОГО СЕНСОРА
ВЫСОКОГО ПРОСТРАНСТВЕННОГО РАЗРЕШЕНИЯ В
АБСОЛЮТНЫХ ЭНЕРГЕТИЧЕСКИХ ЕДИНИЦАХ
(21) Номер заявки: a 20100207
(22) 2010.02.12
(43) 2011.10.30
(71) Заявитель: Государственное научное учреждение "Институт физики
имени Б.И.Степанова Национальной академии наук Беларуси" (BY)
(72) Авторы: Кацев Иосиф Лейбович;
Прихач Александр Сергеевич; Зеге
Элеонора Петровна (BY)
BY 15950 C1 2012.06.30
BY (11) 15950
(13) C1
(19)
(73) Патентообладатель: Государственное
научное учреждение "Институт физики имени Б.И.Степанова Национальной академии наук Беларуси" (BY)
(56) RITCHER R. Int. J. Remote Sensing. 1997. - V. 18. - № 13. -P. 2827-2833.
RU 2082952 C1, 1997.
SU 1578478 A1, 1990.
JP 2006/148205 A.
EP 1059507 A1, 2000.
(57)
Способ определения калибровочного коэффициента для калибровки оптического
спутникового сенсора высокого пространственного разрешения в абсолютных энергетических единицах, в котором посредством указанного сенсора получают сигналы с двух
тестовых участков земной поверхности, размер каждого из которых соответствует пространственному разрешению сенсора, с известными различающимися альбедо и коэффициентами яркости поверхности, осуществляют подспутниковые измерения оптической
толщины атмосферы τатм и освещенности, создаваемой солнцем на двух указанных тестовых участках, а затем рассчитывают искомый калибровочный коэффициент k в соответствии с выражением:
µ0E 0
k=
T0 (µ )E норм (µ 0 )(rp , 2 − rp ,1 ) + T0 (µ 0 )T0 (µ ) (R p, 2 (µ, µ 0 , ϕ) − rp , 2 ) − (R p ,1 )(µ, µ 0 , ϕ) − rp,1 ,
π(N 2 − N1 )
где µ0 – косинус зенитного угла солнца;
µ0E0 – освещенность, создаваемая солнечным излучением на верхней границе земной
атмосферы;
Ni, где i = 1, 2, – полученные сенсором сигналы с i-го тестового участка;
T0(µ) = exp(–τaтm/µ) – коэффициент пропускания атмосферы для прошедшего без рассеяния света;
µ – косинус угла наблюдения;
E норм (µ 0 ) – освещенность тестового участка, нормированная на освещенность, создаваемую солнцем на верхней границе атмосферы;
rp ,i и R p ,i (µ, µ 0 , ϕ) , где i = 1, 2, - средние по поверхности i-го тестового участка альбедо
поверхности и коэффициенты яркости соответственно;
[
[
]]
BY 15950 C1 2012.06.30
T0(µ0) – интегральный коэффициент пропускания атмосферы при зенитном угле солнца, равном arccos µ0;
ϕ – азимут направления наблюдения по отношению к направлению на солнце.
Изобретение относится к области космических технологий, в частности к способам
калибровки во время полета спутниковых сенсоров оптического диапазона в абсолютных
энергетических единицах, и может быть использовано для калибровки спутниковых сенсоров высокого пространственного разрешения.
Известен способ радиационной калибровки спутникового сенсора, включающий измерение сигнала спутниковым сенсором, одновременные подспутниковые наземные измерения параметров атмосферы, коэффициентов яркости и альбедо подстилающей поверхности и расчет спектральной плотности энергетической яркости (СПЭЯ) солнечного
излучения на входе спутникового сенсора (SPOT calibration at the La Crau Test Site (France)
/ R. Santer [et al.] // Remote Sensing of Environment. - 1992. - V.41. - P. 227- 237; Reflectance
and radiance-based methods for the in-flight absolute calibration of multispectral sensors / P. N.
Slater [et al.] // Remote Sensing of Environment. - 1987. - V. 22. - P. 11-37).
Недостатками данного способа являются:
необходимость использования априорных предположений о модели вертикальной
структуры неоднородной молекулярно-газовой и аэрозольной атмосферы;
необходимость проведения расчета переноса излучения в сложной системе атмосфера подстилающая поверхность;
отсутствие учета влияния солнечного излучения, отраженного от неоднородной подстилающей поверхности в соседних пикселях, на сигнал, регистрируемый от тестового
пикселя. Такое влияние имеет место вследствие молекулярного и аэрозольного рассеяния
излучения в атмосфере и распространяется на расстояния порядка 5 км (Dave J. V. Effect
of atmospheric conditions on remote sensing of a surface nonhomogeneity / Photogrammetric
Engineering and Remote Sensing. - 1980. - V.46. - P.1173-1180; Kaufman Y. J. The atmospheric effect on the separability of field classes measured from satellites / Remote Sensing of
Environment. - 1985. - V. 18. - P. 21-34).
Точность калибровки спутникового сенсора по данному способу зависит неконтролируемым образом от точности модели атмосферы, точности кода расчета переноса излучения в системе атмосфера-подстилающая поверхность, однородности подстилающей
поверхности, окружающей тестовый пиксель, и пространственного разрешения спутникового сенсора. Третий из указанных недостатков может приводить к большой погрешности
калибровки спутниковых сенсоров высокого пространственного разрешения.
Наиболее близким техническим решением (прототип) является способ радиационной
калибровки спутникового сенсора, включающий измерение сигнала спутниковым сенсором, одновременные подспутниковые наземные измерения параметров атмосферы, коэффициентов яркости и альбедо подстилающей поверхности и расчет спектральной
плотности энергетической яркости солнечного излучения на входе спутникового сенсора,
причем размер тестового участка в семь раз больше, чем пространственное разрешение
сенсора (Richter R. On the in-flight absolute calibration of high spatial resolution spaceborne
sensors using small ground targets / Int. j. remote sensing. - 1997, V. 18. - № 13. - P. 28272833). Увеличение размера тестового участка в семь раз по сравнению с пространственным разрешением сенсора позволяет уменьшить неконтролируемое влияние неоднородной подстилающей поверхности в соседних пикселях на сигнал, регистрируемый
сенсором от тестового пикселя, и тем самым повысить точность абсолютной калибровки
спутникового сенсора высокого пространственного разрешения.
2
BY 15950 C1 2012.06.30
Недостатками данного способа, как и аналога, являются необходимость использования априорных предположений о модели вертикальной структуры неоднородной молекулярно-газовой и аэрозольной атмосферы, необходимость проведения расчета переноса
излучения в сложной системе атмосфера-подстилающая поверхность, а также неполный
учет влияния неоднородной подстилающей поверхности в соседних пикселях на сигнал,
регистрируемый сенсором от тестового пикселя, что приводит к погрешности калибровки
спутникового сенсора высокого пространственного разрешения.
Задачей изобретения является повышение точности абсолютной энергетической калибровки спутникового сенсора высокого пространственного разрешения и упрощение
процедуры калибровки.
Задача решается за счет того, что в предлагаемом способе определения калибровочного коэффициента для калибровки оптического спутникового сенсора высокого пространственного разрешения в абсолютных энергетических единицах, в котором посредством
указанного сенсора получают сигналы с двух тестовых участков земной поверхности,
размер каждого из которых соответствует пространственному разрешению сенсора, с известными различающимися альбедо и коэффициентами яркости поверхности, осуществляют подспутниковые измерения оптической толщины атмосферы τатм и освещенности,
создаваемой солнцем на двух указанных тестовых участках, а затем рассчитывают искомый калибровочный коэффициент k в соответствии с выражением
T0 (µ )E норм (µ 0 )(rp, 2 − rp ,1 )

µ0 E0

,
k=

π(N 2 − N1 ) + T0 (µ 0 )T0 (µ )[(R p , 2 (µ, µ 0 , ϕ) − rp , 2 ) − (R p ,1 (µ, µ 0 , ϕ) − rp ,1 )]
где Ni (i = 1, 2) - измеряемые спутниковым сенсором сигналы для i-го тестового участка,
R p,i (µ, µ 0 , ϕ) и rp,i (i = 1, 2) - средние по поверхности i-го тестового участка коэффициенты яркости и альбедо поверхности соответственно,
E норм (µ0 ) - освещенность тестового участка, нормированная на освещенность, создаваемую солнцем на верхней границе атмосферы,
µ0E0 - освещенность, создаваемая солнечным излучением на верхней границе земной
атмосферы,
T0(µ) = exp(–τaтm/µ), τатм - оптическая толщина атмосферы,
λ - длина волны излучения,
µ0 - косинус зенитного угла солнца,
µ - косинус угла наблюдения,
ϕ - азимут направления наблюдения по отношению к направлению на солнце.
Сущность изобретения состоит в следующем.
Регистрируемый спутниковым сенсором сигнал N пропорционален СПЭЯ I(λ,µ,µ0,ϕ)
солнечного излучения на входе спутникового сенсора, где λ - длина волны излучения, µ0 косинус зенитного угла солнца, µ - косинус угла наблюдения, ϕ - азимут направления наблюдения по отношению к направлению на солнце. То есть:
(1)
I(λ,µ,µ0,ϕ) = kN,
где k - искомый калибровочный коэффициент.
СПЭЯ солнечного излучения на входе спутникового сенсора равна:
µ E (λ)
(2)
I(λ, µ, µ 0 , ϕ) = 0 0
R (λ, µ, µ 0 , ϕ) ,
π
где µ0Ε0(λ) - спектральная плотность потока солнечного излучения, падающего на верхнюю границу земной атмосферы, R(λ,µ,µ0,ϕ) - спектральный коэффициент яркости на
верхней границе системы атмосфера-подстилающая поверхность. Ниже для сокращения
записи параметр λ будет опускаться.
3
BY 15950 C1 2012.06.30
В общем случае коэффициент яркости на верхней границе атмосферы можно записать
в виде:
T (µ)T(µ 0 )
rp + T0 (µ 0 )T0 (µ) R p (µ, µ 0 , ϕ) − rp + R бок (µ, µ 0 ) ,
R (µ, µ 0 , ϕ) = R атм (µ, µ 0 , ϕ) + 0 *
(3)
1 − rатм rпов
[
]
где Rатм(µ,µ0,ϕ)- коэффициент яркости излучения, отраженного атмосферой, T(µ0) - интегральный (по углам) коэффициент пропускания атмосферы при зенитном угле солнца
arccos µ0, T0(µ) - коэффициент пропускания атмосферы для прямопрошедшего (без рас*
сеяния) света, rатм
- сферическое альбедо атмосферы при освещении ее снизу, rпов - сред-
нее сферическое альбедо окружающей поверхности, R p (µ, µ 0 , ϕ) и rp - средние по
поверхности пикселя коэффициенты яркости и альбедо поверхности. В формуле (3) первый член учитывает вклад излучения, рассеянного в атмосфере без взаимодействия с поверхностью, второй и третий члены в сумме описывают вклад излучения, отраженного от
тестового пикселя, причем третий член учитывает поправку на неламбертовость поверхности. Четвертый член учитывает вклад излучения, отраженного от окружающей поверхности, который часто называют вкладом бокового подсвета.
В частном случае ламбертовской подстилающей поверхности, для которой коэффициент яркости не зависит от угла падения и угла диффузного отражения, т.е.
R p (µ, µ 0 , ϕ) = rp , из формулы (3) имеем:
T0 (µ)T(µ 0 )
rp + R бок (µ, µ 0 ) .
(4)
*
1 − rатм
rпов
С учетом соотношений (1)-(3) измеренный спутниковым сенсором сигнал можно записать в виде:

T0 (µ)T(µ 0 )
1 µ0E0 
r
T
(
)
T
(
)
R
(
,
,
)
r
N=
+
µ
µ
µ
µ
ϕ
−
(5)
α +
p
0 0 0
p
0
p ,
*
k π 
1 − rатм
rпов

где
(6)
α = [ Rатм (µ, µ0, ϕ) + Rбок (µ, µ0)].
Заметим, что здесь величина T0(µ) есть пропускание прямопрошедшего света, которое
определяется оптической толщиной атмосферы τатм, т.е.:
(7)
T0(µ) = exp (–τатм/µ),
а величина:
T(µ 0 )
E(µ 0 ) = µ 0 , E 0
(8)
*
1 − rатм
rпов
- освещенность тестового участка поверхности. Т.е. величина:
E(µ 0 )
T(µ 0 )
=
E норм (µ 0 ) =
(9)
*
µ 0 E 0 1 − rатм
rпов
R (µ, µ 0 , ϕ) = R атм (µ, µ 0 , ϕ) +
(
)
есть освещенность тестового участка, нормированная на освещенность, создаваемую
солнцем на верхней границе атмосферы. Поэтому формулу (5) можно переписать в виде:
1 µ0E 0
N=
α + T0 (µ)E норм (µ 0 ) rp + T0 (µ 0 )T0 (µ)(R p (µ, µ 0 , ϕ) − rp ) .
(10)
k π
Произведя измерения на двух (или более) тестовых участках с различающимися коэффициентами яркости подстилающей поверхности, из формулы (10) получим:
1 µ0E 0
Ni =
α + T0 (µ)E норм rp, i + T0 (µ 0 )T0 (µ)(R p ,i (µ, µ 0 , ϕ) − rp ,i ) , i = 1, 2.
(11)
k π
[
]
[
]
4
BY 15950 C1 2012.06.30
Здесь считается, что тестовые участки находятся близко друг к другу и их освещенности можно считать одинаковыми и равными E(µ0).
Из формулы (11) определяется искомый калибровочный коэффициент:
T0 (µ )E норм (µ 0 )(rp , 2 − rp ,1 )

µ0E0

.
k=
(12)

π(N 2 − N1 ) + T0 (µ 0 )T0 (µ )[(R p , 2 (µ, µ 0 , ϕ) − rp, 2 ) − (R p,1 (µ, µ 0 , ϕ) − rp,1 )]
Все величины в формуле (12) являются либо хорошо известными априори (E0), либо
измеряемыми. Так, величина Ni есть сигнал, регистрируемый спутниковым сенсором, rp, i
и R p, i (µ, µ 0 , ϕ) - средние по поверхности i-го пикселя коэффициент яркости и альбедо по-
верхности, Eнорм (µ0) - нормированная освещенность, связанная формулой (9) с измеряемой реальной освещенностью тестовых площадок E(µ0). Величина T0(µ) рассчитывается
по формуле (7) по измеряемой величине τатм.
В случае ламбертовской поверхности формула (12) упрощается и принимает вид:
µ0E0
[T0 (µ )E норм (µ0 )(rp,2 − rp,1 )] .
k=
(13)
π(N 2 − N1 )
Если с целью повышения точности определения калибровочного коэффициента k измерения производятся на нескольких (более двух) тестовых участках, то искомая величина находится из уравнения (11) каким-либо методом статистической оптимизации,
например методом наименьших квадратов. В случае ламбертовской поверхности формула
(11) имеет вид:
1 µ0E0
(14)
Ni =
α + T0 (µ)E норм (µ 0 ) rp,i
k π
и искомый калибровочный коэффициент k определяется по линейной зависимости
N i ( rp, i ) .
Указанная последовательность выполнения операций позволяет упростить процедуру
и повысить точность радиационной калибровки спутникового сенсора, поскольку при
этом автоматически учитывается вклад в регистрируемый спутниковым сенсором сигнал
как излучения, отраженного атмосферой, так и излучения, отраженного поверхностью,
окружающей тестовый пиксель, т.е. бокового подсвета. Поэтому отпадает необходимость:
а) использования априорных предположений о модели вертикальной структуры неоднородной молекулярно-газовой и аэрозольной атмосферы;
б) учета радиационного влияния отражения солнечного излучения от неоднородной
подстилающей поверхности в соседних пикселях на сигнал, регистрируемый от тестового
пикселя;
в) проведения расчета переноса излучения в сложной системе атмосфераподстилающая поверхность.
Таким образом, радиационная калибровка спутникового сенсора высокого пространственного разрешения предложенным способом производится без неконтролируемых
ошибок, обусловленных априорными предположениями о модели атмосферы и влиянием
бокового подсвета.
[
]
Национальный центр интеллектуальной собственности.
220034, г. Минск, ул. Козлова, 20.
5
Документ
Категория
Без категории
Просмотров
0
Размер файла
95 Кб
Теги
патент, by15950
1/--страниц
Пожаловаться на содержимое документа