JP2001077697A - ディジタル信号処理装置および方法および記憶媒体 - Google Patents

ディジタル信号処理装置および方法および記憶媒体

Info

Publication number
JP2001077697A
JP2001077697A JP25349599A JP25349599A JP2001077697A JP 2001077697 A JP2001077697 A JP 2001077697A JP 25349599 A JP25349599 A JP 25349599A JP 25349599 A JP25349599 A JP 25349599A JP 2001077697 A JP2001077697 A JP 2001077697A
Authority
JP
Japan
Prior art keywords
signal sequence
frequency band
numbered
lifting
output
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
JP25349599A
Other languages
English (en)
Other versions
JP2001077697A5 (ja
JP3710342B2 (ja
Inventor
Hiroshi Kajiwara
浩 梶原
Bokuryo Kaku
睦凌 郭
Shigeo Kato
茂夫 加藤
Juichi Miyaji
壽一 宮道
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Canon Inc
Original Assignee
Canon Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Canon Inc filed Critical Canon Inc
Priority to JP25349599A priority Critical patent/JP3710342B2/ja
Publication of JP2001077697A publication Critical patent/JP2001077697A/ja
Publication of JP2001077697A5 publication Critical patent/JP2001077697A5/ja
Application granted granted Critical
Publication of JP3710342B2 publication Critical patent/JP3710342B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Compression Or Coding Systems Of Tv Signals (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

(57)【要約】 【課題】 ウェーブレットフィルタの適応的な切り替え
を行うことにより、プリスキャンや領域判定情報の伝送
を不要としたディジタル信号処理装置を提案する。 【解決手段】 入力信号系列x(n)を奇数番信号系列x(2
k)と偶数番信号系列x(2k+1)に分離し、この奇数番信号
系列x(2k)と偶数番信号系列x(2k+1)の夫々に対してリフ
ティング処理を施すことにより、第1の低周波帯域係数
信号列r(k)と第1の高周波帯域係数信号列d(k)とを生成
する第1のフィルタ手段(103,105)と、奇数番
信号系列x(2k)と偶数番信号系列x(2k+1)の夫々に対し
て、第2の低周波帯域係数信号列x(2k)と第2の高周波
帯域係数信号列(x(2k+1)-x(2k))とを生成する第2の
フィルタ手段(111)と、第1の低周波出力r(k)と第
2の低周波出力x(2k)のいずれか一方と、第1の高周波
出力d(k)と第2の高周波出力(x(2k+1)-x(2k))のいず
れか一方とを適応的に選択するセレクタ113とを具備
する。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、ディジタル信号処
理装置および方法およびこの方法を記憶した記憶媒体に
関するものである。
【0002】
【従来の技術】画像は非常に多くの情報を含んでおり、
その蓄積・伝送に於いては膨大なデータ量が問題とな
る。そこで蓄積・伝送に際しては、画像の持つ冗長性を
除く、或いはさらに視覚的に認識し難い程度で画素値に
操作を加えること等によってデータ量を削減する高能率
符号化が用いられる。
【0003】近年、エントロピーの低減効果が高いこと
や、画像の段階的伝送に適しているといった理由から、
離散ウェーブレット変換を用いた画像符号化方式の検討
が盛んである。ウェーブレット変換は有効なデータ解析
ツールとみなされ、低ビットレートで高画質を目指すス
ケーラブル画像符号化に広く使われている。
【0004】ウェーブレット変換を簡易・高速に実装す
る方法としてW.Sweldenらにより、リフティング手法(Li
fting Scheme)が提案されている。リフティング手法は
離散ウェーブレット変換を簡単なフィルタリフティング
ステップに分解して実装するものであり、この手法によ
り演算コストとメモリコストとを削減できることが知ら
れている。
【0005】従来の離散ウェーブレット変換の一例とし
てデジタル信号系列x(n)に対する5×3フィルタによる
離散ウェーブレット変換の式を以下に示す。
【0006】[EQ1]
【0007】[EQ2]
【0008】周知のように、r(n)は平滑化の程度を表す
係数であり、従って、低周波帯域出力の係数を表し、d
(n)は詳細度(detail)の程度を表す係数であり、従っ
て、高周波帯域出力の係数を表す。尚、上の[EQ1], [EQ
2]は1次元の離散ウェーブレット変換処理を示したもの
である。2次元離散ウェーブレット変換は、1次元の変
換を水平・垂直方向に順次行うものである。図1乃至図
3は、2次元離散ウェーブレット変換により変換対象画
像をLL,LH,HL,HHの4つの周波数帯域(サブ
バンド)に分解する様子を模式的に示したものである。
その詳細は公知であるので本明細書では説明を省略す
る。
【0009】この5×3フィルタをリフティング方式を
用いて構成した場合のブロック図は図4のようになる。
【0010】同図に於いて、501は偶数番のサンプル
を取り出すサンプリング回路、502は奇数番のサンプ
ルを取り出すサンプリング回路、503,505はフィ
ルタ回路、504,506は加算器、507は信号入力
端子、508,509は信号出力端子である。デジタル
信号系列x(n)は、サンプリング回路501,502によ
り、夫々、偶数アドレスのサンプルxe(= x(2n))と奇数
アドレスのサンプルxo(= x(2n+1))に分けられ、夫々に
対して、2つのフィルタ処理(図中、Lifting Step1と
Lifting Step2)が施される。このようなリフティング
を用いた場合、[EQ1]と[EQ2]の5×3フィルタは次
式のように変形される。
【0011】[EQ3]
【0012】[EQ4]
【0013】上の[EQ3],[EQ4]を前述の[EQ1],[EQ2]と比
較すると、演算回数の削減、必要なデータアクセスの削
減がなされていることが分かる。なお、この逆変換も、
変換処理と対称な処理で実現できる。図5は合成側のブ
ロック図を示したものであり、分解処理と同様に簡易な
処理であることが分かる。尚、合成処理手順の詳細は公
知であるので説明を省略する。なお、通常、前述の5×
3フィルタを実現する場合、除算による丸めの誤差を考
慮して次式のように整数化することにより完全再構成を
可能とすることが一般的である。
【0014】[EQ5]
【0015】[EQ6]
【0016】但し、上式において、floor{X}はXを超え
ない最大の整数を表すものとする。
【0017】
【発明が解決しようとする課題】しかしながら、従来の
離散ウェーブレット変換により、自然画像を主たる符号
化対象とする符号化システムを構成した場合、2値画像
及び混在画像に対しては高い圧縮性能を実現できないと
いう問題がある。これは離散ウェーブレット変換に使用
するフィルタが、変換対象画像の特性に適合していない
場合に発生する。
【0018】即ち、自然画像に対して良好なエントロピ
ー低減効果を持つフィルタを選択した場合には、選択さ
れたフィルタによる変換対象となる画像が、2値画像、
CG画像など、その特性が自然画像と異なるものとなっ
た場合には、良好なエントロピー低減効果を得ることが
できない。
【0019】このようなフィルタ不適合の問題を解決す
る手法として、変換処理の前に、画像毎に、あるいは画
像から分割された領域毎に、その特性を調べて、特性に
適切なフィルタを選択するといった方法が種々提案され
ているが、特性判定のためのプリスキャンを必要とする
ために、メモリコストの増加、処理の複雑化という問題
がある。
【0020】更に、画像毎にフィルタを選択する場合に
は、選択されたフィルタは、そのフィルタを特定する局
所的な特性に対応できないという問題がある。逆に、画
像を分割して各領域毎にフィルタを選択する場合には、
そのフィルタは、画像の局所的性質への対応が可能とな
るものの、画像の細分化による変換効率の低下、フィル
タ選択情報伝送による圧縮性能の低下などの問題を来
す。
【0021】したがって、従来の離散ウェーブレット変
換用の信号処理装置は、ユーザ要求を満たすに十分では
なかった。
【0022】本発明は、前述の問題点に鑑みてなされた
ものであり、リフティング方式により全ての離散ウェー
ブレット変換を簡単なフィルタリフティングステップに
分解できることに着目し、それぞれのリフティングステ
ップにおいて、その時点での低周波帯域信号から最適な
フィルタを選択して適用するディジタル信号処理装置お
よび方式および記憶媒体を提供することを目的とする。
【0023】本発明の第2の目的は、プリスキャンを行
わず、画像の局所的な性質に応じて適切なフィルタを適
用する、効率の良いサブバンド分解を行うディジタル信
号処理装置および方式および記憶媒体を提供することを
目的とする。
【0024】
【課題を解決するための手段】上述の課題を解決するた
めに本発明の、入力信号系列を低周波帯域の信号と高周
波帯域の信号の2つに分解するディジタル信号処理装置
は、入力信号系列を奇数番信号系列と偶数番信号系列に
分離する信号系列分離手段と、前記奇数番信号系列と偶
数番信号系列の夫々に対してリフティング処理を施すこ
とにより、第1の低周波帯域係数信号列と第1の高周波
帯域係数信号列とを生成する第1のフィルタ手段と、前
記2つの信号系列の夫々に対して、前記リフティング処
理と異なるフィルタ処理を施すことにより、第2の低周
波帯域係数信号列と第2の高周波帯域係数信号列とを生
成する第2のフィルタ手段と、前記第1の低周波帯域係
数信号列と前記第2の低周波帯域係数信号列のいずれか
一方と、前記第1の高周波帯域係数信号列と前記第2の
高周波帯域係数信号列のいずれか一方とを適応的に選択
して、夫々、低周波帯域係数信号列出力と高周波帯域係
数信号列出力として出力する適応的選択手段とを具備す
ることを特徴とする。
【0029】所謂ウエーブレットフィルタ出力の適応的
切換が実現される。
【0030】本発明の好適な一態様である請求項2に拠
れば、前記第1のフィルタ手段は、前記奇数番信号系列
と偶数番信号系列の夫々に対して、第1のリフティング
処理と第2のリフティング処理とを直列に行う。
【0031】本発明の好適な一態様である請求項3に拠
れば、前記第2のフィルタ手段は、前記奇数番信号系列
と偶数番信号系列のいずれか一方を取り出して前記第2
の低周波帯域係数信号列とし、前記奇数番信号系列と偶
数番信号系列の差分を前記第2の高周波帯域係数信号列
として前記適応的選択手段に出力する。
【0032】本発明の好適な一態様である請求項4に拠
れば、前記適応的選択手段は、この適応的選択手段によ
る以前の選択結果としての前記低周波帯域係数信号列出
力に基づいて、二値領域か多値領域かを判定する。特
に、処理時以前に生成される低周波係数値に基づいてフ
ィルタ選択を行うので、プリスキャンや領域判定情報の
伝送を必要としないという利点が生まれる。
【0033】本発明の好適な一態様である請求項5に拠
れば、前記第1のフィルタ手段は、前記奇数番信号系列
と偶数番信号系列の夫々に対して、第1のリフティング
処理を行う第1リフティング回路と、この第1リフティ
ング回路の後段に設けられ第2のリフティング処理を行
う第2リフティング回路とを有し、前記適応的選択手段
は、第1と第2の適応的選択回路を有し、前記第1の適
応的選択回路は、前記第1のリフティング回路の出力に
対して適応的選択動作を行い、前記第2のリフティング
回路は、前記第1の適応的選択回路によって選択された
結果に対して前記第2のリフティング処理を行い、前記
第2の適応的選択回路は、前記第2のリフティング回路
によりリフティング処理された結果に対して適応的選択
動作を行うことを特徴とする。
【0038】上記課題は、以下の構成のデジタル信号処
理装置によっても達成される。即ち、請求項6に記載
の、入力信号系列から所定の周波数帯域の信号を得るデ
ィジタル信号処理装置は、入力信号系列を奇数番信号系
列と偶数番信号系列に分離する信号系列分離手段と、前
記奇数番信号系列或いは偶数番信号系列にリフティング
処理を用いたフィルタ処理を施すことにより、所定の周
波数帯域の係数信号列を生成する第1のフィルタ手段
と、前記リフティング処理される信号系列に対して、前
記第1のフィルタ手段が用いるリフティング処理とは異
なるフィルタ処理を施すことにより、前記所定の周波数
帯域の係数信号列を生成する第2のフィルタ手段と前記
第1、第2フィルタ手段で得られる複数の係数信号列の
いずれか一方を適応的に選択して出力する適応的選択手
段とを具備することを特徴とする。
【0043】尚、請求項11のように、前記所定の周波
数帯域は例えば高周波帯域である。
【0044】上記の課題は、上述の信号処理装置に組み
込まれた信号処理方法、更には、その方法を実現するコ
ンピュータプログラムを記憶する記憶媒体によっても達
成される。
【0045】
【発明の実施形態】〈第1実施形態〉以下、本発明を代
表する一実施形態について図面を用いて説明する。
【0046】図6は本発明の第1の実施形態を実行する
ディジタル信号処理装置のブロック図を示したものであ
る。
【0047】同図に於いて101は偶数番サンプルを取
出すサンプリング回路、102は奇数番サンプルを取出
すサンプリング回路、103,105,111はフィル
タ回路、104,106,112は加算器、107は2
値/多値判定回路、108は信号入力端子、109,1
10は信号出力端子、113はセレクタである。
【0048】図6に示した本実施形態を実行するディジ
タル信号処理装置は、信号入力端子108から入力され
る信号x(n)(nは0から2m-1までとする)を、適応的に離
散ウェーブレット変換し、低周波帯域と高周波帯域にお
ける2つの係数(以下、「低周波帯域出力係数」、「高
周波帯域出力係数」と呼ぶ)を生成して、夫々を、端子
109と110から出力するものである。通常有限の信
号系列を離散ウェーブレット変換する場合には信号系列
の両端で信号の折り返しが必要になる。ここでは説明を
簡略化するため、信号系列両端で必要となる特殊処理は
装置外部で行われており、そのように入力信号系列が構
成されているものと仮定する。
【0049】本装置の変換対象画像は、自然画像、2値
/多値混在画像、さらに2値画像とする。即ち、自然画
像、2値/多値混在画像、さらに2値画像が、順不同に
端子108から入力される。尚、2値領域は輝度レベル
0か255で構成されるものとする。
【0050】以下、図6のブロック図を用いて、本実施
形態における各部の動作を詳細に説明する。
【0051】まず、信号入力端子108から入力される
信号x(n)は、偶数番サンプリング回路101と奇数番サ
ンプリング回路102により2つの系列に分けられる。
偶数番サンプリング回路101は入力信号x(n)の偶数番
目のサンプル、即ち、x(2k)(k=0〜m-1)を取り出す。
また、奇数番サンプリング回路102は入力信号x(n)の
奇数番目のサンプル、即ち、x(2k+1)(k=0〜m-1)を取
り出す。
【0052】フィルタ回路103は、
【0053】[EQ7]
【0054】のフィルタ演算を行う、即ち、偶数番のサ
ンプル系列からのx(2k)とx(2k+2)とを用いて、[EQ8]
【0055】の演算を行い、その演算結果を加算器10
4に出力する。加算器104は奇数番サンプリング回路
102の出力x(2k+1)とフィルタ回路103の出力([EQ
8]の出力)の和を求め出力する。従って、この加算器1
04の出力は[EQ6]に示すd(k)に等価である。一方、フ
ィルタ回路105は、
【0056】[EQ9]
【0057】の演算を、即ち、加算器104の出力系列
からd(k-1)とd(k)を用いて、[EQ10]
【0058】を求めて加算器106に出力する。加算器
106は、フィルタ回路105の出力と、偶数番サンプ
リング回路101の出力x(2k)との和を求め、出力す
る。この加算器106の出力は[EQ5]に示すr(k)に相当
する。図4と比較しても明らかなように、上述のフィル
タ103,105と加算器104,106による処理は
5×3フィルタの演算そのものである。これにより、5
×3フィルタによる低周波帯域の係数r(k)がセレクタ1
13の端子aに、高周波帯域の係数d(k)が同じく端子b
に入力される。一方、フィルタ回路111は、偶数番サ
ンプル系列から−x(2k)を求めて、これを加算器112
に出力する。加算器112は、奇数番サンプリング回路
102の出力x(2k+1)とフィルタ回路111の出力−x(2
k)の和を求め、セレクタ113に出力する。これによ
り、偶数番サンプルと奇数番サンプルの差分がセレクタ
113の入力端子dに入力される。一方、サンプリング
回路101の出力x(2k)はセレクタの入力端子cにも供
給されている。かくして、セレクタ113が各端子にお
いて入力する信号は、
【0059】端子a:フィルタ処理結果の低周波帯域の
係数、即ち、r(k) 端子b:フィルタ処理結果の高周波帯域の係数、即ち、
d(k) 端子c:偶数番サンプル、即ち、x(2k) 端子d:偶数番サンプルと奇数番サンプルの差分、即
ち、x(2k+1)−x(2k) となる。
【0063】セレクタ113は、所謂ホールドタイプの
セレクタである。即ち、出力は前回の入力信号を選択し
た結果を維持することができる。このようなホールドタ
イプのセレクタは例えばクロック同期のセレクタにより
簡単に実現できる。
【0064】図6の信号処理装置において、セレクタ1
13の出力端子配列を、109端子に低周波帯域係数出
力が、110端子に高周波帯域係数出力が現れるように
配置されている。2値/多値判定回路107は、この1
09端子の低周波帯域係数出力を入力して、信号系列x
(n)の像域判定を行う。判定回路107の出力である像
域判定信号B/Mの値が、“1”であれば2値領域である
と判断したことを、また、“0”であれば多値領域であ
ると判断したことを示している。像域判定信号B/Mはセ
レクタ113の動作を制御する。
【0065】図7,図8に、セレクタ113の一般的な
動作を示す。即ち、図7に拠れば、判定回路107が二
値判定(B/M=1)を出力すれば、入力端子cの信号が出力
端子109に現れ、入力端子dの信号が出力端子110
に現れる。また、図8に拠れば、判定回路107が多値
判定(B/M=0)を出力すれば、入力端子aの信号が出力端
子109に現れ、入力端子bの信号が出力端子110に
現れる。
【0066】判定回路107の動作を更に説明する。
【0067】上述のように、セレクタ113の出力は2
値/多値判定回路107の出力する像域判定信号B/Mに
応じて切り替わる。しかしながら、像域判定信号B/Mが
0または1いずれの場合であっても、入力系列が2値領
域、即ち、0か255の値しか出現しない場合には、低
周波帯域の係数として出力端子109に現れる信号は、
次の集合Vの13個の要素のいずれかである。
【0068】[EQ11] V = {-64,-32,0,32,64,96,128,159,191,223,255,287,31
9}
【0069】そこで、判定回路107を、低周波帯域の
係数値(即ち、端子109の信号)が集合Vに属する値
であることが所定の回数だけ連続して起こった場合に
は、判定信号B/Mを、
【0070】B/M = 1
【0071】とし、上記事象が認識できていない場合に
は B/M =0
【0072】を出力するように設定する。即ち、判定回
路107が2値領域を認識すると、セレクタ113は、
入力端子c,dの接続を、それぞれ低周波・高周波帯域
の係数として選択し、また、多値領域であると判断した
場合には、入力端子a,bの接続を、それぞれ低周波・
高周波帯域の係数として選択して出力する。以上の動作
により、画像中の2値領域/多値領域に応じてフィルタ
を選択して作られた低周波帯域係数出力、高周波帯域係
数出力が信号出力端子109、信号出力端子110から
それぞれ出力される。
【0073】〈第2実施形態〉第1実施形態は、5×3
フィルタにより生成した低周波帯域の係数の値に基づい
て2値領域・多値領域を判定する方法を採用した。ま
た、ここでの2値領域判定条件は低周波帯域の係数値が
所定回数連続してVに属するということを、2値領域判
定の条件として採用していた。しかしながら、この判定
条件のため、変換処理が多値領域から2値領域に移って
も直ぐには2値領域と判定されず、判定の遅れが生じ
る。
【0074】判定遅れの例として、図9のように変換対
象となる系列x(n)が、ある時点n=iで多値領域から2
値領域へ、また、時点n=i+8で、2値領域から多値
領域へと推移する場合を考える。
【0075】5×3フィルタの場合、ローパスフィルタ
は5タップであるため、2値領域に入った段階では、x
(i)を中心としてフィルタ処理することを考えると、図
9中で、ローパスフィルタ1が系列x(n)を処理すること
となり、したがって、端子109の低周波帯域出力係数
は集合Vに属するとは限らない。フィルタすべてが2値
領域に属するようになる時点(図中、ローパスフィルタ
2が系列x(n)を処理する時点)では、低周波帯域出力係
数がVに属するようになる。この状態が一定回数連続し
て発生した場合に、第1実施形態の判定回路107は、
2値領域と判定しているので、実際は2値領域であって
も多値領域と判定したままとなり、判定の遅れが生じる
のである。
【0076】また、図9中において、ローパスフィルタ
3のように多値領域の信号系列データがかかってしまえ
ば、端子109の低周波帯域出力係数はVに属さなくな
る。したがって、2値領域が十分に長く連続する場合で
なければ、即ち、上記所定回数以上二値領域が連続して
発生しなければ、2値領域判定が有効にならない。
【0077】この第2実施形態においては、2値/多値
判定を各リフティングステップ毎に設け、判定遅れの改
善を図ったものである。
【0078】図10は、本発明の第3の実施形態を実行
するディジタル信号処理装置のブロック図を示したもの
である。
【0079】同図に於いて、301は偶数番サンプリン
グ回路、302は奇数番サンプリング回路、303,3
05,311はフィルタ回路、304,306,312
は加算器、307,316は2値/多値判定回路、30
8は信号入力端子、309,310は信号出力端子、3
13,318はマッピング回路、314,317はセレ
クタ、315は遅延回路である。尚、サンプリング回路
301,302、フィルタ回路303,305,311
の夫々の単独の動作は、前述の第1実施形態のサンプリ
ング回路101,102、フィルタ回路103,10
5,111と同じである。
【0080】図10に示した第2実施形態に係るディジ
タル信号処理装置は信号入力端子308から入力される
信号x(n)(n=0〜2m-1)を、適応的に離散ウェーブレッ
ト変換し、低周波帯域と高周波帯域の2つの係数を生成
して、端子309,310から出力するものである。通
常、有限の信号系列を離散ウェーブレット変換する場合
には信号系列の両端で信号の折り返しが必要になる。第
2実施形態では、説明を簡略化するため、第1実施形態
と同じように、信号系列両端で必要な特殊処理は装置外
部で行われて、入力信号系列が構成されているものと仮
定する。
【0081】本装置の変換対象画像は、第一の実施形態
と同じく、自然画像、2値/多値混在画像、2値画像と
し、2値領域は輝度レベル0か255で構成されるもの
とする。
【0082】以下、図3のブロック図を用いて、第2実
施形態における各部の動作を詳細に説明する。
【0083】まず、信号入力端子308から入力される
信号x(n)は、偶数番サンプリング回路301と奇数番サ
ンプリング回路302により2つの系列に分けられる。
偶数番サンプリング回路301は入力信号x(n)の偶数番
目のサンプルx(2k)(k=0〜m-1)を取り出す。また、
奇数番サンプリング回路302は入力信号x(n)の奇数番
目のサンプルx(2k+1)(k=0〜m-1)を取り出す。
【0084】フィルタ回路303は、偶数番のサンプル
系列からx(2k)とx(2k+2)を用いて、
【0085】[EQ12]
【0086】の演算を行い、演算結果を加算器304に
出力する。加算器304は奇数番サンプリング回路30
2の出力x(2k+1)とフィルタ回路303の出力の和を求
め出力する。この加算器304の出力は前述の[EQ6]に
示すd(k)に相当する。一方、フィルタ回路311は偶数
番のサンプル系列から−x(2k)を求めて加算器312に
出力する。加算器312はフィルタ回路311の出力値
と奇数番サンプリング回路302の出力x(2k+1)の和
(即ち、x(2k+1)−x(2k))を求めて、マッピング回路3
13に出力する。
【0087】加算器312の出力値x(2k+1)−x(2k)が、
集合V'
【0088】[EQ13] V'={-255, 0, 255}
【0089】に属するのは、x(2k),x(2k+1)が共に0ま
たは255である場合のみである。ここで、出力値x(2k
+1)−x(2k)の正負の符号を取り除いても、合成処理側で
判断できるため、マッピング回路313で、“−25
5”を“255”に置き換えることにより、完全再構成
(即ち、ロスレス)の性質を残したまま冗長性を除くこ
とができる。2値/多値判定回路307は、サンプリン
グ回路310出力の偶数番サンプルx(2k)を調べ、2値
領域か多値領域かを判断し、判定結果に応じてセレクタ
314の動作を制御する。入力系列が2値領域の場合、
x(2k)は0か255の値しか出現しない。そこで、x(2k)
が0か255である場合には、判定回路307は2値領
域であると判断して、入力端子dに入力される信号を選
択するように、制御信号B/Mをセレクタ314に出力す
る。また、上記条件が満たされない場合には、多値領域
であると判断して、判定回路307は、入力端子bから
入力される信号を選択して出力するように制御信号B/M
をセレクタ113に送る。
【0090】セレクタ314は、2値/多値判定回路3
07からの制御信号B/Mに応じて、入力端子bから入力
される信号か、入力端子dから入力される信号のいずれ
か一方を選択して、マッピング回路318並びにフィル
タ回路305に出力する。
【0091】マッピング回路318は、第2実施形態を
実行するディジタル信号処理装置を、2次元のウエーブ
レット変換に用いる場合に必要となるものであって、水
平または垂直方向のいずれか一方の方向に既に1次元変
換を適用して得られた係数(即ち、セレクタ314の出
力係数)に、他の方向に再適用する場合のみに機能する
もので、セレクタ314の出力係数において、“12
8”と“255”の入れ替えを行う。この操作によりM
SBプレーンに2値情報が集められるため、伝送の早期
段階での2値領域の再生画質向上を図ることができる。
【0092】フィルタ回路305は、第1実施形態と同
じように、セレクタ314の出力系列からd(k−1)とd
(k)を用いて、
【0093】[EQ14]
【0094】を求めて加算器306に出力する。加算器
306は、フィルタ回路305の出力と偶数番サンプリ
ング回路301の出力x(2k)の和を求め、セレクタ31
7の入力端子cに出力する。一方、遅延回路315は偶
数番サンプリング回路301の出力x(2k)を遅延させ、x
(2k-2)を2値/多値判定回路316に出力する。2値/
多値判定回路316は遅延回路315の出力するx(2k-
2)の値を調べ、x(2k-2)の値が“0”か“255”であ
る場合には、2値領域であると判断して、入力端子aで
入力される信号を選択して出力するよう制御信号B/Mを
セレクタ317に送る。また、この条件が満たされない
場合には多値領域であると判断して、入力端子cから入
力される信号を選択して出力するように制御信号B/Mを
セレクタ317に送る。セレクタ317は2値/多値判
定回路316からの制御信号に応じて、入力端子aから
入力される信号か、入力端子cから入力される信号のい
ずれか一方を選択して出力する。
【0095】以上の動作により、画像中の2値領域/多
値領域に応じてフィルタを選択して作られた低周波帯域
係数出力と高周波成分帯域係数出力とが、夫々、信号出
力端子309、信号出力端子310から出力される。
【0096】上記第1実施形態及び第2実施形態では、
その信号処理回路の出力を2つとし、一方を低周波帯域
係数出力と、他方を高周波帯域係数出力とした。この結
果、第1実施形態と第2実施形態の信号処理回路を、従
来の図4の信号処理回路を用いている回路基板におい
て、その従来の信号処理回路の代わりに置き換えること
ができる。
【0097】なお、第1実施形態と第2実施形態による
分解処理に対応する合成処理は上述の処理を逆順に行う
ものであり、図5に示した従来の合成処理からも容易に
わかるので説明を省略する。
【0098】〈その他の実施形態〉本発明は上述した実
施形態に限定されるものではない。
【0099】上述の第1,第2実施形態では、低周波帯
域出力係数と高周波帯域出力係数の両方を選択する構成
を説明したが、本発明はこれに限られず、高周波帯域出
力係数のみ、或いは低周波帯域出力係数のみを選択出力
するように変更することも可能である。
【0100】例えば第1、第2の実施形態に於いては、
5×3フィルタをベースとした適応的離散ウェーブレッ
ト変換を示したが、2×10,2×6,9×7フィルタ
といった、その他のフィルタをベースとして適応的な離
散ウェーブレット変換を構成しても全く構わない。
【0101】また、変換対象の画像は1画素8ビットの
画像であるとして説明したが、8ビット以上の画像にも
適用しても構わないことは言うまでもない。同様にカラ
ー画像への適応、静止画像以外でも音声信号、動画像へ
の適用も本発明の範疇に含まれる。
【0102】また、上述の実施形態においては2値/多
値領域判定をベースとしてフィルタ処理の切り替えを行
ったが、本発明はこれに限定されるものでなく、その
他、例えば、低周波帯域の係数の分散によってフィルタ
を切り替えたり、係数値の連続性に着目してフィルタの
切り替えを行っても構わない。
【0103】なお、本発明は複数の機器(例えばホスト
コンピュータ、インターフェース機器、リーダ、プリン
タ等)から構成されるシステムの一部として適用して
も、1つの機器(例えば複写機、ファクシミリ装置、デ
ジタルカメラ等)からなる装置の1部に適用してもよ
い。
【0104】また、本発明は上記実施形態を実現するた
めの装置及び方法のみに限定されるものではなく、上記
システム又は装置内のコンピュータ(CPUあるいはM
PU)に、上記実施形態を実現するためのソフトウエア
のプログラムコードを供給し、このプログラムコードに
従って上記システムあるいは装置のコンピュータが上記
各種デバイスを動作させることにより上記実施形態を実
現する場合も本発明の範疇に含まれる。
【0105】またこの場合、前記ソフトウエアのプログ
ラムコード自体が上記実施形態の機能を実現することに
なり、そのプログラムコード自体、及びそのプログラム
コードをコンピュータに供給するための手段、具体的に
は上記プログラムコードを格納した記憶媒体は本発明の
範疇に含まれる。
【0106】この様なプログラムコードを格納する記憶
媒体としては、例えばフロッピーディスク、ハードディ
スク、光ディスク、光磁気ディスク、CD−ROM、磁
気テープ、不揮発性のメモリカード、ROM等を用いる
ことができる。
【0107】また、上記コンピュータが、供給されたプ
ログラムコードのみに従って各種デバイスを制御するこ
とにより、上記実施形態の機能が実現される場合だけで
はなく、上記プログラムコードがコンピュータ上で稼動
しているOS(オペレーティングシステム)、あるいは
他のアプリケーションソフト等と共同して上記実施形態
が実現される場合にもかかるプログラムコードは本発明
の範疇に含まれる。
【0108】更に、この供給されたプログラムコード
が、コンピュータの機能拡張ボードやコンピュータに接
続された機能拡張ユニットに備わるメモリに格納された
後、そのプログラムコードの指示に基づいてその機能拡
張ボードや機能拡張ユニットに備わるCPU等が実際の
処理の一部または全部を行い、その処理によって上記実
施形態が実現される場合も本発明の範疇に含まれる。
【0109】
【発明の効果】以上説明した様に本発明によれば、離散
ウェーブレット変換を簡単なフィルタリフティングステ
ップに分解して実装し、それぞれのリフティングステッ
プにおいて、各時点での低周波帯域信号から最適なフィ
ルタを選択して適用することにより、プリスキャンを行
わず、画像の局所的な性質に応じた効率の良いサブバン
ド分解を行うディジタル信号処理装置および方式および
記憶媒体を実現することができる。
【図面の簡単な説明】
【図1】 一般的な2次元離散ウェーブレット変換の対
象となる画像の一例を示す図である。
【図2】 図1の画像を低周波帯域と高周波帯域との分
解した様子を模式的に示す図である。
【図3】 図2の分解された画像を更に分解する様子を
模式的に示す図である。リフティングを用いた5×3フ
ィルタによる分解処理のブロック図である。
【図4】 従来技術に係る、リフティングを用いた5×
3フィルタによる合成処理のブロック図である。
【図5】 図の処理回路により変換された信号系列を復
元する処理のブロック図である。
【図6】 本発明に係わる第1の実施形態の信号処理装
置のブロック図である。
【図7】 第1の実施形態の信号処理装置に用いられて
いるセレクタの動作を示す図である。
【図8】 第1の実施形態の信号処理装置に用いられて
いるセレクタの動作を示す図である。
【図9】 第1実施形態の領域判定の遅れを示す図であ
る。
【図10】 本発明に係わる第2の実施形態の信号処理
装置のブロック図である。
【符号の説明】
108,507,308 信号入力端子 101,501,301 偶数番サンプリング回路 102,502,302 奇数番サンプリング回路 103,105,111,503,505,303,3
05,311 フィルタ回路 104,106,112,504,506,304,3
06,312 加算器 107,307,316 2値/多値判定回路 113,314,317 セレクタ 315 遅延回路 313,318 マッピング回路、 109,110,508,509,309,310 信
号出力端子
───────────────────────────────────────────────────── フロントページの続き (72)発明者 宮道 壽一 埼玉県久喜市中央1−15−4 Fターム(参考) 5C059 MA24 PP01 PP02 SS20 TC00 TD12 UA12 UA14 UA15 5J064 AA02 BA16 BC08 BC12 BC25 BD03

Claims (17)

    【特許請求の範囲】
  1. 【請求項1】 入力信号系列を低周波帯域の信号と高周
    波帯域の信号の2つに分解するディジタル信号処理装置
    であって、 入力信号系列を奇数番信号系列と偶数番信号系列に分離
    する信号系列分離手段と、 前記奇数番信号系列と偶数番信号系列の夫々に対してリ
    フティング処理を施すことにより、第1の低周波帯域係
    数信号列と第1の高周波帯域係数信号列とを生成する第
    1のフィルタ手段と、 前記2つの信号系列の夫々に対して、前記リフティング
    処理と異なるフィルタ処理を施すことにより、第2の低
    周波帯域係数信号列と第2の高周波帯域係数信号列とを
    生成する第2のフィルタ手段と、 前記第1の低周波帯域係数信号列と前記第2の低周波帯
    域係数信号列のいずれか一方と、前記第1の高周波帯域
    係数信号列と前記第2の高周波帯域係数信号列のいずれ
    か一方とを適応的に選択して、夫々、低周波帯域係数信
    号列出力と高周波帯域係数信号列出力として出力する適
    応的選択手段とを具備することを特徴とするディジタル
    信号処理装置。
  2. 【請求項2】 前記第1のフィルタ手段は、前記奇数番
    信号系列と偶数番信号系列の夫々に対して、第1のリフ
    ティング処理と第2のリフティング処理とを直列に行う
    ことを特徴とする請求項1に記載のディジタル信号処理
    装置。
  3. 【請求項3】 前記第2のフィルタ手段は、前記奇数番
    信号系列と偶数番信号系列のいずれか一方を取り出して
    前記第2の低周波帯域係数信号列とし、前記奇数番信号
    系列と偶数番信号系列の差分を前記第2の高周波帯域係
    数信号列として前記適応的選択手段に出力することを特
    徴とする請求項2に記載のディジタル信号処理装置。
  4. 【請求項4】 前記適応的選択手段は、この適応的選択
    手段による以前の選択結果としての前記低周波帯域係数
    信号列出力に基づいて、二値領域か多値領域かを判定す
    ることを特徴とする請求項1乃至3のいずれかに記載の
    ディジタル信号処理装置。
  5. 【請求項5】 前記第1のフィルタ手段は、前記奇数番
    信号系列と偶数番信号系列の夫々に対して、第1のリフ
    ティング処理を行う第1リフティング回路と、この第1
    リフティング回路の後段に設けられ第2のリフティング
    処理を行う第2リフティング回路とを有し、 前記適応的選択手段は、第1と第2の適応的選択回路を
    有し、 前記第1の適応的選択回路は、前記第1のリフティング
    回路の出力に対して適応的選択動作を行い、 前記第2のリフティング回路は、前記第1の適応的選択
    回路によって選択された結果に対して前記第2のリフテ
    ィング処理を行い、 前記第2の適応的選択回路は、前記第2のリフティング
    回路によりリフティング処理された結果に対して適応的
    選択動作を行うことを特徴とする請求項1記載のディジ
    タル信号処理装置。
  6. 【請求項6】 入力信号系列から所定の周波数帯域の信
    号を得るディジタル信号処理装置であって、 入力信号系列を奇数番信号系列と偶数番信号系列に分離
    する信号系列分離手段と、 前記奇数番信号系列或いは偶数番信号系列にリフティン
    グ処理を用いたフィルタ処理を施すことにより、所定の
    周波数帯域の係数信号列を生成する第1のフィルタ手段
    と、 前記リフティング処理される信号系列に対して、前記第
    1のフィルタ手段が用いるリフティング処理とは異なる
    フィルタ処理を施すことにより、前記所定の周波数帯域
    の係数信号列を生成する第2のフィルタ手段と前記第
    1、第2フィルタ手段で得られる複数の係数信号列のい
    ずれか一方を適応的に選択して出力する適応的選択手段
    とを具備することを特徴とするディジタル信号処理装
    置。
  7. 【請求項7】 前記第1のフィルタ手段は、前記奇数番
    信号系列と偶数番信号系列の夫々に対して、第1のリフ
    ティング処理と第2のリフティング処理とを直列に行う
    ことにより、第1の低周波帯域係数信号列と第1の高周
    波帯域係数信号列とを出力することを特徴とする請求項
    6に記載のディジタル信号処理装置。
  8. 【請求項8】 前記第2のフィルタ手段は、前記奇数番
    信号系列と偶数番信号系列のいずれか一方を取り出して
    第2の低周波帯域係数信号列とし、前記奇数番信号系列
    と偶数番信号系列の差分を第2の高周波帯域係数信号列
    として前記適応的選択手段に出力することを特徴とする
    請求項6に記載のディジタル信号処理装置。
  9. 【請求項9】 前記適応的選択手段は、この適応的選択
    手段による以前の選択結果としての前記低周波帯域係数
    信号列出力に基づいて、二値領域か多値領域かを判定す
    ることを特徴とする請求項6乃至8のいずれかに記載の
    ディジタル信号処理装置。
  10. 【請求項10】 前記第1のフィルタ手段は、前記奇数
    番信号系列と偶数番信号系列の夫々に対して、第1のリ
    フティング処理を行う第1リフティング回路と、この第
    1リフティング回路の後段に設けられ第2のリフティン
    グ処理を行う第2リフティング回路とを有し、 前記適応的選択手段は、第1と第2の適応的選択回路を
    有し、 前記第1の適応的選択回路は、前記第1のリフティング
    回路の出力に対して適応的選択動作を行い、 前記第2のリフティング回路は、前記第1の適応的選択
    回路によって選択された結果に対して前記第2のリフテ
    ィング処理を行い、 前記第2の適応的選択回路は、前記第2のリフティング
    回路によりリフティング処理された結果に対して適応的
    選択動作を行うことを特徴とする請求項6記載のディジ
    タル信号処理装置。
  11. 【請求項11】 前記所定の周波数帯域は高周波帯域で
    あることを特徴とする請求項6に記載のデジタル信号処
    理装置。
  12. 【請求項12】 入力信号系列から所定の周波数帯域の
    信号を得るディジタル信号処理方法であって、 入力信号系列を奇数番信号系列と偶数番信号系列に分離
    する信号系列分離工程と、 前記奇数番信号系列或いは偶数番信号系列にリフティン
    グ処理を用いたフィルタ処理を施すことにより、所定の
    周波数帯域の係数信号列を生成する第1のフィルタ工程
    と、 前記リフティング処理される信号系列に対して、前記第
    1のフィルタ手段が用いるリフティング処理とは異なる
    フィルタ処理を施すことにより、前記所定の周波数帯域
    の係数信号列を生成する第2のフィルタ工程と前記第
    1、第2フィルタ工程で得られる複数の係数信号列のい
    ずれか一方を適応的に選択して出力する適応的選択工程
    とを具備することを特徴とするディジタル信号処理方
    法。
  13. 【請求項13】 前記第1のフィルタ工程は、前記奇数
    番信号系列と偶数番信号系列の夫々に対して、第1のリ
    フティング処理と第2のリフティング処理とを直列に行
    うことにより、第1の低周波帯域係数信号列と第1の高
    周波帯域係数信号列とを出力することを特徴とする請求
    項12に記載のディジタル信号処理方法。
  14. 【請求項14】 前記第2のフィルタ工程は、前記奇数
    番信号系列と偶数番信号系列のいずれか一方を取り出し
    て第2の低周波帯域係数信号列とし、前記奇数番信号系
    列と偶数番信号系列の差分を第2の高周波帯域係数信号
    列として前記適応的選択工程に出力することを特徴とす
    る請求項12に記載のディジタル信号処理方法。
  15. 【請求項15】 前記適応的選択工程は、この適応的選
    択工程による以前の選択結果としての前記低周波帯域係
    数信号列出力に基づいて、二値領域か多値領域かを判定
    することを特徴とする請求項12乃至14のいずれかに
    記載のディジタル信号処理方法。
  16. 【請求項16】 前記所定の周波数帯域は高周波帯域で
    あることを特徴とする請求項12に記載のデジタル信号
    処理方法。
  17. 【請求項17】 請求項2乃至16のいずれかに記載の
    デジタル信号処理方法を実施するコンピュータプログラ
    ムを記憶するプログラム記憶媒体。
JP25349599A 1999-09-07 1999-09-07 ディジタル信号処理装置および方法および記憶媒体 Expired - Fee Related JP3710342B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP25349599A JP3710342B2 (ja) 1999-09-07 1999-09-07 ディジタル信号処理装置および方法および記憶媒体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP25349599A JP3710342B2 (ja) 1999-09-07 1999-09-07 ディジタル信号処理装置および方法および記憶媒体

Publications (3)

Publication Number Publication Date
JP2001077697A true JP2001077697A (ja) 2001-03-23
JP2001077697A5 JP2001077697A5 (ja) 2004-10-28
JP3710342B2 JP3710342B2 (ja) 2005-10-26

Family

ID=17252181

Family Applications (1)

Application Number Title Priority Date Filing Date
JP25349599A Expired - Fee Related JP3710342B2 (ja) 1999-09-07 1999-09-07 ディジタル信号処理装置および方法および記憶媒体

Country Status (1)

Country Link
JP (1) JP3710342B2 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006511175A (ja) * 2002-12-20 2006-03-30 フランス テレコム ソシエテ アノニム ウェーブレット画像の符号化方法及び対応する復号化方法
WO2007058294A1 (ja) * 2005-11-18 2007-05-24 Sony Corporation 符号化装置および方法、復号装置および方法、並びに、伝送システム
US7236997B2 (en) 2001-10-31 2007-06-26 Canon Kabushiki Kaisha Filter processing apparatus and method
JP2008187574A (ja) * 2007-01-31 2008-08-14 Sony Corp 情報処理装置および方法、プログラム、並びに記録媒体
US7512539B2 (en) 2001-06-18 2009-03-31 Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. Method and device for processing time-discrete audio sampled values

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7512539B2 (en) 2001-06-18 2009-03-31 Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. Method and device for processing time-discrete audio sampled values
US7236997B2 (en) 2001-10-31 2007-06-26 Canon Kabushiki Kaisha Filter processing apparatus and method
JP2006511175A (ja) * 2002-12-20 2006-03-30 フランス テレコム ソシエテ アノニム ウェーブレット画像の符号化方法及び対応する復号化方法
WO2007058294A1 (ja) * 2005-11-18 2007-05-24 Sony Corporation 符号化装置および方法、復号装置および方法、並びに、伝送システム
JPWO2007058294A1 (ja) * 2005-11-18 2009-05-07 ソニー株式会社 符号化装置および方法、復号装置および方法、並びに、伝送システム
JP4900720B2 (ja) * 2005-11-18 2012-03-21 ソニー株式会社 符号化装置および方法、並びに、復号装置および方法
US8320693B2 (en) 2005-11-18 2012-11-27 Sony Corporation Encoding device and method, decoding device and method, and transmission system
JP2008187574A (ja) * 2007-01-31 2008-08-14 Sony Corp 情報処理装置および方法、プログラム、並びに記録媒体

Also Published As

Publication number Publication date
JP3710342B2 (ja) 2005-10-26

Similar Documents

Publication Publication Date Title
US6665444B1 (en) Image processing apparatus and method, and storage medium
US6567562B1 (en) Encoding apparatus and method
JP3339054B2 (ja) データ符号化方法及び装置、データ復号方法及び装置、並びに記録媒体
JP2000069292A (ja) 画像処理装置及び方法及び記憶媒体
WO1998005159A1 (en) Picture encoder, picture decoder, picture encoding method, picture decoding method, and medium
JP3710342B2 (ja) ディジタル信号処理装置および方法および記憶媒体
KR20090095979A (ko) 리프팅 기반 방향성 웨이블렛 변환의 적응적 갱신 방법 및갱신 장치와 이를 위한 기록 매체
JP4592035B2 (ja) デジタル信号の変換装置及びその方法
JP4127465B2 (ja) 画像伸張装置及び方法、並びに記録媒体
US20020081036A1 (en) Image encoding device, image decoding device, electronic camera and recording medium
JP4083670B2 (ja) 画像符号化装置及び画像符号化方法
JP2001285643A (ja) 画像変換装置及び方法
JP2003283839A (ja) 画像変換方法および装置
JP4816328B2 (ja) 画像処理システム、画像圧縮システム、画像編集システム、画像処理プログラムおよび画像処理装置
JP3718583B2 (ja) 画像処理装置および画像処理方法
JP5351094B2 (ja) 画像符号化方法,画像符号化装置および画像符号化プログラム
JP3470326B2 (ja) 画像符号化方法
JP4070846B2 (ja) 画像符号化装置、画像復号化装置、画像符号化方法、及び画像復号化方法
JPH11331612A (ja) 画像処理装置及び方法及び記憶媒体
JPH1188183A (ja) ウェーブレット変換装置およびその方法並びにウェーブレット逆変換装置およびその方法並びに画像符号化装置およびその方法並びに画像復号化装置およびその方法
JP2001197498A (ja) 画像処理装置及び方法及び記憶媒体
JP2001169118A (ja) 画像処理装置及びその方法、コンピュータ可読メモリ
JP2002077903A (ja) 画像符号化装置及び画像符号化方法並びに記憶媒体
JP2943925B2 (ja) 画像符号化方法
CN116614632A (zh) 一种应用于led显示控制系统的图像压缩方法、系统及装置

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20050512

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050516

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050712

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050809

R150 Certificate of patent (=grant) or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20080819

Year of fee payment: 3

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

Free format text: PAYMENT UNTIL: 20090819

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20090819

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20100819

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20110819

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20120819

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20120819

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20130819

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees