῍ ,**1

Journal of The Remote Sensing Society of Japan Vol. ,1 No. / (,**1) pp. .1.ῌ.2,

RSSJ

 

¹REº»¼RL7 M=()Ey½y.@&ABCDEFG?HI¢ *+ῌ£m!¤*,ῌ"0 *+ῌ¥#¦$*,ῌ§%¨©*-ῌª&«'*+ῌ'¬§*+ῌ(­¤*+ῌm ®'*+ῌ¯° ±*+ῌ²)!*.ῌ
³0*©+*.ῌ´&µ©*,ῌ¶·¸'+*,
Mapping Global Precipitation with Satellite Borne Microwave Radiometerand Infrared Radiometer Using Kalman Filter
Shunsuke NODA*+, Kazushi SASASHIGE*,, Daisuke KATAGAMI*+, Tomoo USHIO*,, Takuji KUBOTA*-, Ken’ichi OKAMOTO*+,
Yasuhisa IIDA*+, Satoshi KIDA*+, Shouichi SHIGE*+, Suguru SHIMOMURA*+, Kazumasa AONASHI*.,Toshiro INOUE*., Takeshi MORIMOTO*, and Zen-Ichiro KAWASAKI*,
Abstract
Estimates of precipitation at a high time and space resolution are required for many important applications. In thispaper, a new global precipitation map with high spatial (*.+ degree) and temporal (+ hour) resolution using Kalmanﬁlter technique is presented and evaluated. Infrared radiometer data, which are available globally nearly everywhereand nearly all the time from geostationary orbit, are used with the several microwave radiometers aboard the LEO
satellites. IR data is used as a means to move the precipitation estimates from microwave observation during periods
when microwave data are not available at a given location. Moving vector is produced by computing correlations
on successive images of IR data. When precipitation is moved, the Kalman ﬁlter is applied for improving the moving
technique in this research. The new approach showed a better score than the technique without Kalman ﬁlter. The
correlation coe$cient was *.+ better than without the Kalman ﬁlter about 0 hours after the last microwave
overpasses, and the RMS error was improved about *.+ mm/h with the Kalman ﬁlter technique. This approach is
unique in that +) the precipitation estimates from the microwave radiometer is mainly used, ,)the IR temperaturein every hour is also used for the precipitation estimates based on the Kalman ﬁlter theory.
Keywords : Kalman Filter, Infrared Radiometer, Precipitation Map, Microwave Radiometer$%&" '  ()*+$, ' -.,/0*+.
+ῌ ῏ ῍ ῐ ῎

1&2,&3.456*+789:2;$<&=> :?
@96AB?CDEFG?HI., ()JKL7 M=/B-N?CDEFG, DOPQJR?S TU

VWXY,/0JKL2?ZV;[X\M&"]P^?_` aTb cde?f0 gDhi
k.2l\m$%&n3o$6Y ().@&CDj?HI, WCRP World Climate

.pM\qResearch Program ?r? GPCP Global Precipitation Cli-sY Mtu\M&" v?HI.wX\ /0xy,

matology Project .z{u&@9. v?m|3t}~.B.X\po /B? 1 0?

.8 u\M&-" GPCP $,()ῌy7>&0$,.6M=> /0xy$AEy6t./0xy?JKL; u\

BN?CD?7.+&:2,M&;

,**1. +. ,.  ,**1. +*. - 
*+      +ῌ+*,     ,ῌ+*-  *. 
!

*+ Graduate School of Engineering, Osaka Prefecture University, +ῌ+Gakuencho, Sakai, Osaka /33ῌ2/-+, Japan*, Graduate School of Engineering, Osaka University, ,ῌ+ Yamadaoka, Suita, Osaka /0/ῌ*21+, Japan# 474 #

¡ ()?ῌyJKL7r

 
$%&' ()*+

!"#
!,-#



Table +

./ Vol. ,1 No. / ,**1

Characteristics of geostationary infrared
data used in this study.

-0" 12 34%5.6"78 9
# .:!;< <52=>": ?/- GPM Global Precipi?@ABC./-tation Measurement

D E( '  F GH87IJK !
".#$%84L-0" M ?2 #NO4
.=>"

@&P3L7'

(Q):.*+/ ,9R-ST"=%Table , Characteristics of microwave radiometers
used in this study.

%U@S!V - %J".2W" R/ 34 + %5# 0.921"78 + %214L"34X5:&P

?6.+0- 78>"

% , 9?6O4 YZ:[\.;</ #N
=L7

!"#N]Z:.^_>?`@GJ3L-0". M CMORPH The Climate Prediction Center Morphing Technique abL">?2 

!"N]Z:.;<>"Ac+0-S; ML!;,ῌ ῐῑ῏῍῎!"=%.B%/ <52<J#N@de92fJ,-0" V/- IPWG International
Precipitation Working Group/ JghijkJg.

lm2 :&P+07,pC!,- VD+E@<qFr3

&P

L-0"

" 5J+07jk.tuK>

,934Q):&P

jkO4 ;3L"

ῌ

no34X5:

jk.G+/- # *s.HI>
34:&P NOAA  CPC Climate Prediction Center  Table

Z:[\.+0">? <5# *sHI
LMDv2W"M@K3L-w-0"@ NO-G

+ W"34:O4

QR>" M!_J>?2 /7jk. + %2Sz/-0" T03L-0"

A\yjk.P1/ xI{)!,-3L7#N  Z:[\!,-# )2 . km 2W" }Y 0*῍Nῌ0*῍S 2.W3Z:3-0"78 #N@X

Z[\@ ]3L-0J0 ML@#N^W"

>" M-.B_78 CMORPH 2 %_2/

|UVjk}ῌ~56 *.*-0-/7 

jkEjAQ>"

!"N .Z:[78`a.-W" ;:buO40c\0-Z:3"*Fs '.de vfv78(qgG4L/-0" /O/J@4 A\k)E@hJL"J" V78 CPC :buO40c2 g-@W;

34:O4jk.kcx-`a/-0" 7.JgiA\k)E@j¡

3L"As¢ £ +l2 E¤Fs ':&P+¥ ¦\ t# IR mn.+0" c2@o4L"M@p/0

:%2O§J,7qE@W"@¨rs

VM2 lm2 M_/7-.B0 E¤F©¥ªh £ ^L-0J0 /O/Mv

;ST"  noEjAv!;OJ;t30

jk.G+/

7 *.+  + %59# *s>?

ῌ

0-u«>"  CMORPH rvJZ:[lm2 Table , W" / ¬(Q):&P

(Q):&P

\.*+/7# *sHI.^_@ Z:[\



!,-#

!,-=3L7jkO4W®3L7N j

.W3Z:3"w2w­7

NO@W"78 jk.+0-0" ML4¨k.e /-+0-0" N W®A\y

O4/7# ¯%W"°±@W"Mx{) ² ³yz´¬µ JST {m|}lm%~A/ C\ E\k>?.+0-¶>"M· CREST lm]  ¸¹j!,- *.+ ῌ + %5#º) lm»O :!"<<59#.o"

*s

*+ '*sHI lm : ¼½R @`@A\y{).+0-0"02 MA\y{)!"R7s 475  

Fig. +

 !"#

Flow chart of the algorithm used for the present
analysis.

Fig. ,

Schematic diagram of the moving vector derivation to calculate the correlation coe$cients of thesuccessive IR images.

-./$ GSMaP %Global Satellite Mapping of Precipitation& 0123$ 4567$ 8! 9  800

:' ; $ 4<: 
-.$ A

B.!G-.

(=> ῌ? ! @4A B. CD EF$ >H3

ῌ?

!4$ 0./῍RS ' 8!TU$ NVWXY!Z!A

Q3$ GSMaP_MVK$ GSMaP_MV 0:'$ ^_`CD80Oab=c 'ῌ[Y-\]-W

de4/fghiOSj [Y-\]-W
r#sY

KLMNOP3

0./῍4I5J

tu!sY0Sj

Fig. , H!kl m:' n_o )t* 4 IR !p q

' fg[Y-0./῍!y0$ z!_`! )t++hr*vwx> ῌ?

/$ {Y |:80O4}O$ H!~

4!y4 IR !p q !!`4`/  b' 8!  \]-W!  

=KLMN CD' !/$ _` )t* j!S$  !" J`#$=# %/ km y mSjU$ %US y &!4m

|4}bS !O$ [Y-\]-W

3j  @'4$ ()$ *+) +  %0&4n'

!I & ,A j )t++hr* 0!KLMN CD' HSj8!KLMN!I5-    }=.=c ($ *+)!A !NO$ A
B.0Sja3' 8!bSjCD A

Fig. + $ de4/S  "#!\B.4$ H!ynPP3¡¢£!k0 m:' de!\ ¡¢£4/$ 12: + _`,!sY!> ῌ?
4$ IR !pq

-, 

0./῍y'!`=KLMN 

A B.!G :

80cj$ Fig. + nA B. 3a:'HSj$ + _`¥! A ' 4/$ 4~S ¤A B. 

A B.cj4¦§A  D$ !/5¨©6Oª«3j

jA $ Fig. + n7`sY "#:' HSj$

= ' 8!!/5 W¨©6 W 38:
O *.+ mm ῌ hr 6®!°r

¬­!p q O ,1* K 6®!_/ ¯0Sj$ O *.+ mmῌhr 6®!_/ 9S  80S ' de4/$ p ±²:' I7$ ¬­!sY

q 0P~S 0!LM

³8S$ ¬­!D$ ¤:;:' Fig. - H!LM m:' ´µ

"#:'-+ 

cj~S ¤ $ ¶µ!p

A B./$ 12:_`!!sY

q!^PCDO$ H!0}!!sY
!`= @ $ de4/> ῌ?

m:' 8!lP«b$ p q <!`/$ -}== b5!!$ >/?

0./῍0·@!LMOA­:80Om3j ' de4/$ 8

S ' A B. CD D!`= @O!>=LM0!sY/ + _`,P3

32$ ¸B!=Z¹!C3/DºP3O$ ¸B!

80:;S$ H3 ES$ A-}=Z¹!C3/DºP3= ' ; `= @O-

B.cjA S ¤}32$ » ( !Z¹!C3O¼=

± ~:'D$ H3

6®!bp q 0

Q3!Z¹! }ODºP3= ' 8! Dde4/$

P~S 0!LM:;:80cj$ ½F, 476 ,



 Vol. ,1 No. / ,**1

yz1:W@5%> A%yzG.;{|<  G}~%FQk : = + 8%>?
>5.ῌ ῏ΐῒ῎ῑ῍ῐ῔$%&'()*+, + -. *.+/ % 0%!"#12 345 Fig. .  67 2 34,**/  1 8 - 9 * - UTC % + -. *.+/ %: 0*;<: 0*=>%!"#?@

A%&' BCDEF%GH1 IR I&5J)&K J34 : 0*;<: 0*L5MN 0%!"#1 + -. *.+/ >2

345

(!O P57 Q

R S'7TU45VW(!O
1P<45%
Fig. -

Geostationary satellite IR brightness temperature(ordinate) relative to precipitation estimated frommicrowave radiometers (abscissa) during July ,**/.

!(Z'"%!H1[45 3< \%!B] #$@57 ^_ %&'(;<`)*+a;Y,b!1P<45 A%,b!K


.;{7G.;{%&'<@5: ;<X: ;YR cῌdc-.> 1 8 -  ,** mmῌday e^5f

b1g/0&KhiG34Qj A%,b!
&Kk <34 k%>l5A71?345 A

xk+Fxkw

"@.;{ykFxkv

G.;{xk : [W%!4xk+ : + -. %!4j 56(-HQ uv>2 5!"v : G}~#1 w78xg9>Rl5k%%56(!bH?

w :

s}~

%m12 Fig. / ?@c dc-. 37@5n

-H> opqrpsb4F1t!b4?QF : ]C&55A71U;5 3<6-. %!"# Fig. 0

yk : AEF&KG345 >R ]C%&5!"# GSMaP_MV

H : AEF  7BCDEF;<B<   !"# GSMaP_!47%yz

MVK R !4%CD[5%1P45

Fig. - E

A%&'2 34 !"#%8x
[-F>R G.U 5HI134( "@.;{ :

GJ%I*+KLM%NRO %]C&K!b35- ]C PQ>R RSBCDEFT?7%&5"#2>R @5U>PV !b4JRX ( ; hi%!b[/>8x7 WQY5opqrpsT?7 8xY¡@5

R -.%¢w7£ b4RC@5 A%1 s!"#%2 )Z>R BCDEF&K

}~(KQj "@.;{?@[j7(K

B<34 !bH]C*+¤5s}~R BCDEF>B< !tBCDEF1[w - %B<?>\]C&K + ¥3 - ¦§B

,¤' A%\,-R 5- ]C&CDEF1[K -B< !t 7%¨

5B<?k:W@5  BCDEF&5b4B

F5 BCDEFR \R>]©1ª<?T?7«^V '¬¨­§(j =  !4®@5q¯°kª(5V%<48x@5A71®^5 3<

1l51 uv>R ±6(_`75 %  BCDEF&5b4B<R-H²³a&5hiR ª(5µd¶Y(%> &j6(8x1· >l5´;<b®345!4RcF

(¨ª1:W@57TU451 %¸R ¹

Fig. 1  BCDEF1[w;<%-.1¢%ºe7 5» ]C&5B<­§ G.;{ :

¼^ V%B<?%BCDEF&5B<?

6. f%&'-DgW%!b47AEF>

@5hyzi1Z%&'J@5%;?@ A%2

G345  %.R cFR½²;(j?345&' -.1¢5» klhyz

 477 

Fig. .

Fig. /

 !"#$%&'The global *.+ degree hourly rainrate map at * UTC on - July ,**/ in mm/h. GSMaP_MVK
stands for Global Satellite Mapping of Precipitaion using Moving Vector with Kalman ﬁlter.

Fig. 0

GSMaP_MV (Moving vector only) and MVK(Moving vector(Kalman ﬁlter) around Japan
at . UTC on - July ,**/ in mm/h together
with Radar-AMeDAS data (bottom panel).

) 478 )

GSMaP_MV (Moving vector only) and MVK(Moving vector(Kalman ﬁlter) near Japan at/ UTC on - July ,**/ in mm/h together with
Radar-AMeDAS data (bottom panel).



Fig. 1

Correlation coe$cient between the GSMaP_MV(Solid line) and the GSMaP_MVK (Dashed line),
and the precipitation estimated from microwave
radiometers for the period of July ,**/ as a function of the past microwave radiometer overpasses.

 Vol. ,1 No. / ,**1

Fig. 2

Correlation coe$cient between the GSMaP_MV(Solid line) and the GSMaP_MVK (Dashed line),and Radar AMeDAS for the period of July ,**/
as a function of the past microwave radiometer
overpasses.

Fig. 3

Rmse between GSMaP_MV (Solid line) and GSMaP_MVK (Dashed line), and Radar AMeDAS for theperiod of July ,**/ as a function of the past
microwave radiometer overpasses.

<=7"# 0>?@ABCD:EF'G:HI
JKLM# + HI;$#+ , HI:NO(+ ./0
12/34PQ8NORSTU? /:V-W$XYZR+ [R\]^_G!`#:a7" . HIb / HIb(+ ^c-`d ./0 12/3:ef
g!"$hi-;$jR&'$
kl mn-op$qrstruv346#" wxy<z{|} , ~^ RMSE :'G
4 8f4 Fig. 2 z{ Fig. 3 - *+./0
12/346#`#STU? /:V-W$9:( )*+./0

12/3456789:(;$ ,!

!wxy< |} , ~^R` "od +STU? /-W "0>?@ABCD-W "XY784
ST78< HI R` "#$  "  * :+0>?@ABCD-W "XY78Rqr
struR:JKR` "#$ Fig. 2 z{ 3 Wd ST
U?

/-W"ST78<P$Rwxy<+c-=7 |} , ~^+`G'-P7"#$
]G:f(9./0 12/345678l!8
f"#$jR&'$8 0+ + HI- !$:( 7

 ZR7"¡ῌ7"#$ + HI¢£¤:0-+

8 + HI:¥:0Rqrstru4JK780>?@ABCD(XY78:¦- ,!4ST§4 Fig. +* - *+STU?

/:V-W$9:

U? /-W " + ST78῍¨ , ST78

( )*+,!-./0 12/3456789:(;$῍©ª!"#$  !80 «¬R

-wxy< + ,**/ ­ 1 ® ,  * H'G ,. H

7" 0>?@ABCD:V:0RJ¯$R °±

(4 j:§'G9STU? /p(+` .

(  !"#$%&'$

/0 12/34²678l!8f"#$jR
&'$ 8+0>?@ABCD:V(78.ῌ ῎ ῍ ῏

0R:JK(;$ 0>?@ABCD+³´µ¶·¸:8¹ HIEF7`# j:Wº`NO+ §» 479 ῎¼(+ *.+  + HI¢£¤:«½04 ῎ !"#$%&!"#%&) $1& JK4&6L M$NOP+ Q
QGR S3HI, H + ῌ>? TUO P+ VW
012 34567 + (89:8;<=

!"#%&) >?@AB$C*+ DEFG

3HI,

ῌῑῐ῍

Fig. +*

Time series of correlation coe$cient between the
GSMaP_MV(Solid line) and the GSMaP_MVK(Dashed line), and the Radar AMeDAS during** UTC through ,. UTC on , July ,**/. The
product only from the microwave radiometer isalso added to this ﬁgure indicated by “-”.

&I X+ YAZ[\]6L^_`a῎

4+ 

YA bcIde fg3 , h34+ %&3 !"#$fijAklmkn

῎

4T7

, o

p+ ^_`a$q ῎ ders4+ 
tuZ[ 0 vcw$rxN+ yz{|} *.+ ~:+RMSE } *.+ mmῌh ~:6LQ}Z7 , w+

" A#n$ X+ !῍0
$῎ 4 , H + ῌΐN+ $hLL$$ῑ+ $"# %&34 }+ |l ¡¢B£$ }O X+ ῎ v$¤n¥B¦¢ῌ¦¢$ῑL fg w§34 ,῏῎ : ῌ>?+ ¨©ª«¬­® .JST' ¯ῐ°±²>?
³´µ .CREST' >?ῒ¶ (!$·¸¹ºY¡»῏῎¤n¥B) $>?¼½ (+' C. Kummerow, Y. Hong, W. S. Olson, S. Yang, R. F. Adler,
J. McCollum, R. Ferraro, G. Petty, D. B. Shin and T. T.
Wilheit : The Evolution of the Goddard Proﬁling Algorithm(GPROF) for Rainfall Estimation form Passive Microwave
Sensors, J. Appl. Meteorol., vol. .*, no. ++, pp. +2*+ῌ+2,*,,**+.
,' Iida, Y., K. Okamoto, T. Ushio and R. Oki : Simulation ofsampling error of average rainfall rates in space and time byﬁve satellites using radar-AMeDAS composites, Geophys.
Res. Lett., --, L*+2+0, doi : +*.+*,3/,**/GL*,.3+*, ,**0.
-' Hu#man, G. J., R. F. Adler, M. Morrissey, D. T. Bolvin, S.
Curtis, R. Joyce, B. McGavock, J. Susskind : Global : Precipitation at One-Degree Daily Resolution form MultiSatellite Observations. J. Hydrometeor., ,, -0ῌ/*, ,**+.
.' R. J. Joyce, J. E. Janowiak, P. A. Arkin and P. Xie, :CMORPH : A Method that Produces Global Precipitation
Estimates from Passive Microwave and Infrared Data at
High Spatial and Temporal Resolution, J. Hydrometer., vol.
/, No. -, pp. .21ῌ/*-, ,**..
/' E. Ebert, J. Janoviak, and C. Kidd : Comparison of NearReal-Time Precipitation Estimates from Satellite Observations and Numerical Models, Bulletin of the American Society, Vol. 22, pp. .1ῌ0., ,**1.
0' K. Aonashi, A. Shibata and G. Liu : An over-ocean precipitation retrieval using SSM/I multi-channel brightness temperature, J. Meteor. Soc. Japan, vol. 1., pp. 0+1ῌ0-1, +330.
1' K. Aonashi and G. Liu, : Passive microwave precipitation
retrievals using TMI during the Baiu period of +333. Part I :Algorithm description and validation, J. Appl. Meteor., vol.
-3, pp. ,*,.ῌ,*-1, ,***.
2' T. Kubota, S. Shige, H. Hashizume, K. Aonashi, N.
Takahashi, S. Seto, M. Hirose, Y. Takayabu, K. Nakagawa,K. Iwanami, T. Ushio, M. Kachi and K. Okamoto : Global
Precipitation Map using Satelliteborne Microwave Radiometers by the GSMaP Project : Production and Validation,
IEEE Trans. Geosci. Remote Sens., in press, ,**1.

89:8;<=

/ 480 /



<>?@A=
    :       +1     +3     $ %& '()
4 56      :    NO7 -.789
S TU>?V W  +2  O7YZ[\  %&  NO7 -.789S TU>?V W/bc01 defR!g@ijklm!g@ijn%opq=rsdtuI
JKLvw

 Vol. ,1 No. / ,**1    
 :       +302
 !
"#  +31   *+, -./01 +31- 23 % 789:;  +33<=>?@AB +333
CDE ,***
FG

   !H <=>?IIJKL IEEE GRS Society M 
O789
PQR M 
M NX
M +33* 
] +33- NX
"^_] +33. ` ] +33/ 2ah +332 b"] ,**, #$]( )*+
 :       +2    %&   ~'-. /bc01 defR!g@ijklm!g@ijn%
opq=rsdtuIJKLvw 
M

 ! "# $%
 :       +313- xy *,   zF
&{| *.
    /bc01$ *1
  /}c01$ / Jwr AIpI u  J+vwJKL NX
M 

M

 /0

 ,

 &'

12 3   :    NO7 -.  /N  1   $/01 $  +*
  /01$
 NASA d [& j National Space Science and Technology Center/Global Hydrology
and Climate Center USRA M  +,       +/ T  +2   NO7 -.!H  / <= A *!IJK
 N
 NX
 NN
L O789
 IEEE <=+,D M
 #6 78 9:) )9;- 
 : :;  ++ +-   $/01$  +0-   /01$- 2§ :+M  +1  ¨ ¡¢:;Mk²³ J³6  +3 ´
:; <= A jM  / J
wrs A +N¼ IJKL  NX
5 IEEE NX
5
<=+,D M

- .
 
 :       +1     +3    ~'-.
/bc01$ %&    /}c01 +  (r A TRMM )r PR { TRMM defR A TMI I 
Aijn%opJKLvw NXM  4  53
 :       
+33/
- 
-    , +331-  
<=. -.$/01$ ,**+
<=. -./01$ /   +333
¡¢ /0M ,**+ L£ % :; ¤1M ,**.     ,**1
! ¥¦©[ª3«¬,I­ ij®¯I°(±(uRPµe[{rs¶·[¸¹ºJ+vwJKL
 4»NX 4»<=+ ,D 
vw NX

M

½ 481 ½  !"#$%

 

 

+, : -./  -0-0120345678920: &+322 G-.-0-01HI20345JKLMNO 'PQ20( ]& ^G50_`a6789_`34+34b) +33, G-.-0HI20j?@AB205kl) +331 G^-01H
I20345kl) ,**, G-./  -0-0120345'678920( xy) ,**1 G^345z{|}n~
& #) +332*+333 G@D -0b34b& ?
@ABPQ}#34pq& }) LMI'Linear Matrix Inequality( HPQ\@DtvPQno#34pq& JK '20( '.-0(& IEEE) AIAA) ?@ABPQ 0s) r 678
90sno#sb& ()*+ ,-.
/0#
+, : =Ea=E34+&  /1 G=E-0 ¡& ¢% 3
G 0JK¥¦ ';<-0(&  0- G~=E34+ 3
4jª& 34:«) tmDuC?vFD#¬± mF²#^³) VW# #tutD´& r =E0s+,&   
+, : ;<-0=>?@AB34CD& =E
FD
'R:STU 'Split Window(( VW#XYZ[D\=E cde #fg) =EFD =>\hi=
Eῌm@Dno#34pq& r =E0s) r tmD
uC?v0s) w[t=E0ssb& 0JK&     +, : -.-0-0120345= 20&,*** G-.-020j D205& ,**/ G^-01JKMNO]& JK '20(& ^G^-0120
345kl& & !  "#$ %&'+, : -.-0-0120345= 20&+31- G-.-020j£¤205& +312 G^-01JKNO§ '20JK( +313 G) ¨©-0kl 734+&+32/ G@­®DF¯­°"Z-0=E34+bxy&+323 G-.-0xy) +33+ G^k{|) ,*** G-.-0-0120345{|& ,**1 G ¯µ-==0¶bs (ICAE) ¶b·& r -==0s0`¸) =0s¹¸) ºD°0s¹¸»¸¼¸&½ 482 ½