JP2001352549A - Image coding and decoding method - Google Patents

Image coding and decoding method

Info

Publication number
JP2001352549A
JP2001352549A JP2001104091A JP2001104091A JP2001352549A JP 2001352549 A JP2001352549 A JP 2001352549A JP 2001104091 A JP2001104091 A JP 2001104091A JP 2001104091 A JP2001104091 A JP 2001104091A JP 2001352549 A JP2001352549 A JP 2001352549A
Authority
JP
Japan
Prior art keywords
integer
horizontal
image
motion vector
integers
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2001104091A
Other languages
Japanese (ja)
Other versions
JP3777998B2 (en
Inventor
Yuichiro Nakaya
雄一郎 中屋
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.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
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 Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2001104091A priority Critical patent/JP3777998B2/en
Publication of JP2001352549A publication Critical patent/JP2001352549A/en
Application granted granted Critical
Publication of JP3777998B2 publication Critical patent/JP3777998B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Abstract

PROBLEM TO BE SOLVED: To provide a method for simplifying arithmetic operations in the processing of global motion compensation where a motion vector field of the entire image is approximated by a few parameters. SOLUTION: Interpolating/extrapolating motion vectors at representative points 602, 603 and 604 having features in a spatial interval obtain motion vectors in global motion compensation. Since division used to composite predicted images in the global motion compensation can be substituted by shift arithmetic operations, the processing by a computer or an exclusive hardware unit can be simplified.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、画像全体に対して
線形内・外挿または共1次内・外挿に基づくグローバル
動き補償を適用する画像符号化および復号化方法に関す
るものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an image encoding and decoding method for applying global motion compensation based on linear interpolation / extrapolation or bilinear interpolation / extrapolation to an entire image.

【0002】[0002]

【従来の技術】動画像の高能率符号化において、時間的
に近接するフレーム間の類似性を活用する動き補償は情
報圧縮に大きな効果を示すことが知られている。現在の
画像符号化技術の主流となっている動き補償方式は、動
画像符号化方式の国際標準であるH.261,MPEG
1,MPEG2に採用されているブロックマッチングで
ある。この方式では、符号化しようとする画像を多数の
ブロックに分割し、ブロックごとにその動きベクトルを
求める。
2. Description of the Related Art It is known that, in high-efficiency coding of moving images, motion compensation utilizing similarity between temporally adjacent frames has a great effect on information compression. The motion compensation method which is the mainstream of the current image coding technology is H.261, MPEG which is an international standard of the moving image coding method.
1, block matching adopted in MPEG2. In this method, an image to be encoded is divided into a number of blocks, and a motion vector is obtained for each block.

【0003】図1にH.261の符号化器の構成例10
0を示す。H.261は、符号化方式として、ブロック
マッチングとDCT(離散コサイン変換)を組み合わせ
たハイブリッド符号化方式(フレーム間/フレーム内適
応符号化方式)を採用している。減算器102は入力画
像(現フレームの原画像)101とフレーム間/フレー
ム内符号化切り換えスイッチ119の出力画像113
(後述)との差を計算し、誤差画像103を出力する。
この誤差画像は、DCT変換器104でDCT係数に変
換された後に量子化器105で量子化され、量子化DC
T係数106となる。この量子化DCT計数は伝送情報
として通信路に出力されると同時に、符号化器内でもフ
レーム間予測画像を合成するために使用される。以下に
予測画像合成の手順を説明する。上述の量子化DCT係
数106は、逆量子化器108と逆DCT変換器109
を経て復号誤差画像110(受信側で再生される誤差画
像と同じ画像)となる。これに、加算器111において
フレーム間/フレーム内符号化切り換えスイッチ119
の出力画像113(後述)が加えられ、現フレームの復
号画像112(受信側で再生される現フレームの復号画
像と同じ画像)を得る。この画像は一旦フレームメモリ
114に蓄えられ、1フレーム分の時間だけ遅延され
る。したがって、現時点では、フレームメモリ114は
前フレームの復号画像115を出力している。この前フ
レームの復号画像と現フレームの入力画像101がブロ
ックマッチング部116に入力され、ブロックマッチン
グの処理が行われる。ブロックマッチングでは、画像を
複数のブロックに分割し、各ブロックごとに現フレーム
の原画像に最も似た部分を前フレームの復号画像から取
り出すことにより、現フレームの予測画像117が合成
される。このときに、各ブロックが前フレームと現フレ
ームの間でどれだけ移動したかを検出する処理(動き推
定処理)を行う必要がある。動き推定処理によって検出
された各ブロックごとの動きベクトルは、動き情報12
0として受信側へ伝送される。受信側は、この動き情報
と前フレームの復号画像から、独自に送信側で得られる
ものと同じ予測画像を合成することができる。予測画像
117は、「0」信号118と共にフレーム間/フレー
ム内符号化切り換えスイッチ119に入力される。この
スイッチは、両入力のいずれかを選択することにより、
フレーム間符号化とフレーム内符号化を切り換える。予
測画像117が選択された場合(図2はこの場合を表し
ている)には、フレーム間符号化が行われる。一方、
「0」信号が選択された場合には、入力画像がそのまま
DCT符号化されて通信路に出力されるため、フレーム
内符号化が行われることになる。受信側が正しく復号化
画像を得るためには、送信側でフレーム間符号化が行わ
れたかフレーム内符号化が行われたかを知る必要があ
る。このため、識別フラグ121が通信路へ出力され
る。最終的なH.261符号化ビットストリーム123
は多重化器122で量子化DCT係数,動きベクトル,
フレーム内/フレーム間識別フラグの情報を多重化する
ことによって得られる。
FIG. 1 shows a configuration example 10 of an H.261 encoder.
Indicates 0. H.261 adopts a hybrid coding method (inter-frame / intra-frame adaptive coding method) that combines block matching and DCT (discrete cosine transform) as a coding method. The subtracter 102 outputs the input image (original image of the current frame) 101 and the output image 113 of the inter-frame / intra-frame coding switch 119.
A difference between the error image 103 and an error image 103 is output.
This error image is converted into a DCT coefficient by a DCT converter 104 and then quantized by a quantizer 105 to be quantized DCT.
It becomes the T coefficient 106. The quantized DCT count is output to the communication path as transmission information, and is also used in the encoder to synthesize an inter-frame predicted image. The procedure of predictive image synthesis will be described below. The above-described quantized DCT coefficient 106 is converted into an inverse quantizer 108 and an inverse DCT transformer 109.
, And becomes a decoded error image 110 (the same image as the error image reproduced on the receiving side). In addition, in the adder 111, an inter-frame / intra-frame coding changeover switch 119 is provided.
Of the current frame (the same image as the decoded image of the current frame reproduced on the receiving side) is obtained. This image is temporarily stored in the frame memory 114 and is delayed by a time corresponding to one frame. Therefore, at present, the frame memory 114 outputs the decoded image 115 of the previous frame. The decoded image of the previous frame and the input image 101 of the current frame are input to the block matching unit 116, and the block matching process is performed. In the block matching, the predicted image 117 of the current frame is synthesized by dividing the image into a plurality of blocks and extracting a portion most similar to the original image of the current frame from the decoded image of the previous frame for each block. At this time, it is necessary to perform processing (motion estimation processing) for detecting how much each block has moved between the previous frame and the current frame. The motion vector for each block detected by the motion estimation processing is the motion information 12
It is transmitted to the receiving side as 0. The receiving side can independently synthesize the same predicted image as that obtained on the transmitting side from the motion information and the decoded image of the previous frame. The predicted image 117 is input to the inter-frame / intra-frame coding changeover switch 119 together with the “0” signal 118. This switch selects one of both inputs,
Switching between inter-frame coding and intra-frame coding. When the predicted image 117 is selected (FIG. 2 shows this case), inter-frame encoding is performed. on the other hand,
When the “0” signal is selected, the input image is DCT-coded as it is and output to the communication channel, so that intra-frame coding is performed. In order for the receiving side to correctly obtain a decoded image, it is necessary to know whether inter-frame coding or intra-frame coding has been performed on the transmitting side. Therefore, the identification flag 121 is output to the communication path. Final H.261 encoded bit stream 123
Is a quantized DCT coefficient, a motion vector,
It is obtained by multiplexing the information of the intra-frame / inter-frame identification flag.

【0004】図2に図1の符号化器が出力した符号化ビ
ットストリームを受信する復号化器200の構成例を示
す。受信したH.261ビットストリーム217は、分
離器216で量子化DCT係数201,動きベクトル2
02,フレーム内/フレーム間識別フラグ203に分離
される。量子化DCT係数201は逆量子化器204と
逆DCT変換器205を経て復号化された誤差画像20
6となる。この誤差画像は加算器207でフレーム間/
フレーム内符号化切り換えスイッチ214の出力画像2
15を加算され、復号化画像208として出力される。
フレーム間/フレーム内符号化切り換えスイッチはフレ
ーム間/フレーム内符号化識別フラグ203に従って、
出力を切り換える。フレーム間符号化を行う場合に用い
る予測画像212は、予測画像合成部211において合
成される。ここでは、フレームメモリ209に蓄えられ
ている前フレームの復号画像210に対して、受信した
動きベクトル202に従ってブロックごとに位置を移動
させる処理が行われる。一方フレーム内符号化の場合、
フレーム間/フレーム内符号化切り換えスイッチは、
「0」信号213をそのまま出力する。
FIG. 2 shows a configuration example of a decoder 200 that receives an encoded bit stream output from the encoder of FIG. The received H.261 bit stream 217 is divided by a separator 216 into a quantized DCT coefficient 201 and a motion vector 2.
02, separated into an intra-frame / inter-frame identification flag 203. The quantized DCT coefficient 201 is converted into an error image 20 decoded through an inverse quantizer 204 and an inverse DCT transformer 205.
It becomes 6. This error image is inter-frame /
Output image 2 of intra-frame encoding changeover switch 214
15 is added and the result is output as a decoded image 208.
The interframe / intraframe coding changeover switch is operated in accordance with the interframe / intraframe coding identification flag 203.
Switch the output. A predicted image 212 used for performing inter-frame encoding is synthesized by a predicted image synthesis unit 211. Here, a process of moving the position for each block in accordance with the received motion vector 202 is performed on the decoded image 210 of the previous frame stored in the frame memory 209. On the other hand, in the case of intra-frame coding,
The inter-frame / intra-frame encoding switch is
The "0" signal 213 is output as it is.

【0005】ブロックマッチングは現在最も広く利用さ
れている動き補償方式であるが、画像全体が拡大・縮小
・回転している場合には、すべてのブロックに対して動
きベクトルを伝送しなければならず、符号化効率が悪く
なる問題が発生する。この問題に対し、画像全体の動き
ベクトル場を少ないパラメータを用いて表現するグロー
バル動き補償(例えば、M.Hotter, &#34Differential e
stimation of the global motion parameters zoom and
pan&#34, Signal Processing, vol. 16, no.3, pp. 24
9-265, Mar. 1989)が提案されている。これは、画像内
の画素(x,y)の動きベクトル(ug(x,y),vg(x,
y))を、
[0005] Block matching is currently the most widely used motion compensation method. However, when the entire image is enlarged, reduced, or rotated, motion vectors must be transmitted for all blocks. In this case, a problem that encoding efficiency is deteriorated occurs. To solve this problem, global motion compensation (for example, M.Hotter, &# 34Differential e) expressing the motion vector field of the whole image using few parameters
stimation of the global motion parameters zoom and
pan &# 34, Signal Processing, vol. 16, no.3, pp. 24
9-265, Mar. 1989). This corresponds to the motion vector (ug (x, y), vg (x, x) of the pixel (x, y) in the image.
y)),

【0006】[0006]

【数1】 (Equation 1)

【0007】や、[0007]

【0008】[0008]

【数2】 (Equation 2)

【0009】の形式で表し、この動きベクトルを利用し
て動き補償を行う方式である。ここでa0〜a5,b0〜
b7 は動きパラメータである。動き補償を行う際には、
送信側と受信側で同じ予測画像が得られなければならな
い。このために、送信側は受信側へ a0〜a5 または
b0〜b7 の値を直接伝送しても良いが、代わりにいく
つかの代表点の動きベクトルを伝送する方法もある。い
ま、画像の左上端,右上端,左下端,右下端の画素の座
標がそれぞれ(0,0),(r,0),(0,s),(r,s)で
表されるとする(ただし、rとsは正の整数)。このと
き、代表点(0,0),(r,0),(0,s)の動きベクト
ルの水平,垂直成分をそれぞれ(ua,va),(ub,v
b),(uc,vc)とすると、(数1) は
In this method, motion compensation is performed using this motion vector. Where a0-a5, b0-
b7 is a motion parameter. When performing motion compensation,
The same predicted image must be obtained on the transmitting side and the receiving side. To this end, the sender sends a0 to a5 or
Although the values of b0 to b7 may be transmitted directly, there is also a method of transmitting the motion vectors of some representative points instead. Now, it is assumed that the coordinates of the upper left, upper right, lower left, and lower right pixels of the image are represented by (0, 0), (r, 0), (0, s), and (r, s), respectively. (However, r and s are positive integers). At this time, the horizontal and vertical components of the motion vector of the representative points (0, 0), (r, 0), (0, s) are (ua, va), (ub, v
b), (uc, vc), (Equation 1) becomes

【0010】[0010]

【数3】 (Equation 3)

【0011】と書き換えることができる。このことは
a0〜a5 を伝送する代わりに ua,va,ub,vb,u
c,vc を伝送しても同様の機能が実現できることを意
味する。この様子を図3に示す。現フレームの原画像3
02と参照画像301の間でグローバル動き補償が行わ
れたとして、動きパラメータの代わりに代表点303,
304,305の動きベクトル306,307,308
(このとき、動きベクトルは現フレームの原画像の点を
出発点として、参照画像内の対応する点を終点とするも
のとして定義する)を伝送しても良い。これと同じよう
に、4個の代表点(0,0),(r,0),(0,s),(r,
s)の動きベクトルの水平,垂直成分(ua,va),(u
b,vb),(uc,vc),(ud,vd)を用いて (数2) は、
It can be rewritten as This means
Instead of transmitting a0 to a5, ua, va, ub, vb, u
It means that the same function can be realized by transmitting c and vc. This is shown in FIG. Original image 3 of current frame
02 and the reference image 301, it is assumed that global motion compensation has been performed.
304, 305 motion vectors 306, 307, 308
(At this time, the motion vector is defined such that a point in the original image of the current frame is set as a starting point and a corresponding point in the reference image is set as an end point.) Similarly, four representative points (0,0), (r, 0), (0, s), (r,
horizontal and vertical components (ua, va), (u) of the motion vector
b, vb), (uc, vc), (ud, vd),

【0012】[0012]

【数4】 (Equation 4)

【0013】と書き換えることができる。したがって、
b0〜b7 を伝送する代わりに ua,va,ub,vb,u
c,vc,ud,vd を伝送しても同様の機能が実現でき
る。本明細書では (数1) を用いる方式を線形内・外挿
に基づくグローバル動き補償,(数2) を用いる方式を
共1次内・外挿に基づくグローバル動き補償とよぶこと
とする。
It can be rewritten as Therefore,
Instead of transmitting b0 to b7, ua, va, ub, vb, u
The same function can be realized by transmitting c, vc, ud, and vd. In this specification, the method using (Equation 1) is referred to as global motion compensation based on linear interpolation / extrapolation, and the method using (Equation 2) is referred to as global motion compensation based on linear interpolation / extrapolation.

【0014】代表点の動きベクトルを伝送する線形内・
外挿に基づくグローバル動き補償方式を採用した画像符
号化器の動き補償処理部401の構成例を図4に示す。
図1と同じ番号は同じものを指すとする。図1のブロッ
クマッチング部116をこの動き補償処理部401に入
れ替えることにより、グローバル動き補償を行う画像符
号化装置を構成することができる。グローバル動き補償
部402で前フレームの復号画像115と現フレームの
原画像101との間でグローバル動き補償に関する動き
推定が行われ、上記 ua,va,ub,vb,uc,vc の
値が推定される。これらの値に関する情報403は動き
情報120の一部として伝送される。グローバル動き補
償の予測画像404は数3を用いて合成され、ブロック
マッチング部405に供給される。ここでは、グローバ
ル動き補償の予測画像と現フレームの原画像との間でブ
ロックマッチングによる動き補償が行われ、ブロックの
動きベクトル情報406と最終的な予測画像117が得
られる。この動きベクトル情報は動きパラメータ情報と
多重化部407において多重化され、動き情報120と
して出力される。
[0014] In the linear form transmitting the motion vector of the representative point,
FIG. 4 shows a configuration example of the motion compensation processing unit 401 of the image encoder employing the global motion compensation method based on extrapolation.
The same numbers as those in FIG. 1 indicate the same things. By replacing the block matching unit 116 in FIG. 1 with the motion compensation processing unit 401, an image coding device that performs global motion compensation can be configured. The global motion compensation unit 402 performs motion estimation regarding global motion compensation between the decoded image 115 of the previous frame and the original image 101 of the current frame, and estimates the values of ua, va, ub, vb, uc, and vc. You. Information 403 on these values is transmitted as part of the motion information 120. The predicted image 404 of the global motion compensation is synthesized using Expression 3, and is supplied to the block matching unit 405. Here, motion compensation is performed by block matching between the predicted image of global motion compensation and the original image of the current frame, and the motion vector information 406 of the block and the final predicted image 117 are obtained. The motion vector information is multiplexed with the motion parameter information in the multiplexing unit 407, and is output as motion information 120.

【0015】図4とは異なる動き補償処理部501の構
成例を図5に示す。図1と同じ番号は同じものを指すと
する。図1のブロックマッチング部116をこの動き補
償処理部501に入れ替えることにより、グローバル動
き補償を行う画像符号化装置を構成することができる。
この例では、グローバル動き補償の予測画像にブロック
マッチングを適用するのではなく、各ブロックに関して
グローバル動き補償かブロックマッチングのいずれかが
適用される。前フレームの復号画像115と現フレーム
の原画像101との間で、グローバル動き補償部502
とブロックマッチング部505でそれぞれグローバル動
き補償とブロックマッチングが並列に行われる。選択ス
イッチ508は、グローバル動き補償による予測画像5
03とブロックマッチングによる予測画像506の間で
ブロックごとに最適な方式を選択する。代表点の動きベ
クトル504,ブロックごとの動きベクトル507,グ
ローバル動き補償/ブロックマッチングの選択情報50
9は多重化部510で多重化され、動き情報120とし
て出力される。
FIG. 5 shows a configuration example of the motion compensation processing unit 501 different from that of FIG. The same numbers as those in FIG. 1 indicate the same things. By replacing the block matching unit 116 in FIG. 1 with the motion compensation processing unit 501, an image coding device that performs global motion compensation can be configured.
In this example, instead of applying block matching to a predicted image of global motion compensation, either global motion compensation or block matching is applied to each block. Between the decoded image 115 of the previous frame and the original image 101 of the current frame, a global motion compensation unit 502
And the block matching unit 505 performs global motion compensation and block matching in parallel. The selection switch 508 controls the predicted image 5 based on global motion compensation.
03 and a prediction image 506 obtained by block matching, an optimum method is selected for each block. Motion vector 504 of representative point, motion vector 507 for each block, selection information 50 of global motion compensation / block matching
9 is multiplexed by a multiplexing unit 510 and output as motion information 120.

【0016】以上述べたグローバル動き補償を導入する
ことにより、画像の大局的な動きを少ないパラメータを
用いて表現することが可能となり、より高い情報圧縮率
が実現できる。しかし、その一方で符号化および復号化
における処理量は従来の方式と比較して増加する。特に
(数3) および (数4) に見られる除算は、処理を複雑
にする大きな要因となってしまう。
By introducing the above-described global motion compensation, the global motion of an image can be expressed using a small number of parameters, and a higher information compression rate can be realized. However, on the other hand, the processing amount in encoding and decoding is increased as compared with the conventional method. In particular
The division shown in (Equation 3) and (Equation 4) is a major factor complicating the processing.

【0017】[0017]

【発明が解決しようとする課題】画像全体の動きベクト
ル場を少ないパラメータによって近似するグローバル動
き補償では、予測画像の合成のための処理量が多くなる
問題が発生する。本発明の目的は、グローバル動き補償
における除算の処理を2進数のシフト演算に置き換える
ことにより、演算量を減少させることにある。
In global motion compensation in which the motion vector field of the entire image is approximated by a small number of parameters, there is a problem that the processing amount for synthesizing the predicted image increases. An object of the present invention is to reduce the amount of calculation by replacing the division process in global motion compensation with a binary shift operation.

【0018】[0018]

【課題を解決するための手段】グローバル動き補償を行
う際の代表点の座標をうまく選択することにより、除算
処理をシフト演算で実現できるようにする。
Means for Solving the Problems By appropriately selecting the coordinates of a representative point when performing global motion compensation, the division processing can be realized by a shift operation.

【0019】[0019]

【発明の実施の形態】以下の議論では、画素のサンプリ
ング間隔を水平,垂直方向共に1とし、画像の左上端,
右上端,左下端,右下端の画素の座標がそれぞれ(0,
0),(r,0),(0,s),(r,s)で表されるとする
(ただし、rとsは正の整数)。
DESCRIPTION OF THE PREFERRED EMBODIMENTS In the following discussion, the pixel sampling interval is set to 1 in both the horizontal and vertical directions,
The coordinates of the upper right, lower left, and lower right pixels are (0,
0), (r, 0), (0, s), (r, s) (where r and s are positive integers).

【0020】線形内・外挿(アフィン変換)または共1
次内・外挿(共1次変換)を用いた動き補償を行う際に
は、画素ごとの動きベクトルに対して量子化を行うと、
ミスマッチの防止や演算の簡略化などの効果を得ること
ができる(特願平06-193970)。以下では、画素の動き
ベクトルの水平成分と垂直成分が1/m(mは正の整
数)の整数倍であるとする。また、「従来の技術」で説
明した代表点の動きベクトルを用いるグローバル動き補
償を行うと仮定し、各代表点の動きベクトルは1/k
(kは正の整数)の整数倍であるとする。なお、本明細
書では、「画素の動きベクトル」はグローバル動き補償
を行う際に、実際に予測画像を合成するために用いる動
きベクトルのことを指す。一方、「代表点の動きベクト
ル」は画素の動きベクトルを計算するために用いるパラ
メータを意味している。したがって、量子化ステップサ
イズの違いなどが原因で、同じ座標上に存在していても
画素の動きベクトルと代表点の動きベクトルが一致しな
い場合も起こり得る。
Linear interpolation / extrapolation (affine transformation) or 1
When performing motion compensation using next-order interpolation / extrapolation (bilinear transformation), if quantization is performed on a motion vector for each pixel,
Effects such as prevention of mismatch and simplification of calculation can be obtained (Japanese Patent Application No. 06-193970). Hereinafter, it is assumed that the horizontal component and the vertical component of the motion vector of a pixel are integer multiples of 1 / m (m is a positive integer). Also, assuming that global motion compensation using the motion vector of the representative point described in “Prior Art” is performed, the motion vector of each representative point is 1 / k
(K is a positive integer). In this specification, “pixel motion vector” refers to a motion vector used to actually synthesize a predicted image when performing global motion compensation. On the other hand, “representative point motion vector” means a parameter used to calculate a pixel motion vector. Therefore, due to a difference in quantization step size or the like, a case may occur where the motion vector of a pixel does not match the motion vector of a representative point even if they are present on the same coordinates.

【0021】まず線形内・外挿を用いる場合について図
6を用いて説明する。このとき、「従来の技術」で述べ
たように代表点を画像601の隅に位置する点とはせ
ず、(i,j),(i+p,j),(i,j+q)に位置する点
602,603,604とする(i,j,p,qは整
数)。このとき、点602,603,604は画像の内
部に存在していても外部に存在していても良い。代表点
の動きベクトルの水平,垂直成分をk倍したものをそれ
ぞれ(u0,v0),(u1,v1),(u2,v2)とすると(u
0,v0,u1,v1,u2,v2 は整数)、画素(x,y)
の動きベクトルの水平,垂直成分をm倍したもの(u
(x,y),v(x,y))は以下の式で表すことができる
(ただし、x,y,u(x,y),v(x,y)は整数)。
First, the case of using linear interpolation / extrapolation will be described with reference to FIG. At this time, the representative point is not set as the point located at the corner of the image 601, as described in the “Prior Art”, and (i, j), (i + p, j), (i, j + q) (Points i, j, p, and q are integers). At this time, the points 602, 603, and 604 may exist inside or outside the image. If the horizontal and vertical components of the motion vector at the representative point are multiplied by k, they are (u0, v0), (u1, v1), and (u2, v2), respectively.
0, v0, u1, v1, u2, v2 are integers), pixel (x, y)
The horizontal and vertical components of the motion vector are multiplied by m (u
(x, y), v (x, y)) can be represented by the following equation (where x, y, u (x, y) and v (x, y) are integers).

【0022】[0022]

【数5】 (Equation 5)

【0023】ただし、「//」は通常の除算による演算結
果が整数ではない場合にこれを近隣の整数に丸め込む除
算で、演算子としての優先順位は乗除算と同等である。
演算誤差を小さくするためには、非整数値は最も近い整
数に丸め込まれることが望ましい。このとき整数に1/
2を加えた値の丸め込み方法は、(1) 0に近づける方向
に丸め込む、(2) 0から遠ざける方向に丸め込む、(3)
被除数が負の場合は0に近づける方向,正の場合は0か
ら遠ざける方向に丸め込む(除数は常に正であるとす
る)、(4) 被除数が負の場合は0から遠ざける方向,正
の場合は0に近づける方向に丸め込む(除数は常に正で
あるとする)、などが考えられる。これらの中で(3)
と(4)は、被除数の正負に関わらず丸め込みの方向が
変化しないため、正負判定が必要ない分だけ処理量の点
で有利である。(3)を用いた高速処理は以下の式によっ
て実現される。
However, “//” is a division in which the result of a normal division operation is not an integer and is rounded to a neighboring integer, and the priority as an operator is equivalent to multiplication / division.
In order to reduce the calculation error, it is desirable that the non-integer value be rounded to the nearest integer. In this case, 1 /
The rounding method of adding 2 is (1) rounding toward 0, (2) rounding away from 0, (3)
If the dividend is negative, it is rounded toward 0, if it is positive, it is rounded away from 0 (the divisor is always positive). (4) If the dividend is negative, it is rounded away from 0. If the dividend is positive, it is rounded. Rounding in a direction approaching 0 (the divisor is always positive) can be considered. Among these (3)
Since (4) and (4) do not change the rounding direction irrespective of the sign of the dividend, it is advantageous in terms of the processing amount by the amount that the sign judgment is not necessary. The high-speed processing using (3) is realized by the following equation.

【0024】[0024]

【数6】 (Equation 6)

【0025】ただし、「#」は小数点以下を0の方向に
切り捨てる整数の除算であり、一般に計算機では最も実
現しやすい形式の除算である。ここで、LとMは除算の
被除数を常に正に保つための数で、十分に大きな正の整
数である。また、(pqk#2)の項は、除算結果を最
も近い整数に丸め込むために用いられる。
However, "#" is an integer division which rounds off the decimal part in the direction of 0, and is a division which is most easily realized by a computer. Here, L and M are numbers for keeping the dividend dividend positive at all times, and are sufficiently large positive integers. The term (pqk # 2) is used to round the division result to the nearest integer.

【0026】処理を整数化することはそれ自体処理量の
低減に貢献するが、ここで p,q,k をそれぞれ2の
α,β,h0 乗(α,βは正の整数、h0は負ではない
整数)とすると、(数5) の除算は α+β+h0 ビット
のシフト演算で実現できるため、計算機や専用ハードウ
ェアにおける処理量を大きく減らすことができる。さら
にmを2のh1乗とすれば(h1は負ではない整数、h1
<α+β+h0)、(数6) は、
The conversion of the processing into integers contributes to the reduction of the processing amount. Here, p, q, k are each raised to the power of 2 to the power of α, β, h0 (α, β are positive integers, h0 is a negative integer). In this case, since the division of (Equation 5) can be realized by a shift operation of α + β + h0 bits, the amount of processing in a computer and dedicated hardware can be greatly reduced. Further, if m is 2 to the power of h1, (h1 is a non-negative integer, h1
<Α + β + h0), (Equation 6) is

【0027】[0027]

【数7】 (Equation 7)

【0028】と書き換えることができ(「x&#60&#60
α」はxをαビット左にシフトして下位αビットに0を
入れる、「x&#62&#62α」はxをαビット右にシフトし
て上位αビットに0または1を入れる(xが2の補数表
現の場合、xの最上位ビットが1のときは1,0のとき
は0を入れる)ことを意味し、これらの演算子の優先順
位は加減算と乗除算の中間であるとする)、さらに演算
を簡略化することができる。
(“X &# 60 &# 60
“α” shifts x to the left by α bits and puts 0 in the lower α bits, and “x &# 62 &# 62α” shifts x to the right by α bits and puts 0 or 1 in the upper α bits (x is 2 In the case of the complement expression of x, it is assumed that when the most significant bit of x is 1, it is 1, and when it is 0, 0 is inserted.) It is assumed that the precedence of these operators is intermediate between addition and subtraction and multiplication and division.) , The operation can be further simplified.

【0029】線形内・外挿を用いた場合、さらに(i+
p,j+q)に位置する代表点の動きベクトルの水平,垂
直成分をk倍したものを(u3,v3)として、(数5) は
(i,j),(i+p,j),(i+p,j+q)を代表点とすれ
ば、
When linear interpolation / extrapolation is used, (i +
Assuming that the horizontal and vertical components of the motion vector of the representative point located at (p, j + q) are multiplied by k, (u3, v3), (Equation 5) is (i, j), (i + p, j) , (I + p, j + q) as a representative point,

【0030】[0030]

【数8】 (Equation 8)

【0031】(i,j),(i,j+q),(i+p,j+q)を
代表点とすれば、
If (i, j), (i, j + q) and (i + p, j + q) are representative points,

【0032】[0032]

【数9】 (Equation 9)

【0033】(i+p,j),(i,j+q),(i+p,j+
q)を代表点とすれば、
(I + p, j), (i, j + q), (i + p, j +
If q) is the representative point,

【0034】[0034]

【数10】 (Equation 10)

【0035】と書き換えられ、p,q,k,m の値を
2の正の整数乗とすることによって同様に処理量を減ら
すことができる。
The processing amount can be similarly reduced by setting the values of p, q, k, and m to positive integer powers of 2.

【0036】共1次内・外挿を用いた場合には、代表点
(i,j),(i+p,j),(i,j+q),(i+p,j+q)そ
れぞれの動きベクトルの水平,垂直成分をk倍したもの
である(u0,v0),(u1,v1),(u2,v2),(u3,v3)
を用いて(u(x,y),v(x,y))は以下の式で表すこ
とができる。
When the first-order interpolation and extrapolation are used, the representative points (i, j), (i + p, j), (i, j + q), (i + p, j + q) The horizontal and vertical components of each motion vector are multiplied by k (u0, v0), (u1, v1), (u2, v2), (u3, v3).
(U (x, y), v (x, y)) can be expressed by the following equation.

【0037】[0037]

【数11】 [Equation 11]

【0038】この式も p,q,k,m の値をそれぞれ
2の α,β,h0,h1 乗とすることによって、
This equation also takes the values of p, q, k, and m as 2 to the power of α, β, h0, and h1, respectively.

【0039】[0039]

【数12】 (Equation 12)

【0040】と書き換えることができ、上と同様に処理
量を減らすことができる。
The processing amount can be reduced in the same manner as above.

【0041】送信側と受信側で同じグローバル動き補償
予測画像を得るためには、代表点の動きベクトルに関す
る情報を何らかの形で受信側に伝える必要がある。代表
点の動きベクトルそのまま伝送する方法もあるが、画像
の隅の点の動きベクトルを伝送し、この値から代表点の
動きベクトルを計算する方法もある。この方法に関し、
以下に説明する。
In order to obtain the same global motion compensation predicted image on the transmitting side and the receiving side, it is necessary to transmit information about the motion vector at the representative point to the receiving side in some form. Although there is a method of transmitting the motion vector of the representative point as it is, there is also a method of transmitting the motion vector of the point at the corner of the image and calculating the motion vector of the representative point from this value. For this method,
This will be described below.

【0042】まずは線形内・外挿が使用される場合につ
いて説明する。画像の隅の3個の点(0,0),(r,
0),(0,s)の動きベクトルが1/n整数倍の値のみ
とれるとして、これらの水平,垂直成分をn倍した(u
00,v00),(u01,v01),(u02,v02)が伝送されると
する。このとき、点(i,j),(i+p,j),(i,j+
q),(i+p,j+q)それぞれの動きベクトルの水平,
垂直成分をk倍したものである(u0,v0),(u1,v
1),(u2,v2),(u3,v3)を、
First, the case where linear interpolation / extrapolation is used will be described. The three points (0,0), (r,
Assuming that the motion vectors of (0) and (0, s) can only take values that are integral multiples of 1 / n, these horizontal and vertical components are multiplied by n (u
(00, v00), (u01, v01), (u02, v02) are transmitted. At this time, the points (i, j), (i + p, j), (i, j +
q), (i + p, j + q) horizontal of each motion vector,
The vertical component is multiplied by k (u0, v0), (u1, v
1), (u2, v2), (u3, v3)

【0043】[0043]

【数13】 (Equation 13)

【0044】と定義する。ただし、u'(x,y),v'
(x,y)は、(数5) を変形して、
Is defined as Where u '(x, y), v'
(x, y) is obtained by transforming (Equation 5),

【0045】[0045]

【数14】 [Equation 14]

【0046】と定義する。このとき、「///」は通常の除
算による演算結果が整数ではない場合にこれを近隣の整
数に丸め込む除算で、演算子としての優先順位は乗除算
と同等である。(u0,v0),(u1,v1),(u2,v2),(u
3,v3)の中から3点を選び、それらを代表点とするグ
ローバル動き補償を行えば、(0,0),(r,0),(0,
s)を代表点とするグローバル動き補償を近似すること
ができる。もちろんこのときにpとqを2の正の整数乗
とすれば、上で述べたように処理を簡略化することが可
能となる。なお、演算誤差を小さくするためには、「//
/」は非整数値を最も近い整数に丸め込むことが望まし
い。このとき整数に1/2を加えた値の丸め込み方法と
しては、上で述べた(1)〜(4)の方法が考えられる。
ただし、(数5)(画素ごとに計算)の場合と比較して、
(数14) (1枚の画像で3回のみ計算)は演算が実行
される回数が少ないため、(1)または(2)の方法を選
んだとしても全体の演算量に大きな影響は与えない。
Is defined as At this time, "///" is a division in which the result of the ordinary division is not an integer and the result is rounded to a neighboring integer, and the precedence as an operator is equivalent to the multiplication / division. (u0, v0), (u1, v1), (u2, v2), (u
3, v3), and performing global motion compensation using these as representative points, (0, 0), (r, 0), (0,
Global motion compensation with s) as a representative point can be approximated. Of course, at this time, if p and q are set to a positive integer power of 2, the processing can be simplified as described above. In order to reduce the calculation error, "//
It is desirable to round the non-integer value to the nearest integer. At this time, as a method of rounding a value obtained by adding 考 え to an integer, the methods (1) to (4) described above can be considered.
However, compared with the case of (Equation 5) (calculated for each pixel),
(Equation 14) (calculation only three times for one image) has a small number of times the operation is executed, so even if the method (1) or (2) is selected, it does not significantly affect the total amount of operation. .

【0047】画像の隅の点として (数13) の例とは異
なる3点が選ばれた場合も数8〜10を変形することに
よって同様の処理を実現することができる。上記の例に
加え、さらに画像の隅の点(r,s)の動きベクトルの
水平,垂直成分をn倍したものを(u03,v03)とすれ
ば、(数14) は、(u00,v00),(u01,v01),(u03,
v03)が伝送される場合には、
Even when three points different from the example of (Expression 13) are selected as the corner points of the image, the same processing can be realized by modifying Expressions 8 to 10. In addition to the above example, if the horizontal and vertical components of the motion vector at the corner point (r, s) of the image are multiplied by n, then (u03, v03), (Equation 14) becomes (u00, v00) ), (U01, v01), (u03,
v03) is transmitted,

【0048】[0048]

【数15】 (Equation 15)

【0049】と、(u00,v00),(u02,v02),(u03,v
03)が伝送される場合には、
And (u00, v00), (u02, v02), (u03, v
03) is transmitted,

【0050】[0050]

【数16】 (Equation 16)

【0051】と、(u01,v01),(u02,v02),(u03,v
03)が伝送される場合には、
And (u01, v01), (u02, v02), (u03, v
03) is transmitted,

【0052】[0052]

【数17】 [Equation 17]

【0053】と書き換えられる。Can be rewritten as

【0054】共1次内・外挿が行われる場合も同様であ
る。上と同様に画像の隅の4個の代表点(0,0),(r,
0),(0,s),(r,s)の動きベクトルが1/n整数倍
の値のみとれるとして、これらの水平,垂直成分をn倍
した(u00,v00),(u01,v01),(u02,v02),(u03,
v03)が伝送されるとする。このときの代表点(i,
j),(i+p,j),(i,j+q),(i+p,j+q)それぞ
れの動きベクトルの水平,垂直成分をk倍したものであ
る(u0,v0),(u1,v1),(u2,v2),(u3,v3)は上
と同様、(数13) で与えられる。ただし、u'(x,
y),v'(x,y)は、(数11) を変形して、
The same applies to the case where primary interpolation and extrapolation are performed. As above, the four representative points (0,0), (r,
Assuming that the motion vectors of (0), (0, s), and (r, s) can only take values that are integral multiples of 1 / n, these horizontal and vertical components are multiplied by n (u00, v00), (u01, v01). , (U02, v02), (u03,
v03) is transmitted. The representative point (i,
j), (i + p, j), (i, j + q), and (i + p, j + q) are obtained by multiplying the horizontal and vertical components by k times (u0, v0), (u1, v1), (u2, v2), (u3, v3) are given by (Equation 13) similarly to the above. Where u '(x,
y), v '(x, y) is obtained by transforming (Equation 11),

【0055】[0055]

【数18】 (Equation 18)

【0056】と定義される。Is defined as

【0057】画像の隅の点の動きベクトルを伝送し、こ
れに対して内・外挿を行うことによって代表点の動きベ
クトルを求める方式の長所は、画素ごとの動きベクトル
の範囲を限定しやすい点である。例えば、(数4) で与
えられる共1次内・外挿では、(x,y)が画像内の点で
あるとき、ug(x,y)の値は ua,ub,uc,udの最
大値を超えることも、最小値を下まわることもない。し
たがって、グローバル動き推定のときに ua,ub,u
c,ud の値がある制限範囲内(例えば±32画素以内
の範囲)に収まるような制約条件を加えれば、全ての画
素に関して ug(x,y)の値を同じ制限範囲内入れてお
くことができる(もちろんこれは vg(x,y)に関して
も成立する)。こうすると演算に必要な桁数を明確にす
ることができ、ソフトウェアまたはハードウェアを設計
する上で便利である。ただし、以上の議論は計算がすべ
て浮動小数点演算で行われた場合の議論なので、実際の
処理では注意が必要である。画像の隅の点の動きベクト
ルから代表点の動きベクトルを求める演算 (数18) に
は整数への丸め込みが存在するため、計算誤差の影響で
(数12) で求まる動きベクトルが上で述べた制限範囲
内から出る可能性を考慮する必要がある。特に代表点が
画像の内側に位置するような場合には注意が必要であ
る。これは、代表点が囲む長方形の外側に位置する画素
に関しては外挿によって動きベクトルが求められるた
め、丸め込み誤差が増幅される可能性があるためであ
る。外挿によって動きベクトルが求まる例を図7に示
す。画像701に対し、代表点702,703,70
4,705を用いてグローバル動き補償を行うと、画像
内の斜線で示された部分は外挿により動きベクトルが計
算されることになる。これは、斜線部が代表点が囲む長
方形706の外にあるためである。
The advantage of the method of transmitting the motion vector at the corner point of the image and performing the interpolation / extrapolation on the motion vector to obtain the motion vector at the representative point is that the range of the motion vector for each pixel is easily limited. Is a point. For example, in the bilinear interpolation / extrapolation given by (Equation 4), when (x, y) is a point in the image, the value of ug (x, y) is the maximum of ua, ub, uc, ud It does not exceed the value nor fall below the minimum value. Therefore, at the time of global motion estimation, ua, ub, u
If constraints are added so that the values of c and ud fall within a certain limit range (for example, within a range of ± 32 pixels), the values of ug (x, y) must be within the same limit range for all pixels. (Of course, this also holds for vg (x, y)). This makes it possible to clarify the number of digits required for the operation, which is convenient for designing software or hardware. However, the above discussion is for the case where all calculations are performed by floating point arithmetic, so care must be taken in actual processing. In the operation for obtaining the motion vector of the representative point from the motion vector of the point at the corner of the image (Equation 18), rounding to an integer exists.
It is necessary to consider the possibility that the motion vector obtained by (Equation 12) will be out of the above-mentioned limited range. In particular, care must be taken when the representative point is located inside the image. This is because, for pixels located outside the rectangle surrounded by the representative point, a motion vector is obtained by extrapolation, so that a rounding error may be amplified. FIG. 7 shows an example in which a motion vector is obtained by extrapolation. For the image 701, the representative points 702, 703, 70
When global motion compensation is performed using 4,705, a motion vector is calculated by extrapolation of a hatched portion in the image. This is because the hatched portion is outside the rectangle 706 surrounded by the representative point.

【0058】この問題への対策としては、4個の代表点
を、代表点が囲む長方形が画像全体を含むように配置す
る方法が有効である。この例を図8に示す。代表点80
2,803,804,805が囲む長方形806は、画
像801を含んでいる。こうすれば、全ての画素の動き
ベクトルは代表点からの内挿により求められるため、代
表点における丸め込み誤差の影響は画像内では増幅され
ない。したがって、代表点の丸め込み誤差より大きな誤
差が画像内で発生することはなく、誤差の上限が明確に
なる。ただし、代表点が囲む長方形を大きくし過ぎる
と、代表点の動きベクトルがとり得る値の範囲が広くな
るため、演算に必要な桁数が増加し、実装する上では不
利となる。
As a countermeasure against this problem, a method of arranging four representative points such that a rectangle surrounded by the representative points includes the entire image is effective. This example is shown in FIG. Representative point 80
A rectangle 806 surrounded by 2, 803, 804, and 805 includes an image 801. In this case, since the motion vectors of all the pixels are obtained by interpolation from the representative point, the influence of the rounding error at the representative point is not amplified in the image. Therefore, an error larger than the rounding error of the representative point does not occur in the image, and the upper limit of the error becomes clear. However, if the rectangle surrounded by the representative point is too large, the range of values that the motion vector of the representative point can take is widened, and the number of digits required for the calculation increases, which is disadvantageous in mounting.

【0059】以上の議論から、丸め込み誤差の影響を小
さくするため、pの値はr以上であり、qの値はs以上
であるが望ましい。p,q はそれぞれ r,s より小
さい場合にも、なるべく大きな値をとることが望まし
い。また、i,j の値は、画像内のできるだけ広い部
分が代表点により囲まれる領域に入るような値とするこ
とが望ましい。
From the above discussion, it is desirable that the value of p is not less than r and the value of q is not less than s in order to reduce the influence of the rounding error. Even when p and q are smaller than r and s, respectively, it is desirable to take values as large as possible. Further, it is desirable that the values of i, j are set so that a portion as large as possible in the image falls within a region surrounded by the representative points.

【0060】上で述べたように、グローバル動き補償に
共1次内・外挿を用いた場合には、4個の代表点が囲む
長方形に含まれる画素の動きベクトルの各成分は、代表
点の動きベクトルの各成分の最大値と最小値の間の値し
かとれないという性質がある。これに対し、線形内・外
挿が使用された場合には、3個の代表点が囲む3角形内
の画素の動きベクトルが同様の性質を持つ。したがっ
て、線形内・外挿を用いるグローバル動き補償を行う場
合には、画像の四隅の点の動きベクトルを伝送し、画像
の対角線によって分割される2個の直角3角形に対して
それぞれ独立にグローバル動き補償を行う方法が有効で
ある。こうすることにより、4隅の点に対する動きベク
トルの範囲に関する制約が、そのまま画像内のすべての
画素の動きベクトルに適用できる。このとき、i,j,
p,q の値は、2個の直角3角形の間で異なっていて
も良い。また、演算誤差の観点から言えば、外挿によっ
て画素の動きベクトルを計算するケースを避けるため、
代表点の囲む3角形が、グローバル動き補償の対象とな
る直角3角形を含むことが望ましい。この例を図9に示
す。画像901の4隅である点909,903,90
8,910の動きベクトルが伝送され、点909,90
3,910によって構成される直角3角形と、点90
9,910,908によって構成される直角3角形それ
ぞれに対し、独立にグローバル動き補償が行われる。し
たがって、頂点の動きベクトルの範囲に関して制約を設
ければ、画像内の全ての画素の動きベクトルもこの制約
範囲内に入ることになる。点909,903,910に
よって構成される直角3角形は代表点として点902,
903,904を使用し、点909,910,908に
よって構成される直角3角形は代表点として点906,
907,908を使用する。代表点によって構成される
3角形は、それぞれグローバル動き補償の対象となる直
角3角形を中に含んでいる。このため、代表点の動きベ
クトルの丸め込み誤差の影響は、画像内の点において増
幅されることはない。なお、この例では代表点が構成す
る2個の三角形は相似となっているが、必ずしもそうで
ある必要はない。
As described above, when the first-order interpolation and extrapolation are used for global motion compensation, each component of the motion vector of the pixel included in the rectangle surrounded by the four representative points is represented by the representative point. Has the property that only a value between the maximum value and the minimum value of each component of the motion vector can be taken. On the other hand, when linear interpolation / extrapolation is used, the motion vectors of the pixels in the triangle surrounded by the three representative points have similar properties. Therefore, when performing global motion compensation using linear interpolation / extrapolation, the motion vectors at the four corners of the image are transmitted, and global motion is independently performed on two right triangles divided by the diagonal of the image. A method of performing motion compensation is effective. In this way, the restrictions on the range of the motion vector for the four corner points can be applied to the motion vectors of all the pixels in the image as they are. At this time, i, j,
The value of p, q may be different between the two right triangles. Also, from the viewpoint of calculation error, in order to avoid the case of calculating a motion vector of a pixel by extrapolation,
It is preferable that the triangle surrounding the representative point includes a right triangle to be subjected to global motion compensation. This example is shown in FIG. Points 909, 903, and 90 at the four corners of image 901
8,910 motion vectors are transmitted and points 909,90
3,910, a right-angled triangle and a point 90
Global motion compensation is independently performed on each of the right triangles constituted by 9, 910 and 908. Therefore, if a restriction is imposed on the range of the motion vector of the vertex, the motion vectors of all the pixels in the image will fall within this restriction range. A right triangle formed by points 909, 903, and 910 is represented as a point 902 as a representative point.
Using 903 and 904, a right triangle composed of points 909, 910 and 908 is represented by points 906 and 906 as representative points.
907 and 908 are used. The triangles formed by the representative points each include a right-angled triangle to be subjected to global motion compensation. Therefore, the influence of the rounding error of the motion vector at the representative point is not amplified at a point in the image. In this example, the two triangles formed by the representative points are similar but need not be so.

【0061】[0061]

【発明の効果】本発明により、グローバル動き補償の予
測画像合成処理における除算の処理をシフト演算で代用
することが可能となり、ソフトウェアや専用ハードウェ
アによる処理を簡略化することが可能となる。
According to the present invention, it becomes possible to substitute the division processing in the prediction image synthesis processing of the global motion compensation with the shift operation, and to simplify the processing by software and dedicated hardware.

【図面の簡単な説明】[Brief description of the drawings]

【図1】H.261の画像符号化器の構成例を示した図
である。
FIG. 1 is a diagram illustrating a configuration example of an H.261 image encoder.

【図2】H.261の画像復号化器の構成例を示した図
である。
FIG. 2 is a diagram illustrating a configuration example of an H.261 image decoder.

【図3】代表点の動きベクトルを伝送するグローバル動
き補償の例を示した図である。
FIG. 3 is a diagram illustrating an example of global motion compensation for transmitting a motion vector at a representative point.

【図4】グローバル動き補償の予測画像に対してブロッ
クマッチングを行う画像符号化器の動き補償処理部を示
した図である。
FIG. 4 is a diagram illustrating a motion compensation processing unit of an image encoder that performs block matching on a predicted image of global motion compensation.

【図5】ブロックごとにグローバル動き補償とブロック
マッチングを選択する画像符号化器の動き補償処理部を
示した図である。
FIG. 5 is a diagram illustrating a motion compensation processing unit of an image encoder that selects global motion compensation and block matching for each block.

【図6】高速な処理を行うための代表点の配置の例を示
した図である。
FIG. 6 is a diagram showing an example of arrangement of representative points for performing high-speed processing.

【図7】画像内において、外挿によって動きベクトルを
求める領域を示した図である。
FIG. 7 is a diagram showing an area in an image for which a motion vector is obtained by extrapolation.

【図8】画像内のすべての画素の動きベクトルが、代表
点の動きベクトルからの内挿によって求められる場合を
示した図である。
FIG. 8 is a diagram showing a case where motion vectors of all pixels in an image are obtained by interpolation from motion vectors of representative points.

【図9】画像を2個の直角3角形に分割し、それぞれに
対して代表点の動きベクトルからの内挿によるグローバ
ル動き補償を適用した例を示した図である。
FIG. 9 is a diagram illustrating an example in which an image is divided into two right-angled triangles, and global motion compensation by interpolation from a motion vector at a representative point is applied to each of the triangles.

【符号の説明】[Explanation of symbols]

100…画像符号化器、101…入力画像、102…減
算器、103…誤差画像、104…DCT変換器、10
5…DCT係数量子化器、106,201…量子化DC
T係数、108,204…DCT係数逆量子化器、10
9,205…逆DCT変換器、110,206…復号誤
差画像、111,207…加算器、112…現フレーム
の復号画像、113,215…フレーム間/フレーム内
符号化切り換えスイッチの出力画像、114,209…
フレームメモリ、115,210…前フレームの復号画
像、116,405,505…ブロックマッチング部、
117,212…現フレームの予測画像、118,21
3…「0」信号、119,214…フレーム間/フレー
ム内符号化切り換えスイッチ、120,202,40
6,507…動きベクトル情報、121,203…フレ
ーム間/フレーム内識別フラグ、122,407,51
0…多重化器、123…伝送ビットストリーム、200
…画像復号化器、208…出力画像、211…予測画像
合成部、216…分離器、301…参照画像、302…
現フレームの原画像、303,304,305,60
2,603,604,702,703,704,70
5,802,803,804,805,902,90
4,906,907…代表点、306,307,308
…代表点の動きベクトル、401,501…グローバル
動き補償を行う動き補償処理部、402,502…グロ
ーバル動き補償部、403,504…動きパラメータ、
404,503…グローバル動き補償の予測画像、50
6…ブロックマッチングによる予測画像、508…ブロ
ックマッチング/グローバル動き補償切り換えスイッ
チ、509…ブロックマッチング/グローバル動き補償
の選択情報、601,701,801,901…グロー
バル動き補償の対象となる画像、706…代表点が囲む
長方形、903,908…画像の隅の点と代表点を兼用
する点、909,910…画像の隅の点。
100: image encoder, 101: input image, 102: subtractor, 103: error image, 104: DCT converter, 10
5. DCT coefficient quantizer, 106, 201 ... Quantized DC
T coefficient, 108, 204... DCT coefficient inverse quantizer, 10
9, 205: inverse DCT converter, 110, 206: decoded error image, 111, 207: adder, 112: decoded image of current frame, 113, 215: output image of inter-frame / intra-frame encoding changeover switch, 114 , 209 ...
Frame memory, 115, 210... Decoded image of previous frame, 116, 405, 505.
117, 212... Predicted image of the current frame, 118, 21
3: "0" signal, 119, 214 ... Inter-frame / intra-frame coding changeover switch, 120, 202, 40
6,507... Motion vector information, 121, 203... Inter-frame / intra-frame identification flags, 122, 407, 51
0: multiplexer, 123: transmission bit stream, 200
... image decoder, 208 ... output image, 211 ... predicted image synthesizing unit, 216 ... separator, 301 ... reference image, 302 ...
Original image of the current frame, 303, 304, 305, 60
2,603,604,702,703,704,70
5,802,803,804,805,902,90
4,906,907 ... representative points, 306,307,308
... motion vectors of representative points, 401, 501 ... motion compensation processing units for performing global motion compensation, 402, 502 ... global motion compensation units, 403, 504 ... motion parameters,
404, 503: predicted image of global motion compensation, 50
6: predicted image by block matching, 508: block matching / global motion compensation changeover switch, 509: selection information of block matching / global motion compensation, 601, 701, 801, 901: image to be subjected to global motion compensation, 706 ... Rectangles surrounded by the representative points, 903, 908... Points that also serve as representative points and image corner points, and 909, 910... Image corner points.

Claims (28)

【特許請求の範囲】[Claims] 【請求項1】3個の代表点における動きベクトルに対
し、線形内・外挿を行うことによって画像内のすべての
画素に対する動きベクトルを計算する際に、 画素のサンプリング間隔が水平,垂直方向共に1であ
り、かつ、サンプリング点は座標の水平,垂直成分が共
に整数である点の上に存在しているとして、 座標が (i,j),(i+p,j),(i,j+q),(i+p,j+
q) で表される4点の中から3点を代表点として使用し
(iとjは整数、pとqは正の整数)、 かつ、pとqがそれぞれ2のα乗と2のβ乗 (α,βは
正の整数) であることを特徴とするグローバル動き補償
方法。
When calculating motion vectors for all pixels in an image by performing linear interpolation and extrapolation on the motion vectors at three representative points, the sampling intervals of the pixels are set in both the horizontal and vertical directions. 1 and the sampling point is on a point where both the horizontal and vertical components of the coordinates are integers, and the coordinates are (i, j), (i + p, j), (i, j + q), (i + p, j +
q) 3 points out of 4 points are used as representative points.
(i and j are integers, p and q are positive integers), and p and q are 2 α powers and 2 β powers (α and β are positive integers), respectively. Compensation method.
【請求項2】4個の代表点における動きベクトルに対
し、共1次内・外挿を行うことによって画像内のすべて
の画素に対する動きベクトルを計算する際に、 画素のサンプリング間隔が水平,垂直方向共に1であ
り、かつ、サンプリング点は座標の水平,垂直成分が共
に整数である点の上に存在しているとして、 上記4個の代表点の座標が (i,j),(i+p,j),(i,
j+q),(i+p,j+q) で表すことができ (iとjは整
数)、 かつ、pとqがそれぞれ2のα乗と2のβ乗 (α,βは
正の整数) であることを特徴とするグローバル動き補償
方法。
2. When calculating motion vectors for all pixels in an image by performing linear interpolation and extrapolation on motion vectors at four representative points, the sampling intervals of pixels are set to be horizontal and vertical. The coordinates of the four representative points are (i, j) and (i +), assuming that the sampling point exists on a point where both directions are 1 and the horizontal and vertical components of the coordinates are both integers. p, j), (i,
j + q) and (i + p, j + q) (i and j are integers), and p and q are 2 to the power of α and 2 to the power of β (α and β are positive A global motion compensation method characterized by being an integer.
【請求項3】画像における画素のサンプリング間隔が水
平,垂直方向共に1であり、 かつ、サンプリング点は座標の水平,垂直成分が共に整
数である点の上に存在しており、 かつ、代表点 (i,j),(i+p,j),(i,j+q) にお
ける動きベクトルの水平・垂直成分が1/kの整数倍の
値のみをとり(ただし、i,jは整数、p,qは正の整
数、kは2のh0乗、かつ、h0は負ではない整数)、 かつ、各画素の動きベクトルの水平・垂直成分が1/m
の整数倍の値のみをとる(ただし、mは2のh1乗、か
つ、h1は負ではない整数)ときに、 代表点 (i,j),(i+p,j),(i,j+q) における動
きベクトルの水平・垂直成分をk倍したものである (u
0,v0),(u1,v1),(u2,v2) を用いて (ただし、u
0,v0,u1,v1,u2,v2は整数)、 画素 (x,y) の動きベクトルの水平・垂直成分をm倍
したものである (u(x,y),v(x,y)) が (ただし、
x,y,u(x,y),v(x,y) は整数)、 u(x,y)=(u0・p・q+(u1-u0)(x-i)・q+(u2-u0)(y-j)・p)m//(p
・q・k) v(x,y)=(v0・p・q+(v1-v0)(x-i)・q+(v2-v0)(y-j)・p)m//(p
・q・k) によって表され (ただし、「//」は通常の除算による演
算結果が整数ではない場合にこれを近隣の整数に丸め込
む除算で、演算子としての優先順位は乗除算と同等)、 かつ、pとqがそれぞれ2のα乗と2のβ乗 (α,βは
正の整数) であることを特徴とするフレーム間予測画像
の合成方法。
3. The sampling interval of pixels in an image is 1 in both the horizontal and vertical directions, and the sampling point exists on a point where both the horizontal and vertical components of the coordinates are integers. The horizontal and vertical components of the motion vector in (i, j), (i + p, j), and (i, j + q) take only integer multiples of 1 / k (where i and j are integers, p and q are positive integers, k is 2 to the power of h0, and h0 is a non-negative integer), and the horizontal and vertical components of the motion vector of each pixel are 1 / m
(I.e., m is 2 to the power of h1 and h1 is a non-negative integer), the representative points (i, j), (i + p, j), (i, j) + u) is obtained by multiplying the horizontal and vertical components of the motion vector by k times (u
0, v0), (u1, v1), (u2, v2) (where u
0, v0, u1, v1, u2, and v2 are integers), and the horizontal and vertical components of the motion vector of the pixel (x, y) are multiplied by m (u (x, y), v (x, y) ) Is (but
x, y, u (x, y) and v (x, y) are integers), u (x, y) = (u0 · p · q + (u1-u0) (xi) · q + (u2-u0) ( yj) ・ p) m // (p
・ Q ・ k) v (x, y) = (v0 ・ p ・ q + (v1-v0) (xi) ・ q + (v2-v0) (yj) ・ p) m // (p
(Q / k) (however, "//" is a division in which the result of a normal division operation is not an integer and is rounded to a nearby integer, and the precedence as an operator is equivalent to multiplication / division) , And p and q are 2 to the power of α and 2 to the power of β (α and β are positive integers), respectively.
【請求項4】画像における画素のサンプリング間隔が水
平,垂直方向共に1であり、 かつ、サンプリング点は座標の水平,垂直成分が共に整
数である点の上に存在しており、 かつ、代表点 (i,j),(i+p,j),(i+p,j+q) に
おける動きベクトルの水平・垂直成分が1/kの整数倍
の値のみをとり (ただし、i,jは整数、p,qは正の
整数、kは2のh0乗、かつ、h0は負ではない整数)、 かつ、各画素の動きベクトルの水平・垂直成分が1/m
の整数倍の値のみをとる (ただし、mは2のh1乗、か
つ、h1は負ではない整数) ときに、 代表点 (i,j),(i+p,j),(i+p,j+q) における
動きベクトルの水平・垂直成分をk倍したものである
(u0,v0),(u1,v1),(u3,v3) を用いて (ただし、
u0,v0,u1,v1,u3,v3は整数)、 画素 (x,y) の動きベクトルの水平・垂直成分をm倍し
たものである (u(x,y),v(x,y)) が (ただし、
x,y,u(x,y),v(x,y) は整数)、 u(x,y)=(u0・p・q+(u1-u0)(x-i)・q+(u3-u1)(y-j)・p)m//(p
・q・k) v(x,y)=(v0・p・q+(v1-v0)(x-i)・q+(v3-v1)(y-j)・p)m//(p
・q・k) によって表され (ただし、「//」は通常の除算による演
算結果が整数ではない場合にこれを近隣の整数に丸め込
む除算で、演算子としての優先順位は乗除算と同等)、 かつ、pとqがそれぞれ2のα乗と2のβ乗 (α,βは
正の整数) であることを特徴とするフレーム間予測画像
の合成方法。
4. A pixel sampling interval in an image is 1 in both the horizontal and vertical directions, and the sampling point is located on a point where both the horizontal and vertical components of the coordinate are integers. The horizontal and vertical components of the motion vector at (i, j), (i + p, j) and (i + p, j + q) take only integer multiples of 1 / k (where i and j are Integers, p and q are positive integers, k is 2 to the power of h0, and h0 is a non-negative integer), and the horizontal and vertical components of the motion vector of each pixel are 1 / m
(Where m is 2 to the power of h1 and h1 is a non-negative integer), the representative points (i, j), (i + p, j), (i + p , j + q) obtained by multiplying the horizontal and vertical components of the motion vector by k times.
Using (u0, v0), (u1, v1), (u3, v3) (where
u0, v0, u1, v1, u3, v3 are integers) and the horizontal and vertical components of the motion vector of pixel (x, y) are multiplied by m (u (x, y), v (x, y) ) Is (but
x, y, u (x, y) and v (x, y) are integers), u (x, y) = (u0 · p · q + (u1-u0) (xi) · q + (u3-u1) ( yj) ・ p) m // (p
・ Q ・ k) v (x, y) = (v0 ・ p ・ q + (v1-v0) (xi) ・ q + (v3-v1) (yj) ・ p) m // (p
(Q / k) (however, "//" is a division in which the result of a normal division operation is not an integer and is rounded to a nearby integer, and the precedence as an operator is equivalent to multiplication / division) , And p and q are 2 to the power of α and 2 to the power of β (α and β are positive integers), respectively.
【請求項5】画像における画素のサンプリング間隔が水
平,垂直方向共に1であり、 かつ、サンプリング点は座標の水平,垂直成分が共に整
数である点の上に存在しており、 かつ、代表点 (i,j),(i,j+q),(i+p,j+q) に
おける動きベクトルの水平・垂直成分が1/kの整数倍
の値のみをとり (ただし、i,jは整数、p,qは正の
整数、kは2のh0乗、かつ、h0は負ではない整数)、 かつ、各画素の動きベクトルの水平・垂直成分が1/m
の整数倍の値のみをとる (ただし、mは2のh1乗、か
つ、h1は負ではない整数) ときに、 代表点 (i,j),(i,j+q),(i+p,j+q) における
動きベクトルの水平・垂直成分をk倍したものである
(u0,v0),(u2,v2),(u3,v3) を用いて (ただし、
u0,v0,u2,v2,u3,v3は整数)、 画素 (x,y) の動きベクトルの水平・垂直成分をm倍し
たものである (u(x,y),v(x,y)) が (ただし、
x,y,u(x,y),v(x,y) は整数)、 u(x,y)=(u0・p・q+(u3-u2)(x-i)・q+(u2-u0)(y-j)・p)m//(p
・q・k) v(x,y)=(v0・p・q+(v3-v2)(x-i)・q+(v2-v0)(y-j)・p)m//(p
・q・k) によって表され (ただし、「//」は通常の除算による演
算結果が整数ではない場合にこれを近隣の整数に丸め込
む除算で、演算子としての優先順位は乗除算と同等)、 かつ、pとqがそれぞれ2のα乗と2のβ乗 (α,βは
正の整数) であることを特徴とするフレーム間予測画像
の合成方法。
5. A sampling interval of pixels in an image is 1 in both horizontal and vertical directions, and a sampling point is located on a point where both horizontal and vertical components of coordinates are integers. The horizontal and vertical components of the motion vector at (i, j), (i, j + q), (i + p, j + q) take only integer multiples of 1 / k (where i, j are Integers, p and q are positive integers, k is 2 to the power of h0, and h0 is a non-negative integer), and the horizontal and vertical components of the motion vector of each pixel are 1 / m
(I.e., m is 2 to the power of h1 and h1 is a non-negative integer), the representative points (i, j), (i, j + q), (i + p , j + q) obtained by multiplying the horizontal and vertical components of the motion vector by k times.
Using (u0, v0), (u2, v2), (u3, v3) (where
u0, v0, u2, v2, u3, v3 are integers) and the horizontal and vertical components of the motion vector of pixel (x, y) are multiplied by m (u (x, y), v (x, y) ) Is (but
x, y, u (x, y) and v (x, y) are integers), u (x, y) = (u0 · p · q + (u3-u2) (xi) · q + (u2-u0) ( yj) ・ p) m // (p
・ Q ・ k) v (x, y) = (v0 ・ p ・ q + (v3-v2) (xi) ・ q + (v2-v0) (yj) ・ p) m // (p
(Q / k) (however, "//" is a division in which the result of a normal division operation is not an integer and is rounded to a nearby integer, and the precedence as an operator is equivalent to multiplication / division) , And p and q are 2 to the power of α and 2 to the power of β (α and β are positive integers), respectively.
【請求項6】画像における画素のサンプリング間隔が水
平,垂直方向共に1であり、 かつ、サンプリング点は座標の水平,垂直成分が共に整
数である点の上に存在しており、 かつ、代表点 (i+p,j),(i,j+q),(i+p,j+q)
における動きベクトルの水平・垂直成分が1/kの整数
倍の値のみをとり (ただし、i,jは整数、p,qは正
の整数、kは2のh0乗、かつ、h0は負ではない整
数)、 かつ、各画素の動きベクトルの水平・垂直成分が1/m
の整数倍の値のみをとる (ただし、mは2のh1乗、か
つ、h1は負ではない整数) ときに、 代表点 (i+p,j),(i,j+q),(i+p,j+q) にお
ける動きベクトルの水平・垂直成分をk倍したものであ
る (u1,v1),(u2,v2),(u3,v3) を用いて(ただ
し、u1,v1,u2,v2,u3,v3は整数)、 画素 (x,y) の動きベクトルの水平・垂直成分をm倍し
たものである (u(x,y),v(x,y)) が (ただし、
x,y,u(x,y),v(x,y) は整数)、 u(x,y)=(u3・p・q+(u2-u3)(p-x+i)・q+(u1-u3)(q-y+j)・p)m
//(p・q・k) v(x,y)=(v0・p・q+(v2-v3)(p-x+i)・q+(v1-v3)(q-y+j)・p)m
//(p・q・k) によって表され (ただし、「//」は通常の除算による演
算結果が整数ではない場合にこれを近隣の整数に丸め込
む除算で、演算子としての優先順位は乗除算と同等)、 かつ、pとqがそれぞれ2のα乗と2のβ乗 (α,βは
正の整数) であることを特徴とするフレーム間予測画像
の合成方法。
6. A pixel sampling interval in an image is 1 in both the horizontal and vertical directions, and a sampling point exists on a point where both horizontal and vertical components of coordinates are integers. (i + p, j), (i, j + q), (i + p, j + q)
Take only values of the horizontal and vertical components of the motion vector at integer multiples of 1 / k (where i and j are integers, p and q are positive integers, k is 2 to the power of h0, and h0 is a negative And the horizontal and vertical components of the motion vector of each pixel are 1 / m
(Where m is 2 to the power of h1 and h1 is a non-negative integer), the representative points (i + p, j), (i, j + q), (i + p, j + q) by multiplying the horizontal and vertical components of the motion vector by k times (u1, v1), (u2, v2), (u3, v3) (where u1, v1, u2 , V2, u3, and v3 are integers), and the horizontal and vertical components of the motion vector of the pixel (x, y) are multiplied by m (u (x, y), v (x, y)) are
x, y, u (x, y) and v (x, y) are integers), u (x, y) = (u3pq + (u2-u3) (p-x + i) q + (u1 -u3) (q-y + j) ・ p) m
// (p ・ q ・ k) v (x, y) = (v0 ・ p ・ q + (v2-v3) (p-x + i) ・ q + (v1-v3) (q-y + j) ・ p ) m
// Represented by (p, q, k) (however, "//" is a division in which the result of a normal division operation is not an integer and is rounded to a nearby integer. And p and q are 2 to the power of α and 2 to the power of β (α and β are positive integers), respectively.
【請求項7】画像における画素のサンプリング間隔が水
平,垂直方向共に1であり、 かつ、サンプリング点は座標の水平,垂直成分が共に整
数である点の上に存在しており、 かつ、代表点 (i,j),(i+p,j),(i,j+q),(i+
p,j+q) における動きベクトルの水平・垂直成分が1
/kの整数倍の値のみをとり (ただし、i,jは整数、
p,qは正の整数、kは2のh0乗、かつ、h0は負では
ない整数)、 かつ、各画素の動きベクトルの水平・垂直成分が1/m
の整数倍の値のみをとる (ただし、mは2のh1乗、か
つ、h1は負ではない整数) ときに、 代表点 (i,j),(i+p,j),(i,j+q),(i+p,j+
q) における動きベクトルの水平・垂直成分をk倍した
ものである (u0,v0),(u1,v1),(u2,v2),(u3,
v3) を用いて (ただし、u0,v0,u1,v1,u2,v
2,u3,v3は整数)、 画素 (x,y) の動きベクトルの水平・垂直成分をm倍し
たものである (u(x,y),v(x,y)) が (ただし、
x,y,u(x,y),v(x,y) は整数)、 u(x,y)=((u0(p-x+i)+u1(x-i))(q-y+j)+(u2(p-x+i)+u3(x
-i))(y-j))m//(p・q・k) v(x,y)=((v0(p-x+i)+v1(x-i))(q-y+j)+(v2(p-x+i)+v3(x
-i))(y-j))m//(p・q・k) によって表され (ただし、「//」は通常の除算による演
算結果が整数ではない場合にこれを近隣の整数に丸め込
む除算で、演算子としての優先順位は乗除算と同等)、 かつ、pとqがそれぞれ2のα乗と2のβ乗 (α,βは
正の整数) であることを特徴とするフレーム間予測画像
の合成方法。
7. A sampling interval of pixels in an image is 1 in both horizontal and vertical directions, and a sampling point is located on a point where both horizontal and vertical components of coordinates are integers. (i, j), (i + p, j), (i, j + q), (i +
p, j + q), the horizontal and vertical components of the motion vector are 1
/ K only (where i and j are integers,
p and q are positive integers, k is 2 to the power of h0, and h0 is a non-negative integer), and the horizontal and vertical components of the motion vector of each pixel are 1 / m
(Where m is 2 to the power of h1 and h1 is a non-negative integer), the representative points (i, j), (i + p, j), (i, j) + q), (i + p, j +
q) are obtained by multiplying the horizontal and vertical components of the motion vector by k times (u0, v0), (u1, v1), (u2, v2), (u3,
v3) (where u0, v0, u1, v1, u2, v
2, u3 and v3 are integers), and the horizontal and vertical components of the motion vector of pixel (x, y) are multiplied by m (u (x, y), v (x, y))
x, y, u (x, y) and v (x, y) are integers), u (x, y) = ((u0 (p-x + i) + u1 (xi)) (q-y + j ) + (u2 (p-x + i) + u3 (x
-i)) (yj)) m // (pqqk) v (x, y) = ((v0 (p-x + i) + v1 (xi)) (q-y + j) + ( v2 (p-x + i) + v3 (x
-i)) (yj)) m // (p ・ q ・ k) (where `` // '' is a division that rounds the result of a normal division operation to a neighboring integer when the result is not an integer. , And the priority as an operator is equivalent to multiplication / division), and p and q are 2 to the power of α and 2 to the power of β (α and β are positive integers). Synthesis method.
【請求項8】画像の左上端,右上端,左下端,右下端の
画素の座標がそれぞれ (0,0),(r,0),(0,s),
(r,s) で表されるとき (ただし、rとsは正の整
数)、 画像の隅の点 (0,0),(r,0),(0,s) における動
きベクトルの水平・垂直成分が1/nの整数倍の値をと
るとして (ただし、nは正の整数)、これらの動きベク
トルをn倍したものである (u00,v00),(u01,v0
1),(u02,v02) (ただし、u00,v00,u01,v01,
u02,v02は整数) を用いて、 点 (i,j),(i+p,j),(i,j+q),(i+p,j+q)
の動きベクトルの水平・垂直成分をk倍したものである
(u0,v0),(u1,v1),(u2,v2),(u3,v3)が、 u'(x,y)=(u00・r・s+(u01-u00)x・s+(u02-u00)y・r)k///(r・
s・n) v'(x,y)=(v00・r・s+(v01-v00)x・s+(v02-v00)y・r)k///(r・
s・n) u0=u'(i,j) v0=v'(i,j) u1=u'(i+p,j) v1=v'(i+p,j) u2=u'(i,j+q) v2=v'(i,j+q) u3=u'(i+p,j+q) v3=v'(i+p,j+q) によって与えられ (ただし、「///」は通常の除算による
演算結果が整数ではない場合にこれを近隣の整数に丸め
込む除算で、演算子としての優先順位は乗除算と同
等)、これら4点の中の3点とその動きベクトルを、代
表点とその動きベクトルして使用することを特徴とする
請求項1または3ないし6に記載のフレーム間予測画像
の合成方法。
8. The coordinates of pixels at the upper left, upper right, lower left and lower right corners of the image are (0,0), (r, 0), (0, s),
When represented by (r, s) (where r and s are positive integers), the horizontal and vertical positions of the motion vector at the points (0,0), (r, 0), (0, s) at the corners of the image Assuming that the vertical component takes an integer multiple of 1 / n (where n is a positive integer), these motion vectors are multiplied by n (u00, v00), (u01, v0).
1), (u02, v02) (where u00, v00, u01, v01,
u02, v02 are integers), and the points (i, j), (i + p, j), (i, j + q), (i + p, j + q)
The horizontal and vertical components of the motion vector are multiplied by k.
(u0, v0), (u1, v1), (u2, v2), (u3, v3) are given by u '(x, y) = (u00 · rs · s + (u01-u00) x · s + (u02− u00) y ・ r) k /// (r ・
s ・ n) v '(x, y) = (v00 ・ r ・ s + (v01-v00) x ・ s + (v02-v00) y ・ r) k /// (r ・
s ・ n) u0 = u '(i, j) v0 = v' (i, j) u1 = u '(i + p, j) v1 = v' (i + p, j) u2 = u '(i , j + q) v2 = v '(i, j + q) u3 = u' (i + p, j + q) v3 = v '(i + p, j + q) (where "/ //) is a division that rounds the result of a normal division operation to an adjacent integer when the result is not an integer. The precedence as an operator is equivalent to multiplication / division.) 3 points out of these 4 points and their movement 7. The method according to claim 1, wherein the vector is used as a representative point and a motion vector thereof.
【請求項9】画像の左上端,右上端,左下端,右下端の
画素の座標がそれぞれ (0,0),(r,0),(0,s),
(r,s) で表されるとき (ただし、rとsは正の整
数)、 画像の隅の点 (0,0),(r,0),(r,s) における動
きベクトルの水平・垂直成分が1/nの整数倍の値をと
るとして (ただし、nは正の整数)、これらの動きベク
トルをn倍したものである (u00,v00),(u01,v0
1),(u03,v03) (ただし、u00,v00,u01,v01,
u03,v03は整数) を用いて、 点 (i,j),(i+p,j),(i,j+q),(i+p,j+q)
の動きベクトルの水平・垂直成分をk倍したものである
(u0,v0),(u1,v1),(u2,v2),(u3,v3)が、 u'(x,y)=(u00・r・s+(u01-u00)x・s+(u03-u01)y・r)k///(r・
s・n) v'(x,y)=(v00・r・s+(v01-v00)x・s+(v03-v01)y・r)k///(r・
s・n) u0=u'(i,j) v0=v'(i,j) u1=u'(i+p,j) v1=v'(i+p,j) u2=u'(i,j+q) v2=v'(i,j+q) u3=u'(i+p,j+q) v3=v'(i+p,j+q) によって与えられ (ただし、「///」は通常の除算による
演算結果が整数ではない場合にこれを近隣の整数に丸め
込む除算で、演算子としての優先順位は乗除算と同
等)、これら4点の中の3点とその動きベクトルを、代
表点とその動きベクトルして使用することを特徴とする
請求項1または3ないし6に記載のフレーム間予測画像
の合成方法。
9. The coordinates of the upper left, upper right, lower left and lower right pixels of the image are (0,0), (r, 0), (0, s),
When represented by (r, s) (where r and s are positive integers), the horizontal and vertical positions of the motion vector at the points (0,0), (r, 0), (r, s) at the corners of the image Assuming that the vertical component takes an integer multiple of 1 / n (where n is a positive integer), these motion vectors are multiplied by n (u00, v00), (u01, v0).
1), (u03, v03) (where u00, v00, u01, v01,
u03, v03 are integers), and the points (i, j), (i + p, j), (i, j + q), (i + p, j + q)
The horizontal and vertical components of the motion vector are multiplied by k.
(u0, v0), (u1, v1), (u2, v2), (u3, v3) are given by u ′ (x, y) = (u00 · rs · s + (u01−u00) x · s + (u03− u01) y ・ r) k /// (r ・
s ・ n) v '(x, y) = (v00 ・ r ・ s + (v01-v00) x ・ s + (v03-v01) y ・ r) k /// (r ・
s ・ n) u0 = u '(i, j) v0 = v' (i, j) u1 = u '(i + p, j) v1 = v' (i + p, j) u2 = u '(i , j + q) v2 = v '(i, j + q) u3 = u' (i + p, j + q) v3 = v '(i + p, j + q) (where "/ //) is a division that rounds the result of a normal division operation to an adjacent integer when the result is not an integer. The precedence as an operator is equivalent to multiplication / division.) 3 points out of these 4 points and their movement 7. The method according to claim 1, wherein the vector is used as a representative point and a motion vector thereof.
【請求項10】画像の左上端,右上端,左下端,右下端
の画素の座標がそれぞれ (0,0),(r,0),(0,s),
(r,s) で表されるとき (ただし、rとsは正の整
数)、 点 (0,0),(0,s),(r,s) における動きベクトル
の水平・垂直成分が1/nの整数倍の値をとるとして
(ただし、nは正の整数)、これらの動きベクトルをn倍
したものである (u00,v00),(u02,v02),(u03,v0
3) (ただし、u00,v00,u02,v02,u03,v03は整
数) を用いて、 点 (i,j),(i+p,j),(i,j+q),(i+p,j+q)
の動きベクトルの水平・垂直成分をk倍したものである
(u0,v0),(u1,v1),(u2,v2),(u3,v3)が、 u'(x,y)=(u00・r・s+(u03-u02)x・s+(u02-u00)y・r)k///(r・
s・n) v'(x,y)=(v00・r・s+(v03-v02)x・s+(v02-v00)y・r)k///(r・
s・n) u0=u'(i,j) v0=v'(i,j) u1=u'(i+p,j) v1=v'(i+p,j) u2=u'(i,j+q) v2=v'(i,j+q) u3=u'(i+p,j+q) v3=v'(i+p,j+q) によって与えられ (ただし、「///」は通常の除算による
演算結果が整数ではない場合にこれを近隣の整数に丸め
込む除算で、演算子としての優先順位は乗除算と同
等)、これら4点の中の3点とその動きベクトルを、代
表点とその動きベクトルして使用することを特徴とする
請求項1または3ないし6に記載のフレーム間予測画像
の合成方法。
10. The coordinates of the upper left, upper right, lower left and lower right pixels of the image are (0,0), (r, 0), (0, s),
When represented by (r, s) (where r and s are positive integers), the horizontal and vertical components of the motion vector at points (0,0), (0, s), (r, s) are 1 / N
(Where n is a positive integer), and these motion vectors are multiplied by n. (U00, v00), (u02, v02), (u03, v0)
3) (where u00, v00, u02, v02, u03, and v03 are integers), the points (i, j), (i + p, j), (i, j + q), (i + p , j + q)
The horizontal and vertical components of the motion vector are multiplied by k.
(u0, v0), (u1, v1), (u2, v2), (u3, v3) are: u '(x, y) = (u00 · rs · (u03-u02) x · s + (u02− u00) y ・ r) k /// (r ・
s ・ n) v '(x, y) = (v00 ・ r ・ s + (v03-v02) x ・ s + (v02-v00) y ・ r) k /// (r ・
s ・ n) u0 = u '(i, j) v0 = v' (i, j) u1 = u '(i + p, j) v1 = v' (i + p, j) u2 = u '(i , j + q) v2 = v '(i, j + q) u3 = u' (i + p, j + q) v3 = v '(i + p, j + q) (where "/ //) is a division that rounds the result of a normal division operation to an adjacent integer when the result is not an integer. The precedence as an operator is equivalent to multiplication / division.) 3 points out of these 4 points and their movement 7. The method according to claim 1, wherein the vector is used as a representative point and a motion vector thereof.
【請求項11】画像の左上端,右上端,左下端,右下端
の画素の座標がそれぞれ (0,0),(r,0),(0,s),
(r,s) で表されるとき (ただし、rとsは正の整
数)、 点 (r,0),(0,s),(r,s) における動きベクトル
の水平・垂直成分が1/nの整数倍の値をとるとして
(ただし、nは正の整数)、これらの動きベクトルをn倍
したものである (u01,v01),(u02,v02),(u03,v0
3) (ただし、u01,v01,u02,v02,u03,v03は整
数) を用いて、 点 (i,j),(i+p,j),(i,j+q),(i+p,j+q)
の動きベクトルの水平・垂直成分をk倍したものである
(u0,v0),(u1,v1),(u2,v2),(u3,v3)が、 u'(x,y)=(u03・r・s+(u01-u03)(r-x)・s+(u02-u03)(s-y)・
r)k///(r・s・n) v'(x,y)=(v03・r・s+(v01-v03)(r-x)・s+(v02-v03)(s-y)・
r)k///(r・s・n) u0=u'(i,j) v0=v'(i,j) u1=u'(i+p,j) v1=v'(i+p,j) u2=u'(i,j+q) v2=v'(i,j+q) u3=u'(i+p,j+q) v3=v'(i+p,j+q) によって与えられ (ただし、「///」は通常の除算による
演算結果が整数ではない場合にこれを近隣の整数に丸め
込む除算で、演算子としての優先順位は乗除算と同
等)、これら4点の中の3点とその動きベクトルを、代
表点とその動きベクトルして使用することを特徴とする
請求項1または3ないし6に記載のフレーム間予測画像
の合成方法。
11. The coordinates of the upper left, upper right, lower left, and lower right pixels of the image are (0,0), (r, 0), (0, s),
When represented by (r, s) (where r and s are positive integers), the horizontal and vertical components of the motion vector at points (r, 0), (0, s), (r, s) are 1 / N
(Where n is a positive integer), and these motion vectors are multiplied by n. (U01, v01), (u02, v02), (u03, v0)
3) (where u01, v01, u02, v02, u03, and v03 are integers), the points (i, j), (i + p, j), (i, j + q), (i + p , j + q)
The horizontal and vertical components of the motion vector are multiplied by k.
(u0, v0), (u1, v1), (u2, v2), (u3, v3) are u '(x, y) = (u03urs ・ (u01-u03) (rx) ・ s + ( u02-u03) (sy) ・
r) k /// (r ・ s ・ n) v '(x, y) = (v03 ・ rs ・ s + (v01-v03) (rx) ・ s + (v02-v03) (sy) ・
r) k /// (r ・ s ・ n) u0 = u '(i, j) v0 = v' (i, j) u1 = u '(i + p, j) v1 = v' (i + p , j) u2 = u '(i, j + q) v2 = v' (i, j + q) u3 = u '(i + p, j + q) v3 = v' (i + p, j + q (Where "///" is a division in which the result of ordinary division is rounded to the nearest integer when the result of the operation is not an integer, and the precedence as an operator is equivalent to multiplication / division.) 7. The method according to claim 1, wherein three of the points and their motion vectors are used as representative points and their motion vectors.
【請求項12】画像の左上端,右上端,左下端,右下端
の画素の座標がそれぞれ (0,0),(r,0),(0,s),
(r,s) で表されるとき (ただし、rとsは正の整
数)、 点 (0,0),(r,0),(0,s),(r,s) における動き
ベクトルの水平・垂直成分が1/nの整数倍の値をとる
として (ただし、nは正の整数)、これらの動きベクト
ルをn倍したものである (u00,v00),(u01,v01),
(u02,v02),(u03,v03) (ただし、u00,v00,u0
1,v01,u02,v02,u03,v03は整数) を用いて、 u'(x,y)=((s-y)(u00・(r-x)+u01・x)+y(u02・(r-x)+u03・
x))k///(r・s・n) v'(x,y)=((s-y)(v00・(r-x)+v01・x)+y(v02・(r-x)+v03・
x))k///(r・s・n) u0=u'(i,j) v0=v'(i,j) u1=u'(i+p,j) v1=v'(i+p,j) u2=u'(i,j+q) v2=v'(i,j+q) u3=u'(i+p,j+q) v3=v'(i+p,j+q) で表される (u0,v0),(u1,v1),(u2,v2),(u3,
v3) を(ただし、「///」は通常の除算による演算結果が
整数ではない場合にこれを近隣の整数に丸め込む除算
で、演算子としての優先順位は乗除算と同等)、代表点
(i,j),(i+p,j),(i,j+q),(i+p,j+q) の
動きベクトルの水平・垂直成分をk倍したものとして使
用することを特徴とする請求項2または7に記載のフレ
ーム間予測画像の合成方法。
12. The coordinates of the upper left, upper right, lower left and lower right pixels of the image are (0,0), (r, 0), (0, s),
When represented by (r, s) (where r and s are positive integers), the motion vectors at points (0,0), (r, 0), (0, s), (r, s) Assuming that the horizontal and vertical components take an integer multiple of 1 / n (where n is a positive integer), these motion vectors are multiplied by n (u00, v00), (u01, v01),
(u02, v02), (u03, v03) (However, u00, v00, u0
U ′ (x, y) = ((sy) (u00 · (rx) + u01 · x) + y (u02 · (rx) +) using 1, v01, u02, v02, u03, and v03 are integers. u03 ・
x)) k /// (r ・ s ・ n) v '(x, y) = ((sy) (v00 ・ (rx) + v01 ・ x) + y (v02 ・ (rx) + v03 ・
x)) k /// (r ・ s ・ n) u0 = u '(i, j) v0 = v' (i, j) u1 = u '(i + p, j) v1 = v' (i + p, j) u2 = u '(i, j + q) v2 = v' (i, j + q) u3 = u '(i + p, j + q) v3 = v' (i + p, j + (u0, v0), (u1, v1), (u2, v2), (u3,
v3) (where "///" is a division in which the result of normal division is not an integer and is rounded to a nearby integer, and the precedence as an operator is equivalent to multiplication / division).
(i, j), (i + p, j), (i, j + q), (i + p, j + q) are characterized by using the horizontal and vertical components of the motion vector multiplied by k. The method for synthesizing an inter-frame predicted image according to claim 2 or 7.
【請求項13】画像をその対角線によって分割される2
個の直角3角形に分割し、それぞれに含まれる画素に対
して独立に線形内・外挿に基づくグローバル動き補償を
行うことを特徴とする請求項1または3ないし6または
8ないし11に記載のフレーム間予測画像の合成方法。
13. An image which is divided by its diagonal 2
12. The method according to claim 1, wherein the image is divided into a plurality of right-angled triangles, and global motion compensation based on linear interpolation / extrapolation is independently performed on pixels included in each of the triangles. A method for synthesizing an inter-frame prediction image.
【請求項14】「///」が通常の除算による演算結果
を、最も近い整数に丸め込む演算として定義されること
を特徴とする請求項8ないし12に記載のフレーム間予
測画像の合成方法。
14. The method according to claim 8, wherein “///” is defined as an operation for rounding an operation result obtained by ordinary division to the nearest integer.
【請求項15】「///」が通常の除算による演算結果が
整数に1/2を加えた値である場合に、この値を0に近
づける方向に丸め込む演算として定義されることを特徴
とする請求項14に記載のフレーム間予測画像の合成方
法。
15. When "///" is a value obtained by adding 1/2 to an integer as a result of an operation performed by ordinary division, the value is defined as an operation of rounding this value toward zero. The method for synthesizing an inter-frame predicted image according to claim 14.
【請求項16】「///」が通常の除算による演算結果が
整数に1/2を加えた値である場合に、この値を0から
遠ざける方向に丸め込む演算として定義されることを特
徴とする請求項14に記載のフレーム間予測画像の合成
方法。
16. When "///" is a value obtained by adding 1/2 to an integer as a result of an operation by ordinary division, the value is defined as an operation of rounding this value away from 0. The method for synthesizing an inter-frame predicted image according to claim 14.
【請求項17】「//」が通常の除算による演算結果を、
最も近い整数に丸め込む演算として定義されることを特
徴とする請求項1ないし16に記載のフレーム間予測画
像の合成方法。
17. The expression “//” represents an operation result by ordinary division,
17. The method according to claim 1, wherein the method is defined as an operation of rounding to the nearest integer.
【請求項18】「//」が通常の除算による演算結果が整
数に1/2を加えた値である場合に、この値を0に近づ
ける方向に丸め込む演算として定義されることを特徴と
する請求項17に記載のフレーム間予測画像の合成方
法。
18. When "//" is a value obtained by adding 1/2 to an integer as a result of a normal division operation, the value is defined as an operation of rounding this value toward zero. A method for synthesizing an inter-frame prediction image according to claim 17.
【請求項19】「//」が通常の除算による演算結果が整
数に1/2を加えた値である場合に、この値を0から遠
ざける方向に丸め込む演算として定義されることを特徴
とする請求項17に記載のフレーム間予測画像の合成方
法。
19. When "//" is a value obtained by adding 1/2 to an integer as a result of a normal division operation, the value is defined as an operation of rounding this value away from zero. A method for synthesizing an inter-frame prediction image according to claim 17.
【請求項20】「//」が通常の除算による演算結果が整
数に1/2を加えた値である場合に、被除数が負の場合
は0に近づける方向,正の場合は0から遠ざける方向に
丸め込む演算として定義されることを特徴とする請求項
17に記載のフレーム間予測画像の合成方法。
20. When "//" is a value obtained by adding 1/2 to an integer as a result of a normal division operation, the dividend is closer to 0 when the dividend is negative, and away from 0 when the dividend is positive. 18. The method for synthesizing an inter-frame predicted image according to claim 17, wherein the operation is defined as an operation of rounding to an image.
【請求項21】「//」が通常の除算による演算結果が整
数に1/2を加えた値である場合に、被除数が負の場合
は0から遠ざける方向,正の場合は0に近づける方向に
丸め込む演算として定義されることを特徴とする請求項
17に記載のフレーム間予測画像の合成方法。
21. When "//" is a value obtained by adding 1/2 to an integer as a result of the normal division, the direction of the dividend is away from 0 if the dividend is negative, and the direction of approaching 0 if the dividend is positive. 18. The method for synthesizing an inter-frame predicted image according to claim 17, wherein the operation is defined as an operation of rounding to an image.
【請求項22】rを画像の横の長さ,sを画像の縦の長
さとして、pはr以下であり、かつ、2pはrより大き
く、かつ、qはs以下であり、かつ、2qはsより大き
いことを特徴とする請求項1ないし21に記載のフレー
ム間予測画像の合成方法。
22. Assuming that r is the horizontal length of the image and s is the vertical length of the image, p is equal to or less than r, 2p is greater than r, q is equal to or less than s, and 22. The method according to claim 1, wherein 2q is larger than s.
【請求項23】rを画像の横の長さ,sを画像の縦の長
さとして、p/2はrより小さく、かつ、pはr以上
で、かつ、q/2はsより小さく、かつ、qはs以上で
あることを特徴とする請求項1ないし21に記載のフレ
ーム間予測画像の合成方法。
23. Assuming that r is the horizontal length of the image and s is the vertical length of the image, p / 2 is smaller than r, p is not smaller than r, and q / 2 is smaller than s, 22. The method according to claim 1, wherein q is equal to or larger than s.
【請求項24】請求項1ないし23に記載のフレーム間
予測画像の合成方法を用いる画像の符号化方法および復
号化方法。
24. An image coding method and an image decoding method using the inter-frame prediction image synthesizing method according to claim 1. Description:
【請求項25】代表点の動きベクトルに関する情報を直
接符号化することを特徴とする、請求項1ないし7に記
載のフレーム間予測画像の合成方法を用いる画像の符号
化方法。
25. An image coding method using an inter-frame predicted image synthesizing method according to claim 1, wherein information relating to a motion vector at a representative point is directly encoded.
【請求項26】符号化データとして直接符号化されてい
る代表点の動きベクトルに関する情報を再生して用いる
ことを特徴とする、請求項1ないし7に記載のフレーム
間予測画像の合成方法を用いる画像の復号化方法。
26. The method for synthesizing an inter-frame predicted image according to claim 1, wherein information on a motion vector of a representative point directly encoded as encoded data is reproduced and used. Image decoding method.
【請求項27】画像の隅の点の動きベクトルに関する情
報を直接符号化することを特徴とする、請求項8ないし
12または14ないし16に記載のフレーム間予測画像
の合成方法を用いる画像の符号化方法。
27. An image code using the method for synthesizing an inter-frame predicted image according to claim 8, wherein information relating to a motion vector at a corner point of the image is directly encoded. Method.
【請求項28】符号化データとして直接符号化されてい
る画像の隅の点の動きベクトルに関する情報を再生して
用いることを特徴とする、請求項8ないし12または1
4ないし16に記載のフレーム間予測画像の合成方法を
用いる画像の復号化方法。
28. The apparatus according to claim 8, wherein information relating to a motion vector at a corner point of an image directly encoded as encoded data is reproduced and used.
An image decoding method using the method for synthesizing an inter-frame predicted image according to any one of Items 4 to 16.
JP2001104091A 2001-04-03 2001-04-03 Global motion compensation method Expired - Fee Related JP3777998B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001104091A JP3777998B2 (en) 2001-04-03 2001-04-03 Global motion compensation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001104091A JP3777998B2 (en) 2001-04-03 2001-04-03 Global motion compensation method

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2000285235A Division JP3520846B2 (en) 2000-09-14 2000-09-14 Image decoding apparatus and image decoding method

