JPWO2016125206A1 - 合成開口レーダ信号処理装置 - Google Patents
合成開口レーダ信号処理装置 Download PDFInfo
- Publication number
- JPWO2016125206A1 JPWO2016125206A1 JP2016572940A JP2016572940A JPWO2016125206A1 JP WO2016125206 A1 JPWO2016125206 A1 JP WO2016125206A1 JP 2016572940 A JP2016572940 A JP 2016572940A JP 2016572940 A JP2016572940 A JP 2016572940A JP WO2016125206 A1 JPWO2016125206 A1 JP WO2016125206A1
- Authority
- JP
- Japan
- Prior art keywords
- height
- phase
- scatterer
- fringe
- unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000605 extraction Methods 0.000 claims description 123
- 230000015572 biosynthetic process Effects 0.000 claims description 34
- 238000003786 synthesis reaction Methods 0.000 claims description 34
- 239000000284 extract Substances 0.000 claims description 21
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000002123 temporal effect Effects 0.000 claims 1
- 238000000034 method Methods 0.000 abstract description 47
- 238000010586 diagram Methods 0.000 description 28
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000000737 periodic effect Effects 0.000 description 4
- 238000001914 filtration Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
φz=W{(2・π・p・B/λ・R・sinθ)・z} (1)
W{}:2πのラッピング
p:観測方式を表す係数(シングルパスの場合p=1、リピートパスの場合p=2)。
λ:照射電波の波長。
θ:照射電波のオフナディア角。
R:軌道k1、k2の中点から画像中心までの距離。
B:軌道k1、k2の直交ベースラインの長さ。
であり、φzとzは比例する。ただし、位相の値は2πでラッピングされた値である。以降では、直交ベースラインを単にベースラインと呼ぶ。別の散乱体であっても同じ高さzであれば、同じ地形縞の位相φzの値をとるため、観測されたSAR画像の地形縞の位相の値からSAR画像中の散乱体の高さを推定する。さらに、SAR画像の全ての散乱体の高さを推定することでSAR画像を3次元化する。
zmax=(λ・R・sinθ)/(p・B) (2)
である。ここで、式(2)でのBは異なるベースライン長の中で最短のベースライン長である。
本発明は上記のような課題を解決するためになされたもので、SAR画像上の散乱体の高さ方向を推定するにあたって、最短ベースライン長に対応する式(2)の高さzmax以上の高さを持つ散乱体の高さを推定し、その散乱体の画像を抽出することのできる合成開口レーダ信号処理装置を得ることを目的とする。
実施の形態1.
実施の形態1では、異なるベースラインの複数のSAR画像(各画素の地図情報を含む)と、すべてのSAR画像を撮像したセンサの軌道の位置情報(緯度経度または地図座標と高さ)を用いて処理を行う合成開口レーダ信号処理装置について説明する。
Bc=(λ・R・tanθ)/p・r (3)
r:グランドレンジ分解能
φs=W{(2・p・π・(r1−r2))/λ} (4)
r1:プラットフォームk1と散乱体αの距離
r2:プラットフォームk2と散乱体αの距離
である。式(4)で示されるように、反射信号の位相の差(干渉位相)φsはプラットフォームk1、k2それぞれから散乱体αまでの距離の差r1−r2と比例関係にある。ただし、位相の値は2πでラッピングされた値である。
φg=W{(2・p・π・(r’1−r’2))/λ} (5)
r’1:プラットフォームk1と散乱体α’の距離
r’2:プラットフォームk2と散乱体α’の距離
である。軌道座標の位置情報とすべてのSAR画像の地図情報は既知であるため、それらの値からr’1、r’2を算出することができる。したがって、式(5)により、軌道縞φgの値を算出することができる。
地形縞の位相φzの値は、理論的には以下の式(6)のようになる。
φz=W{(2・π・p・B/λ・R・sinθ)・z} (6)
すなわち、地形縞の位相φzと散乱体高さzは比例する。ただし、位相の値は2πでラッピングされた値である。別の散乱体であっても、同じ高さzであれば、同じ補正地形縞の位相φzの値をとる。
φz0=W{(2・π・p・B/λ・R・sinθ)・z0} (7)
m=2・π・p・B/λ・R・sinθ (8)
φz0=W{m・z0} (9)
で求められる。
x0=z0/tanθ (10)
である。したがって、SAR画像を3次元で示す場合、グランドレンジのセンサ側にx0ずれた長さ分を、散乱体の高さz0に応じてx0だけグランドレンジのセンサの反対側に位置を補正する。
実施の形態1ではセンサ軌道の位置情報から、センサと散乱体間の距離r’1、r’2、R、ベースライン長B、及びオフナディア角θを求め、軌道縞φg及び軌道パラメータmを算出する。この際、r’1、r’2、R、Bの精度はセンサ軌道情報の精度に大きく依存するため、センサ軌道の位置精度が必要となる。これに対して、実施の形態2では、センサの軌道の位置情報の精度が十分でない場合にも、r’1、r’2、B、Rの値を用いる代わりにSAR画像を撮像したセンサの照射電波のオフナディア角θを用いて、軌道縞φg、及び軌道パラメータmを高精度に算出する方法を示す。以下、実施の形態1と同じ入力・出力データ、装置、ステップには同一の番号を付し、説明を省く。
vn=exp(j・φs) (11)
vnを空間のフーリエ変換する。これによって、干渉位相φsの空間の変化を周波数領域に変換し、ピークを持つ周期成分のみをBPF(Band Pass Filter)によって取り出す。例えば、周波領域で、ピークをもつ周波数のみを取り出す。さらに、これを逆フ−リエ変換し、この位相の変化成分を軌道縞φgとする。
Δx=(λ・R・cosθ)/(p・B) (12)
と表される。
m=2・π・p・B/λ・R・sinθ
で与えられたが、実施の形態2では式(12)で出力された軌道縞周期Δxを用いて、軌道パラメータmを次式で表す。
m=1/(Δx・tanθ) (13)
ステップST2043(軌道パラメータ算出)では、式(13)の計算により、軌道パラメータmを算出する。
実施の形態1、2では、SAR画像を撮像したセンサの軌道の位置情報や(例えば、衛星の軌道座標)や照射電波のオフナディア角など、センサの情報が必要であるのに対し、実施の形態3では、センサの情報が無い場合にも処理を可能とする方法として、センサの情報の代わりに、異なるベースラインの複数のSAR画像に加えて、指定の高さのGCP(地上基準点、Ground Control Point)データを用いて、SAR画像から指定の高さのGCPの散乱体と同じ高さの散乱体を抽出する形態について説明する。以下、実施の形態1、実施の形態2と同じ入力・出力データ、装置、ステップには同一の番号を付し、説明を省く。
Claims (9)
- 第1のベースライン長を有する2つのセンサそれぞれで生成される第1の2つの合成開口レーダ画像を用いて該第1の2つの合成開口レーダ画像上の同じ散乱体を表す2つの画素に含まれる信号の相対位相を複数の画素に対して表した第1の地形縞及び、第2のベースライン長を有する2つのセンサそれぞれで生成される第2の2つの合成開口レーダ画像を用いて該第2の2つの合成開口レーダ画像上の同じ散乱体を表す2つの画素に含まれる信号の相対位相を複数の画素に対して表した第2の地形縞を算出する干渉位相処理部と、
特定の高さの散乱体に相当する前記第1の地形縞での位相を示す第1の特定位相及び前記特定の高さの散乱体に相当する前記第2の地形縞での位相を示す第2の特定位相を算出する位相算出部、前記干渉位相処理部で算出された前記第1の地形縞及び前記第2の地形縞の中から前記第1の地形縞で前記第1の特定位相となりかつ前記第2の地形縞で前記第2の特定位相となる前記特定の高さに相当する画素を抽出する画素抽出部、を含む抽出処理部と、
を備えたことを特徴とする合成開口レーダ信号処理装置。 - 前記特定の高さは、前記第1のベースライン長と前記第2のベースライン長のうち短いベースライン長を有する2つのセンサのみで測定な可能な高さより高い、
ことを特徴とする請求項1に記載の合成開口レーダ信号処理装置。 - 前記干渉位相処理部は、前記第1または前記第2の2つの合成開口レーダ画像上で高さが既知で同じ高さを有する散乱体を表す少なくとも3つの画素を選定し、前記少なくとも3つの画素に含まれる信号の位相が前記2つの合成開口レーダ画像上で同じとなるように前記2つの合成開口レーダ画像上の画素の位相を補正するバイアス除去部を含む、
ことを特徴とする請求項1または請求項2に記載の合成開口レーダ信号処理装置。 - 前記抽出処理部は、前記第1及び前記第2のベースライン長を有する2つのセンサの軌道情報を用いて、前記第1及び前記第2のベースライン長のそれぞれに対応した軌道パラメータを算出する軌道パラメータ算出部を備え、
前記位相算出部は、前記軌道パラメータ算出部で算出された軌道パラメータを用いて前記第1及び前記第2の特定位相を算出する
ことを特徴とする請求項1乃至3のいずれか1項に記載の合成開口レーダ信号処理装置。 - 前記干渉位相処理部は、
前記第1及び前記第2の2つの合成開口レーダ画像から軌道縞周期を算出する軌道縞周 期算出部を備え、
前記抽出処理部は、
前記軌道縞周期算出部で算出された軌道縞周期と前記合成開口レーダ画像を生成するた めに合成開口レーダから照射された電波のオフナディア角とを用いて、前記第1及び前 記第2のベースライン長のそれぞれに対応した軌道パラメータを算出する軌道パラメー タ算出部を備え、
前記位相算出部は、該軌道パラメータ算出部で算出された軌道パラメータを用いて前記 第1及び前記第2の特定位相を算出する
ことを特徴とする請求項1乃至3のいずれか1項に記載の合成開口レーダ信号処理装置。 - 前記軌道縞周期算出部は、前記第1または前記第2の2つの合成開口レーダ画像上の2つの画素に含まれる信号の相対位相の空間的な変化の度合いを表す周波数スペクトルの中から周波数成分の電力分布に基づき周波数を選択し、該選択された周波数を前記軌道縞周期の周波数として抽出することを特徴とする請求項5に記載の合成開口レーダ信号処理装置。
- 前記抽出処理部は、高さが既知の散乱体を表す少なくとも1つの画素を選択し、前記選択された画素に含まれる前記高さが既知の散乱体と同じ高さの散乱体を抽出する
ことを特徴とする請求項3に記載の合成開口レーダ信号処理装置。 - 前記第1または前記第2の2つの合成開口レーダ画像上の同じ散乱体を表す2つの画素に含まれる位相の差の時間または空間的なばらつきに基づき、前記2つの画素が1種類の反射信号を含むか、複数の種類の反射信号を含むか、
を示すことを特徴とする請求項1乃至3のいずれか1項に記載の合成開口レーダ信号処理装置。 - 前記特定の高さは複数の高さであり、前記画素抽出部は前記複数の高さのそれぞれに相当する画素を抽出し、
前記画素抽出部で前記複数の高さの高さごとに抽出された画素を用いて三次元画像を生成する信号合成処理部
を備えたことを特徴とする請求項1乃至8のいずれか1項に記載の合成開口レーダ信号処理装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2015/000549 WO2016125206A1 (ja) | 2015-02-06 | 2015-02-06 | 合成開口レーダ信号処理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2016125206A1 true JPWO2016125206A1 (ja) | 2017-07-27 |
JP6320577B2 JP6320577B2 (ja) | 2018-05-09 |
Family
ID=56563573
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016572940A Active JP6320577B2 (ja) | 2015-02-06 | 2015-02-06 | 合成開口レーダ信号処理装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20180011187A1 (ja) |
JP (1) | JP6320577B2 (ja) |
DE (1) | DE112015006123T5 (ja) |
WO (1) | WO2016125206A1 (ja) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018056129A1 (ja) * | 2016-09-20 | 2018-03-29 | 日本電気株式会社 | 情報処理装置、情報処理方法、および記憶媒体 |
WO2018123748A1 (ja) * | 2016-12-27 | 2018-07-05 | 日本電気株式会社 | 画像解析装置、画像解析方法及びコンピュータ読み取り可能記録媒体 |
WO2018198212A1 (ja) * | 2017-04-26 | 2018-11-01 | 日本電気株式会社 | 情報処理装置、情報処理方法、およびコンピュータ読み取り可能な記憶媒体 |
US20200166626A1 (en) * | 2017-05-17 | 2020-05-28 | Nec Corporation | Information processing device, information processing method, and storage medium having program stored thereon |
US11029403B2 (en) * | 2017-12-18 | 2021-06-08 | The Government Of The United States Of America, As Represented By The Secretary Of The Navy | Millimeter-wave airborne radar for 3-Dimensional imaging of moving and stationary targets |
WO2019187061A1 (ja) * | 2018-03-30 | 2019-10-03 | 三菱電機株式会社 | レーダ画像処理装置及びレーダ画像処理方法 |
JP7006781B2 (ja) * | 2018-05-16 | 2022-01-24 | 日本電気株式会社 | 合成開口レーダ信号解析装置、合成開口レーダ信号解析方法および合成開口レーダ信号解析プログラム |
CN108627835B (zh) * | 2018-06-29 | 2021-07-27 | 中国科学院电子学研究所 | 全极化差分sar层析的目标重构方法 |
WO2020196308A1 (ja) * | 2019-03-28 | 2020-10-01 | 日本電気株式会社 | 画像処理装置、画像処理方法及び画像処理用コンピュータプログラム |
CN110554377B (zh) * | 2019-09-05 | 2021-04-09 | 中国科学院电子学研究所 | 基于多普勒中心偏移的单通道sar二维流场反演方法及系统 |
WO2021186557A1 (ja) * | 2020-03-17 | 2021-09-23 | 日本電気株式会社 | データ処理装置およびデータ処理方法 |
JP2024034579A (ja) * | 2022-09-01 | 2024-03-13 | 日本電気株式会社 | 信号処理装置および信号処理方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004101299A (ja) * | 2002-09-06 | 2004-04-02 | Mitsubishi Electric Corp | 画像レーダ装置 |
US20080174474A1 (en) * | 2006-11-01 | 2008-07-24 | Imsar | Interferometric switched beam radar apparatus and method |
JP2008249551A (ja) * | 2007-03-30 | 2008-10-16 | Mitsubishi Electric Corp | 情報処理装置及び情報処理方法及びプログラム |
JP2012533051A (ja) * | 2009-07-08 | 2012-12-20 | テレ−リレヴァメント エウローパ−ティ.エルレ.エ. エッセ.エルレ.エルレ. | 同一領域上で取得されたsar画像から得られるインターフェログラムのフィルタリング処理方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4975704A (en) * | 1990-01-26 | 1990-12-04 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Method for detecting surface motions and mapping small terrestrial or planetary surface deformations with synthetic aperture radar |
US5170171A (en) * | 1991-09-19 | 1992-12-08 | Environmental Research Institute Of Michigan | Three dimensional interferometric synthetic aperture radar terrain mapping employing altitude measurement |
US5677693A (en) * | 1994-05-26 | 1997-10-14 | Hughes Aircraft Company | Multi-pass and multi-channel interferometric synthetic aperture radars |
US5608405A (en) * | 1995-10-06 | 1997-03-04 | Lockheed Martin Corporation | Method of generating visual representation of terrain height from SAR data employing multigrid analysis |
US5659318A (en) * | 1996-05-31 | 1997-08-19 | California Institute Of Technology | Interferometric SAR processor for elevation |
US5923278A (en) * | 1996-07-11 | 1999-07-13 | Science Applications International Corporation | Global phase unwrapping of interferograms |
US6011505A (en) * | 1996-07-11 | 2000-01-04 | Science Applications International Corporation | Terrain elevation measurement by interferometric synthetic aperture radar (IFSAR) |
US6046695A (en) * | 1996-07-11 | 2000-04-04 | Science Application International Corporation | Phase gradient auto-focus for SAR images |
US5726656A (en) * | 1996-12-19 | 1998-03-10 | Hughes Electronics | Atmospheric correction method for interferometric synthetic array radar systems operating at long range |
DE19706158C1 (de) * | 1997-02-17 | 1998-05-20 | Deutsch Zentr Luft & Raumfahrt | Verfahren zur Korrektur einer objektabhängigen Spektralverschiebung in Radar-Interferogrammen |
US6911932B1 (en) * | 2003-10-28 | 2005-06-28 | Sandia Corporation | Two antenna, two pass interferometric synthetic aperture radar |
-
2015
- 2015-02-06 WO PCT/JP2015/000549 patent/WO2016125206A1/ja active Application Filing
- 2015-02-06 DE DE112015006123.1T patent/DE112015006123T5/de not_active Withdrawn
- 2015-02-06 US US15/544,278 patent/US20180011187A1/en not_active Abandoned
- 2015-02-06 JP JP2016572940A patent/JP6320577B2/ja active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004101299A (ja) * | 2002-09-06 | 2004-04-02 | Mitsubishi Electric Corp | 画像レーダ装置 |
US20080174474A1 (en) * | 2006-11-01 | 2008-07-24 | Imsar | Interferometric switched beam radar apparatus and method |
JP2008249551A (ja) * | 2007-03-30 | 2008-10-16 | Mitsubishi Electric Corp | 情報処理装置及び情報処理方法及びプログラム |
JP2012533051A (ja) * | 2009-07-08 | 2012-12-20 | テレ−リレヴァメント エウローパ−ティ.エルレ.エ. エッセ.エルレ.エルレ. | 同一領域上で取得されたsar画像から得られるインターフェログラムのフィルタリング処理方法 |
Non-Patent Citations (1)
Title |
---|
C. MAGNARD ET AL.: ""High Resolution Millimeter Wave SAR Interferometry"", GEOSCIENCE AND REMOTE SENSING SYMPOSIUM, 2007. IGARSS 2007. IEEE INTERNATIONAL, JPN6015010062, 23 July 2007 (2007-07-23), US, pages 5061 - 5064, ISSN: 0003750398 * |
Also Published As
Publication number | Publication date |
---|---|
US20180011187A1 (en) | 2018-01-11 |
DE112015006123T5 (de) | 2017-10-26 |
WO2016125206A1 (ja) | 2016-08-11 |
JP6320577B2 (ja) | 2018-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6320577B2 (ja) | 合成開口レーダ信号処理装置 | |
US9417323B2 (en) | SAR point cloud generation system | |
Scheiber et al. | Coregistration of interferometric SAR images using spectral diversity | |
US9523766B2 (en) | Phase error correction in synthetic aperture imaging | |
JP6333196B2 (ja) | エリアの3d画像を生成する方法およびエリアの3d画像を生成するシステム | |
Liang et al. | Measuring azimuth deformation with L-band ALOS-2 ScanSAR interferometry | |
US20150369913A1 (en) | Method and apparatus for correcting ionic distortion of satellite radar interferogram | |
Fornaro et al. | SAR interferometry and tomography: Theory and applications | |
CN102662171A (zh) | 一种sar层析三维成像方法 | |
Ulander et al. | Ultra-wideband SAR interferometry | |
Magnard et al. | Analysis of a maximum likelihood phase estimation method for airborne multibaseline SAR interferometry | |
CN104007439A (zh) | 一种干涉圆迹sar高程估计处理方法 | |
RU2661941C1 (ru) | Способ построения радиолокационного изображения с помощью радиолокационной станции с синтезированной апертурой | |
KR100441590B1 (ko) | 간섭측정용 합성 개구 레이다의 기하학적 특성을 이용하여지형고도를 측정하기 위한 디지털 고도모형 생성방법 | |
US11474236B2 (en) | Radar image processing device and radar image processing method | |
Liang et al. | ScanSAR interferometric processing using existing standard InSAR software for measuring large scale land deformation | |
WO2015151134A1 (ja) | レーダ信号処理装置 | |
WO2000054006A2 (en) | Single-pass interferometric synthetic aperture radar | |
Li et al. | An advanced DSS-SAR InSAR terrain height estimation approach based on baseline decoupling | |
CA2864501A1 (en) | Phase error correction in synthetic aperture imaging | |
Klochko | Algorithms of 3D radio-wave imaging in airborne Doppler radar | |
Meglio et al. | Three dimensional SAR image focusing from non-uniform samples | |
Marston | A correlation-based autofocus algorithm for coherent circular synthetic aperture sonar | |
Tan et al. | Synthetic aperture radar tomography sampling criteria and three-dimensional range migration algorithm with elevation digital spotlighting | |
Zhang et al. | Very high-resolution bistatic SAR imaging with TerraSAR-X as the illuminator in ST mode |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20170410 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20180306 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20180403 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6320577 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |