"撮影画像を読み取る画像読み取り工程S101と,前記画像読み取り工程で読み込んだ撮影画像を周波数解析し,前記撮影画像の波形データを得る周波数解析工程S102と,前記周波数解析工程で解析したピーク周波数のm倍の第一の周波数変換データ及びピーク周波数のn倍の第二の周波数変換データを得る周波数変換工程S103と,前記周波数変換工程で得られた前記第一の周波数変換データ及び前記第二の周波数変換データの共通部分データを抽出する共通部分抽出工程S104と前記共通部分抽出工程で得られた共通部分データを可逆圧縮して第1の圧縮データを得る第1の圧縮工程S105を有し,前記共通部分抽出工程S104においては,前記周波数変換工程で得られた前記第一の周波数変換データ及び前記第二の周波数変換データを波形データとして読み出し,前記第一の周波数データの波形データと前記第二の周波数データの波形データの位相を1ビット分ずつずらしたときに,前記第一の周波数変換データ及び前記第二の周波数変換データが所定の範囲で共通する値を持つ部分を前記共通部分データとして抽出するという操作を,前記撮影画像の波形データの1周期分繰り返す画像圧縮方法。前記撮影画像は医療機関における医療画像である請求項1に記載の画像圧縮方法。前記周波数変換工程S103は2次元フーリエ変換による周波数解析により前記ピーク周波数を解析する工程を含む請求項1に記載の画像圧縮方法。前記m倍は2倍であり前記n倍は3倍である請求項3に記載の画像圧縮方法。前記第1の圧縮工程は,バイナリー可逆圧縮工程である,請求項1に記載の画像圧縮方法。前記周波数変換工程で得られた第一の周波数変換データ及び前記第二の周波数変換データの相違部分データを抽出する相違部分抽出工程S111と前記相違部分抽出工程S111で得られた相違部分データを圧縮して第2の圧縮データを得る第2の圧縮工程S112と,前記第1の圧縮データと前記第2の圧縮データとを合わせて圧縮画像データを得るための圧縮画像データ取得工程S113とを更に含む請求項1に記載の画像圧縮方法。撮影画像を読み取る画像読み取り手段1と,前記画像読み取り手段が読み込んだ撮影画像を周波数変換するための周波数解析手段2と,前記周波数解析手段が解析したピーク周波数のm倍の第一の周波数変換データ及びピーク周波数のn倍の第二の周波数変換データを得る周波数変換手段3と,前記周波数変換手段が得た第一の周波数変換データ及び前記第二の周波数変換データの共通部分データを抽出する共通部分抽出手段4と前記共通部分抽出手段が得た共通部分データを可逆圧縮して第1の圧縮データを得る第1の圧縮手段5と,を有する画像圧縮システム。しかしながら従来の画像圧縮方法では画像圧縮の圧縮量が限られていた。特に病院における医療画像は,連続撮影する場合が多く撮影データの量が病院内のサーバに大きな負荷を与えていた。また医療画像は圧縮する場合であっても可逆圧縮を行う必要があるものも多く存在するため可逆圧縮を行いつつ画像量を軽減する方法が望まれた。本発明はきわめて効率的に画像を圧縮できる画像圧縮方法を提供することを目的とする。本発明は画像撮影装置から発生した元画像と比較して全ての画素値が数学的に全く同一である医療画像の可逆圧縮方法を提供することを目的とする。本発明はさらに従来の方法よりも数倍以上高い圧縮効率の可逆圧縮を行うことができる画像圧縮方法を提供することを目的とする。本発明は基本的には,複数種類の周波数でオーバーサンプリングした画像の共通部分のデータを可逆圧縮することで,きわめて効率的に画像を圧縮できる画像圧縮方法を提供できるという知見に基づくものである。本発明の第1の側面は,画像圧縮方法に関する。この方法は,画像読み取り工程S101と,周波数解析工程S102と,周波数変換工程S103と,共通部分抽出工程S104と第1の圧縮工程S105とを含む。すなわち,この方法は,オリジナルの撮影画像を読み取る。そして,読み込んだ撮影画像を,たとえば,2次元フーリエ変換による周波数解析して,波形データを求め,ピーク周波数を求める。また,本発明は,少なくとも2種類の周波数でオーバーサンプリングを行う。すなわち,本発明では,解析したピーク周波数のm倍の第一の周波数変換データ及びピーク周波数のn倍の第二の周波数変換データを得る。m倍の例は2倍でありn倍の例は3倍である。そして,第一の周波数変換データ及び第二の周波数変換データの共通部分データを抽出する。そして,抽出した共通部分データを,たとえばバイナリー可逆圧縮して第1の圧縮データを得る。このようにすることで,本発明は,きわめて効率的に画像を圧縮できる。この方法は,特に医療機関における医療画像の圧縮に有効である。具体的に説明すると,共通部分抽出工程S104においては,周波数変換工程で得られた第一の周波数変換データ及び第二の周波数変換データを波形データとして読み出し,オリジナルの撮影画像の波形データ,第一の周波数変換データの波形データ,及び第二の周波数変換データの波形データを重ね合わせる。この第一の周波数データの波形データと第二の周波数データの波形データの位相を1ビット分ずつずらしたときに,第一の周波数変換データ及び第二の周波数変換データが所定の範囲で共通する値を持つ部分を共通部分データとして抽出する。このようにして,各共通部分データを抽出する操作を,撮影画像の波形データの1周期分繰り返すことにより,第一の周波数変換データと第二の周波数変換データの共通部分データをすべて抽出する。この側面の好ましい態様は,相違部分抽出工程S111と第2の圧縮工程S112と,圧縮画像データ取得工程S113とを更に含むものである。すなわち,この態様の画像圧縮方法は,第一の周波数変換データ及び第二の周波数変換データの相違部分データを抽出する。相違部分データとは,第一の周波数変換データ及び第二の周波数変換データの共通部分ではない部分のデータである。次に,相違部分データを圧縮して第2の圧縮データを得る。そして,第1の圧縮データと第2の圧縮データとを合わせて圧縮画像データを得る。この方法は,複数種類の周波数でオーバーサンプリングした画像の非共通データをも可逆圧縮することで,可逆圧縮を行うことができる。本発明は,上記の方法を実現できる画像圧縮システム,コンピュータを上記した画像圧縮システムとして実現するためのプログラム,そのようなプログラムを記憶したコンピュータが読み取ることのできる情報記録媒体をも提供する。本発明は複数種類の周波数でオーバーサンプリングした画像の共通データを可逆圧縮することで,きわめて効率的に画像を圧縮できる画像圧縮方法を提供できる。本発明は,複数種類の周波数でオーバーサンプリングした画像の共通データを可逆圧縮することで特に骨の部分正常部分及び罹患部分を明瞭に区別できる医療画像の圧縮方法を提供できる。本発明は複数種類の周波数でオーバーサンプリングした画像の非共通データをも可逆圧縮することで,可逆圧縮を行うことができる画像圧縮方法を提供できる。以下図面を用いて本発明を実施するための最良の形態について説明する。図1は本発明の第1の画像圧縮方法を説明するためのフロー図である。図1に示されるように,この画像圧縮方法は画像読み取り工程S101と,周波数解析工程S102と,周波数変換工程S103と,共通部分抽出工程S104と,第1の圧縮工程S105と,を有する。図2は,本発明の画像圧縮装置のブロック図である。図2に示されるように,この装置は,画像読み取り手段1と,周波数解析手段2と,周波数変換手段3と,共通部分抽出手段4と第1の圧縮手段5とを有する。この画像圧縮装置は,いわゆるコンピュータにより入力画像を圧縮するための装置であってもよい。この場合,コンピュータは,入出力部と,制御部と,演算部と,記憶部とを有する。そして,これらの要素は,バスなどで接続され,情報の授受を行うことができるようにされている。以下に説明する制御工程及び制御手段はハードウェアにより実現されてもよいし,ハードウェアとソフトウェアとの協働により実現されてもよい。後者の場合,主記憶部に制御プログラムが記憶され,各種情報が入力された場合に,主記憶部に記憶された制御プログラムが読みだされる。そして,制御プログラムの指令に従って,各種データが記憶部から読みだされ,演算部で演算処理が行われる。そして,適宜,情報が入出力部から出力される。このような動作を制御部が制御する。画像読み取り工程S101は,オリジナルの撮影画像を読み取るための工程である。撮影画像の例は医療機関における医療画像である。医療機関における医療用の画像は,量が膨大である。このため,医療用の画像を圧縮したいという要望が強い。医療用の医療画像の例は,レントゲン画像やCTスキャン画像である。そして医療画像は一次診断に用いる場合,非圧縮の元画像を用いるか,または圧縮する場合であっても全ての画素値が元画像と数学的に同一である可逆圧縮を用いることが求められる。画像読み取り工程は,撮影されてシステムの記憶部に記憶された画像を読み取る工程であってもよい。また,画像読み取り工程は,レントゲンまたはCTスキャンなどの撮影装置を用いて画像を撮影し,本発明のシステムへ読み取って入力する工程であってもよい。すなわち,レントゲン撮影装置,CTスキャン装置が画像読み取り手段1として機能する。また,コンピュータ内の記憶部又はコンピュータに接続されたデータベースに記憶された撮影画像を読み出す手段も画像読み取り手段1として機能する。周波数解析工程S102は,画像読み取り工程で読み込んだ撮影画像を周波数解析するための工程である。周波数解析工程S102により,撮影画像の波形データを得られる。周波数解析する手段2は,すでに知られている。画像を周波数解析する手段の例は,2次元高速フーリエ変換2次元FFTを用いた解析装置である。周波数解析する手段2は,画像読み取り手段1が記憶部から読みだした撮影画像に対して,周波数解析を行うための演算処理を行う。そして,撮影画像のピーク周波数を解析する。ピーク周波数の例は,サンプリング定理で決定される周波数である。周波数解析する手段2は,求めたピーク周波数を波形データと共に適宜記憶部に記憶する。ピーク周波数は,2次元フーリエ変換による周波数解析によりピーク周波数を解析することで得ることができる。2次元フーリエ変換の例は2次元高速フーリエ変換2次元FFTである。周波数変換工程S103は,周波数解析工程で解析して得られたピーク周波数のm倍の第一の周波数変換データ及びピーク周波数のn倍の第二の周波数変換データを得るための工程である。画像を周波数変換する工程は,たとえば,特開2006-191216号公報に開示される通りすでに知られている。本発明は,ピーク周波数のm倍の周波数及びn倍の周波数を用いてそれぞれサンプリングを行って,第一及び第二の周波数変換データを得る。第一及び第二の周波数変換データは,例えばバイナリデータとして記憶部に記憶される。m倍及びn倍の例は2倍及び3倍である。m及びnの例は正の整数であるが,少数を含む正数であってもよい。n及びmの好ましい例は,n及びmが正の整数であり,かつ素数である。n及びmは,1より大きい数であることが好ましい。n及びmが1より大きい場合,これらの周波数を用いてサンプリングすると,いわゆるオーバーサンプリングを行ったこととなる。周波数変換手段3は,画像読み取り手段が読み込んだ撮影画像を周波数変換し,ピーク周波数のm倍の第一の周波数変換データ及びピーク周波数のn倍の第二の周波数変換データを得るための要素である。周波数変換手段3は,記憶部から周波数解析する手段2が求めたピーク周波数を読みだす。また,周波数変換手段3は,記憶部から読みだされた撮影画像について,ピーク周波数のm倍の周波数及びn倍の周波数を用いてそれぞれサンプリングする演算処理を行い,第一及び第二の周波数変換データを得て,記憶部に記憶する。具体的なサンプリング方法は,たとえば,特開平09-121296号公報,特開2006-191216号公報及び特開2007-192608号公報に開示される方法に従ってサンプリングを行えばよい。本発明は,これらの公知文献を含む公知技術を適宜組み合わせて用いることができる。このサンプリング処理及び記憶部への記憶処理により,撮影画像が周波数変換処理されることとなる。共通部分抽出工程S104は,周波数変換工程で得られた第一の周波数変換データ及び第二の周波数変換データの共通部分データを抽出するための工程である。第一の周波数変換データと第二の周波数変換データとで所定の範囲で共通する値を示す画像データを抽出することで共通部分データを抽出することができる。共通部分抽出手段4は,周波数変換手段が得た第一の周波数変換データ及び第二の周波数変換データの共通部分データを抽出するための要素である。画像読み取り手段1が読み取って記憶部に一時的に記憶した撮影画像の波形データを周波数変換手段2が読み出し,第一の周波数変換データ及び第二の周波数変換データの共通部分データを抽出する演算処理を行う。この工程では,共通部分抽出手段4が,記憶部から第一の周波数変換データ及び第二の周波数変換データを,波形データとして読みだす。そして,例えば,共通部分抽出手段4は,読みだした撮影画像の波形データと,第一の周波数変換波形データ及び第二の周波数変換波形データを重ね合わせる。そして,共通部分抽出手段4は,第一の周波数変換データの波形データと第二の周波数変換データの波形データの位相を1ビット分ずつずらし,第一の周波数変換波形と第二の周波数変換波形が,所定の範囲で共通する値を持つ点を検出することにより,その値を共通部分データとして記憶する。続いて,第一の周波数変換データの波形と第二の周波数変換データの波形の位相を,さらに1ビット分ずつずらし,再び所定の範囲で共通する値を持つ点を検出することにより,その値を共通部分データとして記録する。このような操作をオリジナルの撮影画像の波形データの1周期分繰り返し,その都度共通部分データの値を記録することで,第一の周波数変換波形データ及び第二の周波数変換波形データのすべての共通部分データを抽出できる。この場合の共通する情報は,波形値が全く同じものだけではなく,所定範囲のデータを共通するものとして抽出してもよい。n倍の第一の周波数変換波形データとm倍の第二の周波数変換波形データを用いることにより,オリジナルの撮影画像の波形データの1周期内にはそれぞれn個及びm個のピークが存在するため,共通部分データを与える位相値の部位数はn×m倍となる。この操作を撮影画像の波形データの1周期分繰り返すことにより,第一周波数変換波形データと第二の周波数変換波形データの共通部分データを効率よく抽出することができる。従って,単にオリジナルの画像データのみの中から共通部分を抽出して画像圧縮を行うよりも,数倍以上の効率で圧縮を行うことが可能となる。なお,上記の共通部分に相当しない部分,すなわち後述する相違部分データは,可逆圧縮を行う際には,そのまま保持することが必要となる。第1の圧縮工程S105は,共通部分抽出工程で得られた共通部分データを可逆圧縮して第1の圧縮データを得るための工程である。第1の圧縮工程の例は,バイナリー可逆圧縮工程である。第1の圧縮手段5は,共通部分抽出手段が得た共通部分データを可逆圧縮して第1の圧縮データを得るための要素である。可逆圧縮の方法は,画像処理の分野において公知である。よって,第1の圧縮手段5として公知の装置を適宜用いることができる。第1の圧縮手段5による具体的な圧縮方法は,共通部分データを3バイトずつ1つのブロックとし,1バイトずつずらして順次圧縮するものである。この圧縮に先だって,多値画像グレイデータを所定の閾値を用いて2値データに変換してもよい。このようにして,共通部分データをバイナリー可逆圧縮できることとなる。この画像圧縮システムは,先に説明した画像圧縮方法を実現できる。この方法は,全ての画素値を元画像と同一に保ちしつつ画像容量を飛躍的に圧縮できるため,レントゲンやCTスキャン画像といった医療用画像を圧縮することに有効に利用することができる。従来の圧縮方法では,可逆圧縮の場合最大でも2.5分の1程度の圧縮率が限界であるという問題があったが,本発明の圧縮方法では,このような問題を適切に解消できた。この圧縮したデータを解凍処理することで,原画像を再現できる。なお,上記の方法では,第一の周波数変換データ及び第二の周波数変換データの共通部分以外のデータを捨ててしまってもよい。この場合,捨ててしまったデータを復元できないので,処理が非可逆となる。本発明の好ましい態様は,さらに可逆変換を実行できるものである。すなわち医療過誤事件においては,圧縮画像だけではなく撮影画像そのものを復元し,これを証拠として用いるという要求がある。本発明の方法は,可逆圧縮をも達成できるので,以下この点について説明する。図3は本発明の第2の画像圧縮方法を説明するためのフロー図である。図3に示されるように,この方法は,先に説明した方法に比べ,相違部分抽出工程S111と,第2の圧縮工程S112と,圧縮画像データ取得工程S113とをさらに含む。これを実現するためには,画像圧縮装置が,相違部分抽出部と,第2の圧縮工程部と,圧縮画像データ取得部とをさらに有すればよい。相違部分抽出工程S111は,周波数変換工程で得られた第一の周波数変換データ及び第二の周波数変換データの相違部分データを抽出する工程である。この工程は,共通部分抽出工程と同時に行われてもよい。そして,抽出された相違部分データは,記憶部に記憶されればよい。第2の圧縮工程S112は,相違部分抽出工程S111で得られた相違部分データを圧縮して第2の圧縮データを得るための工程である。この工程における処理は,たとえば,第1の圧縮工程と同様であればよい。圧縮画像データ取得工程S113は,第1の圧縮データと第2の圧縮データとを合わせて圧縮画像データを得るための工程である。圧縮画像データ取得手段が,第1の圧縮データと第2の圧縮データを記憶部に記憶させればよい。この際,第1の圧縮データと第2の圧縮データは,画像データであるから,画像以外のデータと合わせて読みだすことができるようにして記憶部に記憶させればよい。このようにオーバーサンプリングを行ったデータを合わせる方法は,たとえば,特開2007-192608号公報に開示された方法を用いればよい。図3に示されるフローチャートを実現するアルゴリズムに基づき,画像圧縮用のソフトウェアを作成した。作成したソフトウェアをコンピュータにインストールして,画像圧縮システムを得た。図4は,医療画像を圧縮した際の画像を示す。図4Aは,元画像であるレントゲン写真である。図4Bは,JPEGによる可逆圧縮画像を示す。図4Cは,本発明の圧縮方法による可逆圧縮画像を示す。このように本発明の画像圧縮方法は,JPEGによる圧縮画像や元画像と比べて画像精度は数学的に同一であり,肉眼的にも区別できない。さらに,元画像の画像容量は19.6Mバイトであり,JPEGによる可逆圧縮画像の容量は10.2Mバイトであった。これに対し,図4Cに示される本発明の可逆圧縮方法による可逆圧縮画像の容量は,0.94Mバイトであった。すなわち,JPEGによる可逆圧縮では,元データの1/2程度にしか容量を圧縮できない。これに対し,本発明の可逆圧縮方法では,1/20程度まで容量を圧縮できる。このように,本発明は,従来方法と比較して容量をきわめて圧縮できることが分かる。上記と同様の実験を様々な画像データを用いて解析した。その結果を表1及び表2に示す。表中,本発明の可逆とは,基本的には図3に示されるフローチャートに従って可逆圧縮したものを示す。一方,本発明の非可逆とは,基本的には図1に示されるフローチャートに従って非可逆圧縮したものを示す。表1は,静止画の圧縮結果を示す。一方,表2は,動画の圧縮結果を示す。表1及び表2に示されるように,本発明の圧縮方法は,あらゆる医療画像及び医療動画について容量をきわめて圧縮できることが分かる。実施例2では,図5から図8を参照して,本発明における共通部分抽出工程を具体的に説明する。実施例2においては,フーリエ変換によって周波数解析した撮影画像のピーク周波数に対して,第一の周波数変換データと第二の周波数変換データは,それぞれ2倍と3倍の周波数でオーバーサンプリングしている。従って,例えば,50個のデータビットに対する撮影画像のサンプリングレート周波数を30とした場合,第一の周波数変換データと第二の周波数変換データのサンプリングレート周波数は,それぞれ60と90になる。図5aにおいては,撮影画像データの波形データ,第一の周波数変換データの波形データ,及び第二の周波数変換データの波形データを重ねあわせた初期位置を示している。横軸は撮影画像のデータ位置を示し,縦軸は撮影画像の濃淡値を示している。初期位置を示す図5aのグラフ1では,第一の周波数変換データと第二の周波数変換データには共通部分は存在していない。また,図5aに示された状態は,図6に示すデータ表と対応している。図6を参照しても,第一の周波数変換データと第二の周波数変換データには共通部分は存在していないことがわかる。図5bでは,初期データから,第一の周波数変換データの波形データと第二の周波数変換データの波形データを,それぞれ同一方向への位相を1ビットずつずらした状態を示している。図5bのグラフ2のように,第一の周波数変換データの波形データと第二の周波数変換データの波形データを1ビットずつずらしてグラフに表示すると,2倍と3倍のサンプリング周波数のグラフに,所定の範囲で共通する値をとる部分グラフ2の赤い丸で囲った部分が出現することがわかる。図5bに示された状態は,図7に示すデータ表と対応している。図7を参照すると,第一の周波数変換データと第二の周波数変換データは,データ位置4及びデータ位置38において,共通部分薄墨が施された箇所が存在していることがわかる。図5bでは,初期データから,第一の周波数変換データの波形データと第二の周波数変換データの波形データを,それぞれ同一方向への位相を2ビットずつずらした状態を示している。図5cのグラフ3を参照すると,図5bのグラフ2とは異なる位置に,第一の周波数変換データと第二の周波数変換データが所定の範囲で共通する値をとる部分グラフ3の赤い丸で囲った部分が出現することがわかる。図5cに示された状態は,図8に示すデータ表と対応している。図8を参照すると,第一の周波数変換データと第二の周波数変換データは,データ位置8及びデータ位置40おいて,共通部分薄墨が施された箇所が存在していることがわかる。このように1ビットずつ順次サンプリング位置をずらし,共通部分のデータを抽出する操作を,撮影画像データの波形データの1周期分行うことにより,共通部分データと相違部分データを分けて抽出し,それぞれ可逆圧縮することができる。"