JP5787637B2 - 画像処理装置、画像処理方法 - Google Patents

画像処理装置、画像処理方法 Download PDF

Info

Publication number
JP5787637B2
JP5787637B2 JP2011137737A JP2011137737A JP5787637B2 JP 5787637 B2 JP5787637 B2 JP 5787637B2 JP 2011137737 A JP2011137737 A JP 2011137737A JP 2011137737 A JP2011137737 A JP 2011137737A JP 5787637 B2 JP5787637 B2 JP 5787637B2
Authority
JP
Japan
Prior art keywords
coordinate
coordinates
output
pixel
image 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.)
Expired - Fee Related
Application number
JP2011137737A
Other languages
English (en)
Other versions
JP2013004031A5 (ja
JP2013004031A (ja
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.)
Canon Inc
Original Assignee
Canon Inc
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 Canon Inc filed Critical Canon Inc
Priority to JP2011137737A priority Critical patent/JP5787637B2/ja
Priority to US13/490,260 priority patent/US8509568B2/en
Publication of JP2013004031A publication Critical patent/JP2013004031A/ja
Priority to US13/957,375 priority patent/US20130315502A1/en
Publication of JP2013004031A5 publication Critical patent/JP2013004031A5/ja
Application granted granted Critical
Publication of JP5787637B2 publication Critical patent/JP5787637B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N9/00Details of colour television systems
    • H04N9/12Picture reproducers
    • H04N9/31Projection devices for colour picture display, e.g. using electronic spatial light modulators [ESLM]
    • H04N9/3179Video signal processing therefor
    • H04N9/3185Geometric adjustment, e.g. keystone or convergence
    • G06T3/18

Description

本発明は、画像変形技術に関するものである。
従来の画像変形を行う画像処理装置では、まず入力画素を1画面分メモリに記憶する。次に、出力画素の走査順に、出力画素の座標値から入力画素の座標値を座標変換の逆変換によって算出する。さらに、逆変換で得られた入力画素の座標値の整数部を用いて入力画素の近傍の画素をメモリから読み出し、近傍の画素と入力画素の座標値の小数部とを用いて補間計算によって出力画素の値を決定する操作を繰り返す。これらの操作により画像変形を行う画像処理を実現していた。上記方法を以下では、「読み出し時座標変換」と呼ぶ。例えば、特許文献1では、出力画像をタイルに分割し、分割されたタイル毎に読み出し時座標変換を行う画像処理方法が示されている。
ところで、フロントプロジェクタの台形補正やカメラのレンズ補正を行う場合、入力画像から出力画像への拡大率は最小0.6倍〜最大1.25倍程度になる。さらに、映像入出力をリアルタイムで実現しなくてはならないという制約がある。映像入出力を1[pixel/cyc]で行うという制約下で上記に示す画像処理方法を実装すると、必要なメモリ帯域のピーク値が書き込みと読み出しの合計で(1+1/拡大率の最小値)[pixel/cyc]になってしまう。例えば拡大率の最小値が0.6倍であればピーク時に2.67[pixel/cyc]のメモリ帯域が必要になってしまう。
これに対し、以下で「書き込み時座標変換」と呼ぶ方式がある。この方式は、まず、走査順に入力される入力画素に対して座標変換によって出力座標を算出する。次に、出力座標の整数部からメモリの格納先アドレスを算出する。さらに、出力座標の整数部から座標変換の逆変換により入力画像中の座標を求める。逆変換により求めた座標の整数部が元々の入力画素の座標に一致するときに、入力画素とその近傍の画素と、座標変換の逆変換によって求めた座標値の小数部とを用いて出力する画素値を補完計算によって求める。この画素値を先に求めた格納先アドレスに格納する。
書き込み時座標変換方式における必要なメモリ帯域のピーク値は、書き込みと読み出しの合計で(1+拡大率の最大値)[pixel/cyc]になる。拡大率の最大値が1.25倍の時、この値は2.25[pixel/cyc]になり、想定する拡大率の範囲内では読み出し時座標変換よりメモリ帯域のピーク値を小さくすることができる。
しかしながら、上記のままの書き込み時座標変換では、出力されない画素が発生することがあるため、対策が必要である。例えば、書き込み時座標変換を行っている特許文献2では、入力画素の座標をサブピクセル単位で走査することで出力されない画素が発生することを防いでいる。
特開2005−135096号公報 特開平6−149993号公報
特許文献1の方式では、必要なメモリ帯域のピーク値が、想定する拡大率の範囲内で大きくなってしまう。特許文献2の方式では、サブピクセルでの走査をX軸方向とY軸方向で共に1/2ピクセルとした場合でも座標変換で4点分、逆変換で4点分の演算が必要になる。射影変換は座標変換と逆変換の演算量が同じで、どちらの演算にも除算が含まれる。このため、特許文献2の方式では、ハードウェアで実装したときの回路量が大きくなり、ソフトウェア実装であれば処理性能が低下してしまう。
本発明はこのような問題に鑑みてなされたものであり、必要なメモリ帯域のピーク値が小さく、かつ、ハードウェア実装時の回路量が小さい、あるいはソフトウェア実装時の処理性能が高い画像変形技術を提供することを目的とする。
本発明の目的を達成するために、例えば、本発明の画像処理装置は以下の構成を備える。即ち、入力画像が変形処理された出力画像を生成する画像処理装置であって、
前記入力画像を取得する取得手段と、
前記入力画像上の第1座標に対応する前記出力画像上の第2座標の整数部と小数部を、前記変形処理で用いられるパラメータによって特定する特定手段と、
前記第2座標の整数部に対応する前記出力画像上の座標から規定範囲内に位置する前記出力画像上の第1複数座標から、該第1複数座標よりも数が少ない第2複数座標を、該第2座標の小数部に基づいて選択する選択手段と、
前記第2複数座標に対して、前記パラメータに基づく変換の逆変換を行うことで、前記出力画像上の前記第2座標の整数部に対応する前記出力画像の画素データを生成するために参照される前記入力画像の画素データに対応する第3複数座標を判定する判定手段と、
前記判定手段により判定された第3複数座標に対応する前記入力画像の画素データを読み出して、前記第2座標の整数部に対応する前記出力画像の画素データを生成する生成手段と、
前記生成手段により生成された画素データを前記出力画像の画素データとして出力する出力手段と
を有することを特徴とする。
本発明の構成によれば、必要なメモリ帯域のピーク値が小さく、かつ、ハードウェア実装時の回路量が小さい、あるいはソフトウェア実装時の処理性能が高い画像変形技術を提供することができる。
画像処理装置の機能構成例を示すブロック図。 出力座標算出部106のより詳細な機能構成例を示すブロック図。 出力画素周辺座標生成部203のより詳細な機能構成例を示すブロック図。 サブ領域判定部207の動作を説明する為のグラフ。 サブ領域判定部207の動作の変形例を説明するためのグラフ。 画像処理装置の機能構成例を示すブロック図。 出力座標算出部606の機能構成例を示すブロック図。 領域を求める処理を説明するグラフ。 領域を求める処理を説明するグラフ。 領域を求める処理を説明するグラフ。
以下、添付図面を参照し、本発明の好適な実施形態について説明する。なお、以下説明する実施形態は、本発明を具体的に実施した場合の一例を示すもので、特許請求の範囲に記載の構成の具体的な実施例の1つである。
[第1の実施形態]
本実施形態では、座標変換後の出力座標の小数部によって2領域に分割することで座標変換(座標逆変換)が必要な周辺画素数を中心座標も含めて5点まで削減する例と、4領域に分割することで4点まで削減する例とを示す。どの様な領域に分割するかは任意であり、その結果、逆変換(逆写像)が必要な周辺画素数が何点まで削減されるかは領域分割方法によって変わってくる。しかし、どの様な領域に分割するかは容易に変更することができる。
また、本実施形態では、画像変形処理の一例として、射影変換を行う場合について説明するが、画素位置(座標)の座標変換とその逆変換が存在するのであれば、射影変換以外の幾何変換にも適用することができる。
先ず、本実施形態に係る画像処理装置の機能構成例について、図1のブロック図を用いて説明する。画像処理装置には、画像変形処理を行う対象となる画像(入力画像)が入力されるのであるが、本実施形態ではこの入力画像が画素単位で入力されるものとして説明する。しかし、入力画像をどのようにして画像処理装置に入力するのかについては特に限定するものではない。
画像処理装置に入力された入力画素の画素値(若しくは画像処理装置内のメモリから読み出された画素値)は、入力画素値121として入力画素記憶部101に格納されると共に、出力座標算出部106にも入力される。
画素補間処理部102は、入力画素に対応する出力画像上の画素の画素値を求める為に必要な画素(入力画素の周辺画素位置における画素)の画素値を入力周辺画素値122として入力画素記憶部101から読み出す。更に画素補間処理部102は、入力画素について出力座標算出部106が求めた各画素の画素位置(x座標値とy座標値)の小数部127を、出力座標算出部106から取得する。そして画素補間処理部102は、読み出した入力周辺画素値122と、取得した各画素の画素位置の小数部127と、を用いて補間処理を行い、入力画素に対応する出力画像上の画素の画素値を求める。
補間処理には様々な補間技術を適用することができるが、例えば、バイリニア補間を行う場合、入力画素を含む2×2ピクセルからバイリニア補間を行い、補間された画素を出力する。また、バイキュービック補間を行う場合、入力画素を含む4×4ピクセルからバイキュービック補間を行い、補間された画素を出力する。
そして画素補間処理部102は、求めた画素値を入力補間画素値123としてメモリ書き込み部103に対して送出する。
メモリ書き込み部103は、出力座標算出部106から「入力画素に対応する出力画像上の画素の画素位置の整数部128」を受け取り、受け取った整数部128で表される画素位置に対応する記憶部104内のアドレスを特定する。そしてメモリ書き込み部103は、画素補間処理部102から受け取った入力補間画素値123を書き込み画素値124として、この特定したアドレスに格納する。
ランダムなアドレスに対する書き込みは低速である。しかし、連続するメモリアドレスへの書き込みが高速であるDRAMのようなメモリを記憶部104として用いる場合、メモリ書き込み部103内で入力補間画素値123のバッファリングを行う。そして、複数の入力補間画素値123をまとめて記憶部104に書き出すようにする。
メモリ読み出し部105は、記憶部104に格納されている射影変換済みの画像の各画素の画素値を順次読み出し、出力画素値126として出力する。出力する順序については特に限定するものではないが、ラスタ順や縦方向を主走査とする方式や、タイルに分割してタイル毎に読み出しを行う方法などが考えられる。
次に、出力座標算出部106のより詳細な機能構成例について、図2のブロック図を用いて説明する。入力画素座標生成部201は、入力画素値121を受け取ると、この入力画素の入力画像上(入力画像中)における画素位置を求める。例えば、入力画素値121の中にVSYNC、HSYNC信号が含まれているとする。この場合、入力画素座標生成部201はVSYNCを受けると、Y座標を0にリセットする。また入力画素座標生成部201は、HSYNCを受けると、X座標を0にリセットすると共に、前回のHSYNC信号を受けたときから画素値が入力されるとY座標を+1する。画素値が入力された時はX座標を+1する。そして入力画素座標生成部201は、このようにしてカウントしたX座標、Y座標をセットにして、入力画素座標値221として座標変換部202及び入力座標比較部205に対して送出する。なお、入力画素値121に入力画素の画素位置が含まれている場合には、この画素位置をそのまま入力画素座標値221として、座標変換部202及び入力座標比較部205に対して送出する。
座標変換部202は、射影変換用のマトリクスを用いて、入力画素座標値221に対する射影変換を行い、射影変換後の座標値を出力画素座標値222として求める。このマトリクスは、m11〜m33を成分とする3×3の行列で表される。この行列を用いて入力画素座標値221(xi,yi)に対する射影変換を行うことで、この入力画素座標値221(xi,yi)が写像される画素座標位置である出力画素座標値222(xo,yo)を求める。
この時、2×2ピクセルのバイリニア補間や4×4ピクセルのバイキュービックフィルタなどの偶数次のフィルタであればフィルタの中心座標が入力画素座標の中心から(+0.5,+0.5)ずれている。そのため、(xi+0.5,yi+0.5)に対する座標変換を行う。奇数次のフィルタであればフィルタの中心座標が(xi,yi)と同じであるため、(xi,yi)に対する座標変換を行う。以下に示す例は偶数次フィルタの時の出力画素座標値222(xo,yo)算出式である。
xo0=m11・(xi+0.5)+m12・(yi+0.5)+m13
yo0=m21・(xi+0.5)+m22・(yi+0.5)+m23
zo0=m31・(xi+0.5)+m32・(yi+0.5)+m33
xo=xo0/zo0
yo=yo0/zo0
そして座標変換部202は、このようにして求めた出力画素座標値222(xo、yo)を、後段の出力画素周辺座標生成部203に対して送出する。
出力画素周辺座標生成部203は、xoの小数部に0.5を加えた値の小数部が既定値以上であるか否か及び/又はyoの小数部に0.5を加えた値の小数部が既定値以上であるか否かの組み合わせに対して予め設定された複数の参照画素位置を取得する。
出力画素周辺座標生成部203のより詳細な機能構成例について、図3のブロック図を用いて説明する。出力座標量子化部206は、出力画素座標値222(xo、yo)に対して以下の式を適用することで、xoの整数部xoc、yoの整数部yocを算出する。
xoc=floor(xo+0.5)
yoc=floor(yo+0.5)
なお、floor(x)関数は、x以下の最大の整数を返す関数である。そして出力座標量子化部206は、画素位置(xoc、yoc)を、出力画素座標整数部231として、後段のサブ領域判定部207に対して出力する。
一方、出力座標量子化部206は、出力画素座標値222(xo、yo)に対して以下の式を適用することで、xoの小数部xof、yoの小数部yofを算出する。ここでは、後の判定を0≦xof<1,0≦yof<1で行うためにxo,yoからそれぞれの整数部を減じた後に0.5を加えているが、判定の都合によってxof,yofの算出式を変えることもできる。
xof=xo−floor(xo+0.5)+0.5
yof=yo−floor(yo+0.5)+0.5
そして出力座標量子化部206は、このようにして求めたxof、yofを、出力画素座標小数部233として、後段のサブ領域判定部207に対して出力する。
サブ領域判定部207は、xofが既定値以上であるか否か及び/又はyofが既定値以上であるか否かの組み合わせを特定する。ここで、サブ領域判定部207の動作について、図4を用いて説明する。図4のグラフにおいて横軸はxof、縦軸はyofを示す。
図4において領域301は、次のような条件を満たす領域である。先ず、画素位置(xoc、yoc)の右に隣接する画素位置(xoc+1、yoc)に対して上記の射影変換の逆変換を施すことで得られる画素位置が(p1、q1)であったとする。このとき、p1の整数部がxiと一致し、且つq1の整数部がyiと一致するような、xof、yofの範囲が領域301である。然るに、上記の式で求めたxof、yofによって表される位置(xof、yof)が領域301内に位置する場合、画素位置(xoc+1、yoc)に対して上記の射影変換の逆変換を施すことで得られるp1、q1のそれぞれの整数部はxi、yiと一致する。
以下に上記の射影変換の逆変換を行い(po,qo)から(pi,qi)を求める式と、偶数次フィルタで補間係数kを0≦k<1の範囲とする時の(pi,qi)から整数部(pic,qic)を算出する式を示す。なお、奇数次フィルタで補間係数kを−0.5≦k<0.5とするときはpi,qiに0.5を加えて整数部(pic,qic)を求める。
pi0=r11・po+r12・qo+r13
qi0=r21・po+r22・qo+r23
ri0=r31・po+r32・qo+r33
pi=pi0/ri0
qi=qi0/ri0
但し、(r11−r33)は(m11−m33)の逆行列である。射影変換に限定すれば、(r11−r33)は(m11−m33)の逆行列の定数倍であってもよい。
pic=floor(pi)
qic=floor(qi)
同様に、画素位置(xoc、yoc)の下に隣接する画素位置(xoc、yoc+1)に対して上記の射影変換の逆変換を施すことで得られる画素位置が(p2、q2)であったとする。このとき、p2の整数部がxiと一致し、且つq2の整数部がyiと一致するような、xof、yofの範囲が領域302である。然るに、上記の式で求めたxof、yofによって表される位置(xof、yof)が領域302内に位置する場合、画素位置(xoc、yoc+1)に対して上記の射影変換の逆変換を施すことで得られるp2、q2のそれぞれの整数部はxi、yiと一致する。
同様に、画素位置(xoc、yoc)の左下に隣接する画素位置(xoc−1、yoc+1)に対して上記の射影変換の逆変換を施すことで得られる画素位置が(p3、q3)であったとする。このとき、p3の整数部がxiと一致し、且つq3の整数部がyiと一致するような、xof、yofの範囲が領域311である。然るに、上記の式で求めたxof、yofによって表される位置(xof、yof)が領域311内に位置する場合、画素位置(xoc−1、yoc+1)に上記の射影変換の逆変換を施すことで得られるp3、q3のそれぞれの整数部はxi、yiと一致する。
同様に、画素位置(xoc、yoc)の左に隣接する画素位置(xoc−1、yoc)に対して上記の射影変換の逆変換を施すことで得られる画素位置が(p4、q4)であったとする。このとき、p4の整数部がxiと一致し、且つq4の整数部がyiと一致するような、xof、yofの範囲が領域303である。然るに、上記の式で求めたxof、yofによって表される位置(xof、yof)が領域303内に位置する場合、画素位置(xoc−1、yoc)に対して上記の射影変換の逆変換を施すことで得られるp4、q4のそれぞれの整数部はxi、yiと一致する。
同様に、画素位置(xoc、yoc)の上に隣接する画素位置(xoc、yoc−1)に対して上記の射影変換の逆変換を施すことで得られる画素位置が(p5、q5)であったとする。このとき、p5の整数部がxiと一致し、且つq5の整数部がyiと一致するような、xof、yofの範囲が領域304である。然るに、上記の式で求めたxof、yofによって表される位置(xof、yof)が領域304内に位置する場合、画素位置(xoc、yoc−1)に対して上記の射影変換の逆変換を施すことで得られるp5、q5のそれぞれの整数部はxi、yiと一致する。
同様に、画素位置(xoc、yoc)の右上に隣接する画素位置(xoc+1、yoc−1)に対して上記の射影変換の逆変換を施すことで得られる画素位置が(p6、q6)であったとする。このとき、p6の整数部がxiと一致し、且つq6の整数部がyiと一致するような、xof、yofの範囲が領域313である。然るに、上記の式で求めたxof、yofによって表される位置(xof、yof)が領域313内に位置する場合、画素位置(xoc+1、yoc−1)に上記の射影変換の逆変換を施すことで得られるp6、q6のそれぞれの整数部はxi、yiと一致する。
なお、図4では、画素位置(xoc−1、yoc−1)、(xoc+1、yoc+1)に対して上記の射影変換の逆変換を施すことで得られる画素位置は、x座標値の整数部=xi且つy座標値の整数部=yiという条件は満たしていないことになっている。
然るに、図4を例に取り説明すると、サブ領域判定部207は、出力画素座標小数部233としてのxofが、0≦xof<0.5の範囲(サブ領域)に属しているのかそれとも0.5≦xof<1の範囲(サブ領域)に属しているのかを判定する。そしてサブ領域判定部207は、この判定の結果を、後段のサブ領域周辺座標出力部208に対して通知する。
サブ領域周辺座標出力部208は、出力画素座標小数部233としてのxofが、0≦xof<0.5の範囲(サブ領域)に属している旨の通知を受けると、以下の画素位置を参照画素位置として決定する。
(xoc 、yoc+1)
(xoc−1、yoc )
(xoc 、yoc−1)
(xoc−1、yoc+1)
一方、サブ領域周辺座標出力部208は、出力画素座標小数部233としてのxofが、0.5≦xof<1.0の範囲(サブ領域)に属している旨の通知を受けると、以下の画素位置を参照画素位置として決定する。
(xoc+1、yoc )
(xoc 、yoc+1)
(xoc 、yoc−1)
(xoc+1、yoc−1)
そしてサブ領域周辺座標出力部208は、上記のように、xofの属する範囲に応じて決まる参照画素位置と、画素位置(xoc、yoc)と、を出力画素周辺座標値223として出力する。
なお、xofの属する範囲に応じて参照画素位置を決定する方法は特に限定するものではなく、様々な方法が考え得る。例えば、xofの属する範囲毎に、対応する参照画素位置を登録したテーブルを用いて、上記の求めたxocに対応する参照画素位置をこのテーブルから取得しても良い。
座標逆変換部204は、出力画素周辺座標値223としてのそれぞれの画素位置に対して、上記の射影変換の逆変換を施すことで、出力画素周辺座標値223としてのそれぞれの画素位置に対応する逆変換座標値を求める。ここで、xofの属する範囲に応じて決まる参照画素位置に対して、上記の射影変換の逆変換を施すことで得られる画素位置を(xi’、yi’)と表記する。また、画素位置(xoc、yoc)に対して上記の射影変換の逆変換を施す(計算する)ことで得られる画素位置を(xi”、yi”)と表記する。
そして座標逆変換部204は、この求めたそれぞれの逆変換座標値(xi’、yi’)、(xi”、yi”)を入力比較座標値224として、後段の入力座標比較部205に対して送出する。
入力座標比較部205は、画素位置(xi”、yi”)、それぞれの(xi’、yi’)、のうち、x座標値の整数部及びy座標値の整数部がそれぞれxi、yiと同じである画素位置を特定する。そして入力座標比較部205は、特定した画素位置のx座標値の小数部及びy座標値の小数部を小数部127として画素補間処理部102に送出し、画素位置(xoc、yoc)128はメモリ書き込み部103に対して送出する。これにより画素補間処理部102は、上記の補間処理を行うことで、出力画像上の画素位置(xoc、yoc)における画素値を求めることができる。また、メモリ書き込み部103は、画素位置(xoc、yoc)に対応する記憶部104上のアドレスに、画素位置(xoc、yoc)における画素値を格納することができる。
<変形例>
サブ領域判定部207の動作の変形例について、図5を用いて説明する。図5のグラフにおいて横軸はxof、縦軸はyofを示す。図5に示した参照番号が指し示す対象は、形状こそ違うものの、基本的な事項については図4と同じである。
図5において領域312は、次のような条件を満たす領域である。先ず、画素位置(xoc、yoc)の左上に隣接する画素位置(xoc−1、yoc−1)に対して上記の射影変換の逆変換を施すことで得られる画素位置が(p7、q7)であったとする。このとき、p7の整数部がxiと一致し、且つq7の整数部がyiと一致するような、xof、yofの範囲が領域312である。然るに、上記の式で求めたxof、yofによって表される位置(xof、yof)が領域312内に位置する場合、画素位置(xoc−1、yoc−1)に上記の射影変換の逆変換を施すことで得られるp7、q7のそれぞれの整数部はxi、yiと一致する。
同様に、画素位置(xoc、yoc)の右下に隣接する画素位置(xoc+1、yoc+1)に対して上記の射影変換の逆変換を施すことで得られる画素位置が(p8、q8)であったとする。このとき、p8の整数部がxiと一致し、且つq8の整数部がyiと一致するような、xof、yofの範囲が領域310である。然るに、上記の式で求めたxof、yofによって表される位置(xof、yof)が領域310内に位置する場合、画素位置(xoc+1、yoc+1)に上記の射影変換の逆変換を施すことで得られるp8、q8のそれぞれの整数部はxi、yiと一致する。
本変形例の場合、サブ領域周辺座標出力部208は、出力画素座標小数部233としてのxof、yofがそれぞれ、0.0≦xof<0.5、0.0≦yof<0.5の範囲に属している旨の通知を受けると、以下の画素位置を参照画素位置として決定する。
(xoc−1、yoc )
(xoc 、yoc−1)
(xoc−1、yoc−1)
また、サブ領域周辺座標出力部208は、出力画素座標小数部233としてのxof、yofがそれぞれ、0.5≦xof<1.0、0.0≦yof<0.5の範囲に属している旨の通知を受けると、以下の画素位置を参照画素位置として決定する。
(xoc+1、yoc )
(xoc 、yoc−1)
(xoc+1、yoc−1)
また、サブ領域周辺座標出力部208は、出力画素座標小数部233としてのxof、yofがそれぞれ、0.0≦xof<0.5、0.5≦yof<1.0の範囲に属している旨の通知を受けると、以下の画素位置を参照画素位置として決定する。
(xoc−1、yoc )
(xoc 、yoc+1)
(xoc−1、yoc+1)
また、サブ領域周辺座標出力部208は、出力画素座標小数部233としてのxof、yofがそれぞれ、0.5≦xof<1.0、0.5≦yof<1.0の範囲に属している旨の通知を受けると、以下の画素位置を参照画素位置として決定する。
(xoc 、yoc+1)
(xoc+1、yoc )
(xoc+1、yoc+1)
このように、xoの小数部に0.5を加えた値の小数部が既定値以上であるか否か及び/又はyoの小数部に0.5を加えた値の小数部が既定値以上であるか否かの組み合わせについては様々なものがあり、第1の実施形態や本変形例に限定するものではない。そして何れにせよ、それぞれの組み合わせには、予め設定された複数の参照画素位置が対応付けられているので、どのような組み合わせを想定しても、対応する参照画素位置群は一意に特定することができる。
以上の説明により、本変形例を含めた本実施形態によれば、座標変換(座標逆変換)が必要な周辺画素の数を削減することができる。必要となる座標変換数を従来より少なくすることができるため、回路量を従来の書き込み時座標変換方式より少なくすることができる。
また、必要なメモリ帯域のピーク値が小さいまま、ハードウェア実装したときの回路量を小さくすることができる、あるいはソフトウェア実装した時に処理性能が高い射影変換を行うことができる。
[第2の実施形態]
本実施形態に係る画像処理装置の機能構成例について、図6のブロック図を用いて説明する。図6において図1と同じ構成用件については同じ参照番号を付けており、その説明は省略する。本実施形態に係る画像処理装置は、図1に示した構成において、出力座標算出部106を出力座標算出部606に置き換えたものであり、この出力座標算出部606には、上記のマトリクスのデータである変形パラメータ129が入力される。然るに、出力座標算出部606以外の構成用件については第1の実施形態で説明したとおりであるため、以下ではこの出力座標算出部606について説明する。
尚、以下の説明では、第1の実施形態でxo、yo、xoc、yocと表したものをそれぞれxot、yot、xotC、yotCと表記する。即ち、変数名が異なるだけで、指し示すものは同じであり、例えば、xoとxotとは指し示すものは同じである。
出力座標算出部606の機能構成例について、図7のブロック図を用いて説明する。図7において図2と同じ構成用件については同じ参照番号を付けており、その説明は省略する。
座標変換部202は、第1の実施形態で説明したマトリクスのデータとしての変形パラメータ129を受けると、この変形パラメータ129を用いて第1の実施形態と同様にして、入力画素座標値221に対する射影変換を行う。そして座標変換部202は、この射影変換により得られた座標値を出力画素座標値222として、後段の出力画素周辺座標生成部203に対して送出する。
判定条件生成部703は、変形パラメータ129を受けると、この変形パラメータ129を用いて、参照画素位置を決定するための条件を生成する。そして判定条件生成部703は、この生成した条件を、判定条件704として出力画素周辺座標生成部203に対して供給する。判定条件生成部703の動作について、詳しくは後述する。
座標逆変換部204は、第1の実施形態で説明したマトリクスのデータとしての変形パラメータ129を受けると、この変形パラメータ129を用いて第1の実施形態と同様にして、出力画素周辺座標値223に対して射影変換の逆変換を行う。そして座標逆変換部204は、この逆変換により得られた逆変換座標値を入力画素周辺座標値701として、後段の座標出力部702に対して送出する。
座標出力部702は、入力画素周辺座標値701のx座標値の小数部及びy座標値の小数部を小数部127として画素補間処理部102に送出する。更に座標出力部702は、画素位置(xotC、yotC)を整数部128としてメモリ書き込み部103に送出する。
出力座標算出部606は、一例として画像変形を行わない場合は、1つの入力画素値121に対して、1つの整数部128及び1つの小数部127を出力する。具体的には、変形パラメータ129がm11〜m33からなる3x3の単位行列のときがこのような場合に当たる。また、その他の例として、画像変形が拡大部分を含むときは、1つの入力画素値121に対して、複数の整数部128及び小数部127が出力される場合がある。
出力画素周辺座標生成部203は先ず、xoの小数部xotf、yoの小数部yotfを以下の式により求める。
xotf=xot−floor(xot+0.5)+0.5
yotf=yot−floor(yot+0.5)+0.5
そして出力画素周辺座標生成部203は、xotf及びyotfが、判定条件704を満たす場合には、この判定条件704に対応する複数の参照画素位置を取得する。そして出力画素周辺座標生成部203は、この取得した複数の参照画素位置と、画素位置(xotC、yotC)と、を出力画素周辺座標値223として出力する。ここで、出力画素周辺座標生成部203が出力画素周辺座標値223として出力する可能性があるものは、以下の9個の座標値である。
(xotC−1,yotC−1)
(xotC+0,yotC−1)
(xotC+1,yotC−1)
(xotC−1,yotC+0)
(xotC+0,yotC+0)
(xotC+1,yotC+0)
(xotC−1,yotC+1)
(xotC+0,yotC+1)
(xotC+1,yotC+1)
図8のグラフにおいて横軸はxotf、縦軸はyotfを示す。図8において領域801は、次のような条件を満たす領域である。先ず、画素位置(xotC、yotC)の左上に隣接する画素位置(xotC−1、yotC−1)に対して上記の射影変換の逆変換を施すことで得られる画素位置が(p1、q1)であったとする。このとき、p1の整数部がxiと一致し、且つq1の整数部がyiと一致するような、xotf、yotfの範囲が領域801である。然るに、求めたxotf、yotfによって表される位置(xotf、yotf)が領域801内に位置する場合、画素位置(xotC−1、yotC−1)に上記の射影変換の逆変換を施して得られるp1、q1のそれぞれの整数部はxi、yiと一致する。これは、他の領域800,802〜805,808についても同様である。なお、それぞれの領域は重なる場合があり、図8では見易さの都合上、4つに分けて図示している。また、図8では、出力画素周辺座標生成部203が出力画素周辺座標値223として出力する可能性があるものは、以下の7個の座標値としている。
(xotC−1,yotC−1)
(xotC+0,yotC−1)
(xotC−1,yotC+0)
(xotC+0,yotC+0)
(xotC+1,yotC+0)
(xotC+0,yotC+1)
(xotC+1,yotC+1)
判定条件生成部703は、それぞれの領域800〜805,808を規定するパラメータを判定条件704として出力画素周辺座標生成部203に対して送出する。例えば、図8に示す如く、領域801は、直線A806と直線B807とxotf軸とyotf軸とによって囲まれている。直線A806と直線B807は以下のとおり表わすことができる。
y=A1・x+A2 (直線A)
y=B1・x+B2 (直線B)
然るに、判定条件生成部703は、領域801を規定するパラメータとして、以下の条件を示すパラメータを生成する。
0.0 ≦ xotf ≦ 1.0
0.0 ≦ yotf ≦ 1.0
yotf ≧ A1・xotf+A2 (直線A)
yotf ≧ B1・xotf+B2 (直線B)
また判定条件生成部703は、他の領域800,802〜805,808についてもこのような、それぞれの領域を規定するパラメータを生成する。そして判定条件生成部703は、それぞれの領域800〜805,808について生成したパラメータをまとめて判定条件704として出力画素周辺座標生成部203に対して送出する。
ここで、各領域を規定するためのパラメータを求めるための方法について説明する。各領域は、変形パラメータ129に依存しているため、各領域を規定するためのパラメータは、この変形パラメータ129を用いて生成されることになる。具体的には、各隣接画素座標値に対応する領域を囲む直線式の傾き係数と切片係数を求める。
四角形が記されている図9(a)の入力画像(xi−yi平面)に対して変形パラメータ129を用いた画像変形処理を行うことで、変形した四角形が記されている図9(b)の出力画像(xot−yot平面)が得られたとする。図9(a)の四角形は頂点(入力端点)902〜905で構成されており、図9(b)の四角形は頂点(出力端点)906〜909で構成されている。
各領域の求め方について、図10を用いて説明する。図9(a)に示した入力端点902〜905のうち1つを入力端点505として選び、画像変形処理による入力端点505の変換後の端点を出力端点507として求める。出力端点507の画素位置を(Axo,Ayo)とすると、画素位置(floor(Axo)+0.5、floor(Ayo)+0.5)に出力中心点506を設定する。
次に、画像変形処理の逆変換処理による出力中心点506の変換後の点を入力中心点508として求める。入力中心点508の画素位置を(Acxi,Acyi)とすると、以下の4つの画素位置にそれぞれ入力近傍点509〜512を設定する。
(Acxi−0.5,Acyi+0.5)
(Acxi+0.5,Acyi+0.5)
(Acxi+0.5,Acyi+0.5)
(Acxi−0.5,Acyi−0.5)
次に、画像変形処理による入力近傍点509〜512の変換後の点をそれぞれ出力近傍点514〜517として求める。そして出力近傍点514〜517で囲まれている領域のうち、0.0≦xot<1.0且つ0.0≦yot<1.0を満たす領域内を、判定領域518として求める。
以上の判定領域518を求めるための処理を、入力端点902〜905のそれぞれについて行い、それぞれの領域について求めた判定領域518の和集合としての領域を、画素位置(xotC,yotC)に対する領域として出力する。同様の処理を、出力座標系501において出力中心点506を選択するとき、±1.0の8つの隣接画素座標に対して実行し、全ての領域を求めることができる。
出力画素周辺座標生成部203はこの判定条件704を受けると、判定条件704に含まれているそれぞれのパラメータを参照する。そして出力画素周辺座標生成部203は、上記の式で求めたxotf、yotfによって表される位置(xotf、yotf)が領域800〜805,808の何れに属しているのかを判定する。そして出力画素周辺座標生成部203は、上記の式で求めたxotf、yotfによって表される位置(xotf、yotf)が属していると判定した領域に対応する画素位置(参照画素位置)を出力する。
例えば、上記の式で求めたxotf、yotfによって表される位置(xotf、yotf)が、領域801内に位置している場合、画素位置(xotC−1、yotC−1)を参照画素位置として出力する。このようにして出力画素周辺座標生成部203は、上記の式で求めたxotf、yotfによって表される位置(xotf、yotf)が属していると判定した領域に対応する画素位置(参照画素位置)を出力する。
このように、本実施形態では、変形パラメータ129をもとに判定条件704を生成し、判定条件704にもとづき入力画素値121に対応する出力画素周辺座標値223を選択する。拡大率が2倍未満のとき、出力画素周辺座標値223の候補は9点だが、ここまでの説明のとおり、判定条件704を用いることで、出力画素周辺座標値223を4点以下に削減できる。この結果、座標逆変換部204による逆変換回数を削減できる。必要となる座標変換数を従来より少なくすることができるため、回路量を従来の書き込み時座標変換方式より少なくすることができる。また、必要なメモリ帯域のピーク値が小さいまま、ハードウェア実装したときの回路量が小さい、あるいはソフトウェア実装した時に処理性能が高い射影変形を行う画像処理技術を実現することができる。なお、本実施形態は、拡大率が任意の倍率の場合にも近傍点数を増やすことで同様に拡張可能である。
また、図1〜3,6,7に示した各部のうち、メモリとして機能するもの以外については、ソフトウェア(コンピュータプログラム)で実現させても良い。この場合、このメモリを有し、且つこのソフトウェアをインストールしたコンピュータは、このコンピュータが有するCPUがこのソフトウェアを実行することで、上記の画像処理装置として機能することになる。もちろん、図1〜3,6,7に示した全ての構成用件をハードウェアで実現させても良いし、一部をソフトウェアとして実現させても良い。
(その他の実施例)
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア(プログラム)を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ(またはCPUやMPU等)がプログラムを読み出して実行する処理である。

Claims (8)

  1. 入力画像が変形処理された出力画像を生成する画像処理装置であって、
    前記入力画像を取得する取得手段と、
    前記入力画像上の第1座標に対応する前記出力画像上の第2座標の整数部と小数部を、前記変形処理で用いられるパラメータによって特定する特定手段と、
    前記第2座標の整数部に対応する前記出力画像上の座標から規定範囲内に位置する前記出力画像上の第1複数座標から、該第1複数座標よりも数が少ない第2複数座標を、該第2座標の小数部に基づいて選択する選択手段と、
    前記第2複数座標に対して、前記パラメータに基づく変換の逆変換を行うことで、前記出力画像上の前記第2座標の整数部に対応する前記出力画像の画素データを生成するために参照される前記入力画像の画素データに対応する第3複数座標を判定する判定手段と、
    前記判定手段により判定された第3複数座標に対応する前記入力画像の画素データを読み出して、前記第2座標の整数部に対応する前記出力画像の画素データを生成する生成手段と、
    前記生成手段により生成された画素データを前記出力画像の画素データとして出力する出力手段と
    を有することを特徴とする画像処理装置。
  2. 前記生成手段は、前記入力画像の前記第3複数座標に対応する画素データを用いた補間処理を実行することで、前記第2座標の整数部に対応する前記出力画像の画素データを生成することを特徴とする請求項1に記載の画像処理装置。
  3. 前記選択手段は、前記出力画像上の前記第2座標の小数部が第1閾値よりも大きい場合は、前記第2座標より座標値が高い座標を前記第2複数座標の1つとして選択し且つ前記第2座標より座標値が低い座標を前記第2複数座標として選択せず、前記出力画像上の前記第2座標の小数部が前記第1閾値よりも小さい第2閾値よりも小さい場合、前記第2座標より座標値が低い座標を前記第2複数座標の1つとして選択し且つ前記第2座標より座標値が高い座標を前記第2複数座標として選択しないことを特徴とする請求項1又は2に記載の画像処理装置。
  4. 前記特定手段は、前記第2座標を表す第1座標成分値と第2座標成分値のそれぞれについて整数部及び小数部を特定し、
    前記選択手段は、前記第1座標成分値及び前記第2座標成分値のそれぞれの小数部に基づいて、前記第2複数座標を選択することを特徴とする請求項3に記載の画像処理装置。
  5. 入力画像が変形処理された出力画像を生成する画像処理装置が行う画像処理方法であって、
    前記画像処理装置の取得手段が、前記入力画像を取得する取得工程と、
    前記画像処理装置の特定手段が、前記入力画像上の第1座標に対応する前記出力画像上の第2座標の整数部と小数部を、前記変形処理で用いられるパラメータによって特定する特定工程と、
    前記画像処理装置の選択手段が、前記第2座標の整数部に対応する前記出力画像上の座標から規定範囲内に位置する前記出力画像上の第1複数座標から、該第1複数座標よりも数が少ない第2複数座標を、該第2座標の小数部に基づいて選択する選択工程と、
    前記画像処理装置の判定手段が、前記第2複数座標に対して、前記パラメータに基づく変換の逆変換を行うことで、前記出力画像上の前記第2座標の整数部に対応する前記出力画像の画素データを生成するために参照される前記入力画像の画素データに対応する第3複数座標を判定する判定工程と、
    前記画像処理装置の生成手段が、前記判定工程で判定された第3複数座標に対応する前記入力画像の画素データを読み出して、前記第2座標の整数部に対応する前記出力画像の画素データを生成する生成工程と、
    前記画像処理装置の出力手段が、前記生成工程で生成された画素データを前記出力画像の画素データとして出力する出力工程と
    を有することを特徴とする画像処理方法。
  6. 前記生成工程では、前記入力画像の前記第3複数座標に対応する画素データを用いた補間処理を実行することで、前記第2座標の整数部に対応する前記出力画像の画素データを生成することを特徴とする請求項5に記載の画像処理方法。
  7. 前記選択工程では、前記出力画像上の前記第2座標の小数部が第1閾値よりも大きい場合は、前記第2座標より座標値が高い座標を前記第2複数座標の1つとして選択し且つ前記第2座標より座標値が低い座標を前記第2複数座標として選択せず、前記出力画像上の前記第2座標の小数部が前記第1閾値よりも小さい第2閾値よりも小さい場合、前記第2座標より座標値が低い座標を前記第2複数座標の1つとして選択し且つ前記第2座標より座標値が高い座標を前記第2複数座標として選択しないことを特徴とする請求項5又は6に記載の画像処理方法。
  8. コンピュータを、請求項1乃至4の何れか1項に記載の画像処理装置の各手段として機能させるためのコンピュータプログラム。
JP2011137737A 2011-06-21 2011-06-21 画像処理装置、画像処理方法 Expired - Fee Related JP5787637B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2011137737A JP5787637B2 (ja) 2011-06-21 2011-06-21 画像処理装置、画像処理方法
US13/490,260 US8509568B2 (en) 2011-06-21 2012-06-06 Image processing apparatus and image processing method
US13/957,375 US20130315502A1 (en) 2011-06-21 2013-08-01 Image processing apparatus and image processing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011137737A JP5787637B2 (ja) 2011-06-21 2011-06-21 画像処理装置、画像処理方法

Publications (3)

Publication Number Publication Date
JP2013004031A JP2013004031A (ja) 2013-01-07
JP2013004031A5 JP2013004031A5 (ja) 2014-07-31
JP5787637B2 true JP5787637B2 (ja) 2015-09-30

Family

ID=47361919

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011137737A Expired - Fee Related JP5787637B2 (ja) 2011-06-21 2011-06-21 画像処理装置、画像処理方法

Country Status (2)

Country Link
US (2) US8509568B2 (ja)
JP (1) JP5787637B2 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6278716B2 (ja) * 2014-01-22 2018-02-14 キヤノン株式会社 画像処理装置、画像処理方法及びプログラム
WO2018168506A1 (ja) * 2017-03-13 2018-09-20 Necソリューションイノベータ株式会社 画像変換支援装置、画像変換装置、画像変換支援方法、及びコンピュータ読み取り可能な記録媒体
CN116309058B (zh) * 2023-03-15 2024-04-19 长沙观谱红外科技有限公司 一种人体红外图像的放大处理方法

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6256073A (ja) * 1985-09-04 1987-03-11 Canon Inc 画像処理方法
US5204944A (en) * 1989-07-28 1993-04-20 The Trustees Of Columbia University In The City Of New York Separable image warping methods and systems using spatial lookup tables
JP3394551B2 (ja) * 1992-11-09 2003-04-07 松下電器産業株式会社 画像変換処理方法及び画像変換処理装置
DE69533870T2 (de) * 1994-10-19 2005-05-25 Matsushita Electric Industrial Co., Ltd., Kadoma Vorrichtung zur Bilddekodierung
US5870105A (en) * 1996-05-31 1999-02-09 Hewlett-Packard Company System and method for local storage of image data during object to image mapping
US6317158B1 (en) * 1998-10-23 2001-11-13 Avid Technology, Inc. Method and apparatus for positioning an input image into interlaced video
JP3543928B2 (ja) * 1998-10-29 2004-07-21 三菱電機株式会社 画素数変換装置
JP2001266138A (ja) * 2000-03-23 2001-09-28 Seiko Epson Corp 画像変換処理装置および画像変換処理方法ならびに情報記録媒体
US7729563B2 (en) * 2002-08-28 2010-06-01 Fujifilm Corporation Method and device for video image processing, calculating the similarity between video frames, and acquiring a synthesized frame by synthesizing a plurality of contiguous sampled frames
JPWO2004034326A1 (ja) * 2002-10-08 2006-02-09 ソニー株式会社 画像変換装置、画像変換方法、および、画像投射装置
JP4124096B2 (ja) 2003-10-29 2008-07-23 株式会社ニコン 画像処理方法および画像処理装置、並びにプログラム
KR100595254B1 (ko) * 2004-01-17 2006-07-03 엘지전자 주식회사 좌표 변환 장치 및 방법
US7324706B2 (en) * 2004-09-09 2008-01-29 Silicon Optix Inc. System and method for representing a general two dimensional spatial transformation
US20060204125A1 (en) * 2005-03-09 2006-09-14 Kempf Jeffrey M Multi-dimensional keystone correction image projection system and method
US8379066B2 (en) * 2007-06-27 2013-02-19 Christie Digital Systems Usa, Inc. Method and apparatus for scaling an image to produce a scaled image
JP5408906B2 (ja) 2008-05-28 2014-02-05 キヤノン株式会社 画像処理装置
JP2010081024A (ja) * 2008-09-24 2010-04-08 Oki Semiconductor Co Ltd 画像補間処理装置
JP5585885B2 (ja) * 2011-03-18 2014-09-10 株式会社リコー 画像処理装置及び画像処理方法
JP5997480B2 (ja) * 2012-03-30 2016-09-28 キヤノン株式会社 画像処理装置、画像処理方法およびプログラム

Also Published As

Publication number Publication date
US20120328214A1 (en) 2012-12-27
US8509568B2 (en) 2013-08-13
JP2013004031A (ja) 2013-01-07
US20130315502A1 (en) 2013-11-28

Similar Documents

Publication Publication Date Title
JP5914045B2 (ja) 画像処理装置、画像処理方法、及びプログラム
JP2019120749A (ja) 表示制御装置、画像投影システム、制御方法およびプログラム
US8928782B2 (en) Image processing device and image capture device
US9332238B2 (en) Image processing apparatus and image processing method
JP5787637B2 (ja) 画像処理装置、画像処理方法
JP2017191572A (ja) 画像処理装置及びその方法、プログラム
JP2012238171A (ja) 画像処理装置、画像処理方法、及びプログラム
JP2018026007A (ja) 画像処理装置、画像処理方法及びプログラム
JP6487671B2 (ja) 画像処理装置、画像処理方法。およびプログラム
JP6632434B2 (ja) 画像処理装置、画像処理方法、及びプログラム
JP6278716B2 (ja) 画像処理装置、画像処理方法及びプログラム
JP5387276B2 (ja) 画像処理装置及び画像処理方法
JP5387289B2 (ja) 画像処理装置、画像処理方法およびそのプログラム
JP6972089B2 (ja) 画像処理装置、画像処理方法、及びプログラム
US11893706B2 (en) Image correction device
JP6285657B2 (ja) 画像処理装置及び画像処理方法、プログラム
JP2013257665A (ja) 映像処理装置、映像処理装置の制御方法
JP2012019315A (ja) 画像処理装置及び画像処理方法
JP2019106589A (ja) プロジェクター、幾何学歪補正方法およびプログラム
JP5085589B2 (ja) 画像処理装置および方法
JP6562687B2 (ja) 画像変形装置、画像処理方法、およびプログラム
JP2015154209A (ja) 画像処理装置、画像処理方法、及びプログラム
JP6080728B2 (ja) 歪補正装置、歪補正方法及び歪補正プログラム
JP2004227470A (ja) 画像変換方法及び画像変換用コンピュータプログラム
JP2015197862A (ja) 画像処理装置、画像処理方法

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140616

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20140616

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150212

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150306

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150507

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20150728

R151 Written notification of patent or utility model registration

Ref document number: 5787637

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

LAPS Cancellation because of no payment of annual fees