Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP5448757B2 - Radar image processing device - Google Patents
[go: Go Back, main page]

JP5448757B2 - Radar image processing device - Google Patents

Radar image processing device Download PDF

Info

Publication number
JP5448757B2
JP5448757B2 JP2009272323A JP2009272323A JP5448757B2 JP 5448757 B2 JP5448757 B2 JP 5448757B2 JP 2009272323 A JP2009272323 A JP 2009272323A JP 2009272323 A JP2009272323 A JP 2009272323A JP 5448757 B2 JP5448757 B2 JP 5448757B2
Authority
JP
Japan
Prior art keywords
image
unit
coherence
shift
processing
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.)
Active
Application number
JP2009272323A
Other languages
Japanese (ja)
Other versions
JP2011112630A (en
Inventor
昇 大石
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2009272323A priority Critical patent/JP5448757B2/en
Publication of JP2011112630A publication Critical patent/JP2011112630A/en
Application granted granted Critical
Publication of JP5448757B2 publication Critical patent/JP5448757B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Description

この発明は、合成開口レーダ、逆合成開口レーダ及び実開口レーダ等に用いられ、これらのレーダの画像処理を行うレーダ画像処理装置に関する。   The present invention relates to a radar image processing apparatus that is used in a synthetic aperture radar, a reverse synthetic aperture radar, a real aperture radar, and the like and performs image processing of these radars.

一般的に、合成開口レーダ、逆合成開口レーダ及び実開口レーダは、観測対象に対して電波を照射し、その反射波を受信して信号処理をすることによって、レーダ画像を得る。これらの方式のレーダでは、観測対象に電波を照射する際に、観測対象範囲外にも電波の一部が照射されるため、観測対象範囲外に存在する物体によって強い反射波が生じる場合がある。このような場合には、観測対象範囲のレーダ画像を再生すると、観測対象範囲内の像として、観測対象範囲外からの強い反射波に基づく像が映し出されるという問題がある。このような観測範囲内に存在しない物体のレーダ画像中に映り込んだ像を虚像という。   In general, a synthetic aperture radar, a reverse synthetic aperture radar, and a real aperture radar obtain a radar image by irradiating an observation target with a radio wave, receiving the reflected wave, and performing signal processing. In these types of radar, when a radio wave is irradiated to the observation target, a part of the radio wave is also irradiated outside the observation target range, so that a strong reflected wave may be generated by an object outside the observation target range. . In such a case, when a radar image in the observation target range is reproduced, there is a problem that an image based on a strong reflected wave from outside the observation target range is displayed as an image in the observation target range. An image reflected in a radar image of an object that does not exist within the observation range is called a virtual image.

ここで、例えば、非特許文献1に示すような従来のレーダ画像の虚像検出方法では、ビーム右側、ビーム中央及びビーム左側のデータ毎に画像を再生して、ビーム右側又はビーム左側の画像のみに写った像を、虚像であるとして判別可能である。この従来の虚像検出方法によって検出可能となる虚像は、レーダのビームの側方から入り込む強い反射波によって生じる虚像であるアジマスアンビギュイティである。なお、アジマスアンビギュイティ以外の虚像検出方法は、一般的に知られておらず、多くの場合において、レーダ画像利用者の目視によって、再生されたレーダ画像から虚像の検出が行われる。   Here, for example, in a conventional radar image virtual image detection method as shown in Non-Patent Document 1, an image is reproduced for each data on the right side of the beam, the center of the beam, and the left side of the beam, and only the image on the right side or left side of the beam is reproduced. The captured image can be identified as a virtual image. The virtual image that can be detected by this conventional virtual image detection method is azimuth ambiguity, which is a virtual image generated by a strong reflected wave entering from the side of the radar beam. Note that a virtual image detection method other than azimuth ambiguity is not generally known, and in many cases, a virtual image is detected from a reproduced radar image by visual observation of a radar image user.

G. Hajduch, P. Leilde, and V. Kerbaol, “Ship detection on ENVISAT ASAR data : results, limitations and perspectives, ” Proceedings of SEASAR 2006.G. Hajduch, P. Leilde, and V. Kerbaol, “Ship detection on ENVISAT ASAR data: results, limitations and perspectives,” Proceedings of SEASAR 2006.

上記のような従来の虚像検出方法は、アジマスアンビギュイティが側方からの反射波によって生じる特徴を利用したものである。このため、上記のような従来の虚像検出方法では、レンジアンビギュイティと呼ばれるレーダのビームの上方又は下方からの強い反射波による虚像や、干渉波と呼ばれる放送、通信又は他のレーダに利用される電波を受信したことによる虚像等を検出することができないという問題があった。従って、上記のような従来の虚像検出方法では、アジマスアンビギュイティ以外の種類の虚像については検出することはできなかった。   The conventional virtual image detection method as described above uses the feature that the azimuth ambiguity is generated by the reflected wave from the side. For this reason, the conventional virtual image detection method as described above is used for a virtual image due to a strong reflected wave from above or below the radar beam called range ambiguity, broadcasting, communication or other radar called an interference wave. There has been a problem that it is impossible to detect a virtual image or the like caused by receiving a radio wave. Therefore, the conventional virtual image detection method as described above cannot detect a virtual image of a type other than azimuth ambiguity.

この発明は、上記のような課題を解決するためになされたものであり、アジマスアンビギュイティ以外の種類の虚像についても検出対象とすることができるレーダ画像処理装置を得ることを目的とする。   The present invention has been made to solve the above-described problems, and an object of the present invention is to obtain a radar image processing apparatus that can also detect virtual images of types other than azimuth ambiguity.

この発明に係るレーダ画像処理装置は、取得したレーダ画像の一部及び全部のいずれか一方である処理対象画像について、分解能を基本単位とする予め設定されたシフト量に基づいて、前記レーダ画像を表すための2軸の少なくともいずれか一方の軸方向へシフトするシフト処理部と、前記シフト処理部によるシフト後の前記処理対象画像と、そのシフト前の前記処理対象画像とから、前記処理対象画像における各画素と、前記各画素の周囲の画素とについてのコヒーレンスの分布であるコヒーレンスマップを算出するコヒーレンスマップ算出部と、前記処理対象画像内に虚像が存在する場合に、前記コヒーレンスマップに基づいて、前記処理対象領域内に存在する虚像の位置を特定する虚像位置特定部とを備えるものである。 The radar image processing apparatus according to the present invention is configured to process the radar image based on a preset shift amount having a resolution as a basic unit with respect to a processing target image that is one of all or part of the acquired radar image. A shift processing unit that shifts in the axial direction of at least one of two axes for representing, the processing target image after the shift by the shift processing unit, and the processing target image before the shift; A coherence map calculation unit that calculates a coherence map that is a coherence distribution for each pixel in the pixel and pixels around each pixel, and when a virtual image exists in the processing target image, based on the coherence map And a virtual image position specifying unit for specifying the position of the virtual image existing in the processing target area .

この発明に係るレーダ画像処理装置によれば、コヒーレンスマップ算出部が、分解能を基本単位とするシフト量でのシフト後の処理対象画像と、シフト前の処理対象画像とから、処理対象画像における各画素と、各画素の周囲の画素とについてのコヒーレンスの分布であるコヒーレンスマップを算出し、虚像位置特定部が、処理対象画像内に虚像が存在する場合に、コヒーレンスマップに基づいて、処理対象領域内に存在する虚像の位置を特定することにより、アジマスアンビギュイティ以外の種類の虚像についても検出することができる。 According to the radar image processing device of the present invention, the coherence map calculation unit calculates each of the processing target images from the processing target image after the shift with the shift amount having the resolution as a basic unit and the processing target image before the shift. A coherence map, which is a distribution of coherence for pixels and surrounding pixels of each pixel, is calculated, and when a virtual image exists in the processing target image, the virtual image position specifying unit performs processing based on the coherence map. By specifying the position of the virtual image existing inside, it is possible to detect a virtual image of a type other than azimuth ambiguity.

この発明の実施の形態1によるレーダ画像処理装置を示すブロック図である。1 is a block diagram illustrating a radar image processing apparatus according to Embodiment 1 of the present invention. レーダ画像の結像状態を説明するための説明図である。It is explanatory drawing for demonstrating the image formation state of a radar image. 図1のレーダ画像処理装置の動作を示すフローチャートである。2 is a flowchart showing an operation of the radar image processing apparatus of FIG. 1. この発明の実施の形態2によるレーダ画像処理装置を示すブロック図である。It is a block diagram which shows the radar image processing apparatus by Embodiment 2 of this invention. 図4のレーダ画像処理装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the radar image processing apparatus of FIG. この発明の実施の形態3によるレーダ画像処理装置を示すブロック図である。It is a block diagram which shows the radar image processing apparatus by Embodiment 3 of this invention. 図6の画像ぼけ解析部の解析処理に用いられる関数の一例を示すグラフである。It is a graph which shows an example of the function used for the analysis process of the image blur analysis part of FIG. 図6のレーダ画像処理装置の動作を示すフローチャートである。It is a flowchart which shows operation | movement of the radar image processing apparatus of FIG. この発明の実施の形態4によるレーダ画像処理装置を示すブロック図である。It is a block diagram which shows the radar image processing apparatus by Embodiment 4 of this invention. 図9の分解能推定部の推定処理に用いられる関数の一例を示すグラフである。It is a graph which shows an example of the function used for the estimation process of the resolution estimation part of FIG. 図9のレーダ画像処理装置の動作を示すフローチャートである。10 is a flowchart showing the operation of the radar image processing apparatus of FIG.

以下、この発明を実施するための形態について、図面を参照して説明する。
実施の形態1.
図1は、この発明の実施の形態1によるレーダ画像処理装置を示すブロック図である。
図1において、レーダ画像処理装置100は、観測対象へアンテナ(図示せず)から電波を照射することによって得るレーダ画像を処理するものである。また、レーダ画像処理装置100は、レーダ画像格納部101、パラメータ格納部102、画像シフト処理部103、コヒーレンス算出部104、虚像判断部105及び判断結果格納部106を有している。
Hereinafter, embodiments for carrying out the present invention will be described with reference to the drawings.
Embodiment 1 FIG.
FIG. 1 is a block diagram showing a radar image processing apparatus according to Embodiment 1 of the present invention.
In FIG. 1, a radar image processing apparatus 100 processes a radar image obtained by irradiating an observation target with radio waves from an antenna (not shown). The radar image processing apparatus 100 includes a radar image storage unit 101, a parameter storage unit 102, an image shift processing unit 103, a coherence calculation unit 104, a virtual image determination unit 105, and a determination result storage unit 106.

レーダ画像格納部101には、前段の受信処理装置(図示せず)による受信処理を受けたレーダ画像(SAR画像:Synthetic Aperture Radar)の一部又は全部が処理対象画像として格納される。なお、このレーダ画像格納部101に格納される処理対象画像は、合成開口レーダの画像再生処理後のデータでもよく、その画像再生処理の途中のデータでもよい。   In the radar image storage unit 101, a part or all of a radar image (SAR image: Synthetic Aperture Radar) subjected to reception processing by a reception processing device (not shown) in the previous stage is stored as a processing target image. Note that the processing target image stored in the radar image storage unit 101 may be data after the image reproduction processing of the synthetic aperture radar, or data during the image reproduction processing.

パラメータ格納部102には、処理対象画像をずらす大きさであるシフト量や、処理対象画像をずらす方向等を表すシフトパラメータが予め格納されている。このシフトパラメータのシフト量は、1分解能以上が基本である。画像シフト処理部103は、パラメータ格納部102に格納されたシフトパラメータのシフト量に基づいて、レーダ画像格納部101に格納された処理対象画像をシフト(変位)する。コヒーレンス算出部104は、レーダ画像格納部101に格納されたシフト前の処理対象画像と、画像シフト処理部103によってシフトされたシフト後の処理対象画像とについて、相関の度合いを示すコヒーレンス(相関係数)を算出する。   The parameter storage unit 102 stores in advance a shift parameter that indicates the amount by which the processing target image is shifted, the shift direction of the processing target image, and the like. The shift amount of this shift parameter is basically one resolution or more. The image shift processing unit 103 shifts (displaces) the processing target image stored in the radar image storage unit 101 based on the shift amount of the shift parameter stored in the parameter storage unit 102. The coherence calculation unit 104 has a coherence (correlation) indicating a degree of correlation between the pre-shift processing target image stored in the radar image storage unit 101 and the post-shift processing target image shifted by the image shift processing unit 103. Number).

虚像判断部105は、コヒーレンス算出部104によって算出されたコヒーレンスに基づいて、処理対象画像内にぼけ領域が存在しているか否かを判断する。また、虚像判断部105は、処理対象画像内にぼけ領域が存在していると判断した場合に、処理対象画像内に虚像が存在していると判断する。判断結果格納部106には、虚像判断部105によるぼけ領域の有無、又は虚像の有無についての判断結果が格納される。   Based on the coherence calculated by the coherence calculation unit 104, the virtual image determination unit 105 determines whether or not a blur region exists in the processing target image. Further, when the virtual image determination unit 105 determines that a blurred region exists in the processing target image, the virtual image determination unit 105 determines that a virtual image exists in the processing target image. The determination result storage unit 106 stores a determination result regarding the presence / absence of a blur region or the presence / absence of a virtual image by the virtual image determination unit 105.

次に、レーダ画像の結像状態について説明する。レーダ画像処理装置100の検出対象とする虚像は、レーダの想定する観測対象領域の外に存在する物体が観測対象領域に投影されてなる像(即ちアンビギュイティ)や、放送波、通信波又は他のレーダ装置の送信波や妨害波等の干渉波が基となる像である。また、レーダの想定する観測対象領域とは、レーダで画像化したい領域であり、一般的には、アンテナのメインビームが照射された領域である。さらに、レーダ画像においては、像がぼけている場合を除いて、画像中の各地点の観測対象の情報が1分解能以内に収まるという特徴がある。   Next, the imaging state of the radar image will be described. The virtual image to be detected by the radar image processing apparatus 100 is an image (that is, ambiguity) obtained by projecting an object existing outside the observation target region assumed by the radar onto the observation target region, a broadcast wave, a communication wave, It is an image based on interference waves such as transmission waves and interference waves from other radar devices. The observation target region assumed by the radar is a region to be imaged by the radar, and is generally a region irradiated with the main beam of the antenna. Furthermore, the radar image has a feature that the information of the observation target at each point in the image falls within one resolution except when the image is blurred.

図2は、レーダ画像の結像状態を説明するための説明図である。図2の点A,Bは、ぼけずに結像した像であり、図2の領域Cは、ぼけた像、即ちぼけ領域である。また、図2のレーダ画像は、x軸(例えばレンジ方向:距離方向)及びy軸(例えばアジマス方向:角度方向)の2軸の座標によって表されている。また、図2の格子の幅は、画像のx軸及びy軸のそれぞれの分解能の大きさである。   FIG. 2 is an explanatory diagram for explaining the imaging state of the radar image. Points A and B in FIG. 2 are images formed without being blurred, and a region C in FIG. 2 is a blurred image, that is, a blurred region. Further, the radar image of FIG. 2 is represented by two-axis coordinates of an x axis (for example, a range direction: a distance direction) and a y axis (for example, an azimuth direction: an angle direction). Also, the width of the grid in FIG. 2 is the resolution of each of the x-axis and y-axis of the image.

図2において、レーダ画像中で、互いに1分解能以上離れた地点間でコヒーレンスを計算すると、像がぼけていない場合には、コヒーレンスがほぼ0となり、像がぼけている場合にはコヒーレンスが0よりも大きな値となる。また、レーダ画像の再生処理では、観測対象領域からの反射波のみを画像として結像するように(即ちぼけていない像となるように)処理される。このため、観測領域外の電波が基となるアンビギュイティは、ぼけた像としてレーダ画像中に現れる。   In FIG. 2, when coherence is calculated between points separated by 1 resolution or more from each other in the radar image, the coherence is almost 0 when the image is not blurred, and the coherence is 0 when the image is blurred. Is also a large value. In the radar image reproduction process, only the reflected wave from the observation target region is imaged as an image (that is, an image that is not blurred). For this reason, the ambiguity based on radio waves outside the observation region appears as a blurred image in the radar image.

さらに、レーダの送信波は、一定の方式・周波数で変調されており、画像再生の際には、受信波を復調して処理する。これに対して、虚像の基となる干渉波は、自身の送信波とは異なる方式・周波数で変調されている。このため、干渉波は、自身で復調できないため、干渉波による虚像は、アンビギュイティと同様に、レーダ画像中でぼけるという特性がある。なお、レーダ画像処理装置100の前段の受信処理装置において、観測領域内の自身の反射波についての受信処理が不完全な場合にも、レーダ画像中にぼけ領域が生じる。   Further, the transmission wave of the radar is modulated with a fixed method and frequency, and when the image is reproduced, the received wave is demodulated and processed. On the other hand, the interference wave that is the basis of the virtual image is modulated with a method and frequency different from that of its own transmission wave. For this reason, since the interference wave cannot be demodulated by itself, the virtual image due to the interference wave has a characteristic that it is blurred in the radar image like the ambiguity. Note that, in the reception processing device in the preceding stage of the radar image processing device 100, even when the reception processing for the reflected wave in the observation region is incomplete, a blur region occurs in the radar image.

次に、コヒーレンス算出部104によるコヒーレンスの算出方法と、虚像判断部105による虚像又はぼけ領域の検出方法とについてより具体的に説明する。画像上の座標(x,y)の画素値をp(x,y)とし、x軸及びy軸方向のシフト量をそれぞれΔx及びΔyとした場合には、コヒーレンス算出部104は、次の式(1)を用いた演算を実行することによって、コヒーレンスγ(Δx,Δy)を算出する。   Next, a coherence calculation method by the coherence calculation unit 104 and a virtual image or blur region detection method by the virtual image determination unit 105 will be described more specifically. When the pixel value of the coordinates (x, y) on the image is p (x, y) and the shift amounts in the x-axis and y-axis directions are Δx and Δy, respectively, the coherence calculation unit 104 calculates the following equation: The coherence γ (Δx, Δy) is calculated by executing the calculation using (1).

Figure 0005448757
ただし、Aはコヒーレンスの算出範囲(処理対象範囲)であり、処理対象画像の一部の範囲又は全範囲である。*は複素共役である。
Figure 0005448757
However, A is a coherence calculation range (processing target range), which is a partial range or the entire range of the processing target image. * Is a complex conjugate.

そして、虚像判断部105は、コヒーレンス算出部104によって算出されたコヒーレンスに基づいて、処理対象画像内のぼけ領域の有無を判断する。具体的に、虚像判断部105は、コヒーレンスの値がほぼ0であることを確認した場合には、ぼけ領域が処理対象画像内に存在していないと判断する。他方、虚像判断部105は、コヒーレンスの値が0よりも大きいことを確認した場合には、処理対象画像内にぼけ領域が存在すると判断し、そのぼけ領域が虚像であると判断する。即ち、虚像判断部105は、コヒーレンスの値が0よりも大きいことを確認することによって、ぼけ領域及び虚像を検出する。   Then, the virtual image determination unit 105 determines the presence / absence of a blurred region in the processing target image based on the coherence calculated by the coherence calculation unit 104. Specifically, when the virtual image determination unit 105 confirms that the coherence value is approximately 0, the virtual image determination unit 105 determines that the blurred region does not exist in the processing target image. On the other hand, if the virtual image determination unit 105 confirms that the coherence value is greater than 0, the virtual image determination unit 105 determines that a blur area exists in the processing target image, and determines that the blur area is a virtual image. That is, the virtual image determination unit 105 detects a blurred area and a virtual image by confirming that the coherence value is greater than zero.

従って、レーダ画像処理装置100は、コヒーレンスを用いてレーダ画像内のぼけ領域を検出することによって、アジマスアンビギュイティによる虚像のみならず、レンジアンビギュイティによる虚像や干渉波による虚像についても、検出対象とすることができる。   Therefore, the radar image processing apparatus 100 detects not only a virtual image due to azimuth ambiguity but also a virtual image due to range ambiguity and a virtual image due to interference waves by detecting a blur region in the radar image using coherence. Can be targeted.

ここで、ぼけ領域が処理対象画像内に存在していないと判断する(コヒーレンスがほぼ0であると判断する)基準については、レーダ画像利用者が経験に基づいて定めてもよいし、あるいは予め設定した閾値を用いてもよい。この閾値は、レーダ画像利用者が経験に基づいて設定してもよい。また、シミュレーションにより、虚像又はぼけ領域が存在しない場合のコヒーレンスを計算し、そのコヒーレンスの最大値、平均値又は標準偏差等の統計量を指標として用いて、閾値を設定してもよい。   Here, the reference for determining that the blur region does not exist in the processing target image (determining that the coherence is approximately 0) may be determined by the radar image user based on experience or in advance. A set threshold value may be used. This threshold value may be set by a radar image user based on experience. Further, the coherence when a virtual image or a blurred region does not exist may be calculated by simulation, and the threshold value may be set using a statistic such as the maximum value, average value, or standard deviation of the coherence as an index.

次に、動作について説明する。図3は、図1のレーダ画像処理装置100の動作を示すフローチャートである。まず、ステップS101では、画像シフト処理部103が、処理対象画像を、シフト量に合わせてシフトする。この後、ステップS102では、コヒーレンス算出部104が、シフト前の処理対象画像と、シフト後の処理対象画像とのコヒーレンスを算出する。   Next, the operation will be described. FIG. 3 is a flowchart showing the operation of the radar image processing apparatus 100 of FIG. First, in step S101, the image shift processing unit 103 shifts the processing target image according to the shift amount. Thereafter, in step S102, the coherence calculation unit 104 calculates the coherence between the processing target image before the shift and the processing target image after the shift.

続いて、ステップS103では、虚像判断部105が、コヒーレンス算出部104によって算出されたコヒーレンスの値から、処理対象画像内の虚像又はぼけ領域の有無を判断する。そして、ステップS104では、虚像判断部105が、処理対象画像内の虚像又はぼけ領域の有無についての判断結果を判断結果格納部106に格納(保存)する。これにより、レーダ画像処理装置100の動作が終了する。   Subsequently, in step S <b> 103, the virtual image determination unit 105 determines the presence or absence of a virtual image or a blurred region in the processing target image from the coherence value calculated by the coherence calculation unit 104. In step S <b> 104, the virtual image determination unit 105 stores (saves) the determination result regarding the presence or absence of the virtual image or the blurred area in the processing target image in the determination result storage unit 106. Thereby, the operation of the radar image processing apparatus 100 ends.

上記のような実施の形態1のレーダ画像処理装置100によれば、コヒーレンス算出部104が、分解能を基本単位とするシフト量でのシフト後の処理対象画像と、シフト前の処理対象画像とのコヒーレンスを算出する。そして、虚像判断部105が、コヒーレンスに基づいて処理対象画像内にぼけ領域が存在するか否かを判断し、処理対象画像内にぼけ領域が存在すると判断した場合に、処理対象画像内に虚像が存在すると判断する。この構成により、処理対象画像中に虚像が存在する場合に、その虚像についてのぼけた像が存在する画素は分解能単位での周囲の画素と相関を有することから、アジマスアンビギュイティ以外の種類の虚像、即ちレンジアンビギュイティや干渉波による虚像等についても検出対象とすることができる。   According to the radar image processing apparatus 100 of the first embodiment as described above, the coherence calculation unit 104 performs processing between the processing target image after the shift with the shift amount having the resolution as a basic unit and the processing target image before the shift. Calculate coherence. Then, when the virtual image determination unit 105 determines whether or not a blur region exists in the processing target image based on the coherence, and determines that the blur region exists in the processing target image, the virtual image in the processing target image Is determined to exist. With this configuration, when a virtual image exists in the processing target image, the pixel where the blurred image about the virtual image exists correlates with the surrounding pixels in resolution units, and therefore, other types than azimuth ambiguity A virtual image, that is, a range ambiguity, a virtual image due to an interference wave, or the like can also be a detection target.

また、判断結果格納部106に格納された処理対象画像内の虚像又はぼけ領域の有無についての判断結果は、レーダ画像処理装置100の後段の処理装置(図示せず)において、例えば処理対象画像からの虚像又はぼけ領域の除去処理等に用いることができる。   In addition, a determination result about the presence or absence of a virtual image or a blurred area in the processing target image stored in the determination result storage unit 106 is obtained from, for example, the processing target image in a processing device (not shown) at the subsequent stage of the radar image processing device 100. It can be used for the removal process of the virtual image or the blur area.

ここで、非特許文献1に示すような従来の虚像検出方法では、側方から強い反射波を返す特性を持つ観測対象についても、虚像として誤って検出されるという問題があった。
これに対して、実施の形態1のレーダ画像処理装置100では、虚像判断部105が、処理対象画像内におけるぼけ領域の有無の判断結果に基づいて、処理対象画像内における虚像の有無を判断する。この構成により、側方から強い反射波を返す特性を持つ観測対象について、虚像として誤検出がされることを回避できる。
Here, the conventional virtual image detection method as shown in Non-Patent Document 1 has a problem that an observation target having a characteristic of returning a strong reflected wave from the side is erroneously detected as a virtual image.
On the other hand, in the radar image processing apparatus 100 according to the first embodiment, the virtual image determination unit 105 determines the presence / absence of a virtual image in the processing target image based on the determination result of the presence / absence of the blur region in the processing target image. . With this configuration, it is possible to avoid erroneous detection as a virtual image of an observation target having a characteristic of returning a strong reflected wave from the side.

また、非特許文献1に示すような従来の虚像検出方法では、ビーム右側、ビーム中央及びビーム左側の画像をそれぞれ再生するか、又はビーム右側、ビーム中央及びビーム左側で受信した電波を全て用いて1枚の画像を再生した後に、その再生した1枚の画像をビーム右側、ビーム中央及びビーム左側の画像に分解する処理が必要であった。
これに加えて、非特許文献1に示すような従来の虚像検出方法では、ビーム右側、ビーム中央及びビーム左側の画像をそれぞれ作成した後に、それらを1画素以下の誤差で重ね合わせる必要がある。このため、ビーム右側、ビーム中央及びビーム左側の画像を作成しただけでは、例えば航空機の軌道が直線でない等のために1画素以上のずれが生じる場合があり、この場合にはそれぞれの画像を重ね合わせる処理が別途必要であった。
これに対して、実施の形態1のレーダ画像処理装置100では、画像の分解処理や、重ね合わせる処理が不要であることから、従来の虚像検出方法に比べて、処理負担を軽減させることができ、演算時間を短縮させることができる。
Further, in the conventional virtual image detection method as shown in Non-Patent Document 1, images on the right side of the beam, the center of the beam, and the left side of the beam are respectively reproduced, or all radio waves received on the right side of the beam, the center of the beam, and the left side of the beam are used. After reproducing one image, it is necessary to process the reproduced one image into images on the right side, the center of the beam, and the left side of the beam.
In addition, in the conventional virtual image detection method as shown in Non-Patent Document 1, it is necessary to create images on the right side of the beam, the center of the beam, and the left side of the beam, and then superimpose them with an error of one pixel or less. For this reason, if only the images on the right side of the beam, the center of the beam, and the left side of the beam are created, a deviation of one pixel or more may occur due to, for example, the aircraft trajectory not being a straight line. A separate process was required.
On the other hand, the radar image processing apparatus 100 according to the first embodiment does not require image decomposition processing or superimposition processing, so that the processing load can be reduced as compared with the conventional virtual image detection method. The calculation time can be shortened.

さらに、レーダ画像利用者の目視により虚像を検出する方式では、画像上に虚像が比較的弱い信号強度でしか現れない場合があるため、その虚像については、レーダ画像利用者が虚像を目視により検出することが困難であった。
これに対して、実施の形態1のレーダ画像処理装置100では、信号処理によって虚像を検出するので、目視では困難な比較的弱い信号強度の虚像についても、検出対象とすることができる。
Furthermore, in the method of detecting a virtual image by the visual observation of the radar image user, the virtual image may appear only on the image with a relatively weak signal strength. It was difficult to do.
On the other hand, since the radar image processing apparatus 100 according to the first embodiment detects a virtual image by signal processing, a virtual image having a relatively weak signal intensity that is difficult to visually observe can be detected.

なお、実施の形態1における式(1)のコヒーレンスに代えて、画像間の相関を計算する方法、例えば相関関数や変数変換した後に相関係数を計算する方法等を用いてもよい。   Instead of the coherence of the expression (1) in the first embodiment, a method for calculating the correlation between images, for example, a method for calculating a correlation coefficient after performing a correlation function or variable conversion may be used.

また、実施の形態1の図3におけるステップS102を介さずに、座標(x,y)の画素と座標(x+Δx,y+Δy)の画素とのコヒーレンスを直接的に算出してもよい。ただし、虚像又は画像のぼけ領域を検出するためには、シフト量Δx及びΔyのうち、少なくともいずれか一方が分解能以上の値である必要がある。   Further, the coherence between the pixel at coordinates (x, y) and the pixel at coordinates (x + Δx, y + Δy) may be directly calculated without going through step S102 in FIG. 3 of the first embodiment. However, in order to detect a virtual image or a blurred region of an image, at least one of the shift amounts Δx and Δy needs to be a value equal to or higher than the resolution.

実施の形態2.
実施の形態1のレーダ画像処理装置100では、処理対象画像全体のシフト前後のコヒーレンスを用いることによって、処理対象画像内における虚像の有無を判断できるが、虚像又はぼけ領域の位置を特定することはできなかった。これに対して、実施の形態2のレーダ画像処理装置200では、処理対象画像のコヒーレンスに代えて、処理対象画像の各画素のコヒーレンス、即ち処理対象画像の局所的なコヒーレンスの分布であるコヒーレンスマップを作成し、そのコヒーレンスマップに基づいて、虚像又はぼけ領域の位置を特定可能となる。
Embodiment 2. FIG.
In the radar image processing apparatus 100 according to the first embodiment, the presence or absence of a virtual image in the processing target image can be determined by using coherence before and after the entire processing target image is shifted. However, specifying the position of the virtual image or the blur region is not possible. could not. On the other hand, in the radar image processing apparatus 200 according to the second embodiment, instead of the coherence of the processing target image, the coherence map is a coherence of each pixel of the processing target image, that is, a local coherence distribution of the processing target image. And the position of the virtual image or blurred region can be specified based on the coherence map.

図4は、この発明の実施の形態2によるレーダ画像処理装置を示すブロック図である。図4において、実施の形態2のレーダ画像処理装置200の構成の概要は、実施の形態1のブロック構成と同様であるが、実施の形態1のコヒーレンス算出部104、虚像判断部105及び判断結果格納部106に代えて、図4のコヒーレンスマップ算出部204、虚像位置特定部205及び虚像分布格納部206を有している。   FIG. 4 is a block diagram showing a radar image processing apparatus according to Embodiment 2 of the present invention. 4, the outline of the configuration of the radar image processing apparatus 200 according to the second embodiment is the same as the block configuration according to the first embodiment, but the coherence calculation unit 104, the virtual image determination unit 105, and the determination result according to the first embodiment. Instead of the storage unit 106, the coherence map calculation unit 204, the virtual image position specifying unit 205, and the virtual image distribution storage unit 206 of FIG.

コヒーレンスマップ算出部204は、レーダ画像格納部101に格納された処理対象画像と、画像シフト処理部103によるシフト後の処理対象画像とのコヒーレンスマップを算出(生成)する。虚像位置特定部205は、コヒーレンスマップ算出部204によって算出されたコヒーレンスマップに基づいて、虚像の位置を特定する。虚像分布格納部206には、虚像位置特定部205によって特定した虚像の位置の情報が格納される。他の構成は、実施の形態1と同様である。   The coherence map calculation unit 204 calculates (generates) a coherence map between the processing target image stored in the radar image storage unit 101 and the processing target image shifted by the image shift processing unit 103. The virtual image position specifying unit 205 specifies the position of the virtual image based on the coherence map calculated by the coherence map calculating unit 204. The virtual image distribution storage unit 206 stores information on the position of the virtual image specified by the virtual image position specifying unit 205. Other configurations are the same as those in the first embodiment.

次に、コヒーレンスマップ算出部204によるコヒーレンスマップの算出方法についてより具体的に説明する。コヒーレンスマップは、処理対象画像中の各画素において、その画素の周囲の画素とのコヒーレンスに基づいて作成される局所的なコヒーレンスの分布図である。コヒーレンスマップ算出部204は、演算対象画素の座標を(X,Y)とし、シフト量をΔx及びΔyとした場合における局所的なコヒーレンスγ(X,Y,Δx,Δy)について、次の式(2)を用いた演算を実行することによって算出できる。   Next, the calculation method of the coherence map by the coherence map calculation unit 204 will be described more specifically. The coherence map is a local coherence distribution map created based on the coherence with the surrounding pixels of each pixel in the processing target image. The coherence map calculation unit 204 uses the following equation () for local coherence γ (X, Y, Δx, Δy) when the coordinates of the calculation target pixel are (X, Y) and the shift amounts are Δx and Δy. It can be calculated by executing the calculation using 2).

Figure 0005448757
ただし、Aは座標(X,Y)の周囲の画素の領域を示す。
Figure 0005448757
However, A shows the area | region of the pixel around a coordinate (X, Y).

ここで、座標(X,Y)の周囲の画素の設定方法は様々あり、例えば、演算対象画素を中心とした上下左右及び斜め方向1画素ずつの領域である8近傍領域や、演算対象画素の縦方向の複数画素のみの領域や、演算対象画素の横方向の複数画素のみの領域や、演算対象画素を中心とする長方形状に設定した領域等が考えられる。また、座標(X,Y)の周囲の画素の設定方法は、処理対象画像内で一様である必要はなく、処理対象画像内の位置毎にその設定方法を適宜変更してもよい。   Here, there are various methods for setting the pixels around the coordinates (X, Y). For example, there are eight neighboring regions that are one pixel at a time in the vertical and horizontal directions and the diagonal direction around the pixel to be calculated, A region including only a plurality of pixels in the vertical direction, a region including only a plurality of pixels in the horizontal direction of the calculation target pixel, a region set in a rectangular shape centering on the calculation target pixel, and the like can be considered. Further, the setting method of the pixels around the coordinates (X, Y) does not have to be uniform in the processing target image, and the setting method may be appropriately changed for each position in the processing target image.

次に、虚像位置特定部205による虚像又はぼけ領域の位置の特定方法についてより具体的に説明する。虚像位置特定部205は、コヒーレンスマップ算出部204によって算出されたコヒーレンスマップの各画素の画素値(局所的なコヒーレンスの値)を確認する。そして、虚像位置特定部205は、コヒーレンスマップの各画素の画素値がほぼ0であることを確認した場合には、コヒーレンスマップの範囲内に虚像又はぼけ領域が存在しないものとして判断する。   Next, a method for specifying the position of the virtual image or the blurred area by the virtual image position specifying unit 205 will be described more specifically. The virtual image position specifying unit 205 confirms the pixel value (local coherence value) of each pixel of the coherence map calculated by the coherence map calculation unit 204. Then, when the virtual image position specifying unit 205 confirms that the pixel value of each pixel of the coherence map is approximately 0, the virtual image position specifying unit 205 determines that no virtual image or blur area exists within the range of the coherence map.

他方、虚像位置特定部205は、コヒーレンスマップの各画素の画素値が0よりも大きいことを確認した場合には、コヒーレンスマップの範囲内に虚像又はぼけ領域が存在していると判断する。また、虚像位置特定部205に、コヒーレンスマップの各画素の画素値に対する閾値が予め設定されていてもよい。この場合には、虚像位置特定部205は、コヒーレンスマップの各画素の画素値が閾値以上の値となる領域を、虚像又はぼけ領域の存在する位置として特定可能となる。   On the other hand, when the virtual image position specifying unit 205 confirms that the pixel value of each pixel of the coherence map is greater than 0, the virtual image position specifying unit 205 determines that a virtual image or a blur region exists within the range of the coherence map. Further, a threshold value for the pixel value of each pixel of the coherence map may be set in the virtual image position specifying unit 205 in advance. In this case, the virtual image position specifying unit 205 can specify an area where the pixel value of each pixel of the coherence map is a value equal to or larger than the threshold as a position where the virtual image or the blur area exists.

この閾値については、レーダ画像利用者が経験に基づいて設定してもよい。また、コヒーレンスマップの画素値の出現確率分布(ヒストグラム)に基づいて閾値を設定してもよい。さらに、事前情報として虚像又は画像のぼけ領域が存在しない領域が判明している場合には、その判明している領域のコヒーレンスマップの画素値の最大値、平均値又は標準偏差等の統計量に基づいて閾値を設定してもよい。   This threshold value may be set by a radar image user based on experience. Further, the threshold may be set based on the appearance probability distribution (histogram) of the pixel values of the coherence map. In addition, if a virtual image or an area where there is no blurred area of the image is known as prior information, the maximum value, average value, or standard deviation of the pixel value of the coherence map of the known area is added to the statistics. A threshold may be set based on this.

次に、動作について説明する。図5は、図4のレーダ画像処理装置200の動作を示すフローチャートである。図5において、ステップS201では、画像シフト処理部103が、処理対象画像を、シフト量に合わせてシフトする。この後、ステップS202では、コヒーレンスマップ算出部204が、シフト前の処理対象画像と、シフト後の処理対象画像とのコヒーレンスマップを算出する。   Next, the operation will be described. FIG. 5 is a flowchart showing the operation of the radar image processing apparatus 200 of FIG. In FIG. 5, in step S <b> 201, the image shift processing unit 103 shifts the processing target image according to the shift amount. Thereafter, in step S202, the coherence map calculation unit 204 calculates a coherence map between the processing target image before the shift and the processing target image after the shift.

続いて、ステップS203では、虚像位置特定部205が、コヒーレンスマップにおける各画素のコヒーレンスの値から虚像又はぼけ領域の位置を特定する。ステップS204では、虚像位置特定部205が、ステップS203で特定した虚像又はぼけ領域の位置である虚像分布を虚像分布格納部206に格納する。これにより、実施の形態2のレーダ画像処理装置200の動作が終了する。   Subsequently, in step S203, the virtual image position specifying unit 205 specifies the position of the virtual image or the blurred area from the coherence value of each pixel in the coherence map. In step S204, the virtual image position specifying unit 205 stores the virtual image specified in step S203 or the virtual image distribution that is the position of the blurred area in the virtual image distribution storage unit 206. Thereby, the operation of the radar image processing apparatus 200 according to the second embodiment ends.

上記のような実施の形態2のレーダ画像処理装置200によれば、コヒーレンスマップ算出部204が、分解能を基本単位とするシフト量でのシフト後の処理対象画像と、シフト前の処理対象画像とから、コヒーレンスマップを算出する。この構成により、アジマスアンビギュイティ以外の種類の虚像についても検出対象とすることができるとともに、処理対象画像内に存在する虚像又は画像のぼけ領域の位置を特定することができる。   According to the radar image processing apparatus 200 of the second embodiment as described above, the coherence map calculation unit 204 includes the processing target image after the shift with the shift amount having the resolution as a basic unit, the processing target image before the shift, From this, a coherence map is calculated. With this configuration, a virtual image of a type other than the azimuth ambiguity can be set as a detection target, and the position of the virtual image or the blurred region of the image existing in the processing target image can be specified.

また、虚像分布格納部206に格納された虚像分布の情報は、レーダ画像処理装置200の後段の処理装置(図示せず)において、実施の形態1と同様に、例えば処理対象画像からの虚像又はぼけ領域の除去処理等に用いることができる。さらに、虚像やぼけ領域についての位置表示に用いることもできる。また、虚像分布の情報から、虚像の種類を推定することもできる。例えば、レーダ画像のレンジ方向(距離方向)又はアジマス方向(角度方向)の軸に沿って直線状の虚像が存在する場合には、その虚像は干渉波によって生じたものであると推定できる。   Further, the virtual image distribution information stored in the virtual image distribution storage unit 206 is, for example, a virtual image from the processing target image or the like in the processing apparatus (not shown) at the subsequent stage of the radar image processing apparatus 200, as in the first embodiment. It can be used for blur region removal processing or the like. Furthermore, it can also be used to display the position of a virtual image or blurred area. In addition, the type of virtual image can be estimated from the virtual image distribution information. For example, if there is a linear virtual image along the axis in the range direction (distance direction) or azimuth direction (angular direction) of the radar image, it can be estimated that the virtual image is caused by an interference wave.

さらに、このコヒーレンスマップを表示機器に表示する場合には、処理対象画像内に存在する虚像又はぼけ領域について、レーダ画像利用者の目視によって特定することもできる。   Furthermore, when this coherence map is displayed on a display device, a virtual image or a blur region existing in the processing target image can be specified by visual observation of a radar image user.

なお、実施の形態2における式(2)のコヒーレンスに代えて、実施の形態1の場合と同様に、画像間の相関を計算する方法、例えば相関関数や変数変換した後に相関係数を計算する方法等を用いてもよい。   Note that, instead of the coherence of the expression (2) in the second embodiment, a method for calculating correlation between images, for example, a correlation coefficient or a correlation coefficient after variable conversion is calculated, as in the first embodiment. A method or the like may be used.

また、実施の形態2の図5におけるステップS202を介さずに、座標(x,y)の画素と座標(x+Δx,y+Δy)の画素とのコヒーレンスマップを直接的に算出してもよい。ただし、虚像又は画像のぼけ領域を検出するためには、シフト量Δx及びΔyのうち、少なくともいずれか一方が分解能以上の値である必要がある。   Further, the coherence map between the pixel at coordinates (x, y) and the pixel at coordinates (x + Δx, y + Δy) may be directly calculated without going through step S202 in FIG. 5 of the second embodiment. However, in order to detect a virtual image or a blurred region of an image, at least one of the shift amounts Δx and Δy needs to be a value equal to or higher than the resolution.

さらに、実施の形態2におけるコヒーレンスマップ算出部204を実施の形態1におけるレーダ画像処理装置100に追加してもよい。この場合には、コヒーレンスマップ算出部204によって算出されたコヒーレンスマップにおけるコヒーレンスの出現確率分布(ヒストグラム)に基づいて閾値を設定し、その設定した閾値を用いて、虚像判断部105が処理対象画像内に虚像が存在するか否かを判断してもよい。   Furthermore, the coherence map calculation unit 204 in the second embodiment may be added to the radar image processing apparatus 100 in the first embodiment. In this case, a threshold value is set based on the appearance probability distribution (histogram) of the coherence in the coherence map calculated by the coherence map calculation unit 204, and the virtual image determination unit 105 uses the set threshold value in the processing target image. It may be determined whether or not a virtual image exists.

実施の形態3.
実施の形態1のレーダ画像処理装置100は、処理対象画像内における虚像又はぼけ領域の有無を判断し、実施の形態2のレーダ画像処理装置200は、処理対象画像内で虚像又はぼけ領域が存在している領域を特定した。これに対して、実施の形態3のレーダ画像処理装置300は、処理対象画像内に存在する虚像又はぼけ領域の広がり方向及び大きさ(ぼけ具合)を解析する。
Embodiment 3 FIG.
The radar image processing apparatus 100 according to the first embodiment determines the presence or absence of a virtual image or blur area in the processing target image, and the radar image processing apparatus 200 according to the second embodiment has a virtual image or blur area in the processing target image. Identified the area that is. On the other hand, the radar image processing apparatus 300 according to the third embodiment analyzes the spreading direction and size (blurring degree) of a virtual image or a blur area existing in the processing target image.

図6は、この発明の実施の形態3によるレーダ画像処理装置を示すブロック図である。図6において、実施の形態3のレーダ画像処理装置300は、レーダ画像格納部301、パラメータ格納部302、画像シフト処理部303、コヒーレンス算出部304、シフト量対相関対応部305、画像ぼけ解析部306及び解析結果格納部307を有している。レーダ画像格納部301、パラメータ格納部302、画像シフト処理部303及びコヒーレンス算出部304の構成は、実施の形態1の各機能ブロック101〜104の構成と同様である。   FIG. 6 is a block diagram showing a radar image processing apparatus according to Embodiment 3 of the present invention. In FIG. 6, a radar image processing apparatus 300 according to the third embodiment includes a radar image storage unit 301, a parameter storage unit 302, an image shift processing unit 303, a coherence calculation unit 304, a shift amount versus correlation correspondence unit 305, and an image blur analysis unit. 306 and an analysis result storage unit 307. The configurations of the radar image storage unit 301, the parameter storage unit 302, the image shift processing unit 303, and the coherence calculation unit 304 are the same as the configurations of the functional blocks 101 to 104 in the first embodiment.

シフト量対相関対応部305は、コヒーレンス算出部304によって算出されたコヒーレンスの値と、そのコヒーレンスの算出に用いられたシフトパラメータとを、コヒーレンス算出部304及びパラメータ格納部302から受ける。また、シフト量対相関対応部305は、コヒーレンス算出部304によって算出されたコヒーレンスと、そのコヒーレンスの算出の際に用いられたシフト量とを互いに対応付けてシフト量対相関情報として保存する。   The shift amount versus correlation correspondence unit 305 receives from the coherence calculation unit 304 and the parameter storage unit 302 the coherence value calculated by the coherence calculation unit 304 and the shift parameter used to calculate the coherence. Also, the shift amount / correlation correspondence unit 305 stores the coherence calculated by the coherence calculation unit 304 and the shift amount used in calculating the coherence in association with each other and stored as shift amount / correlation information.

さらに、シフト量対相関対応部305は、シフト量対相関情報を保存した後に、未だコヒーレンスを計算していないシフトパラメータ(即ち、シフト量及びシフト方向)にパラメータ格納部302に格納されたシフトパラメータを変更し、処理対象画像についてのシフト量及びシフト方向を変更する。このシフト量の変更する範囲については、予め設定された範囲であり、x軸及びy軸のそれぞれの分解能を最小単位とする範囲である。また、シフト量対相関対応部305は、シフトパラメータを変更した後に、画像シフト処理部303及びコヒーレンス算出部304にそれぞれの処理を繰り返し実行させる。   Furthermore, the shift amount-correlation correspondence unit 305 stores the shift amount-correlation information, and then stores the shift parameter stored in the parameter storage unit 302 in a shift parameter for which coherence has not yet been calculated (that is, the shift amount and the shift direction). To change the shift amount and the shift direction for the processing target image. The range in which the shift amount is changed is a preset range, and is a range in which the resolution of each of the x-axis and the y-axis is the minimum unit. The shift amount correlation correlation unit 305 causes the image shift processing unit 303 and the coherence calculation unit 304 to repeatedly execute the respective processes after changing the shift parameter.

画像ぼけ解析部306は、シフト量対相関対応部305からシフト量対相関情報を受ける。また、画像ぼけ解析部306は、シフト量対相関情報に基づいて、処理対象画像内に存在する虚像又はぼけ領域の広がり方向及びぼけ具合を解析する。解析結果格納部307には、画像ぼけ解析部306による解析結果が格納される。   The image blur analysis unit 306 receives shift amount versus correlation information from the shift amount versus correlation correspondence unit 305. In addition, the image blur analysis unit 306 analyzes the spread direction and the degree of blur of a virtual image or a blur region existing in the processing target image based on the shift amount versus correlation information. The analysis result storage unit 307 stores the analysis result by the image blur analysis unit 306.

次に、画像ぼけ解析部306による解析方法についてより具体的に説明する。画像ぼけ解析部306は、例えば、図7に示すような一次元の関数、シフト量を2軸とした2次元の関数、シフト量若しくはコヒーレンスを変数変換して作成した関数、又はシフト量毎のコヒーレンスの計算結果から算出した統計量等を用いて、処理対象画像内の虚像又はぼけ領域の存在する方向及びぼけ具合を解析する。ここでは、図7の例を用いて説明する。   Next, the analysis method performed by the image blur analysis unit 306 will be described more specifically. The image blur analysis unit 306 includes, for example, a one-dimensional function as illustrated in FIG. 7, a two-dimensional function having two shift axes, a function created by variable conversion of the shift amount or coherence, or the shift amount Using the statistic calculated from the coherence calculation result or the like, the direction in which the virtual image or the blur region exists in the processing target image and the degree of blur are analyzed. Here, a description will be given using the example of FIG.

画像ぼけ解析部306は、ぼけ領域の広がり方向について解析する場合には、画像の軸xのシフト量Δxについて、図7のコヒーレンスの関数を計算した結果と、軸yのシフト量Δyについてコヒーレンスの関数を計算した結果とを作成する。ただし、この際のシフト量の基本単位は各軸方向の画像の分解能とする(例えば分解能の2倍等としてグラフの横軸を標記する)。そして、画像ぼけ解析部306は、Δxに比べてΔyの方がコヒーレンスの低下度合いが小さいことを確認した場合には、y軸方向にぼけ領域が広がっていると判断できる。   When analyzing the spreading direction of the blur region, the image blur analysis unit 306 calculates the coherence function of FIG. 7 for the shift amount Δx of the axis x of the image and the coherence of the shift amount Δy of the axis y. Create the result of calculating the function. However, the basic unit of the shift amount at this time is the resolution of the image in each axis direction (for example, the horizontal axis of the graph is marked as twice the resolution). The image blur analysis unit 306 can determine that the blur region is widened in the y-axis direction when it is confirmed that Δy has a lower degree of reduction in coherence than Δx.

また、画像ぼけ解析部306は、ぼけ領域の大きさについて解析する場合には、図7に示すような閾値を用いる。即ち、画像ぼけ解析部306は、閾値をコヒーレンスが下回ったときのシフト量をぼけ領域の大きさとして解析できる。この閾値については、処理の実施者の経験に基づいても設定してもよい。あるいは、実施の形態2の虚像位置特定結果から虚像が存在していない領域を求め、その領域のコヒーレンスの平均値や分散等の統計量に基づいて設定してもよい。   The image blur analysis unit 306 uses a threshold as shown in FIG. 7 when analyzing the size of the blur region. In other words, the image blur analysis unit 306 can analyze the shift amount when the coherence is below the threshold as the size of the blur region. This threshold value may be set based on the experience of the process operator. Alternatively, a region where a virtual image does not exist may be obtained from the virtual image position specifying result of the second embodiment, and may be set based on a coherence average value or a statistic such as variance.

次に、動作について説明する。図8は、図6のレーダ画像処理装置300の動作を示すフローチャートである。図8において、ステップS301では、画像シフト処理部303が、処理対象画像を、シフト量に合わせてシフトする。この後、ステップS302では、コヒーレンス算出部304が、シフト前の処理対象画像と、シフト後の処理対象画像とのコヒーレンスを算出する。   Next, the operation will be described. FIG. 8 is a flowchart showing the operation of the radar image processing apparatus 300 of FIG. In FIG. 8, in step S301, the image shift processing unit 303 shifts the processing target image according to the shift amount. Thereafter, in step S302, the coherence calculation unit 304 calculates the coherence between the processing target image before the shift and the processing target image after the shift.

続いて、ステップS303では、シフト量対相関対応部305が、コヒーレンス算出部304によって算出されたコヒーレンスと、そのコヒーレンスの算出の際に用いられたシフト量とを互いに対応付けてシフト量対相関情報として保存する。そして、ステップS304では、シフト量対相関対応部305が、シフト量及びシフト方向を更新する。ここで、予め定められた範囲内でシフト量及びシフト方向毎に、ステップS301〜S304の演算処理が繰り返し実行される。   Subsequently, in step S303, the shift amount correlation correlation unit 305 associates the coherence calculated by the coherence calculation unit 304 with the shift amount used in the calculation of the coherence, and shift amount pair correlation information. Save as. In step S304, the shift amount correlation correlation unit 305 updates the shift amount and the shift direction. Here, the arithmetic processing of steps S301 to S304 is repeatedly executed for each shift amount and shift direction within a predetermined range.

ステップS301〜S304の演算処理が終了した後に、ステップS305の演算処理に移行する。ステップS305では、画像ぼけ解析部306が、シフト量対相関情報に基づいて、ぼけ領域の広がり方向及び大きさを解析する。ステップS306では、画像ぼけ解析部306が、ステップS305での解析結果の情報を解析結果格納部307に格納する。これにより、実施の形態3のレーダ画像処理装置300の動作が終了する。   After the calculation process of steps S301 to S304 is completed, the process proceeds to the calculation process of step S305. In step S305, the image blur analysis unit 306 analyzes the spread direction and size of the blur region based on the shift amount versus correlation information. In step S306, the image blur analysis unit 306 stores information on the analysis result in step S305 in the analysis result storage unit 307. Thereby, the operation of the radar image processing apparatus 300 according to the third embodiment is completed.

上記のような実施の形態3のレーダ画像処理装置300によれば、画像ぼけ解析部306が、シフト量対相関対応部305からのシフト量対相関情報に基づいて、処理対象画像内に存在するぼけ領域の広がり方向及び大きさを解析する。この構成により、画像ぼけ解析部306の解析結果の情報を用いて、虚像の分類(類別)をすることが可能となる。例えば、レンジ方向又はアジマス方向に一様なぼけ領域が処理対象画像内に存在する場合には、そのぼけ領域は干渉波によって生じた虚像であると推定することができる。   According to the radar image processing apparatus 300 of the third embodiment as described above, the image blur analysis unit 306 exists in the processing target image based on the shift amount versus correlation information from the shift amount versus correlation correspondence unit 305. Analyze the direction and size of the blur area. With this configuration, it is possible to classify (categorize) virtual images using information on the analysis result of the image blur analysis unit 306. For example, when a blur area that is uniform in the range direction or the azimuth direction exists in the processing target image, it can be estimated that the blur area is a virtual image generated by the interference wave.

また、解析結果格納部307に格納された解析結果の情報は、レーダ画像処理装置300の後段の処理装置(図示せず)において、例えば処理対象画像からの虚像又はぼけ領域の除去処理や、画像補正処理等に用いることができる。さらに、解析結果格納部307に格納された解析結果の情報のうちぼけ領域の大きさの情報は、レーダ画像処理装置300の後段のオートフォーカス処理装置(図示せず)において、ぼけ領域の大きさを動作パラメータとするレーダのオートフォーカス処理に用いることができる。   The analysis result information stored in the analysis result storage unit 307 is, for example, processed in a processing device (not shown) subsequent to the radar image processing device 300 to remove a virtual image or a blur region from the processing target image, It can be used for correction processing and the like. Further, out of the analysis result information stored in the analysis result storage unit 307, the information on the size of the blur region is obtained by the autofocus processing device (not shown) at the subsequent stage of the radar image processing device 300. Can be used for the radar auto-focus processing with the operation parameter as.

実施の形態4.
実施の形態1〜3では、レーダ画像の分解能が既知情報であることが処理の前提条件であった。これに対して、実施の形態4では、レーダ画像の分解能が未知情報である場合に、分解能推定部406によって分解能が推定される。
Embodiment 4 FIG.
In the first to third embodiments, the processing condition is that the resolution of the radar image is known information. On the other hand, in the fourth embodiment, when the resolution of the radar image is unknown information, the resolution estimation unit 406 estimates the resolution.

図9は、この発明の実施の形態4によるレーダ画像処理装置を示すブロック図である。図9において、実施の形態4のレーダ画像処理装置400の構成は、シフト量の基本単位が距離である点と、レーダ画像処理装置300の画像ぼけ解析部306及び解析結果格納部307に代えて分解能推定部406及び分解能推定値格納部407を有する点とを除いて、実施の形態3のレーダ画像処理装置300の構成と同様である。即ち、実施の形態4のレーダ画像処理装置400では、実施の形態3のレーダ画像処理装置300の機能ブロック301〜305が分解能推定用の機能ブロックとされる。   FIG. 9 is a block diagram showing a radar image processing apparatus according to Embodiment 4 of the present invention. In FIG. 9, the configuration of the radar image processing apparatus 400 according to the fourth embodiment is such that the basic unit of the shift amount is a distance, and the image blur analysis unit 306 and the analysis result storage unit 307 of the radar image processing apparatus 300 are replaced. Except for having a resolution estimation unit 406 and a resolution estimation value storage unit 407, the configuration is the same as that of the radar image processing apparatus 300 of the third embodiment. That is, in the radar image processing apparatus 400 of the fourth embodiment, the functional blocks 301 to 305 of the radar image processing apparatus 300 of the third embodiment are function blocks for resolution estimation.

分解能推定部406は、シフト量対相関対応部305からシフト量対相関情報を受ける。また、分解能推定部406は、シフト量対相関情報に基づいて、レーダ画像の分解能を推定する。具体的に、分解能推定部406は、図10に示すようなシフト量とコヒーレンスとの関係を表すグラフを作成する。この図10のグラフと、実施の形態3における図7のグラフとの相違点は、シフト量の単位を分解能ではなく距離(例えばメートル単位)とした点である。   The resolution estimation unit 406 receives the shift amount versus correlation information from the shift amount versus correlation correspondence unit 305. Further, the resolution estimation unit 406 estimates the resolution of the radar image based on the shift amount versus correlation information. Specifically, the resolution estimation unit 406 creates a graph representing the relationship between the shift amount and coherence as shown in FIG. The difference between the graph of FIG. 10 and the graph of FIG. 7 in the third embodiment is that the unit of the shift amount is not a resolution but a distance (for example, a meter unit).

分解能推定部406は、コヒーレンスが予め設定された閾値を下回ったときのシフト量を分解能として推定できる。この閾値の設定方法は、実施の形態3と同様である。分解能推定値格納部407には、分解能推定部406による分解能の推定値が格納される。なお、図11は、図9のレーダ画像処理装置400の動作を示すフローチャートであり、レーダ画像処理装置400の動作の概要としては、実施の形態3のレーダ画像処理装置300の動作と同様である。   The resolution estimation unit 406 can estimate the shift amount when the coherence falls below a preset threshold as the resolution. This threshold setting method is the same as in the third embodiment. The estimated resolution value storage unit 407 stores the estimated resolution value by the resolution estimation unit 406. FIG. 11 is a flowchart showing the operation of the radar image processing apparatus 400 of FIG. 9. The outline of the operation of the radar image processing apparatus 400 is the same as the operation of the radar image processing apparatus 300 of the third embodiment. .

上記のような実施の形態4のレーダ画像処理装置400によれば、分解能推定部406が、シフト量対相関対応部305からのシフト量対相関情報に基づいて、レーダ画像の分解能を推定する。この構成により、レーダ画像の分解能が未知情報である場合に、処理対象画像から分解能を推定することが可能になる。   According to the radar image processing apparatus 400 of the fourth embodiment as described above, the resolution estimation unit 406 estimates the resolution of the radar image based on the shift amount versus correlation information from the shift amount versus correlation correspondence unit 305. With this configuration, when the resolution of the radar image is unknown information, the resolution can be estimated from the processing target image.

なお、実施の形態4のレーダ画像処理装置400の構成と、実施の形態1〜3のレーダ画像処理装置100〜300のそれぞれとの構成を組み合わせてもよい。このように組み合わせた構成では、分解能が未知情報である場合にも、レーダ画像処理装置100〜300が、レーダ画像処理装置400によって算出された分解能を基本単位とするシフト量を用いて、虚像若しくはぼけ領域の検出、位置特定又は解析を行うことができる。   The configuration of the radar image processing apparatus 400 according to the fourth embodiment may be combined with the configuration of each of the radar image processing apparatuses 100 to 300 according to the first to third embodiments. In such a combined configuration, even when the resolution is unknown information, the radar image processing apparatuses 100 to 300 use the shift amount having the resolution calculated by the radar image processing apparatus 400 as a basic unit, Detection, position specification, or analysis of a blurred area can be performed.

100〜400 レーダ画像処理装置、101,301 レーダ画像格納部、102,302 パラメータ格納部、103,303 画像シフト処理部、104,304 コヒーレンス算出部、105 虚像判断部、106 判断結果格納部、204 コヒーレンスマップ算出部、205 虚像位置特定部、206 虚像分布格納部、305 シフト量対相関対応部、306 画像ぼけ解析部、307 解析結果格納部、406 分解能推定部、407 分解能推定値格納部。   100 to 400 Radar image processing apparatus, 101, 301 Radar image storage unit, 102, 302 Parameter storage unit, 103, 303 Image shift processing unit, 104, 304 Coherence calculation unit, 105 Virtual image determination unit, 106 Determination result storage unit, 204 Coherence map calculation unit, 205 virtual image position specifying unit, 206 virtual image distribution storage unit, 305 shift amount correlation correlation unit, 306 image blur analysis unit, 307 analysis result storage unit, 406 resolution estimation unit, 407 resolution estimation value storage unit.

Claims (6)

取得したレーダ画像の一部及び全部のいずれか一方である処理対象画像について、分解能を基本単位とする予め設定されたシフト量に基づいて、前記レーダ画像を表すための2軸の少なくともいずれか一方の軸方向へシフトするシフト処理部と、
前記シフト処理部によるシフト後の前記処理対象画像と、そのシフト前の前記処理対象画像とから、前記処理対象画像における各画素と、前記各画素の周囲の画素とについてのコヒーレンスの分布であるコヒーレンスマップを算出するコヒーレンスマップ算出部と、
前記処理対象画像内に虚像が存在する場合に、前記コヒーレンスマップに基づいて、前記処理対象領域内に存在する虚像の位置を特定する虚像位置特定部と
を備えることを特徴とするレーダ画像処理装置。
At least one of two axes for representing the radar image based on a preset shift amount having a resolution as a basic unit for a processing target image that is one of all or part of the acquired radar image A shift processing unit for shifting in the axial direction of
Coherence that is a coherence distribution of each pixel in the processing target image and pixels around each pixel from the processing target image after the shift by the shift processing unit and the processing target image before the shift. A coherence map calculation unit for calculating a map;
A radar image processing apparatus comprising: a virtual image position specifying unit that specifies a position of a virtual image existing in the processing target area based on the coherence map when a virtual image exists in the processing target image. .
前記虚像位置特定部は、前記コヒーレンスマップの画素値が所定の閾値を超えた画素について、虚像であると判断し、前記処理対象領域内の虚像の位置を特定する
ことを特徴とする請求項記載のレーダ画像処理装置。
The virtual image position specifying unit, the pixel in which the pixel values of the coherence map exceeds a predetermined threshold value, determines that the virtual image, according to claim 1, characterized in that locating the virtual image of the processing target area The radar image processing apparatus described.
前記虚像位置特定部は、前記コヒーレンスマップ算出部によって算出されたコヒーレンスマップにおけるコヒーレンスの出現確率分布に基づいて閾値を定め、その定めた閾値を前記所定の閾値に代えて用いて、前記コヒーレンスマップ内の前記処理対象領域内の虚像の位置を特定する
ことを特徴とする請求項記載のレーダ画像処理装置。
The virtual image position specifying unit determines a threshold based on the appearance probability distribution of the coherence in the coherence map calculated by the coherence map calculation unit, and uses the determined threshold instead of the predetermined threshold in the coherence map. The radar image processing apparatus according to claim 2 , wherein a position of a virtual image in the processing target area is specified.
取得したレーダ画像の一部及び全部のいずれかである処理対象画像について、分解能を基本単位とする予め設定されたシフト量に基づいて、前記レーダ画像を表すための2軸の少なくともいずれか一方の軸方向へシフトするシフト処理部と、
前記シフト処理部によるシフト後の前記処理対象画像と、そのシフト前の前記処理対象画像とのコヒーレンスを算出するコヒーレンス算出部と、
前記コヒーレンス算出部によって算出されたコヒーレンスと、そのコヒーレンスの算出の際に用いられたシフト量とを互いに対応付けてシフト量対相関情報とするとともに、前記処理対象画像についてのシフト量及びシフト方向を変更して、前記シフト処理部及び前記コヒーレンス算出部にそれぞれの処理を繰り返し実行させるシフト量対相関対応部と、
前記シフト量対相関対応部からの前記シフト量対相関情報に基づいて、前記処理対象画像内に存在するぼけ領域の広がり方向及び大きさを解析する画像ぼけ解析部と
を備えることを特徴とするレーダ画像処理装置。
For a processing target image that is a part or all of the acquired radar image, at least one of two axes for representing the radar image based on a preset shift amount with resolution as a basic unit A shift processing unit for shifting in the axial direction;
A coherence calculation unit that calculates coherence between the processing target image after being shifted by the shift processing unit and the processing target image before the shift;
The coherence calculated by the coherence calculation unit and the shift amount used in the calculation of the coherence are associated with each other as shift amount-correlation information, and the shift amount and the shift direction for the processing target image are set. A shift amount-correlation correspondence unit that causes the shift processing unit and the coherence calculation unit to repeatedly execute the respective processes,
An image blur analysis unit that analyzes a spread direction and a size of a blur region existing in the processing target image based on the shift amount pair correlation information from the shift amount pair correlation correspondence unit. Radar image processing device.
前記画像ぼけ解析部は、前記シフト量対相関情報におけるシフト量及びコヒーレンスの関係を示す関数を、前記2軸のそれぞれの軸方向について算出し、その算出した関数の組に基づいて前記処理対象画像に存在するぼけ領域の方向及び大きさを解析する
ことを特徴とする請求項記載のレーダ画像処理装置。
The image blur analysis unit calculates a function indicating a relationship between a shift amount and coherence in the shift amount-correlation information for each axial direction of the two axes, and based on the calculated function set, the processing target image The radar image processing apparatus according to claim 4, wherein a direction and a size of a blur area existing in the area is analyzed.
前記シフト処理部は、画像処理用のシフト処理部であり、
前記処理対象画像について、距離を基本単位とする予め設定されたシフト量に基づいて、前記2軸の少なくともいずれか一方の軸方向へシフトする分解能推定用のシフト処理部と、
前記分解能推定用のシフト処理部によるシフト後の前記処理対象画像と、そのシフト前の前記処理対象画像とのコヒーレンスを算出する分解能推定用のコヒーレンス算出部と、
前記分解能推定用のコヒーレンス算出部によって算出されたコヒーレンスと、そのコヒーレンスの算出の際に用いられたシフト量とを互いに対応付けて分解能推定用のシフト量対相関情報とするとともに、前記処理対象画像についてのシフト量及びシフト方向を変更して、前記分解能推定用のシフト処理部、及び前記分解能推定用のコヒーレンス算出部にそれぞれの処理を繰り返し実行させる分解能推定用のシフト量対相関対応部と、
前記分解能推定用のシフト量対相関対応部からの前記分解能推定用のシフト量対相関情報に基づいて、前記レーダ画像の分解能を推定する分解能推定部と
をさらに備え、
前記画像処理用のシフト処理部は、前記分解能推定部によって推定された分解能を基本単位とするシフト量に基づいて、前記処理対象画像を前記2軸の少なくともいずれか一方の軸方向へシフトする
ことを特徴とする請求項1から請求項までのいずれか1項に記載のレーダ画像処理装置。
The shift processing unit is a shift processing unit for image processing,
A shift processing unit for resolution estimation that shifts in the axial direction of at least one of the two axes based on a shift amount set in advance with the distance as a basic unit for the processing target image;
A resolution estimation coherence calculation unit that calculates coherence between the processing target image after the shift by the resolution estimation shift processing unit and the processing target image before the shift;
The coherence calculated by the coherence calculating unit for resolution estimation and the shift amount used in the calculation of the coherence are associated with each other as shift amount-correlation information for resolution estimation, and the processing target image Changing the shift amount and the shift direction for the resolution estimation, the resolution estimation shift processing unit, and the resolution estimation coherence calculation unit repeatedly executing the respective processing, a resolution estimation shift amount versus correlation correspondence unit,
A resolution estimation unit that estimates the resolution of the radar image based on the resolution estimation shift amount pair correlation information from the resolution estimation shift amount pair correlation correspondence unit;
The shift processing unit for image processing shifts the processing target image in at least one of the two axes based on a shift amount having the resolution estimated by the resolution estimation unit as a basic unit. The radar image processing apparatus according to any one of claims 1 to 5, wherein
JP2009272323A 2009-11-30 2009-11-30 Radar image processing device Active JP5448757B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009272323A JP5448757B2 (en) 2009-11-30 2009-11-30 Radar image processing device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2009272323A JP5448757B2 (en) 2009-11-30 2009-11-30 Radar image processing device

Publications (2)

Publication Number Publication Date
JP2011112630A JP2011112630A (en) 2011-06-09
JP5448757B2 true JP5448757B2 (en) 2014-03-19

Family

ID=44235057

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009272323A Active JP5448757B2 (en) 2009-11-30 2009-11-30 Radar image processing device

Country Status (1)

Country Link
JP (1) JP5448757B2 (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6000074B2 (en) * 2012-11-08 2016-09-28 三菱電機株式会社 Image radar processing apparatus and image radar processing method
JP2015021823A (en) * 2013-07-18 2015-02-02 三菱電機株式会社 Radar image processing apparatus and radar image processing method
CA3028602C (en) * 2016-08-01 2019-04-30 Mitsubishi Electric Corporation Synthetic-aperture radar device
JP7570583B1 (en) * 2022-11-25 2024-10-21 三菱電機株式会社 Radar image processing device, radar image processing method, and radar image processing program

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2594908B2 (en) * 1985-11-15 1997-03-26 株式会社日立製作所 A method for measuring spatial resolution of synthetic aperture radar images.
JP4062771B2 (en) * 1997-05-06 2008-03-19 ソニー株式会社 Image conversion apparatus and method, and recording medium
JP2001208842A (en) * 2000-01-28 2001-08-03 Ntt Data Creation Kk Tie point detection method of sattelite image and its device and recording medium
JP2002214334A (en) * 2001-01-23 2002-07-31 Mitsubishi Heavy Ind Ltd Target detection apparatus
JP3564540B2 (en) * 2001-09-07 2004-09-15 独立行政法人情報通信研究機構 Land use classification processing method using SAR image
JP4443320B2 (en) * 2004-06-21 2010-03-31 三菱電機株式会社 Radar image processing device
JP4818898B2 (en) * 2006-12-25 2011-11-16 三菱電機株式会社 Radar image processing device
JP2009181024A (en) * 2008-01-31 2009-08-13 Nikon Corp Focusing device, optical equipment

Also Published As

Publication number Publication date
JP2011112630A (en) 2011-06-09

Similar Documents

Publication Publication Date Title
TWI767373B (en) A method, a computer program product, an apparatus and a frequency-modulated continuous-wave radar system
US8818077B2 (en) Stereo image matching apparatus and method
KR101248808B1 (en) Apparatus and method for removing noise on edge area
CN107316326B (en) Edge-based disparity map calculation method and device applied to binocular stereo vision
WO2012101967A1 (en) Positioning information forming device, detection device, and positioning information forming method
CN113077459B (en) Image definition detection method and device, electronic equipment and storage medium
US20110038510A1 (en) Image processing apparatus, image processing method, and program
WO2019031086A1 (en) Image processing system, server device, image processing method, and image processing program
JP2008185375A (en) SAR image 3D shape calculation apparatus and SAR image distortion correction apparatus
JP5448757B2 (en) Radar image processing device
JP6673160B2 (en) Distance measuring device
JP7230507B2 (en) Deposit detection device
JP2013002839A (en) Crack detection method
JP6084810B2 (en) Tracking processing apparatus and tracking processing method
US9689668B2 (en) Image processing apparatus and image processing method
JP2013117409A (en) Crack detection method
JP2017083427A5 (en)
CN111630565B (en) Image processing method, edge extraction method, processing equipment and storage medium
JP6473769B2 (en) Correction amount calculation device, correction device, and correction amount calculation method
KR101799143B1 (en) System and method for estimating target size
US8897547B2 (en) Precision improving device for three dimensional topographical data, precision improving method for three dimensional topographical data and recording medium
JP5059855B2 (en) Global motion estimation method
JP6637823B2 (en) Crack detection method
EP2835963A1 (en) Image processing device and method, and image processing program
JP6716769B1 (en) Image inspection method, image inspection apparatus, and image inspection program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120827

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130910

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20130911

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20131111

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: 20131126

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20131224

R150 Certificate of patent or registration of utility model

Ref document number: 5448757

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250