Publications (2)

Publication Number Publication Date
JP2001352549A true JP2001352549A (en) 2001-12-21
JP3777998B2 JP3777998B2 (en) 2006-05-24

Family

ID=18957000

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001104091A Expired - Fee Related JP3777998B2 (en) 2001-04-03 2001-04-03 Global motion compensation method

Country Status (1)

Country Link
JP (1) JP3777998B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009065333A (en) * 2007-09-05 2009-03-26 Sony Corp Image processing apparatus, image processing method, and computer program
JP2009517113A (en) * 2005-11-24 2009-04-30 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Motion compensated CT reconstruction for high contrast objects

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009517113A (en) * 2005-11-24 2009-04-30 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Motion compensated CT reconstruction for high contrast objects
JP2009065333A (en) * 2007-09-05 2009-03-26 Sony Corp Image processing apparatus, image processing method, and computer program
US8126281B2 (en) 2007-09-05 2012-02-28 Sony Corporation Image processing apparatus, method, and computer-readable medium for generating motion compensation images

Also Published As

Publication number Publication date
JP3777998B2 (en) 2006-05-24

Similar Documents

Publication Publication Date Title
JP3183155B2 (en) Image decoding apparatus and image decoding method
US7519117B2 (en) Inter-frame predicted image synthesizing method
JP2001352549A (en) Image coding and decoding method
JP3520800B2 (en) Image decoding apparatus and image decoding method
JP3183277B2 (en) Image encoding device and image encoding method
JP3520845B2 (en) Image encoding apparatus and image encoding method
JP3520846B2 (en) Image decoding apparatus and image decoding method
JP3520858B2 (en) Image encoding apparatus and image encoding method
JPH0865675A (en) Image encoding device and image decoding device
JPH10150665A (en) Method for generating predictive image, and method and device for image encoding
JP3509610B2 (en) Image encoding method and image encoding device
JP3591483B2 (en) Image decoding method and image decoding device
JP3414304B2 (en) Image decoding method and decoding device
JP2003204557A (en) Image encoding/decoding method and image encoding/ decoding apparatus

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A132

Effective date: 20050607

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050805

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060220

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090310

Year of fee payment: 3

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100310

Year of fee payment: 4

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110310

Year of fee payment: 5

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110310

Year of fee payment: 5

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120310

Year of fee payment: 6

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130310

Year of fee payment: 7

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130310

Year of fee payment: 7

FPAY Renewal fee payment (prs date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140310

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees