JP2011243008A - Motion estimation device and program - Google Patents
Motion estimation device and program Download PDFInfo
- Publication number
- JP2011243008A JP2011243008A JP2010114729A JP2010114729A JP2011243008A JP 2011243008 A JP2011243008 A JP 2011243008A JP 2010114729 A JP2010114729 A JP 2010114729A JP 2010114729 A JP2010114729 A JP 2010114729A JP 2011243008 A JP2011243008 A JP 2011243008A
- Authority
- JP
- Japan
- Prior art keywords
- motion
- rank
- motion estimation
- frequency band
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001514 detection method Methods 0.000 claims abstract description 72
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 238000000354 decomposition reaction Methods 0.000 claims description 63
- 239000013598 vector Substances 0.000 claims description 16
- 238000010191 image analysis Methods 0.000 claims description 13
- 238000012545 processing Methods 0.000 description 28
- 238000000034 method Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000012887 quadratic function Methods 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000012615 high-resolution technique Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Compression Or Coding Systems Of Tv Signals (AREA)
- Image Analysis (AREA)
Abstract
Description
本発明は、動画像の時間方向及び/又は空間方向のスペクトルパワーによって動画像の動き量を階層的に分析し、動画像の動きベクトルを高精度に検出する動き推定装置及びプログラムに関する。 The present invention relates to a motion estimation device and a program for hierarchically analyzing the amount of motion of a moving image by spectral power in the time direction and / or spatial direction of the moving image and detecting a motion vector of the moving image with high accuracy.
近年、撮像装置及び表示装置の高精細化が進んでおり、超解像(Super−Resolution)と称される動画像の高解像化技術が研究されている(例えば、特許文献1参照)。いわゆる8Kシステムと呼ばれるスーパーハイビジョン(SHV)のような超高精細映像、又は4Kシステムと呼ばれるデジタルシネマのような高精細映像は、従来のハイビジョン(HV)映像の4倍ないし16倍の高解像度を有するに至っている。 2. Description of the Related Art In recent years, high-definition imaging devices and display devices have been advanced, and a moving image high-resolution technique called super-resolution has been studied (for example, see Patent Document 1). Ultra-high definition video such as Super Hi-Vision (SHV) called 8K system, or high-definition video like Digital Cinema called 4K system has a resolution that is 4 to 16 times higher than that of conventional Hi-Vision (HV) video. It has come to have.
しかしながら、動画像を表示する表示装置の画面が高精細化されるほど、同じ画角で撮影した場合の動領域における1画素あたりの動きボケ量が大きくなる。 However, the higher the definition of the screen of a display device that displays a moving image, the greater the amount of motion blur per pixel in the moving region when shooting at the same angle of view.
例えば、図15(a)に示すように、ハイビジョン(HV:Hi-Vision)画面は水平1920画素×垂直1080ライン、時間60フレーム/秒であり、図15(b)に示すように、スーパーハイビジョン(SHV:Super Hi-Vision)画面は、水平7680画素×垂直4320ライン、時間60フレーム/秒である。ハイビジョン画面用の動画像と同じFOV (FOV: Field Of View)で撮像した動画像をスーパーハイビジョン画面で見ると、水平・垂直解像度ともに4倍となるため、動きのある被写体の動き量は4倍となり、動領域における1画素あたりの動きボケ量も4倍となる。特に、画面全体が大きく変化するスポーツシーン等の高速動きシーンでは、視覚的なボケ感は顕著となる。 For example, as shown in FIG. 15 (a), the Hi-Vision (HV) screen is horizontal 1920 pixels × vertical 1080 lines, time 60 frames / second, and as shown in FIG. The (SHV: Super Hi-Vision) screen is horizontal 7680 pixels × vertical 4320 lines, time 60 frames / second. When a moving image taken with the same FOV (Field Of View) as a moving image for a high-definition screen is viewed on a super high-definition screen, the horizontal and vertical resolutions are four times larger, so the amount of movement of a moving subject is four times larger. Thus, the amount of motion blur per pixel in the motion region is also quadrupled. In particular, in a high-speed motion scene such as a sports scene in which the entire screen changes greatly, the visual blur is remarkable.
尚、或る動画像の1つの画面に対して異なる解像度の画像データを階層的に複数設定して、或る解像度の画像データについて動き量のための評価値を求めるとともに、この解像度とは異なる画像データについて動き量のための評価値を求め、各評価値を加算して得られる値から最終的な動き量を決定する技法が知られている(例えば、特許文献2参照)。 Note that a plurality of image data with different resolutions are set hierarchically for one screen of a certain moving image, and an evaluation value for the amount of motion is obtained for the image data with a certain resolution, which is different from this resolution. A technique is known in which an evaluation value for a motion amount is obtained for image data, and a final motion amount is determined from a value obtained by adding the evaluation values (see, for example, Patent Document 2).
また、或る動画像の1つの画面に対して空間方向に複数階のウェーブレット変換を施して、高周波成分の領域を多く含む階数の優先度を低くする輪郭情報を抽出するとともに、動画像の1つの画面をブロック分割し、輪郭情報で示されるブロックにおけるアクティビティ(画像の局所的性質)が小さいほど優先度を高くなるように、輪郭情報で示されるブロックに同一の優先度を設定して、優先度が低いブロックから順に切り捨て処理を行って符号化データ量を制御する技法が知られている(例えば、特許文献3参照)。 Further, wavelet transform of a plurality of floors is performed in a spatial direction on one screen of a certain moving image to extract contour information that lowers the priority of the floor including a lot of high-frequency component regions. Divide one screen into blocks and set the same priority to the block indicated by the contour information so that the priority becomes higher as the activity (local nature of the image) in the block indicated by the contour information is smaller. A technique is known in which the amount of encoded data is controlled by performing a truncation process in order from the block with the lowest degree (see, for example, Patent Document 3).
前述したように、SHV画面の映像フォーマットは、水平7680画素、垂直4320ライン、時間60フレーム/秒であり、HV画面の映像フォーマットと比較して、水平及び垂直標本化周波数が、時間標本化周波数に対して相対的に増大している。 As described above, the video format of the SHV screen is horizontal 7680 pixels, vertical 4320 lines, time 60 frames / second. Compared with the video format of the HV screen, the horizontal and vertical sampling frequencies are temporal sampling frequencies. It has increased relative to.
従って、SHV画面の動領域は、同じ画角で撮像された動画像で比較した場合、HV画面と比較して大きな動き量(フレーム単位の動きを示す画素数:画素/フレーム)を示し、動領域ではフレーム間の相関が低くなり時間方向の高周波領域のパワーが高くなることが想定されるとともに、動領域のボケ量が大きくなり、空間方向の高周波領域のパワーが低くなることが想定される。符号化処理や超解像処理における複数フレーム間での動き量を推定するには、これらの想定に基づく処理が有効となる。尚、HDTV標準の動画像の動き量は、一般的に数画素/フレーム〜数十画素/フレーム程度であることが知られている。 Therefore, the moving area of the SHV screen shows a larger amount of movement (number of pixels indicating movement in units of frames: pixels / frame) than the HV screen when compared with moving images captured at the same angle of view. In the region, it is assumed that the correlation between frames is low and the power in the high frequency region in the time direction is high, the blur amount in the moving region is large, and the power in the high frequency region in the spatial direction is low. . Processing based on these assumptions is effective in estimating the amount of motion between a plurality of frames in encoding processing and super-resolution processing. Note that it is known that the amount of motion of a moving picture of the HDTV standard is generally about several pixels / frame to several tens of pixels / frame.
一般に、空間高周波成分の多い絵柄は、小さいブロックサイズ(例えば、2×2画素)を用いた動き推定装置が適している。一方、空間高周波成分が少ない絵柄は、小さいブロックサイズでは誤った動き推定装置となる可能性が高くなるため、大きなブロックサイズ(例えば、16×16画素)を用いた動き推定装置が有効である。更に、大きなブロックサイズを用いた動き推定装置では、大きな動きによるボケの影響も考えられるため、大きな動き探索範囲が要求される。 In general, a motion estimation device using a small block size (for example, 2 × 2 pixels) is suitable for a pattern having a large number of spatial high-frequency components. On the other hand, since a pattern with a small number of spatial high-frequency components is likely to be an erroneous motion estimation device with a small block size, a motion estimation device using a large block size (for example, 16 × 16 pixels) is effective. Furthermore, in a motion estimation apparatus using a large block size, the influence of blur due to a large motion can be considered, and thus a large motion search range is required.
一方、このような符号化処理や超解像処理における複数フレーム間での動き量の推定において、特許文献2の技法を適用しても、時空間周波領域のパワーに基づいて階層化するものではなく予め規定した階層数で処理を行うために、処理負担が大きくなり、且つ時間方向のボケの影響が反映された動き量を検出することができない。
On the other hand, even if the technique of
また、このような符号化処理や超解像処理における複数フレーム間での動き量の推定において、特許文献3の技法を適用しても、高周波成分の領域を多く含む階数に依存して符号化データ量を取捨選択するためのブロックの優先度を決定することができるが、時空間周波領域のパワーに基づいて階層化するものではないので、動きベクトルを高精度化させることができない。
In addition, even if the technique of
したがって、大きな動き量は、動領域の動きぼけ量が増し、空間高周波領域のパワー割合を低下させる。大きな動領域面積は、動領域の時間変動面積を増し、時間高周波領域のパワー割合を増加させる。 Therefore, a large amount of motion increases the amount of motion blur in the motion region and decreases the power ratio in the spatial high frequency region. A large dynamic area increases the time variation area of the dynamic area and increases the power ratio of the time high frequency area.
一般に動き量が大きく、動きぼけ量が大きな画像の動き推定は、大きなブロックサイズと探索範囲が適している。SHV画面の動き推定では、大きな動き量から小さな動き量まで、幅広い動き量に対応する動き推定方法が望まれる。 In general, a large block size and search range are suitable for estimating the motion of an image having a large amount of motion and a large amount of motion blur. In the motion estimation of the SHV screen, a motion estimation method that supports a wide range of motion amounts from large motion amounts to small motion amounts is desired.
そこで、本発明の目的は、動画像の時間方向及び/又は空間方向のスペクトルパワーによって動画像の動き量を階層的に分析し、動画像の動き推定情報を高精度に求める動き推定装置及びプログラムを提供することにある。 Therefore, an object of the present invention is to provide a motion estimation apparatus and program for hierarchically analyzing the amount of motion of a moving image based on the spectral power in the time direction and / or spatial direction of the moving image and obtaining motion estimation information of the moving image with high accuracy. Is to provide.
前述のように、動き量と時空間方向の高周波領域のパワーとの間には一定の相関を持つことが多いため、動画像の時間及び空間スペクトルのパワー分析を行って、動き推定装置における適切なブロックサイズを推定するとともに、推定したブロックサイズを空間的に階層化して動きベクトルを検出することが有効である。 As described above, there is often a certain correlation between the amount of motion and the power in the high-frequency region in the space-time direction. It is effective to estimate a block size and detect a motion vector by spatially hierarchizing the estimated block size.
動画像の時空間方向のスペクトルを考察すると、動領域における空間方向の高周波領域のパワーは、動き量の面積が大きくなるにつれて減少する。即ち、大面積の動オブジェクトが大きな動き量を持つ動画像は、空間方向の高周波領域のパワーが小さくなるが、時間方向の高周波領域のパワーは大きくなる傾向がある。これは、画面中で大きな面積のオブジェクトが大きく動く場合は、時間方向の変動が大きくなることに起因する。 Considering the spectrum in the spatio-temporal direction of the moving image, the power of the high frequency region in the moving direction in the moving region decreases as the area of the motion amount increases. That is, a moving image in which a large-area moving object has a large amount of motion has a tendency that the power in the high-frequency region in the spatial direction decreases, but the power in the high-frequency region in the time direction tends to increase. This is because when a large area object moves greatly on the screen, the variation in the time direction becomes large.
そこで、本発明の動き推定装置及びプログラムは、動画像における動領域の面積及び動き量の推定のために、動画像の時間方向及び/又は空間方向のスペクトルパワー分析を行って、画像内の動領域では時間周波数が大きく変動し、動き量が大きい方向に画像は大きく動きぼけする特徴を利用し、時空間周波数帯域のパワーを解析することにより、画像の持つ大凡の動き量と動き方向を推定する。この大凡の動き量と動き方向の推定結果から、画像が持つ周波数帯域ごとのパワーに応じて(階層型動き推定の階数を決定して)、帯域毎にブロックサイズと動き探索範囲の水平方向および垂直方向の大きさを決定する。そして、この階層毎のブロックサイズと動き探索範囲を用いて、階層的動き推定を行う。これにより、大凡の動き量と動き方向の推定により、ブロックサイズと動き探索範囲の水平、垂直方向の大きさが階層的に異なることとなり(等価的に様々なサイズのブロックサイズ及び探索範囲を用いて動き推定を行うこととなり)、従来法よりも高確度に動き推定を行うことが可能となり、さらに、動き方向を水平、垂直だけでなく、斜め方向をも考慮に入れてブロックサイズ及び探索範囲を決定して動き推定を行うことにより、動きぼけ量や方向に応じた動き推定を従来法よりもさらに高確度に行うことが可能となる。 Therefore, the motion estimation apparatus and program of the present invention perform spectral power analysis in the time direction and / or space direction of a moving image to estimate the area of the moving region and the amount of motion in the moving image, and the motion in the image. Estimate the approximate amount of motion and direction of the image by analyzing the power of the spatio-temporal frequency band using the characteristics that the image fluctuates greatly in the direction in which the temporal frequency fluctuates greatly and the amount of motion is large. To do. From the estimation results of the approximate motion amount and motion direction, the block size and the horizontal direction of the motion search range and the horizontal direction of the motion search range for each band are determined according to the power of each frequency band of the image (determining the rank of hierarchical motion estimation). Determine the vertical size. Then, hierarchical motion estimation is performed using the block size and the motion search range for each layer. As a result, the block size and the horizontal and vertical sizes of the motion search range are hierarchically different from each other by estimating the amount of motion and the direction of motion (equivalently, block sizes and search ranges of various sizes are used. It is possible to perform motion estimation with higher accuracy than in the conventional method, and the block size and search range take into account not only horizontal and vertical motion but also diagonal directions. Thus, it is possible to perform motion estimation according to the amount of motion blur and the direction with higher accuracy than the conventional method.
即ち、本発明の動き推定装置は、動画像の動き推定を行う動き推定装置であって、複数フレームのフレーム画像列について時間方向の周波数帯域毎のパワーを算出する算出時間方向周波数分解部と、前記時間方向の周波数帯域別のパワー値と、予め定めた第1閾値とを比較し、前記第1閾値を上回るパワーとなる階数を空間周波数の空間分解階数値として決定する空間分解階数決定部と、前記空間分解階数値に基づいて、動き推定対象のフレームに対してオクターブ分解を行い、前記空間分解階数値に対応する階数毎に、水平、垂直、斜め方向の周波数帯域毎のパワー値を算出する空間方向周波数分解部と、前記空間分解階数値に対応する階数毎に算出した水平、垂直、斜め方向の周波数帯域毎のパワー値と、予め定めた第2閾値とを比較し、当該水平、垂直、斜め方向の周波数帯域毎のパワー値の全てが前記第2閾値を下回るパワーとなる階数を動き検出開始階数として決定し、前記動き検出開始階数の各階層について該水平、垂直、斜め方向の周波数帯域毎のパワー値に応じたブロックサイズ及び動き探索範囲を決定する画像解析部と、前記動き検出開始階数の階数毎に決定されたブロックサイズ及び動き探索範囲に基づいて、前記動き推定対象のフレームに対して階層型の動き推定を行う動き推定部と、を具えることを特徴とする。 That is, the motion estimation device of the present invention is a motion estimation device that performs motion estimation of a moving image, and a calculation time direction frequency decomposition unit that calculates power for each frequency band in the time direction for a frame image sequence of a plurality of frames; A space-resolved rank determining unit that compares the power value for each frequency band in the time direction with a predetermined first threshold and determines a rank that has a power that exceeds the first threshold as a spatially-resolved rank value of the spatial frequency; Based on the spatial decomposition factor, octave decomposition is performed on the frame for motion estimation, and power values are calculated for each frequency band in the horizontal, vertical, and diagonal directions for each rank corresponding to the spatial decomposition factor. A spatial direction frequency resolution unit that compares the power value for each frequency band in the horizontal, vertical, and diagonal directions calculated for each rank corresponding to the spatial resolution rank value and a predetermined second threshold value. The number of floors in which all power values for the horizontal, vertical, and diagonal frequency bands are lower than the second threshold is determined as a motion detection start floor, and the horizontal, vertical, An image analyzer that determines a block size and a motion search range according to a power value for each frequency band in an oblique direction, and the motion based on the block size and the motion search range determined for each rank of the motion detection start floor A motion estimation unit that performs hierarchical motion estimation on the estimation target frame.
また、本発明の動き推定装置において、前記画像解析部は、前記動き検出開始階数における各階層の水平、垂直、斜め方向の周波数帯域毎のパワーの割合を算出して、前記動き検出開始階数における階層毎の基準の動き検出のブロックサイズ及び動き探索範囲に対する倍数を算出し、該倍数を前記基準の動き検出のブロックサイズ及び動き探索範囲に乗じて、当該水平、垂直、斜め方向の周波数帯域毎のパワー値に応じたブロックサイズ及び動き探索範囲を決定する手段を有することを特徴とする。 In the motion estimation apparatus of the present invention, the image analysis unit calculates a ratio of power for each frequency band in the horizontal, vertical, and diagonal directions of each layer in the motion detection start floor, and in the motion detection start floor A multiple of the reference motion detection block size and motion search range for each layer is calculated, and the multiple is multiplied by the reference motion detection block size and motion search range for each horizontal, vertical, and diagonal frequency band. It has a means for determining a block size and a motion search range according to the power value.
また、本発明の動き推定装置において、前記動き推定部は、階層型の動き推定を行う際に、階層間で同位置にブロックがない場合に、階層間のブロックの面積割合に応じた加重平均を行って、前記最終的な動きベクトルを決定する手段を有することを特徴とする。 In the motion estimation device of the present invention, when the motion estimation unit performs hierarchical motion estimation, when there is no block at the same position between layers, the weighted average according to the area ratio of blocks between layers And a means for determining the final motion vector.
さらに、本発明は、動画像の動き推定を行う動き推定装置として構成するコンピュータに、複数フレームのフレーム画像列について時間方向の周波数帯域毎のパワーを算出するステップと、前記時間方向の周波数帯域別のパワー値と、予め定めた第1閾値とを比較し、前記第1閾値を上回るパワーとなる階数を空間周波数の空間分解階数値として決定するステップと、前記空間分解階数値に基づいて、動き推定対象のフレームに対してオクターブ分解を行い、前記空間分解階数値に対応する階数毎に、水平、垂直、斜め方向の周波数帯域毎のパワー値を算出するステップと、前記空間分解階数値に対応する階数毎に算出した水平、垂直、斜め方向の周波数帯域毎のパワー値と、予め定めた第2閾値とを比較し、当該水平、垂直、斜め方向の周波数帯域毎のパワー値の全てが前記第2閾値を下回るパワーとなる階数を動き検出開始階数として決定し、前記動き検出開始階数の各階層について該水平、垂直、斜め方向の周波数帯域毎のパワー値に応じたブロックサイズ及び動き探索範囲を決定するステップと、前記動き検出開始階数の階数毎に決定されたブロックサイズ及び動き探索範囲に基づいて、前記動き推定対象のフレームに対して階層型の動き推定を行うステップと、を実行させるためのプログラムとして構成される。 Further, the present invention provides a computer configured as a motion estimation device that performs motion estimation of a moving image, calculating power for each frequency band in a time direction for a frame image sequence of a plurality of frames, and for each frequency band in the time direction. Comparing the power value of the first frequency with a predetermined first threshold value and determining the rank of the power that exceeds the first threshold value as a spatial resolution factor of the spatial frequency, and based on the spatial resolution factor Performing octave decomposition on the estimation target frame, calculating a power value for each frequency band in the horizontal, vertical, and diagonal directions for each rank corresponding to the spatial decomposition rank, and corresponding to the spatial decomposition rank The power value for each frequency band in the horizontal, vertical, and diagonal directions calculated for each floor is compared with a predetermined second threshold value, and the horizontal, vertical, and diagonal directions are compared. The number of floors in which power values for every several bands are lower than the second threshold is determined as the motion detection start floor, and the power for each frequency band in the horizontal, vertical, and diagonal directions for each layer of the motion detection start floor A step of determining a block size and a motion search range according to a value, and a hierarchical type for the motion estimation target frame based on the block size and the motion search range determined for each rank of the motion detection start floor And a step of performing motion estimation.
本発明によれば、動画像における動き推定にあたって、適切なブロックサイズを推定して動き推定装置を開始することができるだけでなく、時空間方向のスペクトルパワーから画像が持つ大凡の動き量、動き方向(水平・垂直・斜めの動き方向)及び動領域を帯域毎に推定してブロックサイズ及び動き探索範囲の大きさを決定し、帯域毎にブロックサイズ及び動き探索範囲の大きさが異なる階層的に動き推定を行うことで、雑音に強く、且つ高精度の動き推定装置の計算量を削減することができるとともに、動きぼけ量や方向に応じた高確度動き推定を行うことが可能となる。 According to the present invention, when estimating motion in a moving image, not only can the motion estimation device be started by estimating an appropriate block size, but also the approximate amount of motion and motion direction that the image has from the spectral power in the spatio-temporal direction. (Horizontal / vertical / diagonal motion direction) and motion region are estimated for each band to determine the block size and the motion search range size, and the block size and motion search range size differ for each band hierarchically By performing motion estimation, it is possible to reduce the amount of calculation of a motion estimation device that is robust against noise and highly accurate, and it is possible to perform highly accurate motion estimation according to the amount and direction of motion blur.
以下、本発明による一実施例の動き推定装置について説明する。 Hereinafter, a motion estimation apparatus according to an embodiment of the present invention will be described.
一実施例の動き推定装置1として、時間方向及び空間方向の周波数解析にウェーブレット変換によるオクターブ分解処理を用いる場合について説明する。尚、時間方向及び空間方向の周波数解析には、ウェーブレット変換を用いる場合以外に、他の直交変換又はFFT(Fast Fourier transform)を用いることができるが、画像を低解像度化するにしたがってくり返し同じ処理を行って、時間方向及び空間方向の周波数解析を行う点を考慮すれば、ウェーブレット変換によるオクターブ分解処理を用いることが特に処理効率が向上する点で有利である。
As a
[装置構成]
図1に、本発明による一実施例の動き推定装置1を示す。本実施例の動き推定装置1は、時間方向周波数分解部11と、空間分解階数決定部12と、空間方向周波数分解部13と、画像解析部14と、動き推定部15とを備える。尚、各構成要素で処理するのに必要な画像データは、動き推定装置1が備える記憶部(図示せず)に適宜格納して読み出すように構成することができる。
[Device configuration]
FIG. 1 shows a
時間方向周波数分解部11は、動きベクトルの検出を行う基準フレームF(tC)及び動き探索に用いる参照フレームF(tR)を含む、時刻t=t0・・・tmにおける複数フレームのフレーム画像列F(t0),・・・,F(tC),・・・,F(tR),・・・,F(tm)を入力し、基準フレームF(tC)における全画素について、この複数フレームを時間方向に予め規定した最大階数の周波領域に分解した後、全画素における時間方向の周波数帯域毎のパワーを算出し、この複数フレームの全画素に対するパワー値で正規化し、全画素におけるnt階の各時間方向の周波数帯域別のパワー値PT(nt)を算出して空間分解階数決定部12に送出する。
Temporal frequency decomposition unit 11 includes a reference frame F to be used for reference frame F (t C) and the motion search to detect the motion vector (t R), a plurality of frames at time t = t 0 ··· t m F (t 0 ),..., F (t C ),..., F (t R ),..., F (t m ) are input, and the reference frame F (t C ) For all pixels, after decomposing the multiple frames into the frequency domain of the maximum rank specified in advance in the time direction, calculate the power for each frequency band in the time direction for all pixels and normalize the power value for all pixels in the multiple frames. The power value P T ( nt ) for each frequency band in each time direction of the nt floor in all pixels is calculated and sent to the spatial resolution
例えば、図2に示すように、時間方向周波数分解部11は、基準フレームF(tC)及び参照フレームF(tR)を含む、時刻t=t0・・・tmにおけるフレーム画像列F(t0),・・・,F(tC),・・・,F(tR),・・・,F(tm)を入力し、基準フレームF(tC)の全画素について時間方向に予め規定したNmax階(例えば、4階)の離散ウェーブレット分解を行う時間方向1次元Nmax階離散ウェーブレット分解処理部111と、基準フレームF(tC)の全画素における時間方向の周波数帯域毎のパワーを算出し、この複数フレームの全画素に対するパワー値で正規化し、算出した全画素におけるnt階の各時間方向の周波数帯域別のパワー値PT(nt)を算出して空間分解階数決定部12に送出する時間方向周波数帯域別パワー算出部112から構成することができる。
For example, as illustrated in FIG. 2, the time-direction frequency decomposition unit 11 includes a frame image sequence F at time t = t 0 ... T m including a base frame F (t C ) and a reference frame F (t R ). (T 0 ), ..., F (t C ), ..., F (t R ), ..., F (t m ) are input, and the time for all pixels of the reference frame F (t C ) A time-direction one-dimensional Nmax-order discrete wavelet
空間分解階数決定部12は、時間方向周波数分解部11によって算出したnt階の各時間方向の周波数帯域別のパワー値PT(nt)と、時間方向の高周波領域のパワーの割合が大きいほど動領域面積が大きく、且つ動き量が大きいとして判断するための予め定めた閾値Thtとを比較し、閾値Thtを上回るパワーとなる階数(尚、パワーの割合は階数によって区分される)を、時間方向の高周波領域のパワーの割合に応じた空間周波数の分解階数Ns(以下、「空間分解階数値」とも称する)として決定し、決定した空間分解階数値Nsの情報を空間方向周波数分解部13に送出する。時間方向の周波数帯域分割した階数ntの階層が大きいほど、動領域面積が大きく、且つ動き量が大きいと判断することができる。
The spatial resolution
例えば、図3に示すように、空間分解階数決定部12は、時間方向周波数分解部11によって算出したnt階の各時間方向の周波数帯域別のパワー値PT(nt)と閾値Thtとを入力して比較する比較部121と、比較部121の比較によって得られる閾値Thtを上回るパワーとなる階数を、時間方向の高周波領域のパワーの割合に応じた空間分解階数値Nsとして決定し、決定した空間分解階数値Nsの情報を空間方向周波数分解部13に送出する空間分解階数値決定部122から構成することができる。
For example, as illustrated in FIG. 3, the spatial decomposition
空間方向周波数分解部13は、基準フレームF(tC)及び参照フレームF(tR)を入力し、空間分解階数決定部12によって決定した空間分解階数値Nsに基づいて、基準フレームF(tC)及び参照フレームF(tR)の各々の全画素に対して空間Ns階離散ウェーブレット分解を行い、基準フレームF(tC)及び参照フレームF(tR)の各々における空間分解階数値Nsに対応するns階の水平、垂直、斜め方向の周波数帯域毎に平均して正規化したパワー値PH(ns),PV(ns),PD(ns)を算出し、画像解析部14に送出する。
The spatial direction
尚、基準フレームF(tC)及び参照フレームF(tR)の双方に対して空間Ns階離散ウェーブレット分解を実行することは、元画像に対して可変のブロックサイズ及び探索範囲の大きさとする階層型動き推定を行うことができる点で有利であり、特に、動き推定装置を階層的に行うための分解能の決定のためには、基準フレームF(tC)及び参照フレームF(tR)のうちの空間方向の低周波領域のパワーの割合が大きいほうを選定するのが好適となる。以下の説明では、基準フレームF(tC)及び参照フレームF(tR)の双方について空間Ns階離散ウェーブレット分解を行う例を説明する。 Note that performing the spatial Ns-order discrete wavelet decomposition on both the base frame F (t C ) and the reference frame F (t R ) has a variable block size and a search range size for the original image. It is advantageous in that hierarchical motion estimation can be performed, and in particular, for determining resolution for performing the motion estimation device hierarchically, the reference frame F (t C ) and the reference frame F (t R ) Of these, it is preferable to select the one having the larger power ratio in the low frequency region in the spatial direction. In the following description, an example will be described in which spatial Ns-order discrete wavelet decomposition is performed on both the base frame F (t C ) and the reference frame F (t R ).
例えば、図4に示すように、空間方向周波数分解部13は、基準フレームF(tC)及び参照フレームF(tR)を入力し、空間分解階数決定部12によって決定した空間分解階数値Nsに基づいて、基準フレームF(tC)及び参照フレームF(tR)の各々の全画素に対して空間Ns階離散ウェーブレット分解を行う空間方向2次元Ns階離散ウェーブレット分解処理部131と、基準フレームF(tC)及び参照フレームF(tR)の各々における空間分解階数値Nsに対応するns階の水平、垂直、斜め方向の周波数帯域毎に平均して正規化したパワー値PH(ns),PV(ns),PD(ns)を算出して画像解析部14に送出する水平・垂直・斜め方向周波数帯域別パワー算出部132から構成することができる。
For example, as shown in FIG. 4, the spatial direction
画像解析部14は、空間方向周波数分解部13によって算出したns階の水平、垂直、斜め方向の周波数帯域毎のパワー値PH(ns),PV(ns),PD(ns)と、空間方向の低周波領域のパワーの割合が大きいほど(空間方向の高周波領域のパワーの割合が小さいほど)動領域面積が大きく、且つ動き量が大きいと判断するための予め定めた閾値Thsとを入力して比較し、水平、垂直、斜め方向の周波数帯域毎のパワー値PH(ns),PV(ns),PD(ns)の全てが閾値Thsを下回るパワーとなる階数ns(尚、パワーの割合は階数によって区分される)を、動き検出を開始する階数(以下、「動き検出開始階数」と称する)として決定し、動き検出開始階数nsにおける各ns階における水平、垂直、斜め方向の周波数帯域毎のパワーの割合PHR(ns),PVR(ns),PDR(ns)を算出して、動き検出開始階数nsにおける動き検出のブロックサイズ及び動き探索範囲の倍数α,βを算出し、倍数α,βによって動き検出開始階数nsにおける動き検出のブロックサイズ及び動き探索範囲を決定し、それぞれブロックサイズ情報及び動き探索範囲情報として動き検出部15に送出する。動き検出開始階数nsの階層が小さいほど、動領域面積が大きく、且つ動き量が大きいと判断することができる。
The
例えば、図5に示すように、画像解析部14は、空間方向周波数分解部13によって算出したns階の水平、垂直、斜め方向の周波数帯域毎のパワー値PH(ns),PV(ns),PD(ns)と、予め定めた閾値Thsとを入力して比較する比較部141と、比較部141の比較によって水平、垂直、斜め方向の周波数帯域毎のパワー値PH(ns),PV(ns),PD(ns)の全てが閾値Thsを下回るパワーとなる動き検出開始階数nsを決定し、動き検出開始階数nsにおける各ns階における水平、垂直、斜め方向の周波数帯域毎のパワーの割合PHR(ns),PVR(ns),PDR(ns)を算出する水平・垂直・斜め方向パワー割合算出部142と、動き検出開始階数nsにおける動き検出のブロックサイズ及び動き探索範囲の倍数α,βを算出する水平・垂直方向倍数算出部143と、倍数α,βによって動き検出開始階数nsにおける動き検出のブロックサイズ及び動き探索範囲を決定し、それぞれブロックサイズ情報及び動き探索範囲情報として動き検出部15に送出するブロックサイズ・動き探索範囲決定部144として構成することができる。
For example, as shown in FIG. 5, the
動き検出開始階数nsにおける水平、垂直、斜め方向の周波数帯域毎のパワーの割合PHR(ns),PVR(ns),PDR(ns)は、以下の式から得られる。 The power ratios P HR (n s ), P VR (n s ), and P DR (n s ) for each frequency band in the horizontal, vertical, and diagonal directions in the motion detection start rank n s are obtained from the following equations.
水平方向パワーの割合:
PHR(ns)=PH(ns)/(PH(ns)+PV(ns)+PD(ns))
垂直方向パワーの割合:
PVR(ns)=PV(ns)/(PH(ns)+PV(ns)+PD(ns))
斜め方向パワーの割合:
PDR(ns)=PD(ns)/(PH(ns)+PV(ns)+PD(ns))
Horizontal power ratio:
P HR (n s) = P H (n s) / (P H (n s) + P V (n s) + P D (n s))
Percentage of vertical power:
P VR (n s) = P V (n s) / (P H (n s) + P V (n s) + P D (n s))
Diagonal power ratio:
P DR (n s) = P D (n s) / (P H (n s) + P V (n s) + P D (n s))
動き検出開始階数nsにおける動き検出のブロックサイズ及び動き探索範囲の倍数α,βは、以下の式から得られる。尚、aは適宜設定可能な係数であるが、等分配するにはa=1/2とすることができる。 The motion detection block size and the motion search range multiples α and β in the motion detection start rank n s are obtained from the following equations. Note that a is a coefficient that can be set as appropriate, but for equal distribution, a = 1/2.
水平方向の倍数:
α=a×(PVR(ns)+PDR(ns))
/(PHR(ns)+PVR(ns)+PDR(ns))
垂直方向の倍数:
β=(1−a)×(PHR(ns)+PDR(ns))
/(PHR(ns)+PVR(ns)+PDR(ns))
Horizontal multiple:
α = a × (P VR (n s ) + P DR (n s ))
/ (P HR ( ns ) + PVR ( ns ) + PDR ( ns ))
Vertical multiple:
β = (1−a) × (P HR (n s ) + P DR (n s ))
/ (P HR ( ns ) + PVR ( ns ) + PDR ( ns ))
従って、倍数α,βは、斜め方向パワーの割合PDR(ns)を考慮して決定されるので、後述する高確度な動き推定に寄与することになる。 Therefore, the multiples α and β are determined in consideration of the oblique power ratio P DR (n s ), which contributes to highly accurate motion estimation described later.
動き検出開始階数nsにおける動き検出のブロックサイズ情報及び動き探索範囲情報は、以下の式から得られる。動き検出開始階数nsにおける予め定めた基準のブロックサイズをBX(ns)画素、垂直BY(ns)ラインとし、動き検出開始階数nsにおける予め定めた基準の動き探索範囲を水平SX(ns)画素、垂直SSY(ns)ラインとする。 The block size information and motion search range information for motion detection at the motion detection start rank n s are obtained from the following equations. A predetermined reference block size in the motion detection start rank n s is set to B X ( ns ) pixels and a vertical B Y ( ns ) line, and a predetermined reference motion search range in the motion detection start rank n s is horizontal. It is assumed that S X ( ns ) pixels and vertical SS Y ( ns ) lines are used.
ブロックサイズ情報:水平αBX(ns)画素、垂直βBY(ns)ライン
動き探索範囲情報:水平αSX(ns)画素、垂直βSY(ns)ライン
Block size information: horizontal αB X (n s) of pixels, the vertical βB Y (n s) line motion search range information: horizontal αS X (n s) of pixels, the vertical βS Y (n s) line
動き推定部15は、動き検出開始階数nsに応じた空間方向に低周波領域の基準フレームF(tC)及び参照フレームF(tR)の画像を得るために、基準フレームF(tC)及び参照フレームF(tR)の空間方向2次元Ns階離散ウェーブレット分解したデータに対して、動き検出開始階数nsに応じた空間ns階ウェーブレットの再構成を行い、画像解析部14から得られるブロックサイズ情報及び動き探索範囲情報に従う大きさで動き推定を実行し、続いて空間ns−1階ウェーブレットの再構成を行い、当該ブロックサイズ情報及び動き探索範囲情報に従う大きさで動き推定を再度実行し、最上位の階層(即ち、元の画像レベル)にて当該ブロックサイズ情報及び動き探索範囲情報に従う大きさで動き推定を行うまで階数をデクリメントして繰り返す。この動き推定部15の動作は、基準フレームF(tC)及び参照フレームF(tR)を入力し、動き検出開始階数nsに基づいて、基準フレームF(tC)に対して順次ブロックサイズ及び探索範囲の大きさを縮小しながら動き推定を行うことと類似した処理となる。ただし、空間ns階ウェーブレット分解及び再構成を経て順次繰り返すことによる動き推定部15によれば、階層に応じて順次水平、垂直、斜め方向の周波数帯域毎のパワー値PH(ns),PV(ns),PD(ns)に基づいて決定したブロックサイズ情報及び動き探索範囲情報に従って動き推定を行うため、画像シーンに応じた動き検出開始階数nsに応じた動き推定を高確度に行うことができ、高精度化が期待できる。
The
例えば、図6に示すように、動き推定部15は、動き検出開始階数nsに応じた空間方向に低周波領域の基準フレームF(tC)及び参照フレームF(tR)の画像を得るために動き検出開始階数nsに応じた空間ns階ウェーブレットの再構成を行い、画像解析部14から得られるブロックサイズ情報及び動き探索範囲情報に従う大きさで小数画素精度のブロックマッチングによる動き推定を行う階層型動き推定部151と、この動き推定の処理を最上位の階数に対応する元の画像レベルとなるまで階数をデクリメントして繰り返すために、空間方向に1階上位のウェーブレット再構成を実行した画像を階層型動き推定部151に送出する空間1階ウェーブレット再構成部152から構成することができる。従って、階層型動き推定部151は、空間1階ウェーブレット再構成部152から得られる基準フレームF(tC)及び参照フレームF(tR)の再構成画像を用いて、動き推定処理を階層的に繰り返し、最終的な動き推定情報(例えば、動きベクトル)を決定して出力することができる。
For example, as illustrated in FIG. 6, the
以下、本発明による一実施例の動き推定装置1の動作について更に詳細に説明する。
Hereinafter, the operation of the
[装置動作]
図7は、本発明による一実施例の動き推定装置の動作を示す動作フローである。
[Device operation]
FIG. 7 is an operation flow showing the operation of the motion estimation apparatus of one embodiment according to the present invention.
ステップS1にて、動き推定装置1は、基準フレームF(tC)及び参照フレームF(tR)を含む、時刻t=t0・・・tmにおけるフレーム画像列F(t0),・・・,F(tC),・・・,F(tR),・・・,F(tm)を入力して、動き推定装置1が備える記憶部(図示せず)に適宜読み出し可能に格納する。
In step S1, the
ステップS2にて、時間方向周波数分解部11により、フレーム画像列F(t0),・・・,F(tC),・・・,F(tR),・・・,F(tm)を入力し、基準フレームF(tC)における全画素について時間方向に予め規定した最大階数の周波領域に分解した後、全画素における時間方向の周波数帯域毎のパワーを算出する。 In step S2, the time-direction frequency decomposition unit 11 causes the frame image sequence F (t 0 ),..., F (t C ),..., F (t R ) ,. ), And all pixels in the reference frame F (t C ) are decomposed into frequency regions of the maximum rank defined in advance in the time direction, and then the power for each frequency band in the time direction in all pixels is calculated.
例えば、図8に示すように、フレーム画像列F(t0),・・・,F(tC),・・・,F(tR),・・・,F(tm)における或る画素R(k,l)について、時間方向に予め規定した最大階数(Nmax)の周波領域に分解した後、全画素における時間方向周波数帯域毎のパワーを算出することができる。例えば、図9に示すように、16フレームのフレーム画像列F(t)を時間方向にNmax階に分解するとすれば、Nmax=1では、低周波領域L1及び高周波領域H1として分割することができ(図9(a)参照)、Nmax=2では、低周波領域L2及び高周波領域H1,H2として分割することができ(図9(b)参照)、Nmax=3では、低周波領域L3及び高周波領域H1,H2,H3として分割することができ(図9(c)参照)、Nmax=4では、低周波領域L4及び高周波領域H1,H2,H3,H4として分割することができる(図9(d)参照)。 For example, as shown in FIG. 8, the frame image sequence F (t 0), ···, F (t C), ···, F (t R), ···, one at F (t m) After the pixel R (k, l) is decomposed into frequency regions of the maximum rank (Nmax) defined in advance in the time direction, the power for each time direction frequency band in all the pixels can be calculated. For example, as shown in FIG. 9, if a frame image sequence F (t) of 16 frames is decomposed into Nmax floors in the time direction, when Nmax = 1, it is divided into a low frequency region L 1 and a high frequency region H 1. (See FIG. 9A), when Nmax = 2, it can be divided into the low frequency region L 2 and the high frequency regions H 1 and H 2 (see FIG. 9B), and when Nmax = 3, it is low. It can be divided into the frequency region L 3 and the high frequency regions H 1 , H 2 , H 3 (see FIG. 9C), and when Nmax = 4, the low frequency region L 4 and the high frequency regions H 1 , H 2 , H 3 and H 4 (see FIG. 9D).
ステップS3にて、空間分解階数決定部12により、算出した時間方向の周波数帯域毎のパワーPT(nt)と、時間方向の高周波領域のパワーの割合が大きいほど動領域面積が大きく、且つ動き量が大きくなると判断するための予め定めた閾値Thtとを比較し、この閾値Thtを上回るパワーとなる階数を、時間方向の高周波領域のパワーの割合に応じた空間周波数の分解階数(空間分解階数値)Nsとして決定する。
In step S3, the spatial resolution
パワーPT(nt)と閾値Thtとの比較として構成する代わりに、例えば、表1に示すように、時間方向の高周波領域のパワーの割合と空間周波数の分解階数Nsとの間で規定されるテーブルを予め保持しておき、パワーPT(nt)から空間周波数の分解階数Nsを求めるように構成することもできる。 Instead of configuring as a comparison between the power P T (n t ) and the threshold Th t , for example, as shown in Table 1, it is defined between the power ratio in the high-frequency region in the time direction and the decomposition rank Ns of the spatial frequency. It is also possible to store the table in advance and obtain the spatial frequency decomposition rank Ns from the power P T (n t ).
ステップS4にて、空間方向周波数分解部13により、基準フレームF(tC)及び参照フレームF(tR)を入力し、空間分解階数決定部12によって決定した空間周波数の分解階数Nsに基づいて、基準フレームF(tC)及び参照フレームF(tR)に対して空間Ns階離散ウェーブレット分解を実行し、基準フレームF(tC)及び参照フレームF(tR)の各々における空間分解階数値Nsに対応するns階の水平、垂直、斜め方向の周波数帯域毎のパワー値PH(ns),PV(ns),PD(ns)を算出する。尚、動き検出開始階数nsの決定のために、基準フレームF(tC)及び参照フレームF(tR)における空間周波数帯域毎のパワーの割合の大きいほうを選定することや、基準フレームF(tC)又は参照フレームF(tR)における空間周波数帯域毎のパワーの割合のみを算出してもよい。
In step S4, the base frame F (t C ) and the reference frame F (t R ) are input by the spatial
例えば、図10(a)に示すように、基準フレームF(tC)の全画素に対して空間方向に2次元2階離散ウェーブレット分解を実行して、基準フレームF(tC)及び参照フレームF(tR)の各々における空間分解階数値Nsに対応するns階の水平、垂直、斜め方向の周波数帯域毎のパワー値PH(ns),PV(ns),PD(ns)を算出することができる。また、図10(b)に示すように、基準フレームF(tC)の空間方向の低周波領域(例えば、LL2)のみを抽出して基準フレームF(tC)の低周波領域のみの画像を再構成することができる。 For example, as shown in FIG. 10 (a), by performing a two-dimensional second-order discrete wavelet decomposition in the spatial direction for all pixels of the reference frame F (t C), the reference frame F (t C) and the reference frame F n s floor horizontal corresponding to the spatial decomposition floor numerical Ns in each of (t R), a vertical power value P H (n s) of each diagonal frequency bands, P V (n s), P D ( n s ) can be calculated. Further, as shown in FIG. 10 (b), the reference frame F (t C) a low-frequency region (e.g., LL 2) spatial direction only extracted by the reference frame F low-frequency region only of the (t C) Images can be reconstructed.
続いて、画像解析部14により、空間分解階数値Nsに対応するns階の水平、垂直、斜め方向の周波数帯域毎のパワー値PH(ns),PV(ns),PD(ns)と、空間方向の低周波領域のパワーの割合が大きいほど(空間方向の高周波領域のパワーの割合が小さいほど)動領域面積が大きく、且つ動き量が大きいと判断するための予め定めた閾値Thsとを入力して比較し、この比較によって得られる閾値Thsを下回るパワーとなる階数ns(尚、パワーの割合は階数によって区分される)を、動き検出を開始する階数(以下、「動き検出開始階数」と称する)として決定し(ステップS4)、動き検出開始階数nsにおける各ns階における水平、垂直、斜め方向の周波数帯域毎のパワーの割合PHR(ns),PVR(ns),PDR(ns)を算出して(ステップS5)、動き検出開始階数nsにおける動き検出のブロックサイズ及び動き探索範囲の倍数α,βを算出し、倍数α,βによって動き検出開始階数nsにおける動き検出のブロックサイズ及び動き探索範囲を決定する(ステップS6)。ただし、ns≦Nsである。
Subsequently, the
ns階の水平、垂直、斜め方向の周波数帯域毎のパワー値PH(ns),PV(ns),PD(ns)と閾値Thsとの比較で動き検出開始階数nsを求めるように構成する代わりに、例えば、表2に示すように、ns階の水平、垂直、斜め方向の周波数帯域毎のパワー値PH(ns),PV(ns),PD(ns)における低周波領域のパワーの割合と動き検出開始階数nsとの間で規定されるテーブルを予め保持しておき、低周波領域のパワーの割合を算出して動き検出開始階数nsを求めるように構成することもできる。尚、動き検出開始階数nsが大きくなるにつれて、元の画像が低解像度化することを意味しており、元の画像に対して相対的にブロックサイズ及び動き探索範囲の大きさが増大することを意味している。例えば、空間分解能1/16,1/8,1/4,1/2とすれば、それぞれ(ブロックサイズ,動き探索範囲の大きさ)は、(16×16,水平・垂直16画素),(8×8,水平・垂直8画素),(4×4,水平・垂直4画素),(2×2,水平・垂直2画素)などである。ここで、例えば、空間分解能1/16は、元の画像における水平標本化周波数Hs及び垂直標本化周波数Vsにおいて、16画素を1画素として標本化する低解像度化を意味する。 n s floor horizontal, vertical, power value P H (n s) of each diagonal frequency bands, P V (n s), P D (n s) and Comparative motion detection start with the threshold value Th s rank n instead configured to determine a s, for example, as shown in Table 2, n s floor horizontal, vertical, power value P H for each diagonal frequency band (n s), P V ( n s), A table defined between the power ratio of the low frequency region in P D ( ns ) and the motion detection start rank n s is held in advance, and the motion ratio is started by calculating the power ratio of the low frequency region. It can also be configured to determine the rank n s . Note that as the motion detection start rank n s increases, it means that the resolution of the original image is reduced, and the block size and the size of the motion search range increase relative to the original image. Means. For example, if the spatial resolution is 1/16, 1/8, 1/4, 1/2, (block size, size of motion search range) is (16 × 16, horizontal / vertical 16 pixels), ( 8 × 8, horizontal / vertical 8 pixels), (4 × 4, horizontal / vertical 4 pixels), (2 × 2, horizontal / vertical 2 pixels), and the like. Here, for example, the spatial resolution of 1/16 means a reduction in resolution in which 16 pixels are sampled as one pixel at the horizontal sampling frequency Hs and the vertical sampling frequency Vs in the original image.
つまり、図11に示すように、空間方向の低周波領域のパワーの割合によって、動き検出開始階数nsを関連付けることができる。例えば、Ns=4のとき、低周波領域(LL4)及び高周波領域(LL4以外)のそれぞれのパワーを算出して、全体における低周波領域(LL4)の割合が、99.5%以上であれば、動き検出開始階数ns=4として4階層の低周波領域のみの画像を再構成することができる(図11(d)参照)。また、全体における低周波領域(LL4)の割合が、98.0%以上99.5%未満であれば、動き検出開始階数ns=3として3階層の低周波領域(この場合、LL3)のみの画像を再構成することができる(図11(c)参照)。同様に、全体における低周波領域(LL4)の割合が、95.0%以上98.0%未満であれば、動き検出開始階数ns=2として2階層の低周波領域(この場合、LL2)のみの画像を再構成することができ(図11(b)参照)、全体における低周波領域(LL4)の割合が、95.0%未満であれば、動き検出開始階数ns=1として1階層の低周波領域(この場合、LL1)のみの画像を再構成することができる。 That is, as shown in FIG. 11, the motion detection start rank n s can be associated with the ratio of the power in the low frequency region in the spatial direction. For example, when Ns = 4, the respective powers of the low frequency region (LL 4 ) and the high frequency region (other than LL 4 ) are calculated, and the ratio of the low frequency region (LL 4 ) in the whole is 99.5% or more. If so, it is possible to reconstruct an image of only the four-layer low-frequency region with the motion detection start rank n s = 4 (see FIG. 11D). Further, if the ratio of the low frequency region (LL 4 ) in the whole is 98.0% or more and less than 99.5%, the motion detection start rank n s = 3 is set as three layers of low frequency regions (in this case, LL 3 ) Image can be reconstructed (see FIG. 11C). Similarly, if the ratio of the low frequency region (LL 4 ) in the whole is 95.0% or more and less than 98.0%, the motion detection start rank n s = 2 is set as two layers of low frequency regions (in this case, LL 2 ) only images can be reconstructed (see FIG. 11B), and if the ratio of the low frequency region (LL 4 ) in the whole is less than 95.0%, the motion detection start rank n s = 1, it is possible to reconstruct an image of only one layer of a low frequency region (in this case, LL 1 ).
ステップS7にて、動き推定部15により、動き検出開始階数nsに基づいた基準フレームF(tC)及び参照フレームF(tR)の低周波画像に対して、基準フレームF(tC)の低周波画像を所定のブロックサイズに分割し、分割した各ブロックについて、所定の動き探索範囲の大きさで、小数画素精度のブロックマッチングによる動き推定を行う。
In step S7, the
ステップS8にて、動き推定部15により、基準フレームF(tC)及び参照フレームF(tR)に対してブロックサイズ及び動き探索範囲の大きさを縮小しながら動き推定装置を繰り返す効果を得るために、算出していた基準フレームF(tC)及び参照フレームF(tR)の空間ns階離散ウェーブレット分解データに対して動き検出開始階数nsよりも上位の階数の画像となるように空間方向に1階上位のウェーブレット再構成を実行する。
In step S8, the
動き推定部15は、空間1階ウェーブレット再構成部152から得られる基準フレームF(tC)及び参照フレームF(tR)の再構成画像を用いて、最上位の階層(即ち、元の画像レベルにおける動き推定装置)となるまで順次階数をデクリメントして動き推定装置の処理を繰り返し(ステップS9)、最終的な動きベクトルを決定して動き推定情報を出力することができる(ステップS10)。
The
例えば、図12(a)〜(d)に示すように、動き検出開始階数nsが大きいほどブロックサイズが大きくなる様子を示しており、ブロックサイズが大きいほど参照フレームF(tR)における動き探索範囲の大きさも大きくなる。 For example, as shown in FIGS. 12A to 12D, the block size increases as the motion detection start rank n s increases, and the motion in the reference frame F (t R ) increases as the block size increases. The size of the search range is also increased.
つまり、基準フレームF(tC)のns階低周波画像を水平及び垂直のブロックサイズ(Bxns,Byns)の或るブロックBns(上添え字は、階級を示す)に分割し、参照フレームF(tR)の±Sxns,±Synsの範囲(例えば、±2ブロック)で探索し、各ブロックの動きベクトルvnsを算出する。次に、基準フレームF(tC)のns−1階低周波画像を、ブロックサイズ(Bxns,Byns)で分割し、参照フレームF(tR)のns階低周波画像上の同じ位置から2×vnsだけずらした場所を中心位置とする水平及び垂直画素数としてそれぞれ±Sxns,±Synsの範囲で探索し、得られた動きベクトルに2×vnsをベクトル加算して、ns−1階低周波画像における動きベクトルvns−1を算出する。
In other words, the reference frame F (t C) n s Kaihiku frequency image horizontal and vertical block size (Bx ns, By ns) of a certain block B ns (superscript indicates the rank) of the split, the ± Sx ns of the reference frame F (t R), probed with a range of ± Sy ns (e.g., ± 2 blocks), calculates a motion vector v ns of each block. Next, a n s -1 Kaihiku frequency image of the reference frame F (t C), the block size (Bx ns, By ns) is divided by the reference frame F on n s Kaihiku frequency images (t R) each ± Sx ns as the horizontal and vertical number of pixels centered position shifted by
尚、動き推定装置は、2次関数近似による小数画素位置のブロックマッチング法を用いて行うのは、最上位の階数(即ち、1階)でのみ行うのが好適であり、式(1)で与えられる。 Note that the motion estimation apparatus is preferably performed only with the highest rank (that is, the first floor) using the block matching method of decimal pixel positions by quadratic function approximation. Given.
尚、探索位置における画素位置をxとしたとき、SSD(x)は、画素位置におけるSSD値(誤差二乗和)を表し、より具体的には、SSD(0)は中心位置におけるSSD値、SSD(−1)は中心位置から−Sx(Sy)画素の位置におけるSSD値、SSD(1)は中心位置から+Sx(Sy)画素の位置におけるSSD値を表す。式(1)から、水平又は垂直方向の小数画素精度の画素位置(小数画素位置)をそれぞれ算出することができる。例えば、図13に示すように、式(1)から2次関数近似して、小数画素位置として例えば−0.33を得ることができる。 When the pixel position at the search position is x, SSD (x) represents the SSD value (sum of squares of error) at the pixel position. More specifically, SSD (0) is the SSD value at the center position, SSD (−1) represents the SSD value at the position of −Sx (Sy) pixel from the center position, and SSD (1) represents the SSD value at the position of + Sx (Sy) pixel from the center position. From equation (1), pixel positions (decimal pixel positions) with decimal pixel precision in the horizontal or vertical direction can be calculated respectively. For example, as shown in FIG. 13, a quadratic function is approximated from the equation (1) to obtain, for example, −0.33 as the decimal pixel position.
また、ns階におけるブロックサイズは固定でないため、ns階からみてns−1階のブロックが階層間で同位置及び同サイズとならない場合も考えられる。このような場合には、図14に示すように、階層間のブロックの面積割合に応じた加重平均をとるのが好適である。例えば、ns階における面積A1〜A4の4ブロックに対して、ns−1階のブロックがまたがる場合に(面積A1+A2’+A3’+A4’)、ns階における面積A1〜A4の4ブロックにおけるそれぞれ動きベクトルがv1〜v4とすると、ns−1階の動きベクトルは、(A1・2v1+A2’・2v2+A3’・2v3+A4’・2v4)/(A1+A2+A3+A4)として求めることができる。このようにして、最上位の階数(即ち、1階)まで動き推定装置を繰り返すことにより、高精度化を図ることができる。 The block size in n s floor because it is not fixed, n s floor viewed from n s -1 floor blocks can be considered may not be the same position and the same size between tiers. In such a case, as shown in FIG. 14, it is preferable to take a weighted average corresponding to the area ratio of blocks between layers. For example, with respect to four blocks of the area A1~A4 in n s floor, n s -1 floor when the block spans (area A1 + A2 '+ A3' + A4 '), in 4 blocks of the area A1~A4 in n s floor Assuming that the motion vectors are v1 to v4, respectively, the motion vector of the n s −1 floor can be obtained as (A1 · 2v1 + A2 ′ · 2v2 + A3 ′ · 2v3 + A4 ′ · 2v4) / (A1 + A2 + A3 + A4). In this way, high accuracy can be achieved by repeating the motion estimation device up to the highest rank (that is, the first floor).
以上のように、一実施例の動き推定装置によれば、動画像における動き推定にあたって、適切なブロックサイズを推定して動き推定装置を開始することができるだけでなく、時空間方向のスペクトルパワーから画像が持つ大凡の動き量、動き方向(水平・垂直・斜めの動き方向)及び動領域を帯域毎に推定してブロックサイズ及び動き探索範囲の大きさを決定し、帯域毎にブロックサイズ及び動き探索範囲の大きさが異なる階層的に動き推定を行うことで、雑音に強く、且つ高精度の動き推定装置の計算量を削減することができるとともに、動きぼけ量や方向に応じた高確度動き推定を行うことが可能となる。例えば、SHV画面における大きな動きをするオブジェクトを高精度に推定した動き推定装置が可能となる。 As described above, according to the motion estimation device of one embodiment, in motion estimation in a moving image, not only can the motion estimation device be started by estimating an appropriate block size, but also from the spectral power in the spatio-temporal direction. Estimate the approximate amount of motion, motion direction (horizontal / vertical / diagonal motion direction) and motion area of the image for each band to determine the block size and motion search range, and block size and motion for each band. By performing motion estimation hierarchically with different search range sizes, it is possible to reduce the calculation amount of a highly accurate motion estimation device that is resistant to noise and highly accurate motion according to motion blur amount and direction. Estimation can be performed. For example, it is possible to provide a motion estimation device that accurately estimates an object that makes a large motion on the SHV screen.
また、このように高精度で効率的に求めた動きベクトルを既存の符号化装置の符号化処理や超解像処理に適用することで、更なる高品質化が期待できる。 Further, by applying the motion vector obtained with high accuracy and efficiency in this way to encoding processing and super-resolution processing of an existing encoding device, further improvement in quality can be expected.
また、上述の実施例では、各周波領域のパワーの算出のために、ウェーブレット変換を用いる例を説明したが、離散コサイン変換などの既知の直交変換や、フィルターバンク等を用いて各周波領域のパワーを算出することができる。ウェーブレット変換を用いる場合には、画像の周波数変化を高精度に捉えることができる点で優れており、離散コサイン変換を用いる場合には、既存のシステムが離散コサイン変換を用いている場合に装置構成が容易になる。 In the above-described embodiment, an example in which the wavelet transform is used to calculate the power of each frequency region has been described. However, a known orthogonal transform such as a discrete cosine transform, a filter bank, or the like is used. Power can be calculated. When using the wavelet transform, it is superior in that it can capture the frequency change of the image with high accuracy. When using the discrete cosine transform, the device configuration is used when the existing system uses the discrete cosine transform. Becomes easier.
更に、本発明の一態様として、本発明の動き推定装置をコンピュータとして構成させることができる。コンピュータに、前述した本発明の動き推定装置の各構成要素を実現させるためのプログラムは、コンピュータの内部又は外部に備えられる記憶部に記憶される。そのような記憶部は、外付けハードディスクなどの外部記憶装置、或いはROM又はRAMなどの内部記憶装置で実現することができる。コンピュータに備えられる制御部は、中央演算処理装置(CPU)などの制御で実現することができる。即ち、CPUが、各構成要素の機能を実現するための処理内容が記述されたプログラムを、適宜、記憶部から読み込んで、各構成要素の機能をコンピュータ上で実現させることができる。ここで、各構成要素の機能をハードウェアの一部で実現しても良い。 Furthermore, as one aspect of the present invention, the motion estimation apparatus of the present invention can be configured as a computer. A program for causing a computer to realize each component of the motion estimation device of the present invention described above is stored in a storage unit provided inside or outside the computer. Such a storage unit can be realized by an external storage device such as an external hard disk or an internal storage device such as ROM or RAM. The control unit provided in the computer can be realized by controlling a central processing unit (CPU) or the like. In other words, the CPU can appropriately read from the storage unit a program in which the processing content for realizing the function of each component is described, and realize the function of each component on the computer. Here, the function of each component may be realized by a part of hardware.
また、この処理内容を記述したプログラムを、例えばDVD又はCD−ROMなどの可搬型記録媒体の販売、譲渡、貸与等により流通させることができるほか、そのようなプログラムを、例えばネットワーク上にあるサーバの記憶部に記憶しておき、ネットワークを介してサーバから他のコンピュータにそのプログラムを転送することにより、流通させることができる。 In addition, the program describing the processing contents can be distributed by selling, transferring, or lending a portable recording medium such as a DVD or CD-ROM, and such a program can be distributed on a server on a network, for example. Can be distributed by transferring the program from the server to another computer via the network.
また、そのようなプログラムを実行するコンピュータは、例えば、可搬型記録媒体に記録されたプログラム又はサーバから転送されたプログラムを、一旦、自己の記憶部に記憶することができる。また、このプログラムの別の実施態様として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、更に、このコンピュータにサーバからプログラムが転送される度に、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。 In addition, a computer that executes such a program can temporarily store, for example, a program recorded on a portable recording medium or a program transferred from a server in its own storage unit. As another embodiment of the program, the computer may directly read the program from a portable recording medium and execute processing according to the program, and each time the program is transferred from the server to the computer. In addition, the processing according to the received program may be executed sequentially.
以上、具体例を挙げて本発明の実施例を詳細に説明したが、本発明の特許請求の範囲から逸脱しない限りにおいて、あらゆる変形や変更が可能であることは当業者に明らかである。 While the embodiments of the present invention have been described in detail with specific examples, it will be apparent to those skilled in the art that various modifications and changes can be made without departing from the scope of the claims of the present invention.
本発明によれば、動画像の超解像処理及び映像符号化、又はその他の画像処理で用いられる動き推定装置の精度を高めることができる。近年における動画像処理は益々高精細化しており、本発明による動き推定装置は、高精度の動き推定装置が求められる任意の動画像処理の用途に有用である。 ADVANTAGE OF THE INVENTION According to this invention, the precision of the motion estimation apparatus used by the super-resolution process and video encoding of a moving image, or other image processing can be improved. In recent years, moving image processing has become more and more precise, and the motion estimation apparatus according to the present invention is useful for any moving image processing application that requires a highly accurate motion estimation apparatus.
1 動き推定装置
11 時間方向周波数分解部
12 空間分解階数決定部
13 空間方向周波数分解部
14 画像解析部
15 動き推定部
121 比較部
122 空間分解階数値決定部
131 空間方向2次元Ns階離散ウェーブレット分解処理部
132 水平・垂直・斜め方向周波数帯域別パワー算出部
141 比較部
142 水平・垂直・斜め方向パワー割合算出部
143 水平・垂直方向倍数算出部
144 ブロックサイズ・動き探索範囲決定部
151 階層型動き推定部
152 空間1階ウェーブレット再構成部
DESCRIPTION OF
Claims (4)
複数フレームのフレーム画像列について時間方向の周波数帯域毎のパワーを算出する算出時間方向周波数分解部と、
前記時間方向の周波数帯域別のパワー値と、予め定めた第1閾値とを比較し、前記第1閾値を上回るパワーとなる階数を空間周波数の空間分解階数値として決定する空間分解階数決定部と、
前記空間分解階数値に基づいて、動き推定対象のフレームに対してオクターブ分解を行い、前記空間分解階数値に対応する階数毎に、水平、垂直、斜め方向の周波数帯域毎のパワー値を算出する空間方向周波数分解部と、
前記空間分解階数値に対応する階数毎に算出した水平、垂直、斜め方向の周波数帯域毎のパワー値と、予め定めた第2閾値とを比較し、当該水平、垂直、斜め方向の周波数帯域毎のパワー値の全てが前記第2閾値を下回るパワーとなる階数を動き検出開始階数として決定し、前記動き検出開始階数の各階層について該水平、垂直、斜め方向の周波数帯域毎のパワー値に応じたブロックサイズ及び動き探索範囲を決定する画像解析部と、
前記動き検出開始階数の階数毎に決定されたブロックサイズ及び動き探索範囲に基づいて、前記動き推定対象のフレームに対して階層型の動き推定を行う動き推定部と、
を具えることを特徴とする動き推定装置。 A motion estimation device that performs motion estimation of a moving image,
A calculation time direction frequency decomposition unit that calculates power for each frequency band in the time direction for a frame image sequence of a plurality of frames;
A space-resolved rank determining unit that compares the power value for each frequency band in the time direction with a predetermined first threshold and determines a rank that has a power that exceeds the first threshold as a spatially-resolved rank value of the spatial frequency; ,
Based on the spatial resolution factor, octave decomposition is performed on the motion estimation target frame, and a power value is calculated for each horizontal, vertical, and diagonal frequency band for each rank corresponding to the spatial resolution factor. A spatial frequency decomposition unit;
The power value for each frequency band in the horizontal, vertical, and diagonal directions calculated for each rank corresponding to the spatial resolution rank value is compared with a predetermined second threshold value, and the frequency band for each horizontal, vertical, and diagonal direction is compared. The number of floors at which all of the power values are lower than the second threshold is determined as a motion detection start floor, and for each layer of the motion detection start floor, depending on the power value for each horizontal, vertical, and diagonal frequency band An image analysis unit for determining a block size and a motion search range,
A motion estimation unit that performs hierarchical motion estimation on the motion estimation target frame based on a block size and a motion search range determined for each floor of the motion detection start floor;
A motion estimation apparatus comprising:
複数フレームのフレーム画像列について時間方向の周波数帯域毎のパワーを算出するステップと、
前記時間方向の周波数帯域別のパワー値と、予め定めた第1閾値とを比較し、前記第1閾値を上回るパワーとなる階数を空間周波数の空間分解階数値として決定するステップと、
前記空間分解階数値に基づいて、動き推定対象のフレームに対してオクターブ分解を行い、前記空間分解階数値に対応する階数毎に、水平、垂直、斜め方向の周波数帯域毎のパワー値を算出するステップと、
前記空間分解階数値に対応する階数毎に算出した水平、垂直、斜め方向の周波数帯域毎のパワー値と、予め定めた第2閾値とを比較し、当該水平、垂直、斜め方向の周波数帯域毎のパワー値の全てが前記第2閾値を下回るパワーとなる階数を動き検出開始階数として決定し、前記動き検出開始階数の各階層について該水平、垂直、斜め方向の周波数帯域毎のパワー値に応じたブロックサイズ及び動き探索範囲を決定するステップと、
前記動き検出開始階数の階数毎に決定されたブロックサイズ及び動き探索範囲に基づいて、前記動き推定対象のフレームに対して階層型の動き推定を行うステップと、
を実行させるためのプログラム。
In a computer configured as a motion estimation device that performs motion estimation of moving images,
Calculating power for each frequency band in the time direction for a frame image sequence of a plurality of frames;
Comparing the power value for each frequency band in the time direction with a predetermined first threshold value, and determining a rank that is a power that exceeds the first threshold value as a spatially resolved rank value of the spatial frequency;
Based on the spatial resolution factor, octave decomposition is performed on the motion estimation target frame, and a power value is calculated for each horizontal, vertical, and diagonal frequency band for each rank corresponding to the spatial resolution factor. Steps,
The power value for each frequency band in the horizontal, vertical, and diagonal directions calculated for each rank corresponding to the spatial resolution rank value is compared with a predetermined second threshold value, and the frequency band for each horizontal, vertical, and diagonal direction is compared. The number of floors at which all of the power values are lower than the second threshold is determined as a motion detection start floor, and for each layer of the motion detection start floor, depending on the power value for each horizontal, vertical, and diagonal frequency band Determining the determined block size and motion search range;
Performing hierarchical motion estimation on the motion estimation target frame based on the block size and motion search range determined for each floor of the motion detection start floor;
A program for running
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010114729A JP5520687B2 (en) | 2010-05-18 | 2010-05-18 | Motion estimation apparatus and program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010114729A JP5520687B2 (en) | 2010-05-18 | 2010-05-18 | Motion estimation apparatus and program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011243008A true JP2011243008A (en) | 2011-12-01 |
JP5520687B2 JP5520687B2 (en) | 2014-06-11 |
Family
ID=45409586
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010114729A Expired - Fee Related JP5520687B2 (en) | 2010-05-18 | 2010-05-18 | Motion estimation apparatus and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5520687B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015101060A1 (en) * | 2013-12-30 | 2015-07-09 | 华中科技大学 | Decomposition and estimation method for multiple motion parameters in single-arm x-ray angiographic image |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08186815A (en) * | 1994-12-28 | 1996-07-16 | Pioneer Electron Corp | Sub band coding method |
-
2010
- 2010-05-18 JP JP2010114729A patent/JP5520687B2/en not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08186815A (en) * | 1994-12-28 | 1996-07-16 | Pioneer Electron Corp | Sub band coding method |
Non-Patent Citations (2)
Title |
---|
CSNJ201010034053; 松尾康孝 外3名: '超高精細映像の高確度レジストレーション手法の一検討' 映像メディア処理シンポジウム 第14回シンポジウム資料 , 20091007, 第117〜118頁, 電子情報通信学会 画像工学研究専門委員会 * |
JPN6013052812; 松尾康孝 外3名: '超高精細映像の高確度レジストレーション手法の一検討' 映像メディア処理シンポジウム 第14回シンポジウム資料 , 20091007, 第117〜118頁, 電子情報通信学会 画像工学研究専門委員会 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015101060A1 (en) * | 2013-12-30 | 2015-07-09 | 华中科技大学 | Decomposition and estimation method for multiple motion parameters in single-arm x-ray angiographic image |
Also Published As
Publication number | Publication date |
---|---|
JP5520687B2 (en) | 2014-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8958484B2 (en) | Enhanced image and video super-resolution processing | |
JP5727873B2 (en) | Motion vector detection device, encoding device, and program thereof | |
US7881539B2 (en) | Image processing apparatus and method, program recording medium, and program | |
JPH11355778A (en) | Motion field estimation method | |
KR20140135968A (en) | Method and apparatus for performing super-resolution | |
US20090129703A1 (en) | Signal processing method, apparatus, and program | |
EP3637363B1 (en) | Image processing device, image processing method and image processing program | |
JP2004526363A (en) | Apparatus and method for achieving high accuracy and high speed search of intermediate pixels in video coding | |
US20100026897A1 (en) | Method, Apparatus, and Computer Software for Modifying Moving Images Via Motion Compensation Vectors, Degrain/Denoise, and Superresolution | |
JP5669523B2 (en) | Frame interpolation apparatus and method, program, and recording medium | |
JP5717548B2 (en) | Super-resolution auxiliary information generation device, encoding device, decoding device, and programs thereof | |
JP2002539657A (en) | Process, apparatus and use for evaluating an encoded image | |
JP5331643B2 (en) | Motion vector detection apparatus and program | |
JP5829898B2 (en) | Image space super-resolution apparatus and program | |
TW201322769A (en) | Motion vector refining device, motion interpolation apparatus and video refining method thereof | |
JP6202725B2 (en) | Image difference amount calculation device, image encoding device, and program | |
US8830394B2 (en) | System, method, and apparatus for providing improved high definition video from upsampled standard definition video | |
JP5520687B2 (en) | Motion estimation apparatus and program | |
JP5514132B2 (en) | Image reduction device, image enlargement device, and program thereof | |
US20130202199A1 (en) | Using higher order statistics to estimate pixel values in digital image processing to improve accuracy and computation efficiency | |
JP5498970B2 (en) | Image reduction device, image enlargement device, and program thereof | |
Lama et al. | Image interpolation for high-resolution display based on the complex dual-tree wavelet transform and hidden Markov model | |
JP6817784B2 (en) | Super-resolution device and program | |
JP5309097B2 (en) | Motion estimation apparatus and program | |
KR101059649B1 (en) | Encoding device, method and decoding device using adaptive interpolation filter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130110 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131011 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131029 |
|
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: 20140311 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140407 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5520687 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |