JPH02231591A - 合成開口レーダ画像の歪補正方法 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 【産業上の利用分野】

本発明は合成開口レーダ画像処理システムにおける高精
度な画像形状歪補正方法に関する。

【従来の技術１ 衛星または航空機に搭載された合成開口レーダ（以下Ｓ
ＡＲと略す）は全天候で分解能な地表観測画像を得るた
めのセンサである。センサの観測データはホログラムで
あり、一般に地上の処理システムでディジタル処理によ
りホログラムから画像を再生する．この再生処理の詳細
は例えば文献、ヒタチ・レビュー，３３　（１９８４年
）第１６１頁から１６６頁（Ｈｉｔａｃｈｉ　Ｒｅｖｉ
ｅｗ，　３　３（１９８４）ＰＰ．１６１−１６６）に
おいて論じられているが概略は次のとおりである．第２
図に示すごとく、マイクロ波の地表各点からの反射がア
ジマス方向（センサ進行方向）及びレンジ方向（進行方
向に交差するレーダ照射方向）に拡がったものであると
ころのホログラム（２０）を各方向に圧縮する．まずレ
ンジ圧縮（２１）Ｌ，、得られたデータ（２２）に対し
、信号帯域（ルックと言う）に区分しそれぞれについて
アジマス圧縮（２３）する．その結果のデータ（２４）
の足し合わせすなわちルック加算（２５）の結果として
地表画像（２６）が再生される． ところが再生画像には （ｉ）センサを搭載した飛翔体の軌道，姿勢のゆらぎ，
センサが斜め側方を観測していること、地表の曲率，地
球の自転に起伏する低周波的な形状歪 （ｎ）地表各点の標高に起因する倒れ込みである高周波
な形状歪 があるため、画像を地図と対応づけて利用するためには
この形状歪を補正し，歪のない画像に変換する必要があ
る．このための従来技術は、生産研究３４巻９号（１９
８２年９月）の前田らによる論文，ＳＥＡＳＡＴ　ＳＡ
Ｒ画像の地形歪補正、において論じられている．その内
容は概略次のとおりである． （１）補正画像から未補正画像八の写像の決定第３図に
示すごとく，補正画像３１の座標を未補正画像３０の座
標に対応づける写像曹３４を求める．このとき地表の起
伏は無視する。未補正画像の座標系３２をレンジ，アジ
マスに対応して（Ｒ，Ａ）　、補正画像の座標系３３を
所定の地図投影法に対応して（ｕ，ｖ）とすると、この
写像▼は形式的に Ａ　　　　　　　　　　　Ｖ と表わされる．上記従来技術ではこの写像をＵ，■の多
項式とし，その係数を未補正画像及び地図上から読みと
った基準点（Ｇｅｏｄｅｆｉｃ　ＣｏｎｔｒｏｌＰａｉ
ｎｔ，　ＧＣＰと略す》３５の座標から最小２乗法によ
り求めている． （２）標高による補正 与えられた補正画像の画素に対応する地図座標（ｕ，ｖ
）から，数値地形モデル（　Ｄ　ｉｇｉｔａｌＴｅｒｒ
ａｉｎ　Ｍｏｄｅｌ，　ＤＴＭと略す）を用いて該点の
標高ｈを求める．一方飛翔体の軌道データから得られる
該画素撮影時刻の飛翔体すなわちセンサ位置と上記地表
点の３次元座標（ｕ，ｖ，ｈ）から、該点とセンサとの
距離（レンジ）を求めこれをＲ′とする。この結果、座
標を考慮したときの写像は（１）式で求めたＲｔｔＲ’
におきかえたものとなる。 （３）リサンプリング 補正画像の各画素に対し（１），（２）により対応する
未補正画像の座標（Ｒ’　，Ａ）が得られる．第４図に
示すごとくこの座標点４１は、未補注画像の画素４０と
位置が一致しない．そこで点４１の濃淡強度を周囲の点
４０の濃淡強度から補関して求める．これをリサンプリ
ングという．補間にはニアレストネイバ法，キュービッ
クコンポルーシ目ン法等が知られている． 以上により形状歪の補正された画像が得られる．【発明
が解決しようとする課題】 上記従来技術には次の問題があった．すなわち（１）Ｇ
ＣＰの数が（１）式の多項式モデルの次数から定まる所
定の箇数より少ない時には写像が求まらず，歪の補正が
できない。 （２）写像の精度すなわち歪補正の精度がＧＣＰの数，
配置に依存し、かつＧＣＰ以外の場所での精度が予測で
きない。 （３）ＧＣＰの標高がほぼ一定でないと写像の精度すな
わち歪補正の精度が低下する。 本発明は上記従来技術の問題点を解決し、（ｉ）ＧＣＰ
数，配置等の制約がなく （ｉｉ）画像全体について精度の予測が可能な高精度な
ＳＡＲ画像の形状歪補正方法を提供することにある。

【課題を解決するための手段１ 上記目的を達成するために、歪補正係数作成処理とリサ
ンプリング処理からなる画像補正処理において、合成開
口レーダを搭載した飛翔体の軌道，姿勢データおよび地
球形状モデルから合成開口レーダ再生画像の各点を補正
画像中の点に対応づけるモデルを用いて補正画像中の複
数の代表点において該点に対応する再生画像中の位置を
求める処理と、該位置データから補正画像の各画素に対
応する再生画像の処理を求める歪補正係数を作成する処
理と、該係数を用いて得た補正画像の各画素に対応する
再生画像位置の画像強度を再生画像の画素強度のリサン
プリングにより求める処理を設けたものである． さらに、高精度化のために再生画像中の１つ以上の基準
点の該画像ＰＪｉ櫟と、咳点の地図座標とから飛翔体の
軌道，姿勢データを推定する処理を設けだものである。 さらに、地表の起伏による歪の高精度な補正のために歪
補正係数により与えられる補正画像の各画素に対応する
再生画像中の位置を地表の標高データにより修正し、該
修正位置の画像強度を再生画像の画素強度の２次元リサ
ンプリングにより求める処理を設けたものである。 【作用】 飛翔体の軌道，姿勢データおよび地球形状モデルから、
補正画像から未補正画像への写像を決定するので、ＧＣ
Ｐの数が少なくても、特にＧＣＰがなくても歪の補正が
できる。 さらにＧＣＰが与えられた場合には、ＧＣＰを用いて飛
翔体の軌道，姿勢データを推定するので、その推定精度
から写像すなわち歪補正の精度を、ＧＣＰ以外の場所に
ついても予測できる。この際、各ＧＣＰの標高に差があ
っても飛翔体の軌道，姿勢データを推定できるので，Ｇ
ＣＰの配置に制約がなく適用できる。 以上により．ＧＣＰとＤＴＭを用いて高精度なＳＡＲ画
像の歪補正が可能になる。

【実施例】

以下、本発明の一実施例を図を用いて説明する。 第１＠は本発明によるＳＡＲ画像の歪補正処理の全体フ
ロー図である。まず飛翔体の軌道・姿勢パラメータを推
定する。（処理１０）このパラメータを用いて補正画像
から未補正の再生画像への写像Ｗを決定する。（処理１
１）次にこの写像を用いて歪補正係数を決定する。（処
理１２）最後にこの係数を用いて再生画像１４に対しリ
サンプリング処理１３をほどこし、所望の補正画像１５
を得る。以下各処理の詳細を示す。なお処理１０におい
て再生画像１４およびＧＣＰデータ１６を用いる高精度
化の方法について第２の実施例で、また処理１７におい
てＤＴＭ１７を用いる高精度化の方法について第３の実
施例でそれぞれ示す。 処理１０は第５図に示すごとく、軌道パラメータの推定
処理５１と姿勢パラメータの推定処理５３から成る。処
理５１は飛翔体の追跡観測あるいは飛翔体搭載センサか
ら得られる時系列の軌道データ５０を入力とし、それを
時刻の多項式近似したときの係数を求める。すなわち時
刻ｔにおける飛翔体位置を地球固定座標系（ｘ（ｔＬ　
ｙ（ｔ），Ｚ（ｔ））または緯度，経度，高度系で（ｌ
ａｔ　（ｔ）　＝１ｏｎ（ｔ），　ｈ（ｔ））とすると
き各成分をｔの多項式で表わす。例えばｘ　（ｔ）を ？（ｔ）＝ｘ，＋ｘ■ｔ＋ｘ２ｔ”＋−＝＋Ｘｎｔ”＋
ｘ（ｔ）（１）式 で表わす。多項式の次数には３ないし４とすればよい。 またｘ　（ｔ）は多項式近似の残差で高次あるいは高周
波成分である。（１）式のパラメータｘ０，　ｘエ，・
・・ＩＸｎは次の最小２乗法により求める６すなわち評
価 ｆ　｛χ（ｔ）−（ｘ．＋ｘ１ｔ　＋ｘ，　ｔ”＋−・
・・＋　ｘｎ　ｔ　’）｝２ｄ　ｔｔ１ （２）式 を最小にするｘ０，　ｘｌ，・・・，Ｘｎを求める。こ
こに、ｘ　（ｔ）は軌道データ５ｏにより与えられる値
、ｔエとｔ２は処理対象とするシーンの撮影開始および
終了時刻である．いま軌道データ５０が離散時刻で与え
られている場合には（２）式の積分を次の和の形に代え
ればよい。 Σ（ｘ（ｔ）−（ｘ．＋ｘ１ｔ＋ｘ．ｔ”＋−＝＋ｘｎ
ｔ”））”　　　（２）’式他の軌道要素ｙ　（ｔ）　
ｒ　ｚ（ｔ）またはｌａｔ（ｔ），　ｉｏｎ（ｔ），ｈ
（ｔ）についても伺様の処理を行なう．さらに軌道デー
タ５０で与えられる，飛翔体速度の各成分Ｖ　ｘ（ｔ）
　ｒ　Ｖ　ｙ（ｔ）　，　ｖ　ｚ（ｔ）についても同様
に多項式近似係数を求める．処理５３は飛翔体搭載セン
サから得られる時系列の姿勢データ５２を入力とし、そ
れを時刻の多項式近似したときの係数を求める．すなわ
ち時刻ｔにおける飛翔体の姿勢角を軌道座標系あるいは
地球固定座標系で（ｒｏｌｌ（ｔ），ｐｉｔｃｈ（ｔ）
，　ｙａｗ（ｔ））とするとき各成分をｔの多項式で表
わす．例えば、ｒｏｌｌ（ｔ）をｒｏｌｌ（ｔ）＝　ｒ
，＋　ｒ１ｔ　＋　ｒ２ｔ２＋・・・・・・＋ｒｎｔｎ
＋ｒ（ｔ）　　　（３）式で表わす。多項式の次数ｎは
３ないし４とすればよい．またｒ　（ｔ）は多項式近似
の残差で高次および高周波成分である。（３）式のパラ
メータは軌道パラメータと同様最小２乗法により求める
。 すなわち評価 ｆ　　（ｒｏｌｌ（ｔ）−（ｒｌｌ＋ｒｉｔ＋ｒ．ｔ″
＋−−＋　ｒｎ　ｔｆｌ））”　ｄ　ｔ”　　　　　　
　　　　　　　　　　　　　　　　（４）式あるいは Σ（ｒｏｌｌ（ｔ）−（ｒ，，＋ｒ１ｔ＋ｒ，ｔ”＋−
＋ｒｎｔ″）｝”　　（４）’式を最小にするｒｌｌｌ
　ｒエ，ｒ２，・・・Ｈ　　ｒｎを求める．他の姿勢角
成分ｐｉｔｃｈ（ｔ），　ｙａｗ（ｔ）についても同様
である． 処理１１は第６図に示すごとく補正画像６０を分割する
格子の各点６１について，未補正の再生画像６２中の対
応する点６３への写像！６４を求める．この写像▼は直
接得られないので再生画像から補正画像への写像Ｔから
第７図に示す繰返し演算により求める。すなわち未補正
画像の座標系６６における点位置（Ｒ，Ａ）がら幾何モ
デルと地図モデルによる写像Ｔにより補正画像の座標系
６７における座ｍｌ（ｕ，ｖ）を求メル．　　＜処ｙｉ
ｉ７０）これと所望の格子点６１の座標（　ｕ．，　ｖ
，）との誤差を次のように計算する．（処理７１）この
誤差を用いて（Ｒ，Ａ）を（Ｒ＋ΔＲ，Ａ＋ΔＡ）に修
正する．（処理７２）修正量はにより得られる．微係数
は幾何，地図モデルから解析的に、または数値的に求め
る−（ｕｔｖ）の誤差が十分小さくなるまで処理７０〜
７２を繰返す。第８図は写像Ｔの内容を示すものである
．一未補正画像の座ＩｊＩ（Ｒ，Ａ）８０を幾何モデル
８１により地表点の地球固定系の座標（ｘｔ　ｙ＋　ｚ
）８２に変換する．次に地図モデル８３により、（Ｘｙ
３’ｙＺ）を緯度，経度に変換したのち、メルカトール
法等，所定の投影法に従う地図座標（ｕ，ｖ）８４に変
換する．第９図は幾何モデル８１を説明する図である．
まず処理１０で求めた軌道パラメータを用いて飛翔体の
位置と速度ベクトルを求める．着目点のアジマス（Ａ）
座標から時刻ｔが得られる．このｔと軌道パラメータか
ら例えば飛翔体位置のＸ成分を ｘ＝ｘ，＋ｘ，ｔ＋ｘ，ｔ”＋・・・・・・＋ｘｎｔ．
”　　　　　（７）式により求める＊　ｙ，ｚ成分ある
いはｌａｔ，　ｉｏｎ，　ｈ成分も同様に求める。その
結果ｔにおける飛翔体すなわちセンサ位ｆｌｒｓ＝　（
ｘ　　ｙ　　ｚ）”９０を得る。同様に軌道パラメータ
から飛翔体速度のＸ成分を ｌｖ　”　（Ｖｘ　Ｖｙ　ｖｚ）”／　（ｖｘ”＋ｖｙ
”＋ｖｚ”）　　　　　（８）式と得られる。所望の点
ｒ＝　（ｘ　　ｙ　　ｚ）”は、孤９３上にあることか
ら （ｒ　　　ｒｓ）”　（ｒ　−　ｒｓ）　　＝Ｒ２（ｒ
−ｒｓ）ｎ＝０ （１１）式 （１２）式 として得る．ｔにおけるレーダビーム面９２上でセンサ
からの距離がＲである円孤９３と地表面９４との交点ｒ
９５が所望の地表点である。処理１０で求めた姿勢パラ
メータを用いて、軌道と同様にｔにおける姿勢角をｒ０
１１についてを満たす。一方ここでは地表の起伏を無視
し、他表が離心率ｅ，赤道半径Ｒｅの回転楕円体である
と近似すれば、地球固定系９７でｒの座標はｒｏｌｌ　
＝　ｒ，＋　ｒ，　ｔ　十ｒ，　ｔ２＋，，，，，，＋
　ｒｎｔ’　　　　（９）式により、ｐｉｔｃｈ，　ｙ
ａｗについても同様に求める。この姿勢角とセンサの飛
翔体への取付角から、速度ベクトル１ｖからビーム面の
法線ベクトルユ９６への直交変換行列Ｍを得る。このＭ
を用いて工はｎ　”　Ｍ　１　ｖ （１０）式 を満たす。そこで（１１）〜（１３）式を連立させてｒ
が求まる。これが幾何モデル８１である。 次に処理１２の内容を説明する。第１０図に補正画像中
の格子１００の一部である４角形のブロック１０１を示
す。このブロック内の補正画像座標（ｕ，ｖ）から未補
正画像座標（Ｒ，Ａ）への写像▼を ？近似する。この係数Ｒ。，・・・，Ｒ３、Ａ６，・・
・，Ａ３は次のように求まる。すなわち注目ブロックの
四隅■〜■のｕ，ｖ座標と対応するＡ座標からＡ０，・
・・Ａ■を得る。また■〜■のＵ座標と対応するＡ，Ｒ
座標からＲ０，・・・，Ｒ３を得る。以上の処理を補正
画像の全ブロックに対して行なうことにより得られる係
数を歪補正係数と呼ぶ。（１４）式の陪線形モデルはブ
ロックの境界で連続であるという特長がある。 最後にリサンプリング処理１３の内容を第１１図により
説明する。縦リサンプリング処理１１２は座標系（Ｒ，
Ａ）１１０にある未補正の再生画像１１１に対し縦すな
わちＲ方向のリサンプリングを行ない座標系（ｕ，Ａ）
１１３にある中間画像１１４を生成する。すなわち画像
１１４の各画素座標に対し（１４）の第１式を用いて対
応するＲ座標を求め、該座標における画像１１１の強度
を補間により求める。補間では第１２図に示すごとく注
目する画素位置Ｒ（１２０）の両側に隣接する複数の画
素座標１１１・・・ｌ　ｘ，における強度ｆｉ（１”ｌ
ｉ＋・・・，ｉ，）を適当な曲線または直線１２１で補
間して、所望の強度ｆ１２２を求める。一般には ｆ＝ΣＬ−ｗ　（ｉ　−Ｒ）　　　　　　（１５）式の
積和演算となる。重みＷの形によりニアレストネイバ法
，キュービックコンポルーション法が知られている。例
えば文献，土屋等による“全ディジタル方式による高精
度地球観測画像情報処理技術″，日立評論６２．３　（
１９８０年）を参照されたい。次に横リサンプリング処
理１１５で同様は（１４）の第２式を用いて中間画像１
１４から座標系（ｕ，ｖ）１１６にある補正画像１１７
を得る．リサンプリングにおいてブロックに応じた歪補
正係数を用いることは言うまでもない．以上によりＳＡ
Ｒ画像の歪補正が行なえる。 次に第２の実施例としてＧＣＰを用いて歪補正を高精度
化する方法を述べる．第１３図に示すごとく、第５図に
示した方法で得た飛翔体の軌道，姿勢パラメータを高精
度化することにより歪補正の精度を高めるものである．
ａｃｐは地図上の位置すなわち緯，経，高度あるいは（
ｘｔ　ｙｅ　ｚ）が正確に知られており，かつ画像上で
検出可能な特徴点であり、その地図座標はファイルｌ７
に蓄えておく．まず再生画像１４をディスプレイに表示
し、各ＧＣＰの画像座標（Ｒ，Ａ）を読みとる．（１３
０）次にこの座標とファイル１７中の地図座標を用いて
パラメータの精度推定を行なう．（１　３　１）いま推
定対象とする飛翔体の位置，姿勢の多項式係数の集まり
をベクトル五と表わす．第ｉＧＣＰに対し幾何モデル８
ｌを用いて（Ｒ，Ａ》から予測した地表点座標を工ｔ　
（１）と表わす．一方ファイル１７で与えられる該ＧＣ
Ｐの地図座標を−Ｌｌとする．このとき評価 Ｊ＝（ｇ−ｇｏ）”ｗ．（ｇ−＆ｌ）＋Σ（ｙｔ　ｉｔ
　（ｇ））”ｗｓ　（ｙｉ　ｆｉ　（ｇ））（ｌ６）式 を最小にする工が所望のパラメータ値である．ここに１
０は処理５１．５２で得たパラメータ値、Ｗ．はその重
みマトリックス、ｗ１は２１の誤差共分散の逆行列に相
当する重みマトリックスである．この最小２乗法はニュ
ートンラフソン法またはカルマンフィルタにより解くこ
とができる．具体的には先願、特願昭５８−２２７５０
を参照されたい．なお幾何モデル８ｌで用いる（ｌ３）
式のＲｅは、ここでは地球表面の起伏を考慮し、各ＧＣ
Ｐごとに各点の地図座標（ｘ’，ｙ’ｚ’）を用いて で求める． 次に第３の実施例としてＤＴＭを用いて歪補正を高精度
化する方法を述べる．リサンプリング処理１３は第１４
図に示すものに置き換える．すなわち以下の処理を補正
画像の各画素に対し行なう。 まず与えられたＤＴＭｌ７の補間１４０により該画素に
対応する地表点の標高を求める．この標高と軌道パラメ
ータ１４１と歪補正係数１４３から該点観測時のレンジ
距離Ｒ′を算出する．（１４２）このＲ′と歪補正係数
１４３を用いて再生画像１４に対して２次元リサンプリ
ング１４４を施こして所望の補正画像１５を得る．補間
１４０は標高の２次元配列のリサンプリングであり、（
１４）と同等の式を用いて前述の処理１３と同様２段階
の１次元リサンプリングとして実現できる．処理１４２
は第１５図に示すごとく、歪補正係数Ａｓｓ・・・，Ａ
３を用いて（１４）の第２式より対応する再生画像のア
ジマスＡを求める．このＡより撮影時刻ｔを得、軌道パ
ラメータを用いて（７）式より飛翔体位置ｒｓを得る．
一方補正画像の画素（ｕｐ　ｖ）に対応する回転楕円体
１５１上の地表点Ｌ１５２の座標と上記ＤＴＭの補間で
求めた標高とから地表の起伏を考慮した二の直上にある
地表点ｒ’　１５３の座標を得る．ｒａとＬ′からレン
ジ距離Ｒ′を により得る．処理１４４では、第１６図に示すごとく上
述のＲとＡが所望の画素に対応する未補正画像の画素点
１６０であるから，この点の強度を周囲の画素点１６１
の強度から補関する．まず周、囲の各レンジラインにつ
いてアジマス方向に１次元リサンプリングを行ない、該
Ａ位置１６２の強度を求める．次に二九をレンジ方向に
１次元リサンプリングし，所望の点１６０の強度を得る
．ここで起伏を考慮したため，歪の局所的線形性が失な
われたため，第１１図に示すごとき２段階の１次元リサ
ンプリングを適用することはできない．以上により起伏
を考慮した高精度歪補正が可能となる． なお、第２及び第３の実施例を組合せて高精度化を図る
ことが可能なことはもちろんである。 ［発明の効果】 本発明によれば、 （１）飛翔体の軌道，姿勢データおよび地球形状モデル
から、補正画像から未補正画像への写像を決定できるの
で、ＧＣＦの数が少なくても、特にＧＣＰがなくても歪
の補正ができるという効果がある． （２）ＧＣＰを用いて飛翔体の軌道，姿勢データを推定
するので，その推定精度から写像すなわち歪補正の精度
を、ＧＣＰ以外の場所についても予測できるという効果
がある。 （３）各ＧＣＰの標高に差があっても飛翔体の軌道，姿
勢データを推定できるので、ＧＣＰの配置に制約がなく
適用できるという効果がある． （４）ＧＣＰとＤＴＭを用いて、高精度なＳＡＲ画像の
歪補正が，可能になるという効果がある。

【図面の簡単な説明】

第１図は本発明の一実施例の合成開口レーダ画像の歪補
正方法の処理フロー図、第２図は合成開口レーダの画像
再生処理フロー図、第３図は補正画像と未補正画像の間
の写像を示す図、第４図は画像のリサンプリングを説明
するための図、第５図は第１図のフローを構成する軌道
・姿勢パラメータ推定処理のフロー図、第６図は第１図
のフローを構成する写像の決定処理を説明するための図
、第７図は同じく写像の決定処理のフロー図、第８図は
写像の内容を示す図、第９図は幾何モデルを説明するた
めの図、第１０図は画像のブロックを示す図、第１１図
は２段階リサンプリングを示す図、第１２（ｉ！は画像
強度の補間を示す図、第１３図は第２の実施例における
軌道・姿勢パラメータ推定処理のフロー図、第１４図は
第３の実施例におけるリサンプリング処理のフロー図，
第１５図は該処理のレンジ距離計算を説明するための図
、第１６図は同じく２次元リサンプリングを説明するた
めの図である． 第２図 案，１５′目 第乙区 第７目 第７ρ目 第　７２　図 茅／．？目

Claims (1)

1. 【特許請求の範囲】 １、歪補正係数作成処理とリサンプリング処理からなる
   画像補正処理方法において、合成開口レーダを搭載した
   飛翔体の軌道、姿勢データおよび地球形状モデルから合
   成開口レーダ再生画像の各点を補正画像中の点に対応づ
   けるモデルを用いて補正画像中の複数の代表点において
   該点に対する再生画像中の位置を求める処理と、該位置
   データから補正画像の各画素に対応する再生画像の位置
   を求める歪補正係数を作成する処理と、該係数を用いて
   得た補正画像の各画素に対応する再生画像位置の画像強
   度を再生画像の画素強度のリサンプリングにより求める
   処理を設けたことを特徴とする合成開口レーダ画像の歪
   補正方法。 ２、上記再生画像中の１つ以上の基準点の該画像座標と
   、該点の地図座標とから飛翔体の軌道、姿勢データを推
   定することを特徴とする請求項１項記載の合成開口レー
   ダ画像の歪補正方法。 ３、上記歪補正係数として縦方向及び横方向の２種類を
   作成し、縦方向および横方向の１次元リサンプリングを
   順次行なうことを特徴とする請求項１項記載の合成開口
   レーダ画像の歪補正方法。 ４、上記歪補正係数により与えられる補正画像の各画素
   に対応する再生画像中の位置を地表の標高データにより
   修正し、該修正位置の画像強度を再生画像の画素強度の
   ２次元リサンプリングにより求めることを特徴とする請
   求項１項記載の合成開口レーダ画像の歪補正方法。