(57)【要約】      （修正有）【課題】複数の画像の差異をより速く検出することができる技術の提供する。【解決手段】画像差異検出装置は、どの領域に差異が発生するか不明であり、ほぼ同一の空間領域が撮影されている第１の画像及び第２の画像に関して、第１の画像及び第２の画像を符号化した第１及び第２の符号化データにおける差異を、第１の符号化データにおける第１の符号化ブロックと、第２の符号化データにおける第１の符号化ブロックと対応する第２の符号化ブロックごとに、第１の符号化ブロックにおける指標値と第２の符号化ブロックにおける指標値と、に基づき検出し、第１の符号化ブロックと第２の符号化ブロックを、符号化する際に得られる符号化情報に基づいて、平坦領域か複雑領域かに分類する。検出部は、平坦領域と分類された場合、第一の指標値と第二の指標値に基づき、複雑領域と分類された場合、第一の指標値に基づき、差異を検出する。【選択図】図１
2020025290.tif 000002
【特許請求の範囲】
【請求項１】
どの領域に差異が発生するか不明であり、ほぼ同一の空間領域が撮影されている第１の画像及び第２の画像に関して、前記第１の画像及び前記第２の画像を符号化した第１及び第２の符号化データにおける差異を、前記第１の符号化データにおける第１の符号化ブロックと、前記第２の符号化データにおける前記第１の符号化ブロックと対応する第２の符号化ブロックごとに、前記第１の符号化ブロックにおける指標値と前記第２の符号化ブロックにおける指標値と、に基づき検出する検出部と、
前記第１の符号化ブロックと前記第２の符号化ブロックを、符号化する際に得られる符号化情報に基づいて、平坦領域か複雑領域かに分類する分類部と、を有し、
前記検出部は、
平坦領域と分類された場合、第一の指標値と第二の指標値に基づき、
複雑領域と分類された場合、前記第一の指標値に基づき、差異を検出する
画像差異検出装置
【請求項２】
前記第一の指標値は、符号化ブロックの符号量である
請求項1記載の画像差異検出装置
【請求項３】
前記第二の指標値は、前記符号化ブロックの符号量と、前記符号化ブロックにおける周辺の符号化ブロックの符号量と、の分散値である
請求項２記載の画像差異検出装置
【請求項４】
前記分類部は、前記第１の符号化ブロックと、前記第２の符号化ブロックが分類された結果が異なる場合、前記第１の画像と前記第２の画像の撮影日時や画質に基づき定められた優先度に基づき分類する
請求項３記載の画像差異検出装置
【請求項５】
前記第一の指標値は、符号化ブロックの符号量、前記符号化ブロックにおける２値化ビット値幅、前記符号化ブロックにおけるDCT係数のいずれかである
請求項１記載の画像差異検出装置
【請求項６】
どの領域に差異が発生するか不明であり、ほぼ同一の空間領域が撮影されている第１の画像及び第２の画像に関して、前記第１の画像及び前記第２の画像を符号化した第１及び第２の符号化データにおける差異を、前記第１の符号化データにおける第１の符号化ブロックと、前記第２の符号化データにおける前記第１の符号化ブロックと対応する第２の符号化ブロックごとに、前記第１の符号化ブロックにおける指標値と前記第２の符号化ブロックにおける指標値と、に基づき検出する検出ステップと、
前記第１の符号化ブロックと前記第２の符号化ブロックを、符号化する際に得られる符号化情報に基づいて、平坦領域か複雑領域かに分類する分類ステップと、を有し、
前記検出ステップは、
平坦領域と分類された場合、第一の指標値と第二の指標値に基づき、
複雑領域と分類された場合、前記第一の指標値に基づき、差異を検出する
画像差異検出方法
【請求項７】
請求項1から請求項5に記載された検出部と分類部とをコンピュータで実行させるためのコンピュータプログラム
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の画像間の差異を検出する技術に関する。
本願は、２０１５年１２月１５日に日本に出願された特願２０１５−２４４５５６号と、２０１６年６月３０日に日本に出願された特願２０１６−１３１１５５号と、２０１６年６月３０日に日本に出願された特願２０１６−１３１１６２号とに基づき優先権を主張し、その内容をここに援用する。
【背景技術】
【０００２】
近年、監視用動画像からのイベント検出、医用における異常検出、航空・衛星画像からの変化検出などの複数の画像間の差異を検出する技術が注目されている。これに対し、動画像内からのイベントの有無を検出手法が提案されている（例えば、特許文献１参照）。本手法では、対象画像の複雑さ指標値（量子化ステップ幅及び発生符号量から算出）と、時間的に直前にあたる画像の複雑さ指標値との比を基に、対象画像で変化があったか否かを検出する。すなわち、時間軸方向の画像の変化を検出することは可能であり、監視用動画像からのイベント検出には適用することができた。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１４−８６９１３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、画像内のどの領域で変化が起きたかは検出することができず、医用における異常検出や、航空・衛星画像からの変化検出には適用することができない。特に、特許文献１に記載の映像解析装置にあっては、領域の細分化を適用（画像をブロックに分割し、各ブロックの複雑さ指標値を利用）するだけでは、医用の異常検出や航空・衛星画像からの変化検出はできないという問題がある。理由としては、時間軸方向の変化のみに着目しており、空間的な変化は一切考慮されていない。このため、画像全体の特性（例えば明るさや精細感など）が少し異なるだけで、各ブロックの複雑さ指標値も異なってしまう。単純に時間軸方向で直前の画像内の空間的に同一なブロックの複雑さ指標値と対象画像内の対象ブロックの複雑さ指標値比をみるだけでは全てのブロックで一様に変化ありとなってしまうという問題がある。
【０００５】
上記事情に鑑み、本発明は、複数の画像の差異をより速く検出することができる技術の提供を目的とする。
【課題を解決するための手段】
【０００６】
本発明の第１の実施態様における画像差異検出装置は、第１及び第２の画像を符号化した第１及び第２の符号化データにおける符号化ブロックそれぞれの符号量と、前記符号化ブロックそれぞれの符号化情報から得られる指標値とのいずれか一方又は両方に基づいて、前記第１の画像内の第１の領域と、前記第２の画像内において前記第１の領域に対応する第２の領域とにおける差異の有無を検出する検出部、を備える。
【０００７】
本発明の第２の実施態様によれば、上記第１の実施態様の画像差異検出装置において、前記第１及び第２の符号化データを得る際に用いられる符号化では、可変長符号が用いられる。
【０００８】
本発明の第３の実施態様によれば、上記第１又は第２の実施態様の画像差異検出装置において、前記第１及び第２の符号化データから前記符号化ブロックそれぞれの前記符号化情報を取得する符号化情報取得部、を更に備え、前記検出部は、前記符号化情報取得部により取得された前記符号化情報から得られる前記指標値に基づいて、前記第１の領域と前記第２の領域とにおける差異の有無を判定する。
【０００９】
本発明の第４の実施態様によれば、上記第１又は第２の実施態様の画像差異検出装置において、前記第１及び第２の画像において差異の有無が検出される前記第１及び第２の領域をそれぞれ特定する情報を入力する入力部と、前記入力部が入力した情報に基づいて特定される前記第１及び第２の領域に対応する前記符号化ブロックをそれぞれ特定する領域特定部と、を更に備え、前記検出部は、前記領域特定部により特定された前記符号化ブロックの符号量又は前記符号化情報から得られる前記指標値の一方又は両方に基づいて、特定される前記第１及び第２の領域における差異の有無を検出する。
【００１０】
本発明の第５の実施態様によれば、上記第４の実施態様の画像差異検出装置において、前記入力部は、前記第１及び第２の領域を特定する情報として、前記第１及び第２の画像において差異の検出対象となる被写体が存在した場所を示す位置情報を入力する。
【００１１】
本発明の第６の実施態様によれば、上記第１の実施態様の画像差異検出装置において、前記第１及び第２の画像それぞれを複数の前記符号化ブロックに分割し、前記符号化ブロックそれぞれを符号化し、前記第１及び第２の符号化データを生成する符号化部と、前記第１及び第２の画像それぞれを符号化する過程で前記符号化ブロックごとに得られる前記符号化情報に基づいて、前記符号化ブロックそれぞれを複数のカテゴリに分類する分類部と、を更に備え、前記検出部は、前記分類部による分類結果に基づいて、前記第１の領域と前記第２の領域とにおける差異の有無の判定に用いるしきい値を決定し、前記検出部は、前記符号量と前記指標値とのいずれか一方又は両方と、決定した前記しきい値とに基づいて、前記第１の領域と前記第２の領域とにおける差異の有無を判定する。
【００１２】
本発明の第７の実施態様によれば、上記第１の実施態様の画像差異検出装置において、前記検出部は、前記分類部による分類結果に基づいて、前記第１の領域と前記第２の領域とにおける差異の有無の判定に用いる判定条件を決定し、前記検出部は、前記符号量と前記指標値とのいずれか一方又は両方と、決定した前記判定条件とに基づいて、前記第１の領域と前記第２の領域とにおける差異の有無を判定する。
【００１３】
本発明の第８の実施態様によれば、上記第６又は第７の実施態様の画像差異検出装置において、前記分類部は、前記符号化情報から得られる画像の複雑さ度合いに基づいて、前記符号化ブロックを前記複数のカテゴリに分類する。
【００１４】
本発明の第９の実施態様によれば、上記第６又は第７の実施態様の画像差異検出装置において、前記複数のカテゴリは、平坦領域と複雑領域とを含み、前記分類部は、前記符号化情報から得られる画像の複雑さ度合いに基づいて、前記符号化ブロックそれぞれを前記平坦領域と前記複雑領域とに分類する。
【００１５】
本発明の第１０の実施態様によれば、上記第１の実施態様の画像差異検出装置において、前記第１及び第２の画像において検出される差異の特徴を示す種別情報を入力する種別入力部と、前記種別情報に基づいて量子化マトリクスを決定するパラメータ決定部と、前記第１及び第２の画像それぞれを複数の前記符号化ブロックに分割し、前記量子化マトリクスに基づいて前記符号化ブロックそれぞれを符号化し、前記第１及び第２の符号化データを生成する符号化部と、を更に備える。
【００１６】
本発明の第１１の実施態様によれば、上記第１の実施態様の画像差異検出装置において、前記第１及び第２の画像において符号化を開始する複数の開始点を入力する開始点入力部と、前記開始点に基づいて前記第１及び第２の画像それぞれを複数の前記符号化ブロックに分割し、前記符号化ブロックそれぞれを符号化し、前記開始点ごとに前記第１及び第２の符号化データを生成する符号化部と、を更に備え、前記検出部は、前記第１及び第２の画像それぞれを符号化する過程で前記符号化ブロックごとに算出される前記符号化情報から得られる前記指標値に基づいて、前記第１の領域と前記第２の領域とにおける差異の有無を判定する。
【００１７】
本発明の第１２の実施態様によれば、上記第１の実施態様の画像差異検出装置において、前記第１及び第２の符号化データの総符号量の比に基づいて、前記第１の符号化データ又は前記第２の符号化データにおける前記符号化ブロックそれぞれの符号量を調整する調整部、を更に備え、前記検出部は、前記調整部により調整された符号量に基づいて、前記第１の領域と前記第２の領域とにおける差異の有無を判定する。
【００１８】
本発明の第１３の実施態様によれば、上記第１２の実施態様の画像差異検出装置において、前記調整部は、前記第１及び第２の符号化データの総符号量の比が所定値を超える場合、前記第１の画像と前記第２の画像とが全体的に異なる画像である、又は、差異の検出が不可能であると判定する。
【００１９】
本発明の第１４の実施態様によれば、上記第１の実施態様の画像差異検出装置において、前記第１及び第２の画像を符号化したときの前記第１の画像における前記符号化ブロックの位置と前記第２の画像における符号化ブロックの位置との差分に応じた係数を用いて、前記指標値を補正する補正部、を更に備え、前記検出部は、前記補正部により補正された前記指標値に基づいて、前記第１の領域と前記第２の領域とにおける差異の有無を判定する。
【００２０】
本発明の第１５の実施態様によれば、上記第１４の実施態様の画像差異検出装置において、前記補正部は、予め定められた複数の係数から少なくとも一つの係数を前記差分に基づいて選択し、選択した前記係数を用いて前記指標値を補正する。
【００２１】
本発明の第１６の実施態様における画像差異検出方法は、第１及び第２の画像を符号化した第１及び第２の符号化データにおける符号化ブロックそれぞれの符号量と、前記符号化ブロックそれぞれの符号化情報から得られる指標値とのいずれか一方又は両方に基づいて、前記第１の画像内の第１の領域と、前記第２の画像内において前記第１の領域に対応する第２の領域とにおける差異の有無を検出する検出ステップ、を含む。
【００２２】
本発明の第１７の実施態様におけるコンピュータプログラムは、第１及び第２の画像を符号化した第１及び第２の符号化データにおける符号化ブロックそれぞれの符号量と、前記符号化ブロックそれぞれの符号化情報から得られる指標値とのいずれか一方又は両方に基づいて、前記第１の画像内の第１の領域と、前記第２の画像内において前記第１の領域に対応する第２の領域とにおける差異の有無を検出する検出部、としてコンピュータを機能させる。
【発明の効果】
【００２３】
本発明により、複数の画像の差異をより速く検出することが可能となる。
【図面の簡単な説明】
【００２４】
【図１】本発明の第１実施形態による画像差異検出装置の構成を示すブロック図である。
【図２】図１に示す画像差異検出装置の処理動作を示すフローチャートである。
【図３】本発明の第２実施形態による画像差異検出装置の構成を示すブロック図である。
【図４】図３に示す画像差異検出装置の処理動作を示すフローチャートである。
【図５】本発明の第３実施形態による画像差異検出装置の構成を示すブロック図である。
【図６】図５に示す画像差異検出装置の処理動作を示すフローチャートである。
【図７】本発明の第４実施形態による画像差異検出装置の構成を示すブロック図である。
【図８】本発明の第５実施形態による画像差異検出装置の構成を示すブロック図である。
【図９】図８に示す画像差異検出装置の処理動作を示すフローチャートである。
【図１０】本発明の第６実施形態による画像差異検出装置の構成を示すブロック図である。
【図１１】図１０に示す画像差異検出装置の処理動作を示すフローチャートである。
【図１２】本発明の第７実施形態による画像差異検出装置の構成を示すブロック図である。
【図１３】図１２に示す画像差異検出装置の処理動作を示すフローチャートである。
【図１４】本発明の第８実施形態による画像差異検出装置の構成を示すブロック図である。
【図１５】図１４に示す画像差異検出装置の処理動作を示すフローチャートである。
【図１６】図１４に示す動作処理における変化領域か否かを判定する処理の一例を示すフローチャートである。
【図１７】図１４に示す動作処理における変化領域か否かを判定する処理の他の例を示すフローチャートである。
【図１８】本発明の第９実施形態による画像差異検出装置の構成を示すブロック図である。
【図１９】図１８に示す画像差異検出装置の処理動作を示すフローチャートである。
【図２０】本発明の第１０実施形態による画像差異検出装置の構成を示すブロック図である。
【図２１】図２０に示す画像差異検出装置の処理動作を示すフローチャートである。
【図２２】本発明の第１１実施形態による画像差異検出装置の構成を示すブロック図である。
【図２３】図２２に示す画像差異検出装置の処理動作を示すフローチャートである。
【図２４】図２３に示す処理動作においてＨＥＶＣまたはＨ．２６４にて符号化された符号化データから指標値を得る処理動作を示すフローチャートである。
【図２５】図２３に示す処理動作においてＪＰＥＧにて符号化された符号化データから指標値を得る処理動作を示すフローチャートである。
【図２６】図２３に示す処理動作のうち共通領域における複雑さ指標値の補正処理を示すフローチャートである。
【図２７】共通領域における複雑さ指標値の補正処理の説明図である。
【発明を実施するための形態】
【００２５】
以下、図面を参照して、本発明の実施形態による画像差異検出装置を説明する。実施形態による画像差異検出装置は、２枚以上の画像から、画像内の各領域における変化（または差異）を適切に検出するものであり、医用での異常検出や航空・衛星画像からの変化検出にも適用可能である。
【００２６】
＜第１実施形態＞
第１実施形態による画像差異検出装置を説明する。以下では動画像符号化を用いた衛星画像間の変化点検出を例にあげて説明する。図１は第１実施形態による画像差異検出装置１０の構成を示すブロック図である。図２は、図１に示す画像差異検出装置１０の処理動作を示すフローチャートである。図１に示す画像差異検出装置１０は、コンピュータ装置で構成される。画像差異検出装置１０は、変化検出対象サイズ入力部１、符号化パラメータ決定部２、符号化部３、変化検出部４、及び画像入力部５を備える。ここでは、画像差異検出装置１０が、画像Ａに対して画像Ｂではどの領域に変化があるか（２つの画像間において差異がどの空間領域にあるか）を検出する動作について説明する。なお以下では、画像差異検出装置１０が用いる符号化方式はＨＥＶＣであり、符号化パラメータには符号化サイズが含まれる場合を一例として説明する。
【００２７】
図１に示す画像差異検出装置１０の処理動作を説明する。変化検出対象サイズ入力部１は、外部から検出したい変化対象の大きさ（例えば、画素サイズ）を入力し、変化対象の大きさを符号化パラメータ決定部２に対して出力する（ステップＳ１）。変化対象の大きさは、一例として画素サイズが挙げられるが、画像上の範囲を指定する指標であれば画素サイズに限るものではない。続いて、符号化パラメータ決定部２は、入力した変化対象の大きさに応じて符号化パラメータ、例えば各符号化サイズを決定する（ステップＳ２）。各符号化サイズとは、ＭＰＥＧ等のＭＢ（マクロブロック）やＨＥＶＣのＣＵ（符号化ユニット）、ＰＵ（予測ユニット）、ＴＵ（変換ユニット）などのサイズがあげられる。ここでは、ＬＣＵ（Ｌａｒｇｅｓｔ  ＣＵ）サイズを符号化サイズの例として説明する。符号化パラメータ決定部２は、決定したＬＣＵサイズを符号化部３に対して出力する。
【００２８】
一方、画像入力部５は、外部から変化検出対象の少なくとも２枚の画像Ａ（第１の画像）、画像Ｂ（第２の画像）を入力する（ステップＳ３）。画像Ａ、Ｂの撮影日時は異なるが、画像Ａ、Ｂの撮影範囲の空間領域は同じである。画像Ａは、画像Ｂより前の日時に撮影された画像である。この画像Ａ、Ｂは、符号化部３に対して出力される。符号化部３は、入力したＬＣＵサイズで画像Ａおよび画像ＢのＨＥＶＣイントラ符号化を行う（ステップＳ４）。
【００２９】
図１に示す変化検出部４において、各ＬＣＵで変化があるか（差異があるか）否かを検出する処理について説明する。まず、変化検出部４は、検出対象画像における対象ＬＣＵの符号量と、対象ＬＣＵの四方（上下左右）に隣接するＬＣＵの符号量とを比較し、以下の条件Ｘ（（１）〜（８））のいずれかを満たすか否かを確認する（ステップＳ５）。
【００３０】
（１）ｍａｘ（Ｒ（Ｎ−１）／Ｒ（Ｎ），Ｒ（Ｎ）／Ｒ（Ｎ−１））＞Ｒ＿Ｔｈ１
（２）ｍａｘ（Ｒ（Ｎ＋１）／Ｒ（Ｎ），Ｒ（Ｎ）／Ｒ（Ｎ＋１））＞Ｒ＿Ｔｈ１
（３）ｍａｘ（Ｒ（Ｎ−ｘ）／Ｒ（Ｎ），Ｒ（Ｎ）／Ｒ（Ｎ−ｘ））＞Ｒ＿Ｔｈ１
（４）ｍａｘ（Ｒ（Ｎ＋ｘ）／Ｒ（Ｎ），Ｒ（Ｎ）／Ｒ（Ｎ＋ｘ））＞Ｒ＿Ｔｈ１
（５）ｍｉｎ（Ｒ（Ｎ−１）／Ｒ（Ｎ），Ｒ（Ｎ）／Ｒ（Ｎ−１））＜Ｒ＿Ｔｈ２
（６）ｍｉｎ（Ｒ（Ｎ＋１）／Ｒ（Ｎ），Ｒ（Ｎ）／Ｒ（Ｎ＋１））＜Ｒ＿Ｔｈ２
（７）ｍｉｎ（Ｒ（Ｎ−ｘ）／Ｒ（Ｎ），Ｒ（Ｎ）／Ｒ（Ｎ−ｘ））＜Ｒ＿Ｔｈ２
（８）ｍｉｎ（Ｒ（Ｎ＋ｘ）／Ｒ（Ｎ），Ｒ（Ｎ）／Ｒ（Ｎ＋ｘ））＜Ｒ＿Ｔｈ２
【００３１】
ここで、ＲはＬＣＵの符号量を意味する。ｍａｘ（Ｒ（Ｎ−１）、Ｒ（Ｎ））はＲ（Ｎ−１）とＲ（Ｎ）とのうち大きい値を意味する。ｍｉｎ（Ｒ（Ｎ−１），Ｒ（Ｎ））はＲ（Ｎ−１）とＲ（Ｎ）とのうち小さい値を意味する。Ｒ＿Ｔｈ１はＲ＿Ｔｈ１＞１を満たす閾値である。Ｒ＿Ｔｈ２は０＜Ｒ＿Ｔｈ２≦１を満たす閾値である。また、Ｎは画像ＢにおけるＮ番目のＬＣＵを意味し、Ｎ番目のＬＣＵは対象ＬＣＵである。Ｎ−１は対象ＬＣＵの左に隣接するＬＣＵを意味する。Ｎ＋１は対象ＬＣＵの右に隣接するＬＣＵを意味する。Ｎ−ｘは対象ＬＣＵの上に隣接するＬＣＵを意味する。Ｎ＋ｘは対象ＬＣＵの下に隣接するＬＣＵを意味する。このときの符号量比の算出対象は、画像Ｂである。
【００３２】
対象ＬＣＵが条件Ｘを全て満たさない（上下左右の隣接ＬＣＵの符号量と、対象ＬＣＵの符号量とが大きく異ならない）場合（ステップＳ５、ＮＯ）、変化検出部４は、対象ＬＣＵでは変化なしとして、変化検出部４内に対象ＬＣＵの検出結果を蓄積する（ステップＳ７）。
【００３３】
一方、条件Ｘのいずれかを満たす場合（ステップＳ５、ＹＥＳ）、変化検出部４は、画像Ｂにおける対象ＬＣＵと空間的に同位置の画像ＡにおけるＬＣＵの符号量と対象ＬＣＵの符号量との比較を行う。そして、変化検出部４は、画像Ｂの対象ＬＣＵと空間的に同位置の画像ＡのＬＣＵの符号量と画像Ｂの対象ＬＣＵの符号量とに基づいて、以下の条件Ｙ（（９）、（１０））のいずれかを満たすか否かを確認する（ステップＳ６）。
【００３４】
（９）ｍａｘ（Ｒ＿Ａ（Ｎ）／Ｒ＿Ｂ（Ｎ），Ｒ＿Ｂ（Ｎ）／Ｒ＿Ａ（Ｎ））＞Ｒ＿Ｔｈ３
（１０）ｍｉｎ（Ｒ＿Ａ（Ｎ）／Ｒ＿Ｂ（Ｎ），Ｒ＿Ｂ（Ｎ）／Ｒ＿Ａ（Ｎ））＜Ｒ＿Ｔｈ４
【００３５】
ここで、Ｒ＿Ｂは対象ＬＣＵの符号量を意味し、Ｒ＿Ａは対象ＬＣＵと空間的に同位置の画像ＡにおけるＬＣＵの符号量を意味する。Ｒ＿Ｔｈ３はＲ＿Ｔｈ３＞１を満たす閾値であり、Ｒ＿Ｔｈ４は０＜Ｒ＿Ｔｈ４≦１を満たす閾値である。
【００３６】
対象ＬＣＵが条件Ｙを全て満たさない（対象ＬＣＵの符号量と、対象ＬＣＵと空間的に同位置の画像ＡにおけるＬＣＵの符号量とが大きく異ならない）場合（ステップＳ６、ＮＯ）、変化検出部４は、対象ＬＣＵでは変化なしとして、変化検出部４内に対象ＬＣＵの検出結果を蓄積する（ステップＳ７）。
【００３７】
一方、条件Ｙのいずれかを満たす場合（ステップＳ６、ＹＥＳ）、変化検出部４は、対象ＬＣＵでは変化ありとして、検出結果を同様に蓄積する（ステップＳ７）。そして、画像Ｂにおける全てのＬＣＵに対して前述の処理（ステップＳ５，Ｓ６，Ｓ７）を繰り返して完了した時点で、変化検出部４は、対象画像（ここでは、画像Ｂ）の変化検出結果を出力し、処理が終了となる（ステップＳ８）。入力画像が、画像Ａと画像Ｂとの２枚でなく、画像Ａと２枚以上の画像Ｂとである場合は、画像Ｂそれぞれの検出結果が出力されることになる。
【００３８】
＜第２実施形態＞
次に、第２実施形態による画像差異検出装置を説明する。図３は第２実施形態による画像差異検出装置１０Ａの構成を示すブロック図である。図２において、図１に示す画像差異検出装置１０と同一の部分には同一の符号を付し、その説明を省略する。図２に示す画像差異検出装置１０Ａが図１に示す画像差異検出装置１０と異なる点は、変化検出対象種別入力部６が新たに設けられている点である。図４は、図３に示す画像差異検出装置１０Ａの処理動作を示すフローチャートである。図４において、図２に示す処理動作と同一の部分には同一の符号を付し、その説明を省略する。図４に示す処理動作が図２に示す処理動作と異なる点は、ステップＳ４ａ、ステップＳ９、及び、ステップＳ１０が設けられている点である。
【００３９】
第１実施形態では変化検出対象サイズのみが符号化パラメータ決定部２に入力されていた。第２実施形態では、変化検出対象サイズに加え、変化検出対象種別入力部６が外部から入力した変化検出対象とする種別（衛星画像等なら雲、建物、海など）である変化検出対象種別が符号化パラメータ決定部２に入力される（ステップＳ９）。そして、符号化パラメータ決定部２は、変化検出対象種別に応じて符号化時に用いる量子化マトリクス等を決定する（ステップＳ１０）。これを受けて、符号化部３は、符号化パラメータ決定部２により指定された符号化サイズおよび量子化マトリクスで画像Ａおよび画像Ｂを符号化する（ステップＳ４ａ）。この追加された処理動作は、検出の高精度化につながるものである。例えば、エッジが顕著な変化検出対象種別を指定した場合は、符号化パラメータ決定部２は、高周波成分が強調されるような量子化マトリクスを選定する。
【００４０】
＜第３実施形態＞
次に、第３実施形態による画像差異検出装置を説明する。図５は第３実施形態による画像差異検出装置１０Ｂの構成を示すブロック図である。図５において、図１に示す画像差異検出装置１０と同一の部分には同一の符号を付し、その説明を省略する。図５に示す画像差異検出装置１０Ｂが図１に示す画像差異検出装置１０と異なる点は、符号化結果調整部７が新たに設けられている点である。図６は、図５に示す画像差異検出装置１０の処理動作を示すフローチャートである。図６において、図２に示す処理動作と同一の部分には同一の符号を付し、その説明を省略する。図６に示す処理動作が図２に示す処理動作と異なる点は、ステップＳ１１が設けられている点である。
【００４１】
第３実施形態の処理は第１実施形態の処理と同様であるが、変化検出部４で用いる画像Ａ、ＢのＬＣＵ毎の符号量に符号化結果調整部７が重みをつける（調整する）点が異なる。具体的には符号化結果調整部７の処理が追加されている処理であるため、該当部分のみ説明をする。
【００４２】
画像Ａの特性と画像Ｂの特性と（精細感、色味など）がほぼ同じである場合、例えば画像Ａ、Ｂが一般的な３０ｆｐｓなどの映像である場合は、画像Ａと画像Ｂとの符号量は大きく異なることはない。しかし、撮影日時や、撮影環境（屋外だと太陽光の条件や大気の状況など）が異なる場合、同じ部位や地点を撮影した画像であっても、画像Ａの符号量と画像Ｂの符号量とが大きく異なる。これらの画像でＬＣＵ毎の変化検出を行ってもほぼ全てのＬＣＵが変化ありとして検出されてしまい、真の変化を検出することは困難となる。そこで、画像Ａ、Ｂにおいて符号量が異なる場合は、各ＬＣＵの符号量を調整する必要が生じる。本実施形態による画像差異検出装置１０Ｂはその調整を行うものである。
符号化部３は画像Ａ、Ｂの符号化を行い（ステップＳ４）、符号化結果調整部７は画像毎の符号量ｒ＿Ａ、ｒ＿Ｂを算出する。符号化結果調整部７は、符号量ｒ＿Ａ、ｒ＿Ｂを使い、画像Ｂ内ＬＣＵの符号量Ｒ＿Ｂを以下のように調整する（ステップＳ１１）。
Ｒ＿Ｂ’＝（ｒ＿Ａ／ｒ＿Ｂ）×Ｒ＿Ｂ×Ｋ
ここでＲ＿Ｂ’は、フロー後段の処理において、変化検出部４が用いる調整後の符号量となる。また、Ｋは任意で与えられる０より大きい定数である。画像全体やＬＣＵ毎の符号化情報に応じてＫを自動調整してもよい。変化検出部４は、第１実施形態に示した動作と同様に、符号量Ｒ＿Ａと調整された符号量Ｒ＿Ｂ’とを用い、変化検出処理を行う。
【００４３】
＜第４実施形態＞
次に、第４実施形態による画像差異検出装置を説明する。図７は第４実施形態による画像差異検出装置１０Ｃの構成を示すブロック図である。図７において、図５に示す画像差異検出装置１０Ｂと同一の部分には同一の符号を付し、その説明を省略する。図７に示す画像差異検出装置１０Ｃが図５に示す画像差異検出装置１０Ｂと異なる点は、符号化結果調整部７の出力に基づいて、変化検出結果を出力する点である。
【００４４】
第４実施形態による画像差異検出装置１０Ｃは、符号化結果調整部７からの出力を、変化検出部４を介さず変化検出結果として出力する。例えば、符号化結果調整部７は、画像Ａの符号量と画像Ｂの符号量とを比較する。以下（１１）、（１２）が満たされる場合は、画像Ａの特性と画像Ｂの特性とは非常に大きく異なり、全てのＬＣＵに変化がありと検出されるため、符号化結果調整部７は、画像全体として変化あり、または変化検出不能と判断する。符号化結果調整部７は、変化検出部４を介さず、判断結果を画像差異検出装置１０Ｃの出力結果として出力する。
（１１）ｍａｘ（ｒ＿Ａ／ｒ＿Ｂ）＞ｒ＿Ｔｈ１
（１２）ｍｉｎ（ｒ＿Ａ／ｒ＿Ｂ）＜ｒ＿Ｔｈ２
ここで、ｒ＿Ｔｈ１はｒ＿Ｔｈ１＞１を満たす閾値であり、ｒ＿Ｔｈ２は０＜ｒ＿Ｔｈ２≦１を満たす閾値である。
【００４５】
＜第５実施形態＞
次に、第５実施形態による画像差異検出装置を説明する。図８は第５実施形態による画像差異検出装置１０Ｄの構成を示すブロック図である。図８において、図１に示す画像差異検出装置１０と同一の部分には同一の符号を付し、その説明を省略する。図８に示す画像差異検出装置１０Ｄが図１に示す画像差異検出装置１０と異なる点は、画像解析部８が新たに設けられている点である。図９は、図８に示す画像差異検出装置１０Ｄの処理動作を示すフローチャートである。ここでは、画像Ａに対して、画像Ｂではどの領域に変化があるか（差異がどの空間領域にあるか）を画像差異検出装置１０Ｄが検出する処理動作について説明する。
【００４６】
図８に示す画像差異検出装置１０Ｄの処理動作を説明する。画像解析部８は、画像入力部５より画像Ａと画像Ｂとを入力し、変化検出対象サイズ入力部１より所望の変化検出対象サイズを入力し、画像Ａと画像Ｂとそれぞれに対して画像解析を実施する（ステップＳ２１）。画像解析部８は、変化検出対象サイズをもとに分割するブロックサイズを特定し、特定したブロックサイズで画像を分割して得られる検出ブロックそれぞれにおける各色成分の分散値（アクティビティ）を計算する。ここではアクティビティを例にするが、各色成分の平均値や最大値などでも適用可能である。そして、画像解析部８は、画像Ａと画像Ｂとにおける空間的に対応する検出ブロック同士でアクティビティの平均値の比率（差分でもよい）を算出し、算出した比率を符号化パラメータ決定部２に出力する。  
【００４７】
符号化パラメータ決定部２は、入力した変化検出対象サイズにより符号化サイズ等を決定するとともに、画像解析部８より出力されたアクティビティ比率より画像Ａと画像Ｂとに適用する量子化値等のパターンを決定する。後段の符号化にＨ．２６４／ＡＶＣやＨＥＶＣを利用する場合は、符号化パラメータ決定部２は、量子化マトリクスのパターンを決定してもよい。
【００４８】
また、各色成分の解析結果をもとに、符号化パラメータ決定部２は、量子化パラメータ値や量子化マトリクスを色成分ごとに決定してもよい。例えば、検出ブロックｎ（ｎ＝０、１、…、Ｎ）の各色成分のアクティビティをａＡ（ｎ）、ａＢ（ｎ）とすると、画像内のアクティビティ平均値の比率ａｒは、
ａｒ＝ΣａＡ（ｎ）／ΣａＢ（ｎ）
と示される。ここで、符号化パラメータ決定部２は、比率ａｒの値をもとに、画像Ａと画像Ｂとそれぞれに適用する量子化パラメータＱＰＡとＱＰＢとを決定する（ステップＳ２２）。あらかじめ定めた基準量子化パラメータを画像Ａまたは画像Ｂのいずれか一方にそのまま適用し、他方にａｒより算出される量子化パラメータ差分ｅを用いて変換した量子化パラメータを適用する。
【００４９】
例えば、基準量子化パラメータをＱＰＡとすると、変換した量子化パラメータＱＰＢは、
ＱＰＢ＝ＱＰＡ−ｅ
ｅ＝ｆ（ａｒ）
と示される。関数ｆはアクティビティ比率を量子化パラメータ差分に変換する関数である。関数ｆは、アクティビティの比率ａｒが大きい場合に、量子化パラメータ差分ｅが大きくなるような傾向を持つ関数形である。
【００５０】
次に、符号化部３は、符号化パラメータ決定部２により決定された量子化値等を用いて画像Ａと画像Ｂとそれぞれを符号化し、前述した符号化サイズの符号化ブロックごとの符号量を出力する（ステップＳ２３）。符号化の手段には、前述のＨ．２６４／ＡＶＣやＨＥＶＣだけではなく、ＭＰＥＧ−２やＪＰＥＧなど任意の動画像・静止画像の符号化手段が適用可能である。
【００５１】
変化検出部４は、各符号化ブロックの符号量を比較し、所定のしきい値との比較等により、各検出ブロックが変化領域か否かの情報を出力する（ステップＳ２４）。変化検出部４が出力する情報は変化領域か否か、の２判定だけではなく、判定不可なども加えた２つ以上の判定結果でもよい。
【００５２】
上述のように、本実施形態による画像差異検出装置１０Ｄは、アクティビティの比率または差分に基づいて、各画像の符号化時に適用する量子化パラメータを決定している。これにより、撮影環境や撮影日時の影響による変化検出の誤判定を抑制できる。
【００５３】
すなわち、同じ場所を撮影した画像でも、撮影環境や撮影日時の影響で、符号化の結果として出力される符号量が大きく異なる。例えば、晴天のときと曇りのときとでは、全体的な画像の色合いが異なる。このため、符号量をそのまま比較すると、撮影環境や撮影日時の影響により、変化検出の誤判定が生じる。本実施形態による画像差異検出装置１０Ｄは、検出ブロック内の各色成分の分散値をアクティビティとして計算し、このアクティビティの比率または差分に基づいて、量子化パラメータを決定している。量子化パラメータは、符号量を決定する主要素であり、この量子化パラメータをアクティビティの比率または差分に基づいて補正することで、画像差異検出装置１０Ｄは変化検出の誤判定を抑制できる。
【００５４】
＜第６実施形態＞
次に、第６実施形態による画像差異検出装置を説明する。図１０は第６実施形態による画像差異検出装置１０Ｅの構成を示すブロック図である。図１０において、図８に示す画像差異検出装置１０Ｄと同一の部分には同一の符号を付し、その説明を省略する。図１０に示す画像差異検出装置１０Ｅが図８に示す画像差異検出装置１０Ｄと異なる点は、画像解析部８の出力が変化検出部４に入力される点である。図１１は、図１０に示す画像差異検出装置１０Ｅの処理動作を示すフローチャートである。第５実施形態による画像差異検出装置１０Ｄは画像解析部８より出力された画像解析情報を符号化パラメータの決定に利用していた。第６実施形態による画像差異検出装置１０Ｅは符号化情報を用いて変化領域を検出する際にこの画像解析情報を用いる。
【００５５】
本実施形態の画像差異検出装置１０Ｅにおいて、変化検出部４は、画像解析部８より出力された画像Ａと画像Ｂとの検出ブロックごとのアクティビティ比率を入力し、画像Ａと画像Ｂとそれぞれの符号化情報の補正を行う。具体的には、変化検出部４は、画像Ａと画像Ｂとの検出ブロックごとの符号量の情報に対して、アクティビティ比率の情報をもとに符号量のスケーリングを行う。変化検出部４は、スケーリング後の符号量値を用いて所定のしきい値との比較を実施し、各検出ブロックが変化領域か否かを判定し、判定結果を出力する。検出ブロックの符号量に対するスケーリングの計算は、対応する検出ブロックのアクティビティ比率のみを用いて行ってもよいし、指定された周囲領域のアクティビティ比率を考慮して行ってもよい。
【００５６】
例えば、符号化部３が画像Ａ、Ｂの符号化を行い、画像毎の符号量ｒＡ、ｒＢを算出したと仮定する。第５実施形態の画像差異検出装置１０Ｄと同様に、画像解析部８がアクティビティ比率ａｒを求める場合、変化検出部４は、アクティビティ比率ａｒを使い、画像Ｂにおける検出ブロックごとの符号量ＲＢを以下のように調整する。
ＲＢ’＝ＲＢ×Ｋ
Ｋ＝ｇ（ａｒ）
ここで、関数ｇはアクティビティ比率を符号量のスケーリング係数に変換する関数である。関数ｇは、アクティビティの比率ａｒが大きい場合に、Ｋが大きくなるような傾向を持つ関数形である。
【００５７】
この処理動作を図１１を参照して説明する。まず、画像解析部８は、画像Ａおよび画像Ｂの変化抽出対象サイズの情報をもとに、画像Ａ、Ｂに対する画像解析を実施する（ステップＳ２１）。続いて、符号化部３は、指定された符号化パラメータにて画像Ａおよび画像Ｂを符号化する（ステップＳ２３）。
【００５８】
次に、変化検出部４は、画像解析結果をもとに、画像Ａおよび画像Ｂの符号化情報の補正を実施する（ステップＳ２５）。そして、変化検出部４は、補正後の画像Ａおよび画像Ｂの符号化情報をもとに画像Ａに対する画像Ｂの変化領域を特定する（ステップＳ２４）。
【００５９】
上述のように、本実施形態では、アクティビティの比率または差分に基づいて、各画像の符号量を調整している。これにより、撮影環境や撮影日時の影響による変化検出の誤判定を抑制できる。
【００６０】
すなわち、同じ場所を撮影した画像でも、撮影環境や撮影日時の影響で、符号化の結果として出力される符号量が大きく異なる。例えば、晴天のときと曇りのときとでは、全体的な画像の色合いが異なる。このため、符号量をそのまま比較すると、撮影環境や撮影日時の影響により、変化検出の誤判定が生じる。本実施形態による画像差異検出装置１０Ｅは、検出ブロック内の各色成分の分散値をアクティビティとして計算し、このアクティビティの比率または差分に基づいたスケーリングにより、符号量を調整している。これにより、画像差異検出装置１０Ｅは、撮影環境や撮影日時の影響による変化検出の誤判定を抑制できる。
【００６１】
＜第７実施形態＞
次に、第７実施形態による画像差異検出装置を説明する。図１２は第７実施形態による画像差異検出装置１０Ｆの構成を示すブロック図である。図１２において、図８に示す画像差異検出装置１０Ｄと同一の部分には同一の符号を付し、その説明を省略する。図１２に示す画像差異検出装置１０Ｆが図８に示す画像差異検出装置１０Ｄと異なる点は、符号化変化検出対象領域抽出処理部９が設けられている点である。図１３は、図１２に示す画像差異検出装置１０の処理動作を示すフローチャートである。第５実施形態による画像差異検出装置１０Ｄは画像解析部８より出力された画像解析情報を符号化パラメータの決定に利用していた。本実施形態による画像差異検出装置１０Ｆは、符号化パラメータ決定部２、符号化部３、変化検出部４による一連の処理（変化点抽出処理）を行う前に、符号化変化検出対象領域抽出処理部９によりその処理の対象とする画像領域を、画像解析情報を用いて特定する。
【００６２】
符号化変化検出対象領域抽出処理部９は、画像解析部８より出力した画像Ａと画像Ｂとの検出ブロックごとのアクティビティ比率を入力し、アクティビティ比率が所定のしきい値より大きい検出ブロックを後段の符号化による変化点抽出処理の対象外とする。検出ブロックを対象外にする方法としては、
（１）対象外の検出ブロックが省かれるよう画像の切り出しや分割を行う方法、
（２）対象外の検出ブロックを指定の色で塗りつぶす方法、
（３）画像は変化させず、後段の変化検出部にて変化領域と必ず判定されるよう、対象外の検出ブロックにのみ適用される符号化パラメータや変化検出のためのしきい値へと反映させる方法、などを用いることができる。
【００６３】
（１）の場合、分割や切り出しをされた画像ごとに符号化による変化点抽出処理を行う。検出ブロックサイズが小さい場合は、切り出しや分割により細かい画像が多数生成されてしまうことから、対象外と判定される領域がある一定面積以上存在する場合には、対象外でない領域も含めて対象外とみなして、分割数等を削減することも可能である。
（２）の場合、第５実施形態と同様の符号化による変化点抽出処理を行う。
（３）の場合、符号化パラメータ決定部２において、対象外領域においてのみ、符号量差分に対するしきい値を０もしくは非常に小さな値にする。
【００６４】
この処理動作を図１３を参照して説明する。まず、画像解析部８は、画像Ａおよび画像Ｂに対する変化抽出対象サイズの情報をもとに画像解析を実施する（ステップＳ２１）。続いて、符号化変化検出対象領域抽出処理部９は、画像解析結果をもとに、以降の符号化による変化抽出を行う領域を特定する（ステップＳ２６）。
【００６５】
次に、符号化部３は、特定後の画像領域を対象に、領域個々に設定した符号化パラメータにて画像Ａおよび画像Ｂを符号化（必要に応じて、画像の分割などを実施）する（ステップＳ２７）。そして、変化検出部４は、算出された画像Ａと画像Ｂとの符号化情報および画像解析結果による領域抽出の結果をもとに画像Ａに対する画像Ｂの変化領域を特定する（ステップＳ２８）。
【００６６】
上述のように、本実施形態の画像差異検出装置１０Ｆは、アクティビティの比率または差分に基づいて、変化点抽出処理の対象となる領域を決定する。これにより、画像差異検出装置１０Ｆは、撮影環境や撮影日時の影響による変化検出の誤判定を抑制できる。
【００６７】
すなわち、同じ場所を撮影した画像でも、撮影環境や撮影日時の影響で、符号化の結果として出力される符号量が大きく異なる。例えば、晴天のときと曇りのときとでは、全体的な画像の色合いが異なる。更に、雲で多く領域が覆われているような場合には、符号量の差異が大きくなる。そこで、本実施形態による画像差異検出装置１０Ｆは、アクティビティの比率または差分に基づいて、変化点抽出処理の対象となる領域を決定する。これにより、雲がかかっているような画像と雲がかかっていない画像のように、同じ地域を撮影していても、符号量が大きく変化するような画像は、変化検出の対象から除外できる。これにより、画像差異検出装置１０Ｆは、撮影環境や撮影日時の影響による変化検出の誤判定を抑制できる。
【００６８】
＜第８実施形態＞
次に、第８実施形態による画像差異検出装置を説明する。図１４は第８実施形態による画像差異検出装置１０Ｇの構成を示すブロック図である。この図において、図１に示す画像差異検出装置１０と同一の部分には同一の符号を付し、その説明を省略する。図１４に示す画像差異検出装置１０Ｇが図１に示す画像差異検出装置１０と異なる点は、開始点候補入力部１１が新たに設けられている点である。図１５は、図１４に示す画像差異検出装置１０Ｇの処理動作を示すフローチャートである。
【００６９】
本実施形態による画像差異検出装置１０Ｇでは、変化領域特定の際、符号化部３は、異なる開始点で画像Ａ、画像Ｂを符号化し、符号化処理を各画像に対して２回以上実行する。開始点候補入力部１１には、開始点の候補が入力される。変化検出部４は、複数回の符号化処理の結果を総合して、変化領域を特定する。符号化方式としては、ＨＥＶＣ、Ｈ．２６４、ＭＰＥＧ２、ＭＰＥＧ、ＪＰＥＧ等のような、１画面を複数のブロックに分割し、ブロック単位でのイントラ又はインター予測と直交符号により予測符号化し、このブロック毎の符号化データを可変長符号により圧縮符号化する符号化方式が用いられる。符号化部３は、少なくともブロック毎の符号化を行う。
【００７０】
一例としては、開始点を異ならせるずらし幅として、符号化の単位となるブロックサイズの縦横１／２の画素数が適用できる。例えば、８×８ブロックが符号化の一単位とすると、開始点をずれさない場合には位置（０，０）を開始点となり、縦横１／２ブロックサイズずらす場合には、位置（４，４）が開始点となる。この場合、開始点候補入力部１１には、開始点（０，０）と開始点（４，４）とが開始点候補として入力される。そして、符号化部３は、開始点（０，０）からと開始点（４，４）からとの２つのパターンで符号化を行う。そして、変化検出部４は、開始点（０，０）から符号化したときの符号化処理結果と、開始点（４，４）から符号化したときの符号化処理結果とを総合して変化領域を特定する。
【００７１】
次に、図１４に示す画像差異検出装置１０Ｇの処理動作を、図１５のフローチャートを用いて説明する。まず、開始点候補入力部１１は、画像Ａ及び画像Ｂの符号化の開始点（最左上の画素位置）の候補を２つ以上あらかじめ入力し、２つ以上の開始点の候補を設定しておく。
【００７２】
符号化部３は、設定された開始点の全候補について、それぞれを開始点とした画像Ａと画像Ｂとの符号化を実施し（ステップＳ１０１）、符号化により得られた符号化情報（符号量等）を記憶領域に保管する（ステップＳ１０２）。例えば、開始点候補が２つとし、それぞれ（０，０）と（４，４）とする。前者の場合、符号化部３は、画像の最左上の画素位置から符号化を開始する。後者の場合、符号化部３は、画像の最左上の画素位置から右下方向に縦横４画素ずつずらした位置から符号化を開始する。これにより、符号化ブロックの境界が４画素ずつずれることになる。符号化部３は、開始点の全候補について、符号化処理が完了するまで、処理を繰り返す（ステップＳ１０３）。
【００７３】
ステップＳ１０１からステップＳ１０３の処理により、符号化部３で、全ての開始点の候補で画像Ａと画像Ｂとの符号化が実施されたら、続いて、変化検出部４は、各画像の各開始点候補の符号化情報から、当該画像が変化領域か否かを判定する（ステップＳ１０４）。変化検出部４は、全ての画素位置について、変化領域か否かを判定する処理を繰り返す（ステップＳ１０５）。
【００７４】
ステップＳ１０４からステップＳ１０５の処理により、変化領域が検出されたら、最後に、変化検出部４は、画像全体の変化領域判定結果を出力する（ステップＳ１０６）。
【００７５】
ここで、図１５におけるステップＳ１０４では、変化検出部４は、画像Ａと画像Ｂとの各画素位置について、その画素が属する各開始点候補における符号化ブロックの符号化情報を特定し、その符号化情報の差分計算等により、当画素が変化ありか否かを判定する。この処理の詳細を、図１６及び図１７のフローチャートで説明する。
【００７６】
図１６は、ステップＳ１０４における変化領域か否かを判定する処理の一例を示すフローチャートである。ここでは、ある画素ｘについて、変化ありか否かを判定する処理を説明する。
【００７７】
まず、変化検出部４は、開始点候補ｎを用いた符号化にて得られた符号化ブロックのうち当該画素ｘを含む符号化ブロックを画像Ａ及び画像Ｂにおいて特定し（ステップＳ１１１）、特定した符号化ブロックが保有する符号化情報における差分を特定する（ステップＳ１１２）。ここでは符号量Ｒを符号化情報の一例にとり、画像Ａの符号量をＲＡ（ｎ，ｘ）、画像Ｂの符号量をＲＢ（ｎ，ｘ）と表現する。変化検出部４は、符号量ＲＡ（ｎ，ｘ）と符号量ＲＢ（ｎ，ｘ）との差分としきい値ｒＴＨを比較し、符号量差分がしきい値ｒＴＨ以上か否かを判定する（ステップＳ１１３）。すなわち、変化検出部４は、符号量差分としきい値ｒＴＨとを比較し、以下の条件式を満たすかどうかの判定を、すべての開始点候補について行う（ステップＳ１１４）。
｜ＲＡ（ｎ，ｘ）−ＲＢ（ｎ，ｘ）｜≧ｒＴＨ
【００７８】
ただし、図１６に示した処理の一例では、あるひとつの開始点候補にて上式が満たされた場合（ステップＳ１１３：Ｙｅｓ）、変化検出部４は、残りの開始点候補に対する判定をやめ、当該画素を変化ありとする判定結果を出力する（ステップＳ１１５）。すべての開始点候補について、上式が満たされない場合（ステップ１１３：Ｎｏ）、変化検出部４は、当該画素を変化なしとする判定結果を出力する（ステップＳ１１６）。
【００７９】
図１７は、ステップＳ１０４における変化領域か否かを判定する処理の他の例を示すフローチャートである。ここでは、ある画素ｘについて、変化ありか否かを判定する処理を説明する。
【００８０】
ステップＳ１２１からステップＳ１２２までの処理は、図１６に示したステップＳ１１１からステップＳ１１２までの処理と同じである。この例では、ステップＳ１２３で、変化検出部４は、特定した符号化ブロックに関する符号量差分としきい値ｒＴＨとを比較し、符号量差分がしきい値ｒＴＨ未満か否かを判定する。すなわち、変化検出部４は、符号量差分としきい値ｒＴＨとを比較し、以下の条件式を満たすかどうかの判定を、すべての開始点候補について行う。
【００８１】
｜ＲＡ（ｎ，ｘ）−ＲＢ（ｎ，ｘ）｜＜ｒＴＨ
【００８２】
ただし、図１７に示した処理の一例では、あるひとつの開始点候補にて上式が満たされた場合（ステップＳ１２３：Ｙｅｓ）、変化検出部４は、残りの開始点候補に対する判定をやめ、当該画素を変化なしとする判定結果を出力する（ステップＳ１２５）。すべての開始点候補について、上式が満たされない場合（ステップＳ１２３：Ｎｏ）、変化検出部４は、当該画素を変化ありとする判定結果を出力する（ステップＳ１２６）。
【００８３】
＜第９実施形態＞
次に、第９実施形態による画像差異検出装置を説明する。図１８は第９実施形態による画像差異検出装置１０Ｈの構成を示すブロック図である。図１８において、図１に示す画像差異検出装置１０と同一の部分には同一の符号を付し、その説明を省略する。図１８に示す画像差異検出装置１０Ｈが図１に示す画像差異検出装置１０と異なる点は、ブロック特性分類部１２と、領域ごと利用符号化情報・しきい値入力部１３とが新たに設けられている点である。図１９は、図１８に示す画像差異検出装置１０Ｈの処理動作を示すフローチャートである。
【００８４】
図１８におけるブロック特性分類部１２は、平坦領域、複雑領域等、符号化ブロック毎の特性をカテゴリに分類する。変化検出部４は、例えば、平坦領域では符号量絶対値のみを用いて変化を検出し、複雑領域では符号量絶対値と符号量分散値とを用いて変化を検出する。このように、ブロック特性分類部１２は、カテゴリに応じて制約条件を設け、変化検出に用いる符号化情報と判断基準とを変化させる。また、このとき、ブロック特性分類部１２は、平坦領域か複雑領域か等のカテゴリに応じて、しきい値を変化させる。また、分類されたカテゴリには、その度合いに応じた段階が設定されてもよい。例えば、平坦領域に関して、平坦１、平坦２、・・・のように平坦の度合いが強くなる毎に添え字を大きくしたカテゴリを設定する等である。
【００８５】
次に、図１８示す画像差異検出装置１０Ｈの処理動作を、図１９のフローチャートを用いて説明する。まず、符号化部３は、画像Ａおよび画像Ｂを符号化し（ステップＳ２０１）、各画像の符号化情報を記憶領域に保管する（ステップＳ２０２）。符号化方式としては、ＨＥＶＣ、Ｈ．２６４、ＭＰＥＧ２、ＭＰＥＧ、ＪＰＥＧ等のような、１画面を複数の符号化ブロックに分割し、符号化ブロック単位でのイントラ又はインター予測と直交符号により予測符号化し、符号化ブロック毎の符号化データを可変長符号により圧縮符号化する符号化方式が用いられる。符号化部３は、少なくとも符号化ブロック毎の符号化を行う。
【００８６】
次に、ブロック特性分類部１２は、その符号化情報を用いて、画像Ａおよび画像Ｂにおける各符号化ブロックを、その領域特性に応じたカテゴリに分類する。ここでは、符号化情報として符号量を例にとり、この符号量を用いて、符号化ブロックを平坦な領域か複雑な領域か２つのカテゴリに分類する場合について説明する。
【００８７】
まず、ブロック特性分類部１２は、画像Ａの各符号化ブロックの符号量を入力して、１次元データの非階層クラスタリングを実施する。クラスタリングの方法としては、例えば、ｋ−ｍｅａｎｓ法などがあげられる。ブロック特性分類部１２は、ｋ−ｍｅａｎｓ法などにより符号量データ群を２つのクラスタに分類し、重心の符号量値が小さいクラスタ（小クラスタ）と大きいクラスタ（大クラスタ）とに分別する（ステップＳ２０３）。つまり、画像Ａのある符号化ブロックｙの符号量をＲＡ（ｙ）とし、ｋ−ｍｅａｎｓ法で算出された小クラスタの重心符号量をｒＡ１とし、大クラスタの重心符号量をｒＡ２としたとき、ブロック特性分類部１２は、下記のようにクラスタリングを行う。
【００８８】
Ｉｆ  （ＲＡ（ｙ）≦（ｒＡ１＋ｒＡ２）／２） ：符号化ブロックｙは小クラスタである
Ｅｌｓｅ ：符号化ブロックｙは大クラスタである
【００８９】
ブロック特性分類部１２は、上述と同様の処理を、画像Ｂに対しても同様に実施する。
【００９０】
次に、ブロック特性分類部１２は、符号化ブロックｙが平坦領域と複雑領域とのいずれに分類されるかを判定する（ステップＳ２０４）。例えば、ブロック特性分類部１２は、画像Ａと画像Ｂとのクラスタ種類を確認し、以下のようにして、符号化ブロックｙを平坦領域と複雑領域とに分類する。
【００９１】
Ｉｆ  （ＲＡ（ｙ）≦（ｒＡ１＋ｒＡ２）／２）  
＆＆  （ＲＢ（ｙ）≦（ｒＢ１＋ｒＢ２）／２）  ：  符号化ブロックｙは平坦領域である
Ｅｌｓｅ  ：  符号化ブロックｙは複雑領域である
【００９２】
なお、平坦領域か複雑領域かを、優先度の高低をつけて、いずれかに分類する方法としては、上記の方法以外に、画像Ａや画像Ｂがもつ外部情報（撮影日時や画質など）をキーにして、画像Ａと画像Ｂとに優先度を定め、クラスタ種類が異なる場合には、その優先度をもとに領域特性を定める方法が適用できる。
【００９３】
ステップＳ２０４の処理は、全ての符号化ブロックについて完了するまで、繰り返して行われる（ステップＳ２０５）。
【００９４】
続いて、変化検出部４は、符号化ブロックの特性が平坦特性か否かを判定する。（ステップＳ２０６）。そして、変化検出部４は、当該符号化ブロックの特性が平坦特性である場合（ステップＳ２０６：Ｙｅｓ）、平坦特性用の変化判定を実施し（ステップＳ２０７）、当該符号化ブロックの特性が複雑特性である場合（ステップＳ２０６：Ｎｏ）、複雑特性用の変化判定を実施する（ステップＳ２０８）。
【００９５】
ステップＳ２０６からステップＳ２０８の処理は、全ての符号化ブロックについて完了するまで、繰り返して行われる（ステップＳ２０９）。
【００９６】
ここで、ステップＳ２０７の平坦特性用の変化判定、および、ステップＳ２０８の複雑特性用の変化判定は、以下のようにして行うことができる。
【００９７】
（１）平坦領域では符号量値および符号量分散値を指標として用い、複雑領域では符号量値のみを用いる場合
平坦領域における変化判定のしきい値をそれぞれ、符号量しきい値ｒ’ＴＨ、符号量分散しきい値ｖＴＨ、複雑領域における変化判定のしきい値を、符号量しきい値ｒ’’ＴＨとする。このとき、符号化ブロックｙが平坦領域である場合、以下の判定式により、変化検出部４は符号化ブロックｙに変化があるか否かを判定する。
【００９８】
Ｉｆ  （｜ＲＡ（ｙ）−ＲＢ（ｙ）｜≦ｒ’ＴＨ）
＆＆  （ＶＡ（ｙ）≦ｖＴＨ）
＆＆  （ＶＢ（ｙ）≦ｖＴＨ）  ：  符号化ブロックｙに変化なし
Ｅｌｓｅ ：  符号化ブロックｙに変化あり
【００９９】
また、符号化ブロックｙが複雑領域である場合、以下の判定式により、変化検出部４は符号化ブロックｙに変化があるか否かを判定する。
【０１００】
Ｉｆ  （｜ＲＡ（ｙ）−ＲＢ（ｙ）｜≦ｒ’’ＴＨ）  ：  符号化ブロックｙに変化なし
Ｅｌｓｅ  ：  符号化ブロックｙに変化あり
【０１０１】
ここで、ＶＡ（ｙ）、ＶＢ（ｙ）は、画像Ａおよび画像Ｂの符号化ブロックｙにおける周辺の符号化ブロックと符号化ブロックｙとの符号量の分散値である。例えば、符号化ブロックｙを中心とした上下左右ｐ個の隣接符号化ブロック（（２ｐ＋１）＊（２ｐ＋１）個の符号化ブロック）の符号量値から、変化検出部４は分散値を計算する。ｐはあらかじめ外部より指定される。
【０１０２】
（２）平坦領域、複雑領域ともに符号量値のみを指標として用いる場合
平坦領域における変化判定のしきい値を符号量しきい値ｒ’ＴＨ、複雑領域における変化判定のしきい値を、符号量しきい値ｒ’’ＴＨとする。このとき、符号化ブロックｙが平坦領域である場合は、以下の判定式により、変化検出部４は符号化ブロックに変化があるか否かを判定する。
【０１０３】
Ｉｆ  （｜ＲＡ（ｙ）−ＲＢ（ｙ）｜≦ｒ’ＴＨ）  ：  符号化ブロックｙに変化なし
Ｅｌｓｅ  ：  符号化ブロックｙに変化あり
【０１０４】
符号化ブロックｙが複雑領域である場合は、以下の判定式により、変化検出部４は符号化ブロックに変化があるか否かを判定する。
【０１０５】
Ｉｆ  （｜ＲＡ（ｙ）−ＲＢ（ｙ）｜≦ｒ’’ＴＨ）  ：  符号化ブロックｙに変化なし
Ｅｌｓｅ  ：  符号化ブロックｙに変化あり
【０１０６】
一例として、当該符号化ブロックが平坦領域であれば、変化検出部４は、あらかじめ定められた平坦領域向けの符号化情報を変化判定指標とそれに対応するしきい値とを用いて変化判定を行う。また、当該符号化ブロックが複雑領域であれば、変化検出部４は、変化判定指標としきい値とのいずれか、もしくは両方を、平坦領域とは異なるものを用いる。
【０１０７】
ステップＳ２０６からステップＳ２０９の処理により、変化領域が検出されたら、最後に、変化検出部４は、画像全体の変化領域判定結果を出力する（ステップＳ２１０）。
【０１０８】
＜第１０実施形態＞
次に、第１０実施形態による画像差異検出装置を説明する。図２０は第１０実施形態による画像差異検出装置１１０の構成を示すブロック図である。図２０に示すように、第１１実施形態による画像差異検出装置１１０は、画像入力部１０１と、符号化部１０２と、ブロックごと符号化情報特定部１０３と、複雑さ指標値算出部１０４と、変化検出部１０５とを備える。図２１は、図２０に示す画像差異検出装置１１０の処理動作を示すフローチャートである。
【０１０９】
本実施形態による画像差異検出装置１１０では、画像入力部１０１が非圧縮の画像Ａと画像Ｂとを入力し、符号化部１０２が符号化を実施する。画像Ａ、Ｂは、撮影日時は異なる画像であるが、撮影範囲の空間領域が同じ画像である。そして、ブロックごと符号化情報特定部１０３は、符号化途中で算出される各種符号化情報を特定する。複雑さ指標値算出部１０４は、符号化情報から複雑さ指標値を算出する。そして、変化検出部１０５は、この複雑さ指標値を用いて、変化領域を特定する。
【０１１０】
次に、図２０に示す画像差異検出装置１１０の処理動作を、図２１のフローチャートを用いて説明する。まず、画像入力部１０１は、画像Ａと画像Ｂとを入力し、符号化部１０２は、画像Ａおよび画像Ｂを符号化する（ステップＳ３０１）。符号化にはＨＥＶＣ、Ｈ．２６４、ＪＰＥＧなどが用いられる。次に、ブロックごと符号化情報特定部１０３は、符号化途中で算出される各種符号化情報を記録し、ログとしてメモリ領域に保管する。複雑さ指標値を算出するための符号化情報としては、以下のものを用いることができる。
【０１１１】
・符号化ブロック（ＣＵ、ＣＴＵ）における符号量
・符号化ブロック（ＣＵ、ＣＴＵ）における２値化ビット値幅（ビン）
・符号化ブロック（ＣＵ、ＣＴＵ）におけるＤＣＴ係数
【０１１２】
続いて、複雑さ指標値算出部１０４は、これら符号化情報から、複雑さ指標値を算出する（ステップＳ３０２）。複雑さ指標値としては、以下の様な符号化情報、もしくはこれらの加重和などが指標値として適用できる。加重和により得られる指標値は、例えば、符号化情報に含まれるパラメータごとに画像Ａのパラメータと画像Ｂのパラメータとの差分の二乗値又は絶対値を加算した値、又は、符号化情報に含まれるパラメータごとに画像Ａのパラメータと画像Ｂのパラメータとの差分の二乗値又は絶対値を重み付け加算した値である。
【０１１３】
・符号化ブロック（ＣＵ、ＣＴＵ）における総符号量
・符号化ブロック（ＣＵ、ＣＴＵ）におけるＤＣＴ係数に関与する符号量
・符号化ブロック（ＣＵ、ＣＴＵ）におけるＹ／Ｃｂ／Ｃｒそれぞれの符号量
・符号化ブロック（ＣＵ、ＣＴＵ）における２値化ビット値幅（ビン）
・符号化ブロック（ＣＵ、ＣＴＵ）におけるＤＣＴ係数値または有意係数個数
【０１１４】
次に、変化検出部１０５は、この複雑さ指標値を用いて、変化検出を行う（ステップＳ３０３）。なお、変化の検出方法は、前述の実施形態と同様である。
【０１１５】
このように、前述の第１実施形態から第９実施形態の画像差異検出装置は、変化領域を特定するのに、符号量を比較している。これに対して、第１０実施形態の画像差異検出装置は、複雑さ指標値を比較している。
【０１１６】
＜第１１実施形態＞
次に、第１１実施形態による画像差異検出装置を説明する。図２２は第１１実施形態よる画像差異検出装置２１０の構成を示すブロック図である。図２３は、図２２に示す画像差異検出装置２１０の処理動作を示すフローチャートである。図２２に示す画像差異検出装置２１０は、コンピュータ装置で構成されてもよい。画像差異検出装置２１０は、符号化データ入力部２０１、可変長符号復号部２０２、ブロックごと符号化情報特定部２０３、複雑さ指標値算出部２０４、共通領域特定部２０５、複雑さ指標値補正部２０６、変化検出部２０７、位置情報メタデータ入力部２０８、及び、線形和係数入力部２０９を備える。
【０１１７】
前述までの第１から第１０までの各実施形態では、画像Ａおよび画像Ｂとして、非圧縮の画像を入力している。これに対して、本実施形態では、符号化データ入力部２０１は、符号化された画像Ａおよび画像Ｂの符号化データを入力する。画像Ａ、Ｂは、撮影日時が異なる画像であるが、撮影範囲の空間領域が同じ画像である。符号化方式としては、ＨＥＶＣ、Ｈ．２６４、ＭＰＥＧ２、ＭＰＥＧ、ＪＰＥＧ等が用いられる。可変長符号復号部２０２は、符号化された画像Ａ及び画像Ｂの符号化データに対して可変長符号の復号を行う。ブロックごと符号化情報特定部２０３は、復号途中で算出される各種符号化情報を特定する。複雑さ指標値算出部２０４は、復号途中で算出される各種符号化情報から、各符号化ブロックのテクスチャの複雑さを模擬する複雑さ指標値を算出する。すなわち、可変
長符号復号部２０２は、可変長符号の復号を完了することにより、符号化ブロックそれぞれの符号化情報を取得する取得部として動作する。可変長符号の復号は、符号化データから画像を復号する過程において行われる処理の一部である。
【０１１８】
位置情報メタデータ入力部２０８は、外部情報として、グローバルな位置情報をメタデータとして入力する。共通領域特定部２０５は、画像Ａと画像Ｂとにおける同一の被写体を共通領域として特定する。複雑さ指標値補正部２０６は、画像Ａおよび画像Ｂの共通領域上での複雑さ指標値を算出する。複雑さ指標値補正部２０６は、共通領域上での複雑さ指標値を算出するための線形和演算を行うフィルタを有している。線形和係数入力部２０９は、線形和演算の係数を入力する。同一の被写体とは、位置情報が保存される画像ファイル間において、画像ファイル間の位置情報により同一位置と判断される領域である。すなわち、同一の被写体とは、ほぼ同一空間と判断される領域である。また、画像ファイルに保存される位置情報は、別途取得されてもよい。また、画像ファイルは、衛星画像、航空写真、市販されるデジタルカメラや携帯電話機等で撮影される画像であってもよい。
【０１１９】
次に、図２２に示す画像差異検出装置１０の処理動作を、図２３のフローチャートを用いて説明する。
【０１２０】
まず、符号化データ入力部２０１は、符号化された画像Ａおよび画像Ｂの符号化データを入力する。可変長符号復号部２０２は、それらの符号化データに対して可変長符号の復号を行う（ステップＳ４０１）。次に、ブロックごと符号化情報特定部２０３は、各符号化ブロックが保有する符号化情報を特定し（ステップＳ４０２）、複雑さ指標値算出部２０４は、特定された符号化情報をもとにして複雑さ指標値を算出する（ステップＳ４０３）。
【０１２１】
入力する画像Ａおよび画像Ｂの符号化データの符号化方式としては、ＨＥＶＣ、Ｈ．２６４、ＪＰＥＧなどが考えらえる。これらの符号化方式は、１画面を複数の符号化ブロックに分割し、符号化ブロック単位のイントラ又はインター予測とＤＣＴ等の直交符号により予測符号化を行い、符号化ブロック毎の符号を可変長符号により圧縮符号化している。可変長符号として、ＨＥＶＣやＨ．２６４では、二値算術符号が用いられる。ＪＰＥＧでは、ハフマンコードが用いられる。したがって、上述のステップＳ４０１からステップＳ４０３の処理は、符号化方式に応じて、具体的には、以下のような処理となる。
【０１２２】
（１）ＨＥＶＣまたはＨ．２６４にて符号化された符号化データが入力される場合
図２４は、ＨＥＶＣまたはＨ．２６４にて符号化された符号化データから複雑さ指標値を得る処理動作を示すフローチャートである。まず、可変長符号復号部２０２は、算術符号（ＣＡＢＡＣ）を復号し、２値信号を取得する（ステップＳ５０１）。そして、可変長符号復号部２０２は、各シンタクス要素を取得し、この２値信号を多値信号に変換する（ステップＳ５０２）。次に、ブロックごと符号化情報特定部２０３は、各符号化ブロック（ＣＵ、ＣＴＵ）が保有する符号化情報を特定する（ステップＳ５０３）。この符号化情報としては、例えば、以下の様なものが適用できる。
【０１２３】
・符号化ブロック（ＣＵ、ＣＴＵ）における符号量
・符号化ブロック（ＣＵ、ＣＴＵ）における２値化ビット値幅（ビン）
・符号化ブロック（ＣＵ、ＣＴＵ）におけるＤＣＴ係数
【０１２４】
続いて、複雑さ指標値算出部２０４は、これらの符号化情報をもとに、被写体を構成する符号化ブロックの構造を反映する指標値である複雑さ指標値を算出する（ステップＳ５０４）。複雑さ指標値としては、例えば、以下の様なもの、もしくはこれらの加重和などが適用できる。加重和により得られる指標値は、例えば、符号化情報に含まれるパラメータごとに画像Ａのパラメータと画像Ｂのパラメータとの差分の二乗値又は絶対値を加算した値、又は、符号化情報に含まれるパラメータごとに画像Ａのパラメータと画像Ｂのパラメータとの差分の二乗値又は絶対値を重み付け加算した値である。
【０１２５】
・符号化ブロック（ＣＵ、ＣＴＵ）における総符号量
・符号化ブロック（ＣＵ、ＣＴＵ）におけるＤＣＴ係数に関与する符号量
・符号化ブロック（ＣＵ、ＣＴＵ）におけるＹ／Ｃｂ／Ｃｒそれぞれの符号量
・符号化ブロック（ＣＵ、ＣＴＵ）における２値化ビット値幅（ビン）
・符号化ブロック（ＣＵ、ＣＴＵ）におけるＤＣＴ係数値または有意係数個数
【０１２６】
（２）ＪＰＥＧにて符号化された符号化データが入力される場合
図２５は、ＪＰＥＧにて符号化された符号化データから複雑さ指標値を得る処理動作を示すフローチャートである。まず、可変長符号復号部２０２は、ハフマン符号を復号し、多値信号を出力する（ステップＳ６０１）。次に、ブロックごと符号化情報特定部２０３は、各符号化ブロック（８ｘ８）のデータ境界を特定する（ステップＳ６０２）。続いて、複雑さ指標値算出部２０４は、各符号化ブロックが保持しているＤＣＴの係数値、もしくは、有意係数の個数などを複雑さ指標として出力する（ステップＳ６０３）。
【０１２７】
図２３に示す処理において、ステップＳ４０１からステップＳ４０３で、複雑さ指標値が算出されたら、共通領域特定部２０５は、位置情報メタデータ入力部２０８から位置情報メタデータを入力し、両画像間にて位置的に対応する画素位置を位置情報メタデータを用いて特定し、両画像間の共通領域を導出する（ステップＳ４０４）。位置情報メタデータとしては、対象の画像が衛星画像などである場合には、その画像中の各座標が緯度経度もしくはＵＴＭ座標系等で示される地理上の位置情報が用いられる。
【０１２８】
続いて、複雑さ指標値補正部２０６は、対応する画素位置を含む符号化ブロックを各画像にて特定し、符号化ブロック間の位置差分（ずれ）の分だけ、複雑さ指標値を補正する（ステップＳ４０５）。
【０１２９】
図２６は、図２３に示すステップＳ４０５での共通領域における複雑さ指標値の補正処理を示すフローチャートである。また、図２７は、共通領域における複雑さ指標値の補正処理の一例を示す図である。
【０１３０】
まず、複雑さ指標値補正部２０６は、画像Ａおよび画像Ｂの各符号化ブロックにおけるテクスチャの複雑さを示す複雑さ指標値を入力する（ステップＳ７０１）。図２７に示すように、画像Ａの符号化ブロックをＧＡ（ｓ）としたとき、その複雑さ指標値をＣＡ（ＧＡ（ｓ））とする。ここでｓは画像左上からラスタスキャン順に符号化ブロックの順列を定めた際の番号を示す。同様に画像Ｂの複雑さ指標値をＣＢ（ＧＢ（ｓ））と定める。
【０１３１】
次に、複雑さ指標値補正部２０６は、外部情報により特定した各画像間の共通領域（地理的に重複する領域）を特定する（ステップＳ７０２）。共通領域は画素の領域（範囲）として示される。画像Ａと画像Ｂとの共通領域は、画像Ａでは、ｋ１≦ｘ≦ｋ２、ｌ１≦ｙ≦ｌ２、画像Ｂでは、ｍ１≦ｘ≦ｍ２、ｎ１≦ｙ≦ｎ２とする。ここで、（ｘ，ｙ）は画素位置を示す値である。
【０１３２】
続いて、複雑さ指標値補正部２０６は、共通領域のみを画像として符号化すると想定した際の符号化ブロック位置と、画像Ａおよび画像Ｂのそれぞれの符号化データ上で実際に符号化されている符号化ブロック位置との位置差分（差分画素数）を特定する（ステップＳ７０３）。ここでは、図２７に示すように、画像Ａの符号化ブロック位置の差分画素数を（ＰＡ，ＱＡ）とする。
【０１３３】
最後に、複雑さ指標値補正部２０６は、フィルタ演算により共通領域の各符号化ブロック位置における複雑さ指標値を算出し、複雑さ指標値の補正を行う（ステップＳ７０４）。補正後の画像Ａの共通領域における符号化ブロックｔの複雑さ指標値をＣＡ’（ｔ）とすると、ＣＡ’（ｔ）は、以下のように、記述される。
【０１３４】
ＣＡ’（ｔ）＝Σα（ｓ，ＰＡ，ＱＡ）＊ＣＡ（ＧＡ（ｓ））
【０１３５】
ここで、αは、差分画素数に応じた係数である。つまり、この例では、補正後の複雑さ指標値は、もともとの複雑さ指標値群の線形和で算出される。複雑さ指標値補正部２０６は、画像Ｂについても同様の手順によって、補正後の複雑さ指標値を算出する。なお、線形和係数αは、画素位置に応じて０≦α（ｓ，ＰＡ，ＱＡ）≦１を変動するように定め、かつ、Σα＝１となるように外部より与える。フィルタ演算の対象となる符号化ブロックは、例えば、符号化データにおける符号化ブロックのうち、符号化ブロックｔと重複する領域がある符号化ブロックである。
【０１３６】
なお、位置差分に基づいて画像Ａおよび画像Ｂの複雑さ指標値の補正を行う際の線形演算の係数は、前述の位置差分に対応する係数として線形和係数入力部２０９へ入力され、符号化ブロックに起因する位置ずれを補正するための係数として設定される。また、複雑さ指標値補正部２０６は、線形演算を行う係数として予め複数記憶しておき、位置差分に対応する係数をその中から選定してもよい。また、複雑さ指標値ＣＡからＣＡ’へのフィルタ演算には、非線形演算が用いられてもよい。フィルタ係数は、予め画像差異検出装置２１０に記憶されていてもよい。この場合、複雑さ指標値補正部２０６は、このフィルタ係数を用いてフィルタ演算を行ってもよい。このとき、異なるフィルタ演算（例えば、非線形演算）が用いられてもよい。
【０１３７】
図２７に示すように、共通領域における符号化ブロックｔの複雑さ指標値ＣＡ’（ｔ）は、画像Ａの符号化ブロックＧＡ（ｓ）の複雑さ指標値ＣＡ（ＧＡ（ｓ））に対して、差分画素数（ＰＡ，ＱＡ）を基に線形演算を行うことで求められる。同様に、共通領域における符号化ブロックｔの複雑さ指標値ＣＢ’（ｔ）は、画像Ｂの符号化ブロックＧＢ（ｓ）の複雑さ指標値ＣＢ（ＧＢ（ｓ））に対して、差分画素数（ＰＢ，ＱＢ）を基に線形演算を行うことで求められる。
【０１３８】
図２３において、ステップＳ４０５で、画像Ａおよび画像Ｂの複雑さ指標値を補正したら、変化検出部２０７は、補正後の複雑さ指標値を用いて、画像Ａおよび画像Ｂとで同一となる共通領域を特定し、共通領域における変化検出を行う（ステップＳ４０６）。補正後の複雑さ指標値を用いて変化検出を行う処理は、他の実施形態と同様である。
【０１３９】
なお、各実施形態では、撮影した日時が異なる複数の画像における変化を検出する場合を説明したが、各実施形態による画像差異検出装置は、撮影した場所、日時に限定されず複数の画像における差異を検出する用途にも適用可能である。
【０１４０】
以上説明したように、衛星画像、航空画像、医療画像等において、撮影日時が異なる同一空間領域の画像同士を比較することにより、変化領域（差異がある位置）を高速に検出することができるようになる。
【０１４１】
第１１実施形態による画像差異検出装置は、符号化された画像Ａおよび画像Ｂから画像の変化を検出する際に、符号データの復号を完了させず、復号途中で得られる符号化情報から複雑さ指標値を求めて、変化領域を検出している。すなわち、一般的な画像の符号化方式では、１画面を複数の符号化ブロックに分割し、符号化ブロック単位のイントラ又はインター予測と直交符号により予測符号化をする段階と、この符号化ブロック毎の符号化データを可変長符号により圧縮符号化する段階とを有している。したがって、画像を復号するには、可変長符号の復号段階と、直交符号の復号段階とが必要であり、演算量が膨大になる。第１１実施形態による画像差異検出装置は、可変長符号の復号だけを行った段階で得られる符号化情報から複雑さ指標値を求めて、変化領域を検出している。このため、符号化データとして与えられる画像間における差異を検出する際の演算量が大幅に削減され、高速動作が可能になる。
【０１４２】
なお、各実施形態において、符号化ブロックに基づいた差異の検出に代えて、検出ブロックに基づいた差異の検出を行ってもよい。検出ブロックのサイズと符号化ブロックのサイズとが同じ場合、画像差異検出装置は同じ処理を行う。検出ブロックのサイズと符号化ブロックのサイズとが異なる場合、画像差異検出装置は、符号化ブロックの符号化量又は指標値に基づいて、検出ブロック符号量又は指標値を算出する。例えば、符号化ブロックのサイズが検出ブロックのサイズよりも小さい場合、画像差異検出装置は、検出ブロックに含まれる符号化ブロックそれぞれの符号量又は指標値を加算することにより、検出ブロックの符号量又は指標値を算出する。符号化ブロックのサイズが検出ブロックのサイズより大きい場合、画像差異検出装置は、検出ブロックを含む符号化ブロックの符号量又は指標値に、符号化ブロックの面積に対する検出ブロックの面積の比を乗算し、乗算により得られる値を検出ブロックの符号量又は指標値とする。すなわち、符号化ブロックのサイズと検出ブロックのサイズとが異なる場合、画像差異検出装置は、符号化ブロックと検出ブロックとサイズの差又は比に応じた線形演算を符号化ブロックの符号量又は指標値に対して行うことにより、検出ブロックの符号量又は指標値を得る。サイズの差に応じて符号化ブロックの符号量又は指標値を算出する処理は、変化検出部、ブロック特性分類部、又は、複雑さ指標値算出部により行われる。
【０１４３】
また、ＤＣＴ係数が指標値に用いられる場合における指標値の加算は、例えば、符号化ブロックごとに高周波数側のＤＣＴ係数を所定数抽出し、抽出したＤＣＴ係数の総和を算出する演算を意味する。なお、高周波数側のＤＣＴ係数に代えて低周波数側のＤＣＴ係数を用いてもよい。また、差異の有無を顕著に表す周波数成分をＤＣＴ係数から抽出するようにしてもよいし、総和の算出に代えて当該ＤＣＴ係数に対する重みを他のＤＣＴ係数より大きくした重み付け加算を行ってもよい。また、ＤＣＴ係数が指標値に用いられる場合、ＤＣＴ係数の絶対値が予め定められた閾値を超えるＤＣＴ係数の個数に基づいて差異の有無を検出してもよい。すなわち、閾値を超えるパワーを有する有意なＤＣＴ係数の個数により差異の有無を判定してもよい。
【０１４４】
なお、各実施形態で用いられる符号化方式としてはＨＥＶＣ以外、例えばＡＶＣ、ＭＰＥＧ、ＭＰＥＧ−２、ＪＰＥＧ等を用いることも可能である。さらにマクロブロックや符号化サイズより大きいまたは小さいサイズで分割した領域毎に差異の検出を実施することも可能であり、同様の効果を奏する。また、符号化情報としては、符号量の代わりに、各領域のアクティビティ平均値・最大値、ＤＣＴで得られるＤＣ係数および、各符号化モード等、符号化の結果で得られる各情報を用いることも可能であり、同様の効果を奏する。
【０１４５】
また、前述した説明では撮影した日時が異なる複数の画像の変化を検出するものとなっているが、前述した各実施形態による画像差異検出装置は、撮影した場所、日時に限定されず複数の画像間の差異を検出する用途にも適用可能である。
また、前述した説明においては、矩形の小領域（符号化ブロック）に分割する符号化方式を用いて、差異がある小領域を特定する画像差異検出装置の構成を説明したが、領域の特定は必須ではなく、画像差異検出装置は、単に２つの画像に差異があることを符号量の違いによって検出するようにしてもよい。
また、各実施形態による画像差異検出装置は、符号量を比較することで変化の判断を行う。比較される画像が同じ地点から撮影されたものであっても画像の精細感などが異なると全体として情報量が大きく変わってしまうため「全体の符号量を近づける処理」を行うことを前提として行っている。
また、各実施形態による画像差異検出装置は、符号化を行う規格によらず使用することができる。符号化方式が「可変ビットレート」である場合に効果が顕著に表れる。
また、前述した各実施形態による画像差異検出装置は、他の技術と合わせて使用することも可能である。
また、前述した各実施形態における変化の判定は、変化の判定のみならず、通常の符号化とセットで使うことも可能である。したがって、前述の画像差異検出装置は、他の目的を達成するものであってもよい。
【０１４６】
以上説明したように、衛星画像、航空画像、医療画像等において、撮影日時が異なる同一空間領域の画像同士を比較することにより、変化領域（差異がある位置）を高速に検出することができるようになる。
【０１４７】
画像間の変化を検出する従来技術と比較して、符号量を指標として用いる各実施形態の画像差異検出装置は、高速に差異を検出し、かつ検出に要する演算量を削減することができる。従来技術と同様の環境で、従来技術と各実施形態による画像差異検出装置とを比較をした際、結果が出力されるまで要する時間が大幅に短縮される結果が得られた。
【０１４８】
上述した実施形態における画像差異検出装置をコンピュータで実現するようにしてもよい。その場合、この機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現してもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでもよい。また上記プログラムは、前述した機能の一部を実現するためのものであってもよく、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであってもよく、ＦＰＧＡ（Field Programmable Gate Array）等のプログラマブルロジックデバイスを用いて実現されるものであってもよい。
以上、この発明の実施形態について図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。例えば、複数の実施形態の構成を組み合わせることによって一つの実施形態として実現されてもよい。
【産業上の利用可能性】
【０１４９】
本発明は、複数の画像の差異をより速く検出することが不可欠な用途に適用できる。
【符号の説明】
【０１５０】
１０，１０Ａ，１０Ｂ，１０Ｃ，１０Ｄ，１０Ｅ，１０Ｆ，１０Ｇ，１０Ｈ，１１０，２１０  画像差異検出装置
１  変化検出対象サイズ入力部
２  符号化パラメータ決定部
３，１０２  符号化部
４，１０５，２０７  変化検出部
５，１０１  画像入力部
６  変化検出対象種別入力部
７  符号化結果調整部
８  画像解析部
９  符号化変化検出対象領域抽出処理部
１１  開始点候補入力部
１２  ブロック特性分類部
１３  領域ごと利用符号化情報・しきい値入力部
１０３  ブロックごと符号化情報特定部
１０４，２０４  複雑さ指標値算出部
２０１  符号化データ入力部
２０２  可変長符号復号部
２０３  ブロックごと符号化情報特定部
２０５  共通領域特定部
２０６  複雑さ指標値補正部
２０８  位置情報メタデータ入力部
２０９  線形和係数入力部
【図１】
2020025290.tif 000003
【図２】
2020025290.tif 000004
【図３】
2020025290.tif 000005
【図４】
2020025290.tif 000006
【図５】
2020025290.tif 000007
【図６】
2020025290.tif 000008
【図７】
2020025290.tif 000009
【図８】
2020025290.tif 000010
【図９】
2020025290.tif 000011
【図１０】
2020025290.tif 000012
【図１１】
2020025290.tif 000013
【図１２】
2020025290.tif 000014
【図１３】
2020025290.tif 000015
【図１４】
2020025290.tif 000016
【図１５】
2020025290.tif 000017
【図１６】
2020025290.tif 000018
【図１７】
2020025290.tif 000019
【図１８】
2020025290.tif 000020
【図１９】
2020025290.tif 000021
【図２０】
2020025290.tif 000022
【図２１】
2020025290.tif 000023
【図２２】
2020025290.tif 000024
【図２３】
2020025290.tif 000025
【図２４】
2020025290.tif 000026
【図２５】
2020025290.tif 000027
【図２６】
2020025290.tif 000028
【図２７】
2020025290.tif 000029