画像処理装置、画像処理方法、及び画像処理プログラム
【課題】ＪＰＥＧ２０００等の圧縮アルゴリズムにおいて、圧縮処理の対象とされる画像の特徴に応じて必要な情報を的確に残しつつ、従来より高品質の圧縮処理を可能にする。
【解決手段】特性解析部１７７は、離散ウェーブレット変換部１７３によって離散ウェーブレット変換されたタイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とし、量子化部１７４が、特性解析部１７７によって解析された特性に対応する量子化テーブルを記憶部１２から読み出し、当該読み出した量子化テーブルを用いて、上記離散ウェーブレット変換されたタイルの各サブバンドデータを量子化する。
【解決手段】特性解析部１７７は、離散ウェーブレット変換部１７３によって離散ウェーブレット変換されたタイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とし、量子化部１７４が、特性解析部１７７によって解析された特性に対応する量子化テーブルを記憶部１２から読み出し、当該読み出した量子化テーブルを用いて、上記離散ウェーブレット変換されたタイルの各サブバンドデータを量子化する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置、画像処理方法、及び画像処理プログラムに関し、特に、画像データを圧縮する技術に関する。
【背景技術】
【０００２】
近年、デジタル機器の発展に伴って、画像データ等の様々な圧縮方式が提案されているが、圧縮方式の１つとしてＪＰＥＧ（Joint Photographic Experts Group）２０００が知られている。ＪＰＥＧ２０００は画像の圧縮・展開の方法を規定したもので、ＪＰＥＧを発展させたものである。ＪＰＥＧ２０００では、画像データを複数のタイルに分割し、タイル毎に離散ウェーブレット変換、量子化、エントロピー符号化等の処理を施して符号化データを生成する。なお、特許文献１及び２に記載されているように、ＪＰＥＧ圧縮では、圧縮処理の対象とされる画像全体に基づいて量子化処理に用いる量子化テーブルを選択し、当該選択した量子化テーブルを用いて量子化を行う方式が提案されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００２−２４７３７３号公報
【特許文献２】特開２００４−１１０３６８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上記従来一般のＪＰＥＧ２０００の画像圧縮方式では、量子化処理に用いる量子化テーブルは、圧縮率の指定やタイルサイズ等に基づいて１つに定められている。また、圧縮率の制御は、符号化後で各ビットに対する重み付で優先順位を決めて制御するようになっているため、入力される画像の特徴に応じた最適な符号化は行われていない。また、上記特許文献１に記載の発明は、ＪＰＥＧ圧縮に関するものであり、ＪＰＥＧ２０００の画像圧縮方式に関するものではなく、さらに、入力された画像全体に基づいて量子化テーブルを選択するため、画像に応じた細かな圧縮処理ができない。
【０００５】
本発明は、上記の問題を解決するためになされたもので、ＪＰＥＧ２０００等の圧縮アルゴリズムにおいて、圧縮処理の対象とされる画像の特徴に応じて必要な情報を的確に残しつつ、従来よりも高品質の圧縮処理を可能にすることを目的とする。
【課題を解決するための手段】
【０００６】
本発明の請求項１に記載の発明は、画像データを輝度成分と色差成分からなる表色系に変換する表色系変換部と、
前記表色系変換部によって変換された画像データを複数のタイルに分割する分割部と、
前記分割部によって分割されたタイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶する量子化テーブル記憶部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析部と、
前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化部と、
前記量子化部によって量子化された画像データを符号化して圧縮データを出力する符号化部と
を備えた画像処理装置である。
【０００７】
また、請求項２に記載の発明は、請求項１に記載の画像処理装置であって、前記量子化部は、前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを、予め設定されている圧縮率に応じた係数を用いて算出した値に変更し、当該変更した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化するものである。
【０００８】
また、請求項５に記載の発明は、画像データを輝度成分と色差成分からなる表色系に変換する表色系変換ステップと、
前記表色系変換ステップで変換された画像データを複数のタイルに分割する分割ステップと、
前記分割ステップで分割された前記タイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換ステップと、
前記離散ウェーブレット変換ステップで離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析ステップと、
前記離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶しており、前記特性解析ステップで解析された特性に対応する量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化ステップと、
前記量子化ステップによって量子化された画像データを符号化して圧縮データを出力する符号化ステップと
を有する画像処理方法である。
【０００９】
また、請求項６に記載の発明は、コンピュータを、
画像データを輝度成分と色差成分からなる表色系に変換する表色系変換部と、
前記表色系変換部によって変換された画像データを複数のタイルに分割する分割部と、
前記分割部によって分割されたタイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶する量子化テーブル記憶部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析部と、
前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化部と、
前記量子化部によって量子化された画像データを符号化して圧縮データを出力する符号化部と
して機能させる画像処理プログラムである。
【００１０】
これらの発明では、離散ウェーブレット変換された上記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドデータの特性を当該タイル全体の特性とし、量子化処理時には、当該解析されたタイルの特性に対応する量子化テーブルを用いて当該タイルの各サブバンドデータを量子化する。これにより、ＪＰＥＧ２０００等の圧縮アルゴリズムにおいて、圧縮処理の対象とされる画像データをなす当該タイルの各サブバンドデータに対し、画像の特性に応じて必要な情報を的確に残しつつ、従来よりも高品質の圧縮処理を行うことができる。なお、ここでは、当該離散ウェーブレット変換されたタイルを構成するサブバンドデータの一部のサブバンドデータであって、当該タイルの特質が最も現れやすいサブバンドデータ（水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータ）についてその特性（文字画像、網点画像、又は写真画像のうちいずれであるか等）を解析して用いるので、当該タイル全体の特性を比較的短時間で解析可能である。
【００１１】
また、発明の請求項３に記載の発明は、請求項１又は請求項２に記載の画像処理装置であって、前記特性解析部は、前記サブバンドデータの周波数成分と、当該サブバンドデータを構成する各画素値の分布を示すヒストグラムとを用いて、当該サブバンドデータの特性を解析するものである。
【００１２】
この発明によれば、当該サブバンドデータの上記特性解析処理が、サブバンドデータの周波数成分と、当該サブバンドデータを構成する各画素値の分布を示すヒストグラムとを用いて行われるので、当該サブバンドデータの特性解析処理、例えば、当該サブバンドデータの特性が文字画像、網点画像、又は写真画像のうちいずれであるかを判定する処理を的確に行うことができる。
【００１３】
また、発明の請求項４に記載の発明は、請求項１乃至請求項３のいずれかに記載の画像処理装置であって、前記特性解析部は、前記サブバンドデータの前記特性解析に加えて、前記タイルを構成する他の予め定められたサブバンドデータについても前記特性解析を行い、前記量子化部は、当該特性解析により得られた各サブバンドデータの特性に応じた量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、当該タイルの各サブバンドデータを量子化するものである。
【００１４】
この発明では、上記特性解析が、上記低域通過フィルタリングが行われたサブバンドデータに加えて、タイルを構成する他の予め定められたサブバンドデータについても更に行われ、上記量子化処理は、当該特性解析により得られた各サブバンドデータの特性に応じた各サブバンドデータ毎の量子化テーブルを用いて、当該各サブバンドデータを量子化するので、当該サブバンドデータの特性解析処理、例えば、当該サブバンドデータの特性が文字画像、網点画像、又は写真画像のうちいずれであるかを判定する処理をより的確に行うことができる。
【発明の効果】
【００１５】
本発明によれば、ＪＰＥＧ２０００等の圧縮アルゴリズムにおいて、圧縮処理の対象とされる画像データをなす当該タイルの各サブバンドデータに対し、画像の特性に応じて必要な情報を的確に残しつつ、従来よりも高品質の圧縮処理を行うことができる。
【図面の簡単な説明】
【００１６】
【図１】本実施の形態における画像処理装置の電気的構成を示すブロック図である。
【図２】画像圧縮部の概略構成を示す図である。
【図３】タイル分割を示す図である。
【図４】離散ウェーブレット変換について示す模式図である。
【図５】画像圧縮部が画像圧縮プログラムに従って実行する画像圧縮処理の実施形態を示したフローチャートである。
【図６】（ａ）は文字画像及びそのヒストグラムの例を、（ｂ）は網点画像及びそのヒストグラムの例、（ｃ）は写真画像及びそのヒストグラムの例を示す。
【図７】（ａ）は文字画像に用いる量子化テーブルの例を示す図、（ｂ）は網点画像に用いる量子化テーブルの例を示す図、（ａ）は写真画像に用いる量子化テーブルの例を示す図である。
【発明を実施するための形態】
【００１７】
本発明における画像処理装置、画像処理方法、及び画像処理プログラムの実施の形態について図面を参照して説明する。図１は、本実施の形態における画像処理装置１の電気的構成を示すブロック図である。なお、画像処理装置１は、パーソナルコンピュータ（パソコン）やワークステーション、携帯情報端末等の情報処理装置で実現され、図１のブロック図は画像処理装置１がこれらの情報処理装置によって実現された場合を例に図示している。この他、画像処理装置１としての情報処理装置は、コピー機やスキャナ等の画像読取装置や、デジタルカメラやデジタルビデオカメラ、デジタルカメラを搭載した携帯電話等の撮像装置に組み込まれて用いられてもよい。
【００１８】
図１に示すように、画像処理装置１は、制御部１１、記憶部１２、入力操作部１３、表示部１４、Ｉ／Ｆ部１５、ネットワークＩ／Ｆ部１６、画像圧縮部１７及びＲＡＭ１８等を備えて構成される。制御部１１は、ＣＰＵ（Central Processing Unit：中央処理装置）等によって構成され、入力された指示信号等に応じて記憶部１２に記憶された画像処理プログラムを読み出して処理を実行し、各機能部への指示信号の出力、データ転送等を行って画像処理装置１を統括的に制御する。
【００１９】
記憶部１２は、画像処理装置１の備える種々の機能を実現するためのプログラムやデータ等を記憶する。本実施の形態では、記憶部１２は画像圧縮プログラム（特許請求の範囲でいう画像処理プログラムの一例）１２１を記憶する。画像圧縮プログラム１２１は、Ｉ／Ｆ部１５又はネットワークＩ／Ｆ部１６が入力した画像データをＪＰＥＧ２０００方式で圧縮して符号化データを生成するためのプログラムである。例えば、ＣＰＵ等でなる制御部１１が、画像圧縮プログラム１２１に従って動作することで画像圧縮部１７として機能する。この画像圧縮プログラム１２１は、ＣＤ−ＲＯＭ又はＤＶＤからの読込、或いは、Ｉ／Ｆ部１５又はネットワークＩ／Ｆ部１６によるインターネット上のサーバからダウンロードにより記憶部１２に記憶される。さらに、記憶部１２は、量子化処理に用いられる後述の量子化テーブルを記憶している。
【００２０】
入力操作部１３は、各種操作ボタンやマウス等のポインティングデバイスを備え、ユーザによって操作がなされると、操作信号を制御部１１へ出力する。表示部１４は、液晶ディスプレイ等の表示画面であり、入力操作部１３から入力された内容に応じた表示を行ったり、制御部１１による処理内容や処理結果を表示する。
【００２１】
Ｉ／Ｆ部１５はＩＥＥＥ１３９４やＵＳＢ等のインターフェイスであり、外部装置と直接データの送受信を行うことができる。ネットワークＩ／Ｆ部１６は、ＬＡＮボード等の通信モジュールから構成され、ネットワークＩ／Ｆ部１６と接続されたネットワーク（不図示）を介して外部装置と種々のデータの送受信を行う。
【００２２】
画像圧縮部１７は、記憶部１２に記憶された画像圧縮プログラム１２１を実行することによって、Ｉ／Ｆ部１５やネットワークＩ／Ｆ部１６を介して入力された画像データをＪＰＥＧ２０００方式で圧縮して符号化データを生成する。ＲＡＭ（Random Access Memory）１８は、制御部１１が実行するための制御プログラムや処理を実行する際のデータ等を一時的に保存する。また、画像圧縮部１７が画像圧縮プログラム１２１に従って処理を実行する際に生成される符号化データや圧縮率調整データ等を一時的に保存する。
【００２３】
図２は、画像圧縮部１７の概略構成を示す図である。図３はタイル分割を示す図である。図４は離散ウェーブレット変換について示す模式図である。
【００２４】
画像圧縮部１７は、表色系変換部１７１と、タイル分割部（分割部）１７２と、離散ウェーブレット変換部１７３と、特性解析部１７７と、量子化部１７４と、符号化部１７５として機能する。
【００２５】
まず、表色系変換部１７１が、入力された画像データの表色系を例えばＲＧＢ方式からＹＣｂＣｒ方式へコンポーネント変換する。
【００２６】
続いて、タイル分割部１７２が、表色系が変換された画像データを図３に示すように複数のタイルに分割する。ここで、タイル分割部１７２は、主走査方向に配列したタイルを１つのバンドとして区分する。離散ウェーブレット変換部１７３から符号化部１７５による処理は当該バンド毎に行う。図３を用いて説明すると、タイルＡ〜Ｅ、タイルＦ〜Ｋ、タイルＬ〜Ｑ、タイルＲ〜Ｖをそれぞれ１つのバンドとして区分し、バンド別に離散ウェーブレット変換部１７３以降の処理を行う。従って、以下で説明する処理は全てバンド単位で行われるものとして説明する。
【００２７】
図２に戻り、離散ウェーブレット変換部１７３は、１つのバンドの各タイルに対して離散ウェーブレット変換を施して係数データを出力する。図４は、離散ウェーブレット変換について示した模式図である。離散ウェーブレット変換部１７３は、図４（ａ）のタイル画像（０ＬＬ）に対して、離散ウェーブレット変換を施し、図４（ｂ）に示すようにサブバンド１ＬＬ、１ＨＬ、１ＬＨ及び１ＨＨに分解する。続いて、低周波成分のサブバンド１ＬＬに対して離散ウェーブレット変換を施し、図４（ｃ）に示すようにサブバンド２ＬＬ、２ＨＬ、２ＬＨ及び２ＨＨに分解する。尚、図４及び本実施形態では２レベル変換（２段階の変換）について図示しているが、離散ウェーブレット変換の回数は特に制限しない。
【００２８】
特性解析部１７７は、離散ウェーブレット変換部１７３によって離散ウェーブレット変換された上記タイルの各サブバンドのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンド２ＬＬの特性を、当該サブハンドデータの周波数成分に基づいて解析する。この特性とは、本実施形態では、当該ＪＰＥＧ２０００方式による圧縮処理の対象となっている画像データが、文字画像（文字の画像を含む画像）、網点画像（網点を示す画像群を含む画像）、又は写真画像（写真を示す画像を含む画像）のうちいずれであるかを示す。但し、当該特性をこの例に限定する趣旨ではない。本実施形態では、特性解析部１７７は、上記離散ウェーブレット変換されたサブバンド２ＬＬが、文字画像、網点画像、又は写真画像のうちいずれであるかを解析するものとして説明する。特性解析部１７７は、当該解析したサブバンド２ＬＬの特性を、当該サブバンド２ＬＬを含む上記タイル全体の特性とする。
【００２９】
量子化部１７４は、離散ウェーブレット変換された係数データに対して線形量子化を行う。量子化部１７４は、当該量子化処理を行うとき、(1)特性解析部１７７によって解析された上記タイルの特性に対応する量子化テーブルを記憶部１２から読み出し、(2)当該読み出した量子化テーブルを、当該ＪＰＥＧ２０００方式による圧縮処理の圧縮率として予め設定されている値に応じた係数を用いて算出した値に変更し、(3)当該変更した量子化テーブルを用いて、離散ウェーブレット変換部１７３による離散ウェーブレット変換済みの上記タイルの各サブバンドデータを量子化する。
【００３０】
なお、記憶部１２には、上記離散ウェーブレット変換されたタイルの各サブバンドデータの量子化に用いる量子化テーブルが、当該タイルの特性毎にそれぞれ記憶されている。
【００３１】
更に、符号化部１７５が、量子化部１７４により量子化された量子化データに対してエントロピー符号化を行う。このエントロピー符号化は、 EBCOT (Embedded Block Coding with Optimized Truncation) と呼ばれるアルゴリズムにより行われる。エントロピー符号化が終わった後、符号化部１７５は、ファイルにデータを書き込むための符号ストリームの形成を行う。この符号ストリーム処理については、各文献により周知のものであるため説明を省略する。
【００３２】
図５は、画像圧縮部１７が画像圧縮プログラム１２１に従って実行する画像圧縮処理の第１実施形態を示したフローチャートである。なお、以下の各実施形態について各フローチャートで示す画像処理は、特許請求の範囲における画像処理方法の一例である。
【００３３】
まず、表色系変換部１７１及びタイル分割部１７２は、画像データの色空間変換、タイル分割を行い、分割されたタイルをバンド別に区分する（Ｓ１）。続いて、離散ウェーブレット変換部１７３は、当該画像圧縮処理の対象とするバンドを選択する（Ｓ２）。さらに、離散ウェーブレット変換部１７３は、当該選択したバンドの中から離散ウェーブレット変換の対象とするタイルを選択する（Ｓ３）。
【００３４】
この後、離散ウェーブレット変換部１７３は、当該選択されたタイル（０ＬＬ）に対して、離散ウェーブレット変換を施してサブバンド１ＬＬ、１ＨＬ、１ＬＨ及び１ＨＨに分解し、更に、低周波成分のサブバンド１ＬＬに対して離散ウェーブレット変換を施し、サブバンド２ＬＬ、２ＨＬ、２ＬＨ及び２ＨＨに分解する（Ｓ４）。
【００３５】
続いて、特性解析部１７７は、サブバンド２ＬＬの特性（パターン）分析（上記文字画像、網点画像、又は写真画像のいずれであるかの判定）を開始する。まず、特性解析部１７７は、サブバンド２ＬＬを構成する各画素について画素値の分布を示すヒストグラムを作成する（Ｓ５）。このヒストグラムは、横軸に画素値（明度）、縦軸に画素値（明るさ）毎の画素数を積み上げて示したものである。さらに、特性解析部１７７は、サブバンド２ＬＬのデータについてDCT（離散コサイン変換）による周波数の測定を行う（Ｓ６）。
【００３６】
ここで、文字画像、網点画像、及び写真画像は、例えば、以下のような傾向を示す。文字画像及びそのヒストグラムの例を図６（ａ）、網点画像及びそのヒストグラムの例を図６（ｂ）、写真画像及びそのヒストグラムの例を図６（ｃ）に示す。
〈文字画像〉
ヒストグラム：地肌（文字以外を示す背景の画像）の部分を示す大きな山と、文字を構成する部分を示す小さな山とが存在
周波数：高周波成分がランダムに発生
〈網点画像〉
ヒストグラム：地肌部分と網点で構成された濃度部分が一様に存在
周波数：特定の高周波成分が均一に発生
〈写真画像〉
ヒストグラム：地肌部分は少なく、特定の濃度領域に広く濃度が分布
周波数：低周波成分が多く、高周波成分がほとんど発生しない（文字原稿より高周波成が少ない）
【００３７】
この傾向に基づいて、特性解析部１７７は、サブバンド２ＬＬの上記ヒストグラム及び周波数から、当該サブバンド２ＬＬのデータが、文字画像、網点画像、又は写真画像のいずれであるかを判定する（Ｓ７）。すなわち、特性解析部１７７は、上記文字画像、網点画像、及び写真画像のそれぞれのヒストグラム及び周波数が示す傾向を定義して記憶しておき、Ｓ５及びＳ６で得られたヒストグラム及び周波数を解析し、当該ヒストグラム及び周波数の両方が、上記各画像のいずれかの傾向に当てはまると判断した場合に、当てはまった画像種を当該サブバンド２ＬＬのデータの特性とし、当該特性を、当該サブバンド２ＬＬを含むタイル全体の特性とする。
【００３８】
そして、量子化部１７４は、特性解析部１７７によって解析された上記特性に対応する量子化テーブルを記憶部１２から読み出す（Ｓ８）。すなわち、量子化部１７４は、特性解析部１７７によって解析された上記特性に対応する量子化テーブルを決定する。例えば、特性が文字画像の場合には、当該サブバンド２ＬＬを含むタイル全体の成分を均等に量子化してエッジ成分を残すため、図７（ａ）に例を示すような量子化テーブルが用意されている。網点画像の場合には、当該サブバンド２ＬＬ〜２ＨＨ成分を適度に残し、１ＨＬ〜１ＨＨ成分を削除するため、図７（ｂ）に例を示すような量子化テーブルが用意されている。写真画像の場合には、当該サブバンド２ＬＬ成分を優先し、１ＨＨ成分を大きく削除するため、図７（ｃ）に例を示すような量子化テーブルが用意されている。なお、サブバンド２ＬＬのデータが、文字画像、網点画像、又は写真画像のいずれでもない（サブバンド２ＬＬのヒストグラム及び周波数の両方が上記各画像のいずれの傾向にも当てはまらない）場合、サブバンド２ＬＬのデータは、文字と網点、文字と写真が混在している複合画像であると判断する。この複合画像には、量子化部１７４は、量子化テーブルの変更は行わず、初期設定が用いられる。図７（ｄ）に、初期設定の量子化テーブルの例を示す。
【００３９】
続いて、量子化部１７４は、当該読み出した量子化テーブルを、当該画像圧縮処理で予め定めている圧縮率に応じた係数を用いて算出した値に変更し、当該算出した量子化テーブルを量子化に用いる量子化テーブルとする（Ｓ９）。
【００４０】
例えば、量子化部１７４は、
圧縮率が５０％以上の場合は、求められた量子化テーブルをそのまま使用、
４０％以上５０％未満の場合は、量子化テーブルを１．１倍
３０％以上４０％未満の場合は、量子化テーブルを１．２倍
２０％以上３０％未満の場合は、量子化テーブルを１．３倍
１０％以上２０％未満の場合は、量子化テーブルを１．４倍
１０％未満の場合は、量子化テーブルを１．５倍、とする。
【００４１】
なお、ここでは、当該画像圧縮処理で予め定めている圧縮率に応じて、量子化処理に用いる量子化テーブルを算出により求めているが、当該画像圧縮処理で予め定めている圧縮率が高い（符号量を小さくする）場合は、記憶部１２が、当該高い圧縮率として細かく設定された複数の圧縮率のそれぞれに対応する量子化テーブルを記憶しておき、量子化部１７４が、当該各圧縮率に対応する量子化テーブルを記憶部１２から選択して読み出すことで、量子化処理に用いる量子化テーブルを求めてもよい。これにより、当該画像圧縮処理で予め定めている圧縮率が高い場合に、当該圧縮率に応じた量子化テーブルを詳細に設定することが可能である。
【００４２】
そして、量子化部１７４は、当該算出した量子化テーブルを用いて上記タイル全体の全サブバンド（２ＬＬ〜１ＨＨ）を量子化する（Ｓ１０）。この後、符号化部１７５が符号化を行う（Ｓ１１）。さらに、符号化部１７５により、公知の圧縮率制御が行われる（Ｓ１２）。
【００４３】
なお、当該ＪＰＥＧ２０００の圧縮率制御は、例えば、符号量と、画質への影響度の２つのパラメータとを用いて行われる。ここで、符号量は、上記符号化処理で符号化処理された結果そのものが用いられる。画質への影響度は、各画素位置における画素値によって求められ、符号化処理時に量子化後のデータを用いて求められる。この影響度は、上位bitほど影響度が高いと判断されるように係数が設定されており、該当bit位置ごとに求められる。当該影響度は、具体的には、サブバンド1LLの最上位bit (10bit目) から2HHの最下位bit(1bit目）まで求められる。尚、上記符号化処理は、当該影響度の計算と同じ処理単位で行われる。当該影響度の処理単位は、サブバンド数×bit数で分割されたものとなるため、入力10bitで1LL〜2HHに分割する場合、7×10bitの70個の符号データと当該符号データについての影響度とがタイルごとに算出される。そして、各タイル毎に、この求められた影響度の高い符号データから順番に出力を行い、指定された符号量に達した時点で符号データの出力を停止させる制御を行うことが、当該ＪＰＥＧ２０００の圧縮率制御の一例となる。
【００４４】
そして、離散ウェーブレット変換部１７３は、この時点で画像処理の対象としているバンドについて未処理のタイルが残っているかを判断し（Ｓ１３）、未処理のタイルが残っている場合には（Ｓ１３でＹＥＳ）、Ｓ３乃至Ｓ１２までの処理を次のタイルについて繰り返す。
【００４５】
一方、未処理のタイルが残っていない場合には（Ｓ１３でＮＯ）、上記分割したバンドのうち、離散ウェーブレット変換未処理のバンドが残っているか否かを判断し（Ｓ１４）、残っている場合には（Ｓ１４でＹＥＳ）、処理はＳ２に戻り、残りのバンドについてＳ２乃至Ｓ１３の処理を行う。
【００４６】
この後、離散ウェーブレット変換処理を全てのバンドに対して終えると（Ｓ１４でＮＯ）、符号化部１７５は、上記符号処理及び圧縮率制御後のデータを用いてパケットデータを作成し（Ｓ１５）、当該符号化処理及び圧縮率制御後の符号データを出力する（Ｓ１６）。
【００４７】
なお、上記一連の処理においては、上記タイルを構成するサブハンドのうち、サブバンド２ＬＬのみについてＳ５乃至Ｓ９の処理を行っているが、当該Ｓ５乃至Ｓ９の処理を変形して、上記タイルを構成するサブバンドのうちの特定のサブバンドまで画像の特性を判定するものとしてもよい。例えば、サブバンド２ＬＬについてＳ５乃至Ｓ９と同様の処理を行って上記画像の特性を判断した後に、特性解析部１７７が、サブバンド２ＨＬ、２ＬＨ、２ＨＨについてそれぞれヒストグラムを作成し、高周波成分であるサブバンド２ＨＬ、２ＬＨ、２ＨＨの画像の特性を、当該ヒストグラムを用いて分析する。入力画像からサブバンド２ＬＬで求められた低周波成分を除いた画像の成分がサブバンド２ＨＬ、２ＬＨ、２ＨＨであるため、サブバンド２ＨＬ、２ＬＨ、２ＨＨの部分は、サブバンド２ＬＬとの差と考えることができ、サブバンド２ＨＬ、２ＬＨ、２ＨＨの画像特性を示すヒストグラムとなる。このヒストグラムにおいて、高濃度や低濃度の部分に存在する画素の数により、特性解析部１７７は、サブバンド２ＨＬ、２ＬＨ、２ＨＨの画像特性が、文字画像、網点画像、写真画像のいずれであるかを判定する。
【００４８】
これらサブバンド２ＨＬ、２ＬＨ、２ＨＨのヒストグラムでは、(a)文字画像であれば、高濃度と低濃度に文字エッジの量に従った画素数の分布が存在し、(b)網点画像であれば、高濃度と低濃度に網点の線数に対応した一定数以上の画素数が存在し（網点が画像一面に存在するために多くの高周波成分が存在する）、(c)写真画像であれば、高濃度と低濃度の画素数は少ない（エッジとなる部分がほとんど存在しない）。
【００４９】
このため、高周波成分であるサブバンド２ＨＬ、２ＬＨ、２ＨＨにおける高濃度及び低濃度の画素数は、網点画像＞文字画像＞写真画像の関係となる。特性解析部１７７は、各画像種における高濃度及び低濃度の画素数の各境界部分を閾値として設定し、サブバンド２ＨＬ、２ＬＨ、２ＨＨのそれぞれのヒストグラムが示す画素分布が、網点画像、文字画像、又は写真画像のいずれに属するかを判定することで、サブバンド２ＨＬ、２ＬＨ、２ＨＨの画像特性を判定する。
【００５０】
このように判定された各サブバンド２ＨＬ、２ＬＨ、２ＨＨと、既に判定されているサブバンド２ＬＬの各画像特性に基づいて、例えば、サブバンド２ＬＬ、２ＨＬ、２ＬＨ、２ＨＨのうちで最も多くのサブバンドが示す画像特性を、これらサブバンドを含むタイル全体の画像特性とする（画像特性が同数の場合はサブバンド２ＬＬの画像特性とする）。量子化部１７４は、当該特定した画像特性に対応する量子化テーブルを記憶部１２から読み出して、当該読み出した量子化テーブルを上記のように圧縮率を変更し、当該圧縮率変更後の量子化テーブルを、上記タイル全体の全サブバンド（２ＬＬ〜１ＨＨ）の量子化に用いる量子化テーブルとする。
【００５１】
また、さらには、各サブバンド１ＬＬ〜２ＨＨのそれぞれについて、上記Ｓ５〜Ｓ７の処理と同様にして画像特性を求めると共に、各画像特性のそれぞれに対応する量子化テーブルを当該各サブバンド毎に記憶部１２が記憶しておき、量子化部１７４が、これら各サブバンド１ＬＬ〜２ＨＨの各画像特性に応じた量子化テーブルを各サブバンド毎に読み出し、当該読み出した量子化テーブルを用いて、各サブバンド別に量子化を行うようにしてもよい。
【００５２】
なお、上記図１乃至図７に示した実施形態に係る構成及び処理は、本発明に係る画像処理装置、画像処理方法、及び画像処理プログラムの構成及び処理の単なる一例に過ぎず、本発明に係る画像処理装置、画像処理方法、及び画像処理プログラムの構成及び処理を上記に示した内容に限定するものではない。
【符号の説明】
【００５３】
１ 画像処理装置
１２ 記憶部
１２１ 画像圧縮プログラム
１７ 画像圧縮部
１７１ 表色系変換部
１７２ タイル分割部
１７３ 離散ウェーブレット変換部
１７４ 量子化部
１７５ 符号化部
１７７ 特性解析部
【技術分野】
【０００１】
本発明は、画像処理装置、画像処理方法、及び画像処理プログラムに関し、特に、画像データを圧縮する技術に関する。
【背景技術】
【０００２】
近年、デジタル機器の発展に伴って、画像データ等の様々な圧縮方式が提案されているが、圧縮方式の１つとしてＪＰＥＧ（Joint Photographic Experts Group）２０００が知られている。ＪＰＥＧ２０００は画像の圧縮・展開の方法を規定したもので、ＪＰＥＧを発展させたものである。ＪＰＥＧ２０００では、画像データを複数のタイルに分割し、タイル毎に離散ウェーブレット変換、量子化、エントロピー符号化等の処理を施して符号化データを生成する。なお、特許文献１及び２に記載されているように、ＪＰＥＧ圧縮では、圧縮処理の対象とされる画像全体に基づいて量子化処理に用いる量子化テーブルを選択し、当該選択した量子化テーブルを用いて量子化を行う方式が提案されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００２−２４７３７３号公報
【特許文献２】特開２００４−１１０３６８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上記従来一般のＪＰＥＧ２０００の画像圧縮方式では、量子化処理に用いる量子化テーブルは、圧縮率の指定やタイルサイズ等に基づいて１つに定められている。また、圧縮率の制御は、符号化後で各ビットに対する重み付で優先順位を決めて制御するようになっているため、入力される画像の特徴に応じた最適な符号化は行われていない。また、上記特許文献１に記載の発明は、ＪＰＥＧ圧縮に関するものであり、ＪＰＥＧ２０００の画像圧縮方式に関するものではなく、さらに、入力された画像全体に基づいて量子化テーブルを選択するため、画像に応じた細かな圧縮処理ができない。
【０００５】
本発明は、上記の問題を解決するためになされたもので、ＪＰＥＧ２０００等の圧縮アルゴリズムにおいて、圧縮処理の対象とされる画像の特徴に応じて必要な情報を的確に残しつつ、従来よりも高品質の圧縮処理を可能にすることを目的とする。
【課題を解決するための手段】
【０００６】
本発明の請求項１に記載の発明は、画像データを輝度成分と色差成分からなる表色系に変換する表色系変換部と、
前記表色系変換部によって変換された画像データを複数のタイルに分割する分割部と、
前記分割部によって分割されたタイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶する量子化テーブル記憶部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析部と、
前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化部と、
前記量子化部によって量子化された画像データを符号化して圧縮データを出力する符号化部と
を備えた画像処理装置である。
【０００７】
また、請求項２に記載の発明は、請求項１に記載の画像処理装置であって、前記量子化部は、前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを、予め設定されている圧縮率に応じた係数を用いて算出した値に変更し、当該変更した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化するものである。
【０００８】
また、請求項５に記載の発明は、画像データを輝度成分と色差成分からなる表色系に変換する表色系変換ステップと、
前記表色系変換ステップで変換された画像データを複数のタイルに分割する分割ステップと、
前記分割ステップで分割された前記タイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換ステップと、
前記離散ウェーブレット変換ステップで離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析ステップと、
前記離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶しており、前記特性解析ステップで解析された特性に対応する量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化ステップと、
前記量子化ステップによって量子化された画像データを符号化して圧縮データを出力する符号化ステップと
を有する画像処理方法である。
【０００９】
また、請求項６に記載の発明は、コンピュータを、
画像データを輝度成分と色差成分からなる表色系に変換する表色系変換部と、
前記表色系変換部によって変換された画像データを複数のタイルに分割する分割部と、
前記分割部によって分割されたタイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶する量子化テーブル記憶部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析部と、
前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化部と、
前記量子化部によって量子化された画像データを符号化して圧縮データを出力する符号化部と
して機能させる画像処理プログラムである。
【００１０】
これらの発明では、離散ウェーブレット変換された上記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドデータの特性を当該タイル全体の特性とし、量子化処理時には、当該解析されたタイルの特性に対応する量子化テーブルを用いて当該タイルの各サブバンドデータを量子化する。これにより、ＪＰＥＧ２０００等の圧縮アルゴリズムにおいて、圧縮処理の対象とされる画像データをなす当該タイルの各サブバンドデータに対し、画像の特性に応じて必要な情報を的確に残しつつ、従来よりも高品質の圧縮処理を行うことができる。なお、ここでは、当該離散ウェーブレット変換されたタイルを構成するサブバンドデータの一部のサブバンドデータであって、当該タイルの特質が最も現れやすいサブバンドデータ（水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータ）についてその特性（文字画像、網点画像、又は写真画像のうちいずれであるか等）を解析して用いるので、当該タイル全体の特性を比較的短時間で解析可能である。
【００１１】
また、発明の請求項３に記載の発明は、請求項１又は請求項２に記載の画像処理装置であって、前記特性解析部は、前記サブバンドデータの周波数成分と、当該サブバンドデータを構成する各画素値の分布を示すヒストグラムとを用いて、当該サブバンドデータの特性を解析するものである。
【００１２】
この発明によれば、当該サブバンドデータの上記特性解析処理が、サブバンドデータの周波数成分と、当該サブバンドデータを構成する各画素値の分布を示すヒストグラムとを用いて行われるので、当該サブバンドデータの特性解析処理、例えば、当該サブバンドデータの特性が文字画像、網点画像、又は写真画像のうちいずれであるかを判定する処理を的確に行うことができる。
【００１３】
また、発明の請求項４に記載の発明は、請求項１乃至請求項３のいずれかに記載の画像処理装置であって、前記特性解析部は、前記サブバンドデータの前記特性解析に加えて、前記タイルを構成する他の予め定められたサブバンドデータについても前記特性解析を行い、前記量子化部は、当該特性解析により得られた各サブバンドデータの特性に応じた量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、当該タイルの各サブバンドデータを量子化するものである。
【００１４】
この発明では、上記特性解析が、上記低域通過フィルタリングが行われたサブバンドデータに加えて、タイルを構成する他の予め定められたサブバンドデータについても更に行われ、上記量子化処理は、当該特性解析により得られた各サブバンドデータの特性に応じた各サブバンドデータ毎の量子化テーブルを用いて、当該各サブバンドデータを量子化するので、当該サブバンドデータの特性解析処理、例えば、当該サブバンドデータの特性が文字画像、網点画像、又は写真画像のうちいずれであるかを判定する処理をより的確に行うことができる。
【発明の効果】
【００１５】
本発明によれば、ＪＰＥＧ２０００等の圧縮アルゴリズムにおいて、圧縮処理の対象とされる画像データをなす当該タイルの各サブバンドデータに対し、画像の特性に応じて必要な情報を的確に残しつつ、従来よりも高品質の圧縮処理を行うことができる。
【図面の簡単な説明】
【００１６】
【図１】本実施の形態における画像処理装置の電気的構成を示すブロック図である。
【図２】画像圧縮部の概略構成を示す図である。
【図３】タイル分割を示す図である。
【図４】離散ウェーブレット変換について示す模式図である。
【図５】画像圧縮部が画像圧縮プログラムに従って実行する画像圧縮処理の実施形態を示したフローチャートである。
【図６】（ａ）は文字画像及びそのヒストグラムの例を、（ｂ）は網点画像及びそのヒストグラムの例、（ｃ）は写真画像及びそのヒストグラムの例を示す。
【図７】（ａ）は文字画像に用いる量子化テーブルの例を示す図、（ｂ）は網点画像に用いる量子化テーブルの例を示す図、（ａ）は写真画像に用いる量子化テーブルの例を示す図である。
【発明を実施するための形態】
【００１７】
本発明における画像処理装置、画像処理方法、及び画像処理プログラムの実施の形態について図面を参照して説明する。図１は、本実施の形態における画像処理装置１の電気的構成を示すブロック図である。なお、画像処理装置１は、パーソナルコンピュータ（パソコン）やワークステーション、携帯情報端末等の情報処理装置で実現され、図１のブロック図は画像処理装置１がこれらの情報処理装置によって実現された場合を例に図示している。この他、画像処理装置１としての情報処理装置は、コピー機やスキャナ等の画像読取装置や、デジタルカメラやデジタルビデオカメラ、デジタルカメラを搭載した携帯電話等の撮像装置に組み込まれて用いられてもよい。
【００１８】
図１に示すように、画像処理装置１は、制御部１１、記憶部１２、入力操作部１３、表示部１４、Ｉ／Ｆ部１５、ネットワークＩ／Ｆ部１６、画像圧縮部１７及びＲＡＭ１８等を備えて構成される。制御部１１は、ＣＰＵ（Central Processing Unit：中央処理装置）等によって構成され、入力された指示信号等に応じて記憶部１２に記憶された画像処理プログラムを読み出して処理を実行し、各機能部への指示信号の出力、データ転送等を行って画像処理装置１を統括的に制御する。
【００１９】
記憶部１２は、画像処理装置１の備える種々の機能を実現するためのプログラムやデータ等を記憶する。本実施の形態では、記憶部１２は画像圧縮プログラム（特許請求の範囲でいう画像処理プログラムの一例）１２１を記憶する。画像圧縮プログラム１２１は、Ｉ／Ｆ部１５又はネットワークＩ／Ｆ部１６が入力した画像データをＪＰＥＧ２０００方式で圧縮して符号化データを生成するためのプログラムである。例えば、ＣＰＵ等でなる制御部１１が、画像圧縮プログラム１２１に従って動作することで画像圧縮部１７として機能する。この画像圧縮プログラム１２１は、ＣＤ−ＲＯＭ又はＤＶＤからの読込、或いは、Ｉ／Ｆ部１５又はネットワークＩ／Ｆ部１６によるインターネット上のサーバからダウンロードにより記憶部１２に記憶される。さらに、記憶部１２は、量子化処理に用いられる後述の量子化テーブルを記憶している。
【００２０】
入力操作部１３は、各種操作ボタンやマウス等のポインティングデバイスを備え、ユーザによって操作がなされると、操作信号を制御部１１へ出力する。表示部１４は、液晶ディスプレイ等の表示画面であり、入力操作部１３から入力された内容に応じた表示を行ったり、制御部１１による処理内容や処理結果を表示する。
【００２１】
Ｉ／Ｆ部１５はＩＥＥＥ１３９４やＵＳＢ等のインターフェイスであり、外部装置と直接データの送受信を行うことができる。ネットワークＩ／Ｆ部１６は、ＬＡＮボード等の通信モジュールから構成され、ネットワークＩ／Ｆ部１６と接続されたネットワーク（不図示）を介して外部装置と種々のデータの送受信を行う。
【００２２】
画像圧縮部１７は、記憶部１２に記憶された画像圧縮プログラム１２１を実行することによって、Ｉ／Ｆ部１５やネットワークＩ／Ｆ部１６を介して入力された画像データをＪＰＥＧ２０００方式で圧縮して符号化データを生成する。ＲＡＭ（Random Access Memory）１８は、制御部１１が実行するための制御プログラムや処理を実行する際のデータ等を一時的に保存する。また、画像圧縮部１７が画像圧縮プログラム１２１に従って処理を実行する際に生成される符号化データや圧縮率調整データ等を一時的に保存する。
【００２３】
図２は、画像圧縮部１７の概略構成を示す図である。図３はタイル分割を示す図である。図４は離散ウェーブレット変換について示す模式図である。
【００２４】
画像圧縮部１７は、表色系変換部１７１と、タイル分割部（分割部）１７２と、離散ウェーブレット変換部１７３と、特性解析部１７７と、量子化部１７４と、符号化部１７５として機能する。
【００２５】
まず、表色系変換部１７１が、入力された画像データの表色系を例えばＲＧＢ方式からＹＣｂＣｒ方式へコンポーネント変換する。
【００２６】
続いて、タイル分割部１７２が、表色系が変換された画像データを図３に示すように複数のタイルに分割する。ここで、タイル分割部１７２は、主走査方向に配列したタイルを１つのバンドとして区分する。離散ウェーブレット変換部１７３から符号化部１７５による処理は当該バンド毎に行う。図３を用いて説明すると、タイルＡ〜Ｅ、タイルＦ〜Ｋ、タイルＬ〜Ｑ、タイルＲ〜Ｖをそれぞれ１つのバンドとして区分し、バンド別に離散ウェーブレット変換部１７３以降の処理を行う。従って、以下で説明する処理は全てバンド単位で行われるものとして説明する。
【００２７】
図２に戻り、離散ウェーブレット変換部１７３は、１つのバンドの各タイルに対して離散ウェーブレット変換を施して係数データを出力する。図４は、離散ウェーブレット変換について示した模式図である。離散ウェーブレット変換部１７３は、図４（ａ）のタイル画像（０ＬＬ）に対して、離散ウェーブレット変換を施し、図４（ｂ）に示すようにサブバンド１ＬＬ、１ＨＬ、１ＬＨ及び１ＨＨに分解する。続いて、低周波成分のサブバンド１ＬＬに対して離散ウェーブレット変換を施し、図４（ｃ）に示すようにサブバンド２ＬＬ、２ＨＬ、２ＬＨ及び２ＨＨに分解する。尚、図４及び本実施形態では２レベル変換（２段階の変換）について図示しているが、離散ウェーブレット変換の回数は特に制限しない。
【００２８】
特性解析部１７７は、離散ウェーブレット変換部１７３によって離散ウェーブレット変換された上記タイルの各サブバンドのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンド２ＬＬの特性を、当該サブハンドデータの周波数成分に基づいて解析する。この特性とは、本実施形態では、当該ＪＰＥＧ２０００方式による圧縮処理の対象となっている画像データが、文字画像（文字の画像を含む画像）、網点画像（網点を示す画像群を含む画像）、又は写真画像（写真を示す画像を含む画像）のうちいずれであるかを示す。但し、当該特性をこの例に限定する趣旨ではない。本実施形態では、特性解析部１７７は、上記離散ウェーブレット変換されたサブバンド２ＬＬが、文字画像、網点画像、又は写真画像のうちいずれであるかを解析するものとして説明する。特性解析部１７７は、当該解析したサブバンド２ＬＬの特性を、当該サブバンド２ＬＬを含む上記タイル全体の特性とする。
【００２９】
量子化部１７４は、離散ウェーブレット変換された係数データに対して線形量子化を行う。量子化部１７４は、当該量子化処理を行うとき、(1)特性解析部１７７によって解析された上記タイルの特性に対応する量子化テーブルを記憶部１２から読み出し、(2)当該読み出した量子化テーブルを、当該ＪＰＥＧ２０００方式による圧縮処理の圧縮率として予め設定されている値に応じた係数を用いて算出した値に変更し、(3)当該変更した量子化テーブルを用いて、離散ウェーブレット変換部１７３による離散ウェーブレット変換済みの上記タイルの各サブバンドデータを量子化する。
【００３０】
なお、記憶部１２には、上記離散ウェーブレット変換されたタイルの各サブバンドデータの量子化に用いる量子化テーブルが、当該タイルの特性毎にそれぞれ記憶されている。
【００３１】
更に、符号化部１７５が、量子化部１７４により量子化された量子化データに対してエントロピー符号化を行う。このエントロピー符号化は、 EBCOT (Embedded Block Coding with Optimized Truncation) と呼ばれるアルゴリズムにより行われる。エントロピー符号化が終わった後、符号化部１７５は、ファイルにデータを書き込むための符号ストリームの形成を行う。この符号ストリーム処理については、各文献により周知のものであるため説明を省略する。
【００３２】
図５は、画像圧縮部１７が画像圧縮プログラム１２１に従って実行する画像圧縮処理の第１実施形態を示したフローチャートである。なお、以下の各実施形態について各フローチャートで示す画像処理は、特許請求の範囲における画像処理方法の一例である。
【００３３】
まず、表色系変換部１７１及びタイル分割部１７２は、画像データの色空間変換、タイル分割を行い、分割されたタイルをバンド別に区分する（Ｓ１）。続いて、離散ウェーブレット変換部１７３は、当該画像圧縮処理の対象とするバンドを選択する（Ｓ２）。さらに、離散ウェーブレット変換部１７３は、当該選択したバンドの中から離散ウェーブレット変換の対象とするタイルを選択する（Ｓ３）。
【００３４】
この後、離散ウェーブレット変換部１７３は、当該選択されたタイル（０ＬＬ）に対して、離散ウェーブレット変換を施してサブバンド１ＬＬ、１ＨＬ、１ＬＨ及び１ＨＨに分解し、更に、低周波成分のサブバンド１ＬＬに対して離散ウェーブレット変換を施し、サブバンド２ＬＬ、２ＨＬ、２ＬＨ及び２ＨＨに分解する（Ｓ４）。
【００３５】
続いて、特性解析部１７７は、サブバンド２ＬＬの特性（パターン）分析（上記文字画像、網点画像、又は写真画像のいずれであるかの判定）を開始する。まず、特性解析部１７７は、サブバンド２ＬＬを構成する各画素について画素値の分布を示すヒストグラムを作成する（Ｓ５）。このヒストグラムは、横軸に画素値（明度）、縦軸に画素値（明るさ）毎の画素数を積み上げて示したものである。さらに、特性解析部１７７は、サブバンド２ＬＬのデータについてDCT（離散コサイン変換）による周波数の測定を行う（Ｓ６）。
【００３６】
ここで、文字画像、網点画像、及び写真画像は、例えば、以下のような傾向を示す。文字画像及びそのヒストグラムの例を図６（ａ）、網点画像及びそのヒストグラムの例を図６（ｂ）、写真画像及びそのヒストグラムの例を図６（ｃ）に示す。
〈文字画像〉
ヒストグラム：地肌（文字以外を示す背景の画像）の部分を示す大きな山と、文字を構成する部分を示す小さな山とが存在
周波数：高周波成分がランダムに発生
〈網点画像〉
ヒストグラム：地肌部分と網点で構成された濃度部分が一様に存在
周波数：特定の高周波成分が均一に発生
〈写真画像〉
ヒストグラム：地肌部分は少なく、特定の濃度領域に広く濃度が分布
周波数：低周波成分が多く、高周波成分がほとんど発生しない（文字原稿より高周波成が少ない）
【００３７】
この傾向に基づいて、特性解析部１７７は、サブバンド２ＬＬの上記ヒストグラム及び周波数から、当該サブバンド２ＬＬのデータが、文字画像、網点画像、又は写真画像のいずれであるかを判定する（Ｓ７）。すなわち、特性解析部１７７は、上記文字画像、網点画像、及び写真画像のそれぞれのヒストグラム及び周波数が示す傾向を定義して記憶しておき、Ｓ５及びＳ６で得られたヒストグラム及び周波数を解析し、当該ヒストグラム及び周波数の両方が、上記各画像のいずれかの傾向に当てはまると判断した場合に、当てはまった画像種を当該サブバンド２ＬＬのデータの特性とし、当該特性を、当該サブバンド２ＬＬを含むタイル全体の特性とする。
【００３８】
そして、量子化部１７４は、特性解析部１７７によって解析された上記特性に対応する量子化テーブルを記憶部１２から読み出す（Ｓ８）。すなわち、量子化部１７４は、特性解析部１７７によって解析された上記特性に対応する量子化テーブルを決定する。例えば、特性が文字画像の場合には、当該サブバンド２ＬＬを含むタイル全体の成分を均等に量子化してエッジ成分を残すため、図７（ａ）に例を示すような量子化テーブルが用意されている。網点画像の場合には、当該サブバンド２ＬＬ〜２ＨＨ成分を適度に残し、１ＨＬ〜１ＨＨ成分を削除するため、図７（ｂ）に例を示すような量子化テーブルが用意されている。写真画像の場合には、当該サブバンド２ＬＬ成分を優先し、１ＨＨ成分を大きく削除するため、図７（ｃ）に例を示すような量子化テーブルが用意されている。なお、サブバンド２ＬＬのデータが、文字画像、網点画像、又は写真画像のいずれでもない（サブバンド２ＬＬのヒストグラム及び周波数の両方が上記各画像のいずれの傾向にも当てはまらない）場合、サブバンド２ＬＬのデータは、文字と網点、文字と写真が混在している複合画像であると判断する。この複合画像には、量子化部１７４は、量子化テーブルの変更は行わず、初期設定が用いられる。図７（ｄ）に、初期設定の量子化テーブルの例を示す。
【００３９】
続いて、量子化部１７４は、当該読み出した量子化テーブルを、当該画像圧縮処理で予め定めている圧縮率に応じた係数を用いて算出した値に変更し、当該算出した量子化テーブルを量子化に用いる量子化テーブルとする（Ｓ９）。
【００４０】
例えば、量子化部１７４は、
圧縮率が５０％以上の場合は、求められた量子化テーブルをそのまま使用、
４０％以上５０％未満の場合は、量子化テーブルを１．１倍
３０％以上４０％未満の場合は、量子化テーブルを１．２倍
２０％以上３０％未満の場合は、量子化テーブルを１．３倍
１０％以上２０％未満の場合は、量子化テーブルを１．４倍
１０％未満の場合は、量子化テーブルを１．５倍、とする。
【００４１】
なお、ここでは、当該画像圧縮処理で予め定めている圧縮率に応じて、量子化処理に用いる量子化テーブルを算出により求めているが、当該画像圧縮処理で予め定めている圧縮率が高い（符号量を小さくする）場合は、記憶部１２が、当該高い圧縮率として細かく設定された複数の圧縮率のそれぞれに対応する量子化テーブルを記憶しておき、量子化部１７４が、当該各圧縮率に対応する量子化テーブルを記憶部１２から選択して読み出すことで、量子化処理に用いる量子化テーブルを求めてもよい。これにより、当該画像圧縮処理で予め定めている圧縮率が高い場合に、当該圧縮率に応じた量子化テーブルを詳細に設定することが可能である。
【００４２】
そして、量子化部１７４は、当該算出した量子化テーブルを用いて上記タイル全体の全サブバンド（２ＬＬ〜１ＨＨ）を量子化する（Ｓ１０）。この後、符号化部１７５が符号化を行う（Ｓ１１）。さらに、符号化部１７５により、公知の圧縮率制御が行われる（Ｓ１２）。
【００４３】
なお、当該ＪＰＥＧ２０００の圧縮率制御は、例えば、符号量と、画質への影響度の２つのパラメータとを用いて行われる。ここで、符号量は、上記符号化処理で符号化処理された結果そのものが用いられる。画質への影響度は、各画素位置における画素値によって求められ、符号化処理時に量子化後のデータを用いて求められる。この影響度は、上位bitほど影響度が高いと判断されるように係数が設定されており、該当bit位置ごとに求められる。当該影響度は、具体的には、サブバンド1LLの最上位bit (10bit目) から2HHの最下位bit(1bit目）まで求められる。尚、上記符号化処理は、当該影響度の計算と同じ処理単位で行われる。当該影響度の処理単位は、サブバンド数×bit数で分割されたものとなるため、入力10bitで1LL〜2HHに分割する場合、7×10bitの70個の符号データと当該符号データについての影響度とがタイルごとに算出される。そして、各タイル毎に、この求められた影響度の高い符号データから順番に出力を行い、指定された符号量に達した時点で符号データの出力を停止させる制御を行うことが、当該ＪＰＥＧ２０００の圧縮率制御の一例となる。
【００４４】
そして、離散ウェーブレット変換部１７３は、この時点で画像処理の対象としているバンドについて未処理のタイルが残っているかを判断し（Ｓ１３）、未処理のタイルが残っている場合には（Ｓ１３でＹＥＳ）、Ｓ３乃至Ｓ１２までの処理を次のタイルについて繰り返す。
【００４５】
一方、未処理のタイルが残っていない場合には（Ｓ１３でＮＯ）、上記分割したバンドのうち、離散ウェーブレット変換未処理のバンドが残っているか否かを判断し（Ｓ１４）、残っている場合には（Ｓ１４でＹＥＳ）、処理はＳ２に戻り、残りのバンドについてＳ２乃至Ｓ１３の処理を行う。
【００４６】
この後、離散ウェーブレット変換処理を全てのバンドに対して終えると（Ｓ１４でＮＯ）、符号化部１７５は、上記符号処理及び圧縮率制御後のデータを用いてパケットデータを作成し（Ｓ１５）、当該符号化処理及び圧縮率制御後の符号データを出力する（Ｓ１６）。
【００４７】
なお、上記一連の処理においては、上記タイルを構成するサブハンドのうち、サブバンド２ＬＬのみについてＳ５乃至Ｓ９の処理を行っているが、当該Ｓ５乃至Ｓ９の処理を変形して、上記タイルを構成するサブバンドのうちの特定のサブバンドまで画像の特性を判定するものとしてもよい。例えば、サブバンド２ＬＬについてＳ５乃至Ｓ９と同様の処理を行って上記画像の特性を判断した後に、特性解析部１７７が、サブバンド２ＨＬ、２ＬＨ、２ＨＨについてそれぞれヒストグラムを作成し、高周波成分であるサブバンド２ＨＬ、２ＬＨ、２ＨＨの画像の特性を、当該ヒストグラムを用いて分析する。入力画像からサブバンド２ＬＬで求められた低周波成分を除いた画像の成分がサブバンド２ＨＬ、２ＬＨ、２ＨＨであるため、サブバンド２ＨＬ、２ＬＨ、２ＨＨの部分は、サブバンド２ＬＬとの差と考えることができ、サブバンド２ＨＬ、２ＬＨ、２ＨＨの画像特性を示すヒストグラムとなる。このヒストグラムにおいて、高濃度や低濃度の部分に存在する画素の数により、特性解析部１７７は、サブバンド２ＨＬ、２ＬＨ、２ＨＨの画像特性が、文字画像、網点画像、写真画像のいずれであるかを判定する。
【００４８】
これらサブバンド２ＨＬ、２ＬＨ、２ＨＨのヒストグラムでは、(a)文字画像であれば、高濃度と低濃度に文字エッジの量に従った画素数の分布が存在し、(b)網点画像であれば、高濃度と低濃度に網点の線数に対応した一定数以上の画素数が存在し（網点が画像一面に存在するために多くの高周波成分が存在する）、(c)写真画像であれば、高濃度と低濃度の画素数は少ない（エッジとなる部分がほとんど存在しない）。
【００４９】
このため、高周波成分であるサブバンド２ＨＬ、２ＬＨ、２ＨＨにおける高濃度及び低濃度の画素数は、網点画像＞文字画像＞写真画像の関係となる。特性解析部１７７は、各画像種における高濃度及び低濃度の画素数の各境界部分を閾値として設定し、サブバンド２ＨＬ、２ＬＨ、２ＨＨのそれぞれのヒストグラムが示す画素分布が、網点画像、文字画像、又は写真画像のいずれに属するかを判定することで、サブバンド２ＨＬ、２ＬＨ、２ＨＨの画像特性を判定する。
【００５０】
このように判定された各サブバンド２ＨＬ、２ＬＨ、２ＨＨと、既に判定されているサブバンド２ＬＬの各画像特性に基づいて、例えば、サブバンド２ＬＬ、２ＨＬ、２ＬＨ、２ＨＨのうちで最も多くのサブバンドが示す画像特性を、これらサブバンドを含むタイル全体の画像特性とする（画像特性が同数の場合はサブバンド２ＬＬの画像特性とする）。量子化部１７４は、当該特定した画像特性に対応する量子化テーブルを記憶部１２から読み出して、当該読み出した量子化テーブルを上記のように圧縮率を変更し、当該圧縮率変更後の量子化テーブルを、上記タイル全体の全サブバンド（２ＬＬ〜１ＨＨ）の量子化に用いる量子化テーブルとする。
【００５１】
また、さらには、各サブバンド１ＬＬ〜２ＨＨのそれぞれについて、上記Ｓ５〜Ｓ７の処理と同様にして画像特性を求めると共に、各画像特性のそれぞれに対応する量子化テーブルを当該各サブバンド毎に記憶部１２が記憶しておき、量子化部１７４が、これら各サブバンド１ＬＬ〜２ＨＨの各画像特性に応じた量子化テーブルを各サブバンド毎に読み出し、当該読み出した量子化テーブルを用いて、各サブバンド別に量子化を行うようにしてもよい。
【００５２】
なお、上記図１乃至図７に示した実施形態に係る構成及び処理は、本発明に係る画像処理装置、画像処理方法、及び画像処理プログラムの構成及び処理の単なる一例に過ぎず、本発明に係る画像処理装置、画像処理方法、及び画像処理プログラムの構成及び処理を上記に示した内容に限定するものではない。
【符号の説明】
【００５３】
１ 画像処理装置
１２ 記憶部
１２１ 画像圧縮プログラム
１７ 画像圧縮部
１７１ 表色系変換部
１７２ タイル分割部
１７３ 離散ウェーブレット変換部
１７４ 量子化部
１７５ 符号化部
１７７ 特性解析部
【特許請求の範囲】
【請求項１】
画像データを輝度成分と色差成分からなる表色系に変換する表色系変換部と、
前記表色系変換部によって変換された画像データを複数のタイルに分割する分割部と、
前記分割部によって分割されたタイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶する量子化テーブル記憶部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析部と、
前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化部と、
前記量子化部によって量子化された画像データを符号化して圧縮データを出力する符号化部と
を備えた画像処理装置。
【請求項２】
前記量子化部は、前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを、予め設定されている圧縮率に応じた係数を用いて算出した値に変更し、当該変更した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する請求項１に記載の画像処理装置。
【請求項３】
前記特性解析部は、前記サブバンドデータの周波数成分と、当該サブバンドデータを構成する各画素値の分布を示すヒストグラムとを用いて、当該サブバンドデータの特性を解析する請求項１又は請求項２に記載の画像処理装置。
【請求項４】
前記特性解析部は、前記サブバンドデータの前記特性解析に加えて、前記タイルを構成する他の予め定められたサブバンドデータについても前記特性解析を行い、前記量子化部は、当該特性解析により得られた各サブバンドデータの特性に応じた量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、当該タイルの各サブバンドデータを量子化する請求項１乃至請求項３のいずれかに記載の画像処理装置。
【請求項５】
画像データを輝度成分と色差成分からなる表色系に変換する表色系変換ステップと、
前記表色系変換ステップで変換された画像データを複数のタイルに分割する分割ステップと、
前記分割ステップで分割された前記タイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換ステップと、
前記離散ウェーブレット変換ステップで離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析ステップと、
前記離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶しており、前記特性解析ステップで解析された特性に対応する量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化ステップと、
前記量子化ステップによって量子化された画像データを符号化して圧縮データを出力する符号化ステップと
を有する画像処理方法。
【請求項６】
コンピュータを、
画像データを輝度成分と色差成分からなる表色系に変換する表色系変換部と、
前記表色系変換部によって変換された画像データを複数のタイルに分割する分割部と、
前記分割部によって分割されたタイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶する量子化テーブル記憶部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析部と、
前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化部と、
前記量子化部によって量子化された画像データを符号化して圧縮データを出力する符号化部と
して機能させる画像処理プログラム。
【請求項１】
画像データを輝度成分と色差成分からなる表色系に変換する表色系変換部と、
前記表色系変換部によって変換された画像データを複数のタイルに分割する分割部と、
前記分割部によって分割されたタイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶する量子化テーブル記憶部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析部と、
前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化部と、
前記量子化部によって量子化された画像データを符号化して圧縮データを出力する符号化部と
を備えた画像処理装置。
【請求項２】
前記量子化部は、前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを、予め設定されている圧縮率に応じた係数を用いて算出した値に変更し、当該変更した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する請求項１に記載の画像処理装置。
【請求項３】
前記特性解析部は、前記サブバンドデータの周波数成分と、当該サブバンドデータを構成する各画素値の分布を示すヒストグラムとを用いて、当該サブバンドデータの特性を解析する請求項１又は請求項２に記載の画像処理装置。
【請求項４】
前記特性解析部は、前記サブバンドデータの前記特性解析に加えて、前記タイルを構成する他の予め定められたサブバンドデータについても前記特性解析を行い、前記量子化部は、当該特性解析により得られた各サブバンドデータの特性に応じた量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、当該タイルの各サブバンドデータを量子化する請求項１乃至請求項３のいずれかに記載の画像処理装置。
【請求項５】
画像データを輝度成分と色差成分からなる表色系に変換する表色系変換ステップと、
前記表色系変換ステップで変換された画像データを複数のタイルに分割する分割ステップと、
前記分割ステップで分割された前記タイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換ステップと、
前記離散ウェーブレット変換ステップで離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析ステップと、
前記離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶しており、前記特性解析ステップで解析された特性に対応する量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化ステップと、
前記量子化ステップによって量子化された画像データを符号化して圧縮データを出力する符号化ステップと
を有する画像処理方法。
【請求項６】
コンピュータを、
画像データを輝度成分と色差成分からなる表色系に変換する表色系変換部と、
前記表色系変換部によって変換された画像データを複数のタイルに分割する分割部と、
前記分割部によって分割されたタイルの画像データに対して離散ウェーブレット変換を施す離散ウェーブレット変換部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルのサブバンドデータの量子化に用いる量子化テーブルを、当該タイルの特性毎にそれぞれ記憶する量子化テーブル記憶部と、
前記離散ウェーブレット変換部によって離散ウェーブレット変換された前記タイルの各サブバンドデータのうち、少なくとも水平及び垂直方向に低域通過フィルタリングが行われたサブバンドデータの特性を解析し、当該サブバンドの特性を当該タイルの特性とする特性解析部と、
前記特性解析部によって解析された特性に対応する量子化テーブルを前記量子化テーブル記憶部から読み出し、当該読み出した量子化テーブルを用いて、前記離散ウェーブレット変換された前記タイルの各サブバンドデータを量子化する量子化部と、
前記量子化部によって量子化された画像データを符号化して圧縮データを出力する符号化部と
して機能させる画像処理プログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【公開番号】特開２０１２−１１４７８６（Ｐ２０１２−１１４７８６Ａ）
【公開日】平成２４年６月１４日（２０１２．６．１４）
【国際特許分類】
【出願番号】特願２０１０−２６３４５９（Ｐ２０１０−２６３４５９）
【出願日】平成２２年１１月２６日（２０１０．１１．２６）
【出願人】（０００００６１５０）京セラドキュメントソリュ−ションズ株式会社 (13,173)
【Ｆターム（参考）】
【公開日】平成２４年６月１４日（２０１２．６．１４）
【国際特許分類】
【出願日】平成２２年１１月２６日（２０１０．１１．２６）
【出願人】（０００００６１５０）京セラドキュメントソリュ−ションズ株式会社 (13,173)
【Ｆターム（参考）】
[ Back to top ]