(57)【特許請求の範囲】
【請求項１】
圧縮センシングを利用して対象信号について複数回の観測を行う信号観測装置であって、
対象信号の光量を単一の検出器で観測する観測部と、
行列状に配列された複数の素子であって前記観測部に伝送される前記対象信号の光量を個別に制御可能な複数の素子を有するフィルタ部と、
各行列要素の値が予め定められた第一の観測行列に基づいて前記フィルタ部の各素子を制御する第一の制御と、各行列要素の値が乱数に基づく値を有する第二の観測行列に基づいて前記フィルタ部の各素子を制御する第二の制御と、を選択的に用いて、前記フィルタ部を介して伝送された前記対象信号の光量を前記観測部に複数回の観測をさせ、前記第一の制御を用いて観測された対象信号の光量と、前記第二の制御を用いて観測された対象信号の光量との差分値を算出させる制御部と、
算出された前記差分値を用いて対象信号を復元する信号復元装置へ前記差分値を送信する通信部と、を備え、
前記第一の観測行列の全ての行列要素の値は、同じ値であり、
前記第一の制御では、前記フィルタ部の各素子は、前記第一の観測行列に基づいて、前記観測部に伝送される前記対象信号の光量を制限する制限状態に制御される、
信号観測装置。
【請求項２】
前記複数回の観測のうち、最初の観測で、前記第一の制御が適用される、
請求項１に記載の信号観測装置。
【請求項３】
前記複数回の観測のうち、真ん中の観測で、前記第一の制御が適用される、
請求項１に記載の信号観測装置。
【請求項４】
前記複数回の観測のうち、所定の回数毎の観測で、前記第一の制御が適用される、
請求項１に記載の信号観測装置。
【請求項５】
前記複数回の観測において、前記第一の制御と前記第二の制御とが交互に適用される、
請求項１に記載の信号観測装置。
【請求項６】
圧縮センシングを利用して対象信号について複数回の観測を行う信号観測方法であって、
行列状に配列された複数の素子であって観測部に伝送する前記対象信号の光量を個別に制御可能な複数の素子の各々を、各行列要素の値が予め定められた第一の観測行列に基づいて制御する第一の制御を実行し、
各行列要素の値が乱数に基づく値を有する第二の観測行列に基づいて前記複数の素子の各々を制御する第二の制御を実行し、
前記第一の制御の実行又は前記第二の制御の実行で制御された前記複数の素子を介して前記観測部に伝送された前記対象信号を観測し、
前記観測は、複数回行われ、
複数回の前記観測のうち少なくとも１回は、前記第一の制御の実行によって制御された前記複数の素子を介して前記観測部に伝送された対象信号を観測し、
前記第一の制御を用いて観測された対象信号の光量と、前記第二の制御を用いて観測された対象信号の光量との差分値を算出し、
算出された前記差分値を用いて対象信号を復元する信号復元装置へ前記差分値を送信し、
前記第一の観測行列の全ての行列要素の値は、同じ値であり、
前記第一の制御では、前記複数の素子の各々は、前記第一の観測行列に基づいて、前記観測部に伝送される前記対象信号の光量を制限する制限状態に制御する、
信号観測方法。
【請求項７】
圧縮センシングを利用して対象信号について複数回の観測を行うコンピュータ読み取り可能なプログラムであって、
前記プログラムは、コンピュータに
行列状に配列された複数の素子であって観測部に伝送する前記対象信号の光量を個別に制御可能な複数の素子の各々を、各行列要素の値が予め定められた第一の観測行列に基づいて制御する第一の制御を実行し、
各行列要素の値が乱数に基づく値を有する第二の観測行列に基づいて前記複数の素子の各々を制御する第二の制御を実行し、
前記第一の制御の実行又は前記第二の制御の実行で制御された前記複数の素子を介して前記観測部に伝送された前記対象信号を観測する
処理を実行させ、
前記観測は、複数回行われ、
複数回の前記観測のうち少なくとも１回は、前記第一の制御の実行によって制御された前記複数の素子を介して前記観測部に伝送された対象信号を観測し、
前記第一の制御を用いて観測された対象信号の光量と、前記第二の制御を用いて観測された対象信号の光量との差分値を算出し、
算出された前記差分値を用いて対象信号を復元する信号復元装置へ前記差分値を送信し、
前記第一の観測行列の全ての行列要素の値は、同じ値であり、
前記第一の制御では、前記複数の素子の各々は、前記第一の観測行列に基づいて、前記観測部に伝送される前記対象信号の光量を制限する制限状態に制御する、
プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、圧縮センシングを用いた信号観測装置に関する。
【背景技術】
【０００２】
近年、圧縮センシングを用いた信号観測技術が注目されている。この圧縮センシングは、圧縮サンプリングあるいはスパースサンプリングとも呼ばれる。圧縮センシングは、信号がスパース（疎）であるという条件下で、標本化定理以下のサンプル数で取得した信号から、元信号を復元することを可能にする技術であり、信号の取得と圧縮を同時に行う技術である。
【０００３】
信号がスパースであるとは、信号をウェーブレット空間や離散コサイン空間などへ射影した際、多くの係数がほぼ０になるという事象を意味する。信号のスパース性を利用した信号の復元手法としては、Ｌ１ノルム最小化などが利用される（例えば、非特許文献１を参照）。
【０００４】
圧縮センシングを用いたレンズのない画像取得が特許文献１に開示されている。この特許文献１の画像取得では、シャッターアレイを用いて、被写体から反射される光を単一の検出器上に選択的に向かわせる。この光の検出を複数回行うことで、単一の検出器で高解像の画像を取得する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特表２０１５−５１０３５６号公報
【非特許文献】
【０００６】
【非特許文献１】田中利幸著「圧縮センシングの数理」  電子情報通信学会  基礎・境界ソサイエティ  Fundamentals Review Vol.4 No1、２０１０年７月（Ｐ３９−Ｐ４７）
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、高解像の画像を取得するための複数回の観測の間に、外環境の変化や検出器の状態変化などが発生する場合がある。このような場合に、上記従来の技術では、外環境の変化や検出器の状態変化などにより生じるノイズによって、圧縮信号から復元される画像の画質が劣化したり、圧縮信号から画像を復元できなくなったりする場合がある。
【０００８】
そこで、本発明は、圧縮センシングを用いた信号観測において、ノイズによる信号の劣化を抑制することができる信号観測装置及び信号観測方法を提供する。
【課題を解決するための手段】
【０００９】
本発明の一態様に係る信号観測装置は、圧縮センシングを利用して対象信号について複数回の観測を行う信号観測装置であって、対象信号の光量を単一の検出器で観測する観測部と、行列状に配列された複数の素子であって前記観測部に伝送される前記対象信号の量を個別に制御可能な複数の素子を有するフィルタ部と、各行列要素の値が予め定められた第一の観測行列に基づいて前記フィルタ部の各素子を制御する第一の制御と、各行列要素の値が乱数に基づく値を有する第二の観測行列に基づいて前記フィルタ部の各素子を制御する第二の制御と、を選択的に用いて、前記フィルタ部を介して伝送された前記対象信号の光量を前記観測部に複数回の観測をさせ、前記第一の制御を用いて観測された対象信号の光量と、前記第二の制御を用いて観測された対象信号の光量との差分値を算出させる制御部と、算出された前記差分値を用いて対象信号を復元する信号復元装置へ前記差分値を送信する通信部と、を備え、前記第一の観測行列の全ての行列要素の値は、同じ値であり、前記第一の制御では、前記フィルタ部の各素子は、前記第一の観測行列に基づいて、前記観測部に伝送される前記対象信号の光量を制限する制限状態に制御される。
【００１０】
なお、これらの包括的又は具体的な態様は、システム、方法、集積回路、コンピュータプログラム又はコンピュータ読み取り可能なＣＤ−ＲＯＭなどの記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラム及び記録媒体の任意な組み合わせで実現されてもよい。
【発明の効果】
【００１１】
本発明の一態様に係る信号観測装置によれば、圧縮センシングを用いた信号観測において、ノイズによる信号の劣化を抑制することができる。
【図面の簡単な説明】
【００１２】
【図１】図１は、実施の形態１における信号観測システムの外観図である。
【図２】図２は、実施の形態１における信号観測システムの構成の一例を示すブロック図である。
【図３】図３は、実施の形態１における観測行列の一例を示す図である。
【図４Ａ】図４Ａは、図３で示した２行２列の観測行列に基づいて２ｘ２画素の観測行列表示部を制御する例を示す図である。
【図４Ｂ】図４Ｂは、図３で示した２行２列の観測行列に基づいて３ｘ３画素の観測行列表示部を制御した結果の一例を示す図である。
【図４Ｃ】図４Ｃは、図３で示した２行２列の観測行列に基づいて４ｘ４画素の観測行列表示部を制御した結果の一例を示す図である。
【図４Ｄ】図４Ｄは、図３で示した２行２列の観測行列に基づいて５ｘ５画素の観測行列表示部を制御した結果の一例を示す図である。
【図４Ｅ】図４Ｅは、図３で示した２行２列の観測行列に基づいて６ｘ６画素の観測行列表示部を制御した結果の一例を示す図である。
【図５】図５は、実施の形態１における信号観測装置による信号観測処理の一例を示すフローチャートである。
【図６】図６は、実施の形態１における信号観測装置による初期設定処理の一例を示すフローチャートである。
【図７】図７は、実施の形態１における観測行列生成部による観測行列生成処理の一例を示すフローチャートである。
【図８】図８は、実施の形態１における観測行列表示部の制御処理の一例を示すフローチャートである。
【図９Ａ】図９Ａは、実施の形態１における第二の観測行列に基づく観測行列表示部の制御結果の一例を示す図である。
【図９Ｂ】図９Ｂは、実施の形態１における第一の観測行列に基づく観測行列表示部の制御結果の一例を示す図である。
【図９Ｃ】図９Ｃは、実施の形態１における第二の観測行列に基づく観測行列表示部の制御結果の他の一例を示す図である。
【図１０】図１０は、実施の形態１における信号復元装置による信号復元処理の一例を示すフローチャートである。
【図１１】図１１は、実施の形態１における信号復元装置の通信部によるデータセット受信処理の一例を示すフローチャートである。
【図１２】図１２は、実施の形態１における信号復元装置の復元行列生成部による復元行列生成処理を示すフローチャートである。
【図１３Ａ】図１３Ａは、実験で用いた被写体を示す図である。
【図１３Ｂ】図１３Ｂは、実験で復元された画像を示す図である。
【図１４】図１４は、実施の形態２における信号観測装置による信号観測処理の一例を示すフローチャートである。
【図１５】図１５は、実施の形態３における信号観測装置による信号観測処理の一例を示すフローチャートである。
【図１６】図１６は、実施の形態４における信号観測装置による信号観測処理の一例を示すフローチャートである。
【図１７】図１７は、実施の形態５における信号観測システムの構成の一例を示すブロック図である。
【図１８】図１８は、実施の形態５における識別番号とアドレスとの対応関係を示す情報の一例を示す図である。
【図１９】図１９は、実施の形態５における信号観測装置による信号観測処理の一例を示すフローチャートである。
【図２０】図２０は、実施の形態５における観測行列生成部による観測行列生成処理の一例を示すフローチャートである。
【図２１】図２１は、実施の形態５における復元行列生成部による復元行列生成処理を示すフローチャートである。
【発明を実施するための形態】
【００１３】
本発明の一態様に係る信号観測装置は、圧縮センシングを利用して１つの対象信号について複数回の観測を行う信号観測装置であって、対象信号の量を観測する観測部と、行列状に配列された複数の素子であって前記観測部に伝送される前記対象信号の量を個別に制限可能な複数の素子を有するフィルタ部と、各行列要素の値が予め定められた第一の観測行列に基づいて前記フィルタ部の各素子を制御する第一の制御と、各行列要素の値が乱数に基づく値を有する第二の観測行列に基づいて前記フィルタ部の各素子を制御する第二の制御と、を選択的に用いて、前記フィルタ部を介して伝送された前記対象信号の量を前記観測部に複数回の観測をさせる制御部と、を備える。
【００１４】
この構成によれば、行列要素の値が予め定められた値を有する第一の観測行列を用いた第一の制御と、行列要素の値が乱数に基づく値を有する第二の観測行列を用いた第二の制御とを選択的に用いて、複数回の観測を行うことができる。したがって、第一の制御によって得られた対象信号の観測値と第二の制御によって得られた対象信号の観測値との差分値を求めることでき、その差分値を用いて対象信号を復元することが可能となる。その結果、複数回の観測における観測行列による観測値の変化を強調して信号の復元を行うことができ、ノイズによる信号の劣化を抑制することができる。
【００１５】
また例えば、前記第一の観測行列の行列要素の値は、すべて同じであるとしてもよい。
【００１６】
この構成によれば、第一の観測行列として、行列要素の値がすべて同じである観測行列を用いることができる。したがって、フィルタ部の複数の素子を一様に制御することができ、複数回の観測における観測行列による観測値の変化を適切に強調することができる。
【００１７】
また例えば、前記第一の制御では、前記フィルタ部の各素子は、前記第一の観測行列に基づいて、前記観測部に伝送される前記対象信号の量を制限する制限状態に制御されてもよい。
【００１８】
この構成によれば、第一の制御においてフィルタ部の各素子を制限状態に制御することができる。したがって、対象信号の観測部への伝送が制限された状態における観測値を用いて差分値を求めることができるので、第二の制御によって得られた対象信号の観測値からノイズを減少させることができ、ノイズによる信号の劣化を抑制することができる。
【００１９】
また例えば、前記複数回の観測のうち、最初の観測で、前記第一の制御が適用されてもよい。
【００２０】
この構成によれば、最初の観測で第一の制御を適用することができる。したがって、後続する第二の制御を用いた観測では、既に第一の制御を用いた観測値が得られているので、差分値をすぐに求めることができ、観測値のためのバッファを減少させることができる。
【００２１】
また例えば、前記複数回の観測のうち、真ん中の観測で、前記第一の制御が適用されてもよい。
【００２２】
この構成によれば、真ん中の観測で第一の制御を適用することができる。したがって、複数回の観測において、比較的時間的に近い観測値同士で差分値を求めることができ、経時変化に伴うノイズの影響を抑制することが可能となる。
【００２３】
また例えば、前記複数回の観測のうち、所定の回数毎の観測で、前記第一の制御が適用されてもよい。
【００２４】
この構成によれば、所定回数毎の観測で第一の制御を適用することができる。したがって、複数回の観測において、比較的時間的に近い観測値同士で差分値を求めることができ、経時変化に伴うノイズの影響を抑制することが可能となる。
【００２５】
また例えば、前記複数回の観測において、前記第一の制御と前記第二の制御とが交互に適用されてもよい。
【００２６】
この構成によれば、第一の制御及び第二の制御を交互に適用することができる。したがって、複数回の観測において、時間的に隣り合う観測値同士で差分値を求めることができ、経時変化に伴うノイズの影響を抑制することが可能となる。
【００２７】
また例えば、前記制御部は、さらに、前記第一の制御を用いて観測された対象信号の量と、前記第二の制御を用いて観測された対象信号の量との差分値を算出し、前記信号観測装置は、さらに、算出された前記差分値を用いて対象信号を復元する信号復元装置へ前記差分値を送信する通信部を備えてもよい。
【００２８】
この構成によれば、信号観測装置において差分値を算出することができる。したがって、信号復元装置において差分値を算出する必要がなく、対象信号の復元を容易にすることができる。さらに、信号観測装置１１０から信号復元装置１２０への通信量を削減することもできる。
【００２９】
なお、これらの包括的又は具体的な態様は、システム、方法、集積回路、コンピュータプログラム又はコンピュータ読み取り可能なＣＤ−ＲＯＭなどの記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラム及び記録媒体の任意な組み合わせで実現されてもよい。
【００３０】
以下、実施の形態について、図面を用いて詳細に説明する。
【００３１】
なお、以下で説明する実施の形態は、いずれも包括的又は具体的な例を示す。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、ステップ、ステップの順序などは、一例であり、請求の範囲を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、最上位概念を示す独立請求項に記載されていない構成要素については、任意の構成要素として説明される。
【００３２】
（実施の形態１）
［信号観測システムの構成］
図１は、実施の形態１における信号観測システムの外観図である。図２は、実施の形態１における信号観測システムの構成の一例を示すブロック図である。図１において、破線は、見えない部分を表す。図１及び図２に示すように、信号観測システム１００は、信号観測装置１１０、及び、信号復元装置１２０を備える。
【００３３】
信号観測装置１１０は、圧縮センシングを用いて観測の対象となる対象信号を複数回観測し、観測結果として、後述する圧縮信号を出力する。本実施の形態では、対象信号が、被写体９９から放射又は反射された可視光信号である場合について説明する。また、信号観測装置１１０は、同じ被写体９９から放射又は反射する可視光信号を複数回観測するものとする。なお、対象信号は、可視光信号に限定されない。例えば、対象信号は、赤外線信号、テラヘルツ波信号、又は紫外線信号であってもよい。
【００３４】
信号観測装置１１０は、入力部１１１、制御部１１２、乱数シード生成部１１３、観測行列生成部１１４、観測行列表示部１１５、観測部１１６、及び、通信部１１７を備える。
【００３５】
入力部１１１は、ユーザからの入力を受け付け、制御部１１２に入力を伝達する。例えば、入力部１１１はキーボードやタッチパネルなどである。具体的には、入力部１１１は、対象信号（画像）の画素数の入力を受ける。ここでは、対象信号の画素数とは、例えば、対象信号が示す画像の縦方向の画素数（つまり第一の行数）及び横方向の画素数（つまり第一の列数）で表される。
【００３６】
制御部１１２は、入力部１１１からの情報を得て、乱数シード生成部１１３、観測行列生成部１１４、観測行列表示部１１５、観測部１１６、及び、通信部１１７を制御する。
【００３７】
なお、制御部１１２は、制御機能を備えるものであればよく、どのように実現されてもよい。例えば、制御部１１２は、専用のハードウェアで構成されてもよい。また例えば、制御部１１２は、各構成要素に適したソフトウェアプログラムを実行することによって実現されてもよい。この場合、制御部１１２は、例えば、演算処理部（図示せず）と制御プログラムを記憶する記憶部（図示せず）とを備えてもよい。演算処理部としては、ＭＰＵ（Micro Processing Unit）、ＣＰＵ（Central Processing Unit）が例示される。記憶部としては、半導体メモリが例示される。なお、制御部１１２は、集中制御を行う単独の制御部で構成されていてもよく、互いに協働して分散制御を行う複数の制御部で構成されていてもよい。
【００３８】
乱数シード生成部１１３は、観測行列生成部１１４で利用するための複数の乱数シードを生成する。
【００３９】
観測行列生成部１１４は、入力部１１１から得られる対象信号の画素数（第一の行数及び第一の列数）から決定されるサイズを有する複数の観測行列を生成する。観測行列生成部１１４は、例えば、入力部１１１より入力された対象信号の縦方向の画素数を生成すべき観測行列の行数（第一の行数）として決定する。また、観測行列生成部１１４は、例えば、入力部１１１より入力された対象信号の横方向の画素数を生成すべき観測行列の列数（第一の列数）として決定する。これにより生成すべき観測行列のサイズとして、第一の行数および第一の列数が、決定される。
【００４０】
この場合、生成される観測行列に含まれる行列要素の数は、第一の行数と第一の列数との積で示される。
【００４１】
また、観測行列に含まれる各行列要素の値をどのようにして決定するのかについては、後述する。
【００４２】
観測行列とは、例えば、後述する第一の観測行列および第二の観測行列を含む。
【００４３】
観測行列生成部１１４は、入力部１１１から得られる第一の行数及び第一の列数を有する第一の観測行列を生成する。
【００４４】
生成される第一の観測行列は、予め定められたパターンを有する。第一の観測行列が予め定められたパターンを有するとは、例えば、第一の観測行列に含まれる全てまたは、ほぼ全ての行列要素の値が予め決められていることを意味する。
【００４５】
予め定められたパターンの第一の観測行列は、例えば、第一の観測行列に含まれる複数の行列要素の値が一様な値を有する。第一の観測行列に含まれる複数の行列要素の値が一様であるとは、例えば、第一の観測行列に含まれる全て又はほぼ全ての行列要素が実質的に同じ値（例えば、０）を有することを意味する。本実施の形態では、第一の観測行列に含まれる全ての行列要素の値を０としたものを例に説明をする。
【００４６】
さらに、観測行列生成部１１４は、入力部１１１から得られる第一の行数及び第一の列数を有する第二の観測行列を生成する。観測行列生成部１１４は、例えば、第二の観測行列の各行列要素の値を、乱数シード生成部１１３が生成した乱数シードに基づき算出する。
【００４７】
観測行列表示部１１５は、フィルタ部の一例である。観測行列表示部１１５は、対象信号が観測部１６０へと伝送される伝送経路上に配置されている。例えば、観測行列表示部１１５は、行列状に配列された複数の素子を有する。
【００４８】
複数の素子は、観測部１１６に伝送される対象信号の量を個別に制限可能である。例えば、観測行列表示部１１５は、透過率を個別に制御可能な複数の画素を有する透過型液晶ディスプレイである。
【００４９】
例えば、透過型液晶ディスプレイに含まれる複数の画素のそれぞれは、行列状に配列された複数の素子の何れかと対応する。観測行列表示部１１５の複数の素子のそれぞれは、その状態を制御することにより、対応する画素の透過率を制御する。観測行列表示部１１５は、例えば、各素子の状態を制御することで、各画素の透過率を制御する。結果として、観測行列表示部１１５は、観測部１６０へ伝送される対象信号の量を制御することができる。
【００５０】
観測行列表示部１１５は、図１に示すように、被写体９９及び観測部１１６の間に配置される。
【００５１】
制御部１１２は、観測行列生成部１１４が生成した観測行列に基づいて観測行列表示部１１５の各素子の状態を変更して、観測行列表示部１１５を介して伝送された対象信号の量の観測を観測部１１６に行わせる。
【００５２】
例えば、制御部１１２は、観測行列生成部１１４が生成した観測行列に含まれる各行列要素を観測行列表示部１１５の複数の素子のいずれかに対応付ける。
【００５３】
制御部１１２は、観測行列に含まれる一の行列要素と観測行列表示部１１５の一の素子とを一対一に対応付けてもよい。
【００５４】
または、制御部１１２は、観測行列に含まれる一の行列要素と観測行列表示部１１５の複数の素子とを対応付けてもよい。但し、制御部１１２は、観測行列に含まれる複数の行列要素を観測行列表示部１１５の中の同一の素子に対応付けることはしない。
【００５５】
制御部１１２は、観測行列の行列要素と対応付けた観測行列表示部１１５の素子の状態を当該行列要素の値に基づいて制御する。なお、制御部１１２は、観測行列表示部１１５の複数の素子のうち、観測行列の行列要素と対応付けていない素子がある場合、この素子の状態を予め定めた状態に制御する。
【００５６】
このように観測行列の各行列要素を観測行列表示部１１５の複数の素子のいずれと対応付け、各行列要素の値に基づいて対応する素子の状態を制御するとともに、観測行列の行列要素と対応付けていない観測行列表示部１１５の素子の状態を予め定めた状態に制御することを、観測行列に基づいて観測行列表示部１１５の各素子（または各画素）を制御するという。
【００５７】
上述したように観測行列生成部１１４は、第一の観測行列および第二の観測行列を生成する。制御部１１２は、第一の観測行列に基づいて観測行列表示部１１５の各素子を制御する第一の制御と、第二の観測行列に基づいて観測行列表示部１１５の各素子を制御する第二の制御とを選択的に用いて、観測部１１６に対象信号の観測を複数回、行わせる。
【００５８】
観測部１１６に対象信号の観測を行わせるというのは、観測行列表示部１１５を介して観測部１１６が備える単一の検出器に伝送される対象信号の量を測定（または検出、または出力ともいう）させることである。
【００５９】
制御部１１２が、第一の制御を用いて、観測部１１６に対象信号の観測を行わせるとは、例えば、制御部１１２が、第一の制御を行っている間に観測部１１６に対象信号の観測を行わせることである。
【００６０】
制御部１１２が、第二の制御を用いて、観測部１１６に対象信号の観測を行わせるとは、例えば、制御部１１２が第二の制御を行っている間に観測部１１６に対象信号の観測を行わせることである。
【００６１】
制御部１１２が、第一の制御と第二の制御とを選択的に用いて、観測部１１６に対象信号の観測を複数回、行わせるというのは、例えば、制御部１１２が、第一の制御および第二の制御のうちの一方を選択し、選択した制御を行っている間に観測部１１６に観測を行わせる動作を複数回行うことである。
【００６２】
上述の動作を行う複数回の回数は、予め決められていてもよい。
【００６３】
この複数回の動作には、制御部１１２が、第一の制御を選択し、観測部１１６に観測を行わせる動作が少なくとも１回含まれる。
【００６４】
また、複数回の動作には、制御部１１２が、第二の制御を選択し、観測部１１６に観測を行わせる動作が少なくとも２回含まれる。この場合、制御部１１２が、第二の制御を選択する毎に、乱数シード生成部１１３が異なる乱数シードを生成し、観測行列生成部１１４が、生成した乱数シード基づき第二の観測行列を生成するのであれば、第二の制御を選択する毎に生成される第二の観測行列は、互いに異なる行列となる。
【００６５】
また複数回の動作の内、１回目の動作（つまり最初の動作、または観測部１１６による最初の観測ともいう）が、制御部１１２が、第一の制御を選択し、観測部１１６に観測を行わせる動作であってもよい。
【００６６】
または、複数回の動作の内、真ん中の回数目の動作（例えば、第二の制御を選択する回数を２で割ったときの商の値に対応する回数）が、制御部１１２が、第一の制御を選択し、観測部１１６に観測を行わせる動作であってもよい。
【００６７】
または、制御部１１２が、第一の制御を選択し、観測部１１６に観測を行わせる動作および制御部１１２が、第二の制御を選択し、観測部１１６に観測を行わせる動作を交互に行ってもよい。この場合、制御部１１２が、第一の制御を選択し、観測部１１６に観測を行わせる動作および制御部１１２が、第二の制御を選択し、観測部１１６に観測を行わせる動作を１組の動作とし、この組にした動作を複数回行うとしてもよい。
【００６８】
次に、制御部１１２が、観測行列を用いて観測行列表示部１１５の各素子を制御する具体例を説明する。
【００６９】
制御部１１２は、観測行列に含まれる各行列要素の値に応じて、観測行列表示部１１５の対応する素子（または画素ともいう）の状態を制御する。例えば、観測行列の行列要素の値が０であれば、制御部１１２は、当該行列要素に対応する観測行列表示部１１５の素子の状態を制限状態に制御する。一方、観測行列の行列要素の値が１であれば、当該行列要素に対応する観測行列表示部１１５の素子の状態を非制限状態に制御する。
【００７０】
制限状態とは、観測部１１６に伝送する対象信号の量を制限する状態であり、例えば非透過状態である。非制限状態とは、観測部１１６に伝送する対象信号の量を制限しない状態であり、例えば透過状態である。
【００７１】
ここで、図３及び図４Ａ〜図４Ｅを参照して、観測行列の行列要素と、観測行列表示部１１５の素子との対応関係について説明する。以降の図では、ハッチング領域は、素子の状態を制限状態（非透過状態）としたときの画素を表し、白領域は、素子の状態を非制限状態（透過状態）としたときの画素を表す。また、０の値を有する行列要素に対応する素子の状態は制限状態に制御され、１の値を有する行列要素に対応する素子の状態は非制限状態に制御される場合について説明する。
【００７２】
図３は、観測行列の一例を示す。図３に示す観測行列は、２行２列の観測行列である。この例では、入力部１１１より、対象信号の縦方向の画素数として２、および対象信号の横方向の画素数として２が入力されたときに生成される観測行列（第二の観測行列）の一例を示している。
【００７３】
観測行列生成部１１４は、決定されたサイズ（この例では２行２列）の第一の観測行列および第二の観測行列をそれぞれ生成する。
【００７４】
観測行列生成部１１４は、第一の観測行列を生成する場合、各行列要素の値として、対応する素子を制限状態するために用いる値（この例では、０）とする。
【００７５】
また、観測行列生成部１１４は、第二の観測行列を生成する場合、各行列要素の値は、所定の確率分布に従う乱数値とする。具体的には、例えば、乱数シード生成部１１３が生成した乱数シードに基づいて、所定の確率分布に従う乱数を第二の観測行列の行列要素の個数分（この例では４個）、算出し、算出した値を第二の観測行列の行列要素の値とする。
【００７６】
図３では、観測行列の各行列要素は、０又は１の値を有する。
【００７７】
図４Ａは、図３で示した２行２列の観測行列に基づいて２ｘ２画素の観測行列表示部１１５を制御した結果の一例を示す図である。図４Ａの例では、例えば、観測行列表示部１１５を構成する、４個の素子が２行２列の行列状に配列されている。また、４個の素子は、それぞれ、状態を制御することにより素子に対応する画素の透過率が制御される。この図４Ａの例では、観測行列の第一の行数及び第一の列数と観測行列表示部１１５の第二の行数及び第二の列数とが一致しているため、観測行列の４つの行列要素のいずれかと観測行列表示部１１５の４つの素子（画素）のいずれかとを一対一で対応付けることができる。この例では制御部１１２は、行番号及び列番号が同じ行列要素と素子とを対応付けている。
【００７８】
制御部１１２は、第二の観測行列の各行列要素と、観測行列表示部１１５の各素子との対応付けを行い、第二の観測行列の行列要素の値に基づいて、対応する観測行列表示部１１５の素子の状態を制御する。
【００７９】
図４Ｂは、図３で示した２行２列の観測行列に基づいて３ｘ３画素の観測行列表示部１１５を制御した結果の一例を示す図である。図４Ｂの例では、例えば、観測行列表示部１１５を構成する９個の素子が３行３列の行列状に配列されている。また、９個の素子は、それぞれ、状態を制御することにより素子に対応する画素の透過率が制御される。
【００８０】
この図４Ｂの例では、観測行列の第一の行数及び第一の列数と観測行列表示部１１５の第二の行数及び第二の列数とが一致していない。そのため、制御部１１２は、観測行列表示部１１５の左上の素子を基準とする２行２列に配列された４つの素子（２×２の素子ともいう）のいずれかと、観測行列の４つの行列要素のいずれかとを一対一で対応付ける。その際、制御部１１２は、対応する行列要素が存在しない、観測行列表示部１１５の右端及び下端の複数の素子の状態を、予め定めた状態に制御する。予め定めた状態とは、制限状態（非透過状態）である。
【００８１】
制御部１１２は、第二の観測行列の行列要素の値に基づいて、対応する観測行列表示部１１５の素子の状態を制御する。
【００８２】
図４Ｃは、図３で示した２行２列の観測行列に基づいて４ｘ４画素の観測行列表示部を制御した結果の一例を示す図である。図４Ｃの例では、例えば、観測行列表示部１１５を構成する１６個の素子が４行４列の行列状に配列されている。また、１６個の素子は、それぞれ、状態を制御することにより素子に対応する画素の透過率が制御される。
【００８３】
この図４Ｃの例では、観測行列表示部１１５の第二の行数及び第二の列数は、観測行列の第一の行数及び第一の列数の２倍である。このとき、観測行列表示部１１５の２×２の素子を１つのグループとしてまとめてみれば、観測行列表示部１１５は、４つのグループが２行２列の行列状に配列されていると考えられる。
【００８４】
このように考えれば、観測行列表示部１１５の４つのグループのいずれかと、観測行列の４つの行列要素のいずれかとを一対一で対応付けることができる。
【００８５】
つまり、制御部１１２は、観測行列の各行列要素を対応するグループに含まれる複数の素子（この例では、２ｘ２の素子）と対応付ける。
【００８６】
そして、制御部１１２は、観測行列の各行列要素の値に応じて、当該行列要素に対応するグループに含まれる複数の素子の状態を制御する。
【００８７】
つまり、制御部１１２は、観測行列に含まれる行列要素の数と観測行列表示部１１５の複数の素子の数とに基づいて複数の素子をグループ化し、グループ単位で複数の素子の状態を変更する。
【００８８】
観測行列表示部１１５の複数の素子をグループ化する具体的な処理を以下に説明する。
【００８９】
具体的には、第一の行数より第二の行数が大きい場合に、制御部１１２は、観測行列表示部１１５の複数の素子の行列の行を、第一の行数及び第二の行数に基づく数の行からなるグループにグループ化し、同一の行のグループに含まれる同一の列の素子の状態を同一に制御する。第一の行数及び第二の行数に基づく数は、例えば、第一の行数で第二の行数を除算したときの商が示す数（ここでは２（＝４／２））である。
【００９０】
また、第一の列数より第二の列数が大きい場合に、制御部１１２は、観測行列表示部１１５の複数の素子の行列の列を、第一の列数及び第二の列数に基づく数の列からなるグループにグループ化し、同一の列のグループに含まれる同一の行の素子の状態を同一の状態に制御する。第一の列数及び第二の列数に基づく数は、例えば、第一の列数で第二の列数を除算したときの商が示す数（ここでは２（＝４／２））である。
【００９１】
図４Ｄは、図３で示した２行２列の観測行列に基づいて５ｘ５画素の観測行列表示部１１５を制御した結果の一例を示す図である。図４Ｄの例では、例えば、観測行列表示部１１５を構成する２５個の素子が５行５列の行列状に配列されていることを意味する。また、２５個の素子は、それぞれ、状態を制御することにより素子に対応する画素の透過率が制御される。
【００９２】
この図４Ｄの例では、観測行列表示部１１５の第二の行数及び第二の列数は、観測行列の第一の行数及び第一の列数の２倍にそれぞれ１を足した数である。例えば、図４Ｄに示す観測行列表示部１１５の左端及び上端に配列した複数の素子以外の素子は、４×４の素子で構成される。この場合、４×４の素子で構成される部分において、２×２の素子を１つのグループとしてまとめてみれば、観測行列表示部１１５は、４つのグループが２行２列の行列状に配列されていると考えられる。
【００９３】
このように考えれば、観測行列表示部１１５の４つのグループのいずれかと観測行列の４つの行列要素のいずれかとを一対一で対応付けることができる。
【００９４】
つまり、制御部１１２は、観測行列の各行列要素を対応するグループに含まれる複数の素子（この例では、２ｘ２の素子）と対応付ける。その際、制御部１１２は、対応する行列要素が存在しない、観測行列表示部１１５の左端及び上端の複数の素子の状態を、予め定めた状態に制御する。予め定めた状態とは、制限状態（非透過状態）である。
【００９５】
そして、制御部１１２は、観測行列の各行列要素の値に応じて、当該行列要素に対応するグループに含まれる複数の素子の状態を制御する。
【００９６】
観測行列表示部１１５の複数の素子をグループ化する具体的な処理を以下に説明する。      
【００９７】
つまり、第一の行数で第二の行数を除算したときに余りがある場合、制御部１１２は、観測行列表示部１１５の複数の素子の行列の行から余りが示す数（ここでは、１（＝５％２）ここで％は剰余演算子）の行を除外して残った行を、商（ここでは、２（＝５／２））が示す数の行からなるグループにグループ化し、除外された行に含まれる素子の状態を、観測部１１６に伝送する対象信号の量を制限する制限状態に制御する。また、第一の列数で第二の列数を除算したときに余りがある場合、制御部１１２は、観測行列表示部１１５の複数の素子の行列の列から余りが示す数（ここでは、１（＝５％２））の列を除外して残った列を、商（ここでは、２（＝５／２））が示す数の列からなるグループにグループ化し、除外された列に含まれる素子の状態を、観測部１１６に伝送する対象信号の量を制限する制限状態に制御する。
【００９８】
図４Ｅは、図３で示した２行２列の観測行列に基づいて６ｘ６画素の観測行列表示部１１５を制御した結果の一例を示す図である。図４Ｅの例では、例えば、観測行列表示部１１５を構成する３６個の素子が６行６列の行列状に配列されている。また、３６個の素子は、それぞれ、状態を制御することにより素子に対応する画素の透過率が制御される。
【００９９】
この図４Ｅの例では、観測行列表示部１１５の第二の行数及び第二の列数は、観測行列の第一の行数及び第一の列数の３倍である。観測行列表示部１１５の３×３の素子を１つのグループとしてまとめてみれば、観測行列表示部１１５は、４つのグループが２行２列の行列状に配列されていると考えられる。
【０１００】
このように考えれば、観測行列表示部１１５の４つのグループのいずれかと、観測行列の４つの行列要素のいずれかとを一対一で対応付けることができる。
【０１０１】
つまり、制御部１１２は、観測行列の各行列要素を対応するグループに含まれる複数の素子（この例では３ｘ３の素子）と対応付ける。
【０１０２】
そして、制御部１１２は、観測行列の各行列要素の値に応じて、対応するグループに含まれる複数の素子を制御する。
【０１０３】
なお、図４Ａ〜図４Ｅは、観測行列表示部１１５の制御結果の一例であり、これに限定されない。観測行列の第一の行数及び第一の列数、並びに観測行列表示部１１５の第二の行数、第二の列数が特定できれば、上述の考えを適用することにより、観測行列に基づいて、観測行列表示部１１５を構成する素子の状態を制御することができる。
【０１０４】
また、例えば、図４Ｂ又は図４Ｄにおいて、必ずしも左上又は右下の素子を基準に観測行列の行列要素と観測行列表示部１１５の素子とが対応付けられなくてもよい。例えば、左上又は右下の素子の代わりに、右上、左下又は中央の素子を基準に観測行列の行列要素と観測行列表示部１１５の素子とが対応付けられてもよい。また、ユーザから指示された位置の素子を基準に、観測行列の行列要素と観測行列表示部１１５の素子とが対応付けられてもよい。
【０１０５】
なお、第一の観測行列を用いて、観測行列表示部１１５の各素子を制御するには、上述の説明を適用すればよい。
【０１０６】
観測部１１６は、例えば、単一の検出器で構成される。検出器とは、例えば、光センサ（例えば、フォトダイオード）で構成される。観測部１１６は、観測行列表示部１１５を介して観測部１１６に伝送された対象信号の量を複数回観測する。具体的には、観測部１１６は、観測行列表示部１１５が第一の制御又は第二の制御によって制御されるたびに、観測行列表示部１１５を介して伝送された対象信号の量を測定する（または観測するともいう）。また、観測部１１６で、測定される対象信号の量を対象信号の観測値ともいう。
【０１０７】
通信部１１７は、第一の制御を用いて観測された対象信号の観測値と、第二の制御を用いて観測された対象信号の観測値との差分値を圧縮信号として、信号復元装置１２０へ送信する。具体的には、通信部１１７は、第一の制御で観測行列表示部１１５が制御されたときに観測部１１６で観測された対象信号の観測値と、第二の制御で観測行列表示部１１５が制御されたときに観測部１１６で観測された対象信号の観測値との差分値を圧縮信号として信号復元装置１２０へ送信する。
【０１０８】
また、通信部１１７は、当該第二の制御で用いられた第二の観測行列の生成に利用した乱数シードを差分値とともに信号復元装置１２０へ送信する。例えば、通信部１１７は、差分値と乱数シードとをセットにしたデータセットを信号復元装置１２０へ送信する。さらに、通信部１１７は、観測行列の第一の行数及び第一の列数（または、対象信号の縦の画素数および横の画素数）と、観測回数とを信号復元装置１２０へ送信する。
【０１０９】
信号復元装置１２０は、例えば、信号観測装置１１０から受信した圧縮信号及び乱数シードに基づいて、対象信号を復元する。より具体的には、信号復元装置１２０は、受信した複数のデータセットのそれぞれに含まれる圧縮信号及び乱数シードに基づいて、対象信号を復元する。
【０１１０】
信号復元装置１２０は、入力部１２１、制御部１２２、通信部１２３、復元行列生成部１２４、復元部１２５、及び、出力部１２６を備える。
【０１１１】
入力部１２１は、ユーザからの入力を受け、当該入力を制御部１２２に伝達する。例えば、入力部１２１はキーボードやタッチパネルである。
【０１１２】
制御部１２２は、入力部１２１から情報を得て、通信部１２３、復元行列生成部１２４、復元部１２５、及び、出力部１２６を制御する。なお、制御部１２２は、制御部１１２と同様に、制御機能を備えるものであればよく、どのように実現されてもよい。
【０１１３】
通信部１２３は、信号観測装置１１０から圧縮信号を含む情報を受信する。具体的には、通信部１２３は、信号観測装置１１０から、対象信号の画素数（観測行列の第一の行数及び第一の列数）と、観測回数とを受信する。また、通信部１２３は、信号観測装置１１０から、乱数シード及び差分値を含むデータセットを複数受信する。
【０１１４】
復元行列生成部１２４は、対象信号の画素数、観測回数、及び、複数のデータセットのそれぞれに含まれる乱数シードに基づき、複数の第二の観測行列を生成する。復元行列生成部１２４は、信号観測装置１１０の観測行列生成部１１４が生成した複数の第二の観測行列のそれぞれと同一の第二の観測行列を生成する。つまり、復元行列生成部１２４は、同一の乱数シードに対しては、観測行列生成部１１４と同一の第二の観測行列を生成することができる。そして、復元行列生成部１２４は、生成した複数の第二の観測行列を用いて、圧縮信号から対象信号を復元するための復元行列を生成する。
【０１１５】
復元部１２５は、復元行列生成部１２４が生成した復元行列と、観測行列の第一の行数及び第一の列数と、観測回数と、複数のデータセットに含まれる圧縮信号としての差分値と、基底行列と、に基づき、対象信号を復元する。ここで、基底行列は、信号観測装置１１０で観測した対象信号がスパース性を示す空間の基底行列である。なお、ここでの復元処理は、圧縮センシングで広くつかわれているシンプレックス法や内点法などの公知の手法が利用されればよい（例えば非特許文献１を参照）。
【０１１６】
出力部１２６は、復元部１２５が復元した対象信号を出力する。
【０１１７】
以上が、本実施の形態に係る圧縮信号観測システムの構成についての説明である。
【０１１８】
［信号観測システムの動作］
次に、本実施の形態に係る信号観測システムの動作を、図５〜図１２を参照して具体的に説明する。
【０１１９】
［信号観測処理］
まず、信号観測装置１１０による信号観測処理について図５〜図９Ｃを参照して具体的に説明する。
【０１２０】
図５は、信号観測装置１１０による信号観測処理の一例を示すフローチャートである。図６は、信号観測装置１１０の初期設定処理の一例を示すフローチャートである。具体的には、図６は、図５のステップＳ１０１の詳細を示す。
【０１２１】
図５に示すように、信号観測装置１１０は、信号観測のための初期設定を行う（Ｓ１０１）。具体的には、図６に示すように、制御部１１２は、入力部１１１を介してユーザから、観測行列表示部１１５の横の画素数（以後ＤｅｖＷともいう）と、縦の画素数（以後ＤｅｖＨともいう）の入力を受け、ＤｅｖＷ及びＤｅｖＨの設定を行う（Ｓ２０１）。ＤｅｖＷは、観測行列表示部１１５を構成する行列状に配列された複数の素子の列の数（第二の列数）に相当する。また、ＤｅｖＨは、観測行列表示部１１５を構成する行列状に配列された複数の素子の行の数（第二の行数）に相当する。
【０１２２】
また、制御部１１２は、入力部１１１を介してユーザから、対象信号の縦の画素数（以後ＰｉｃＨともいう）と、対象信号の横の画素数（以後ＰｉｃＷともいう）との入力を受け、ＰｉｃＨ及びＰｉｃＷの設定を行う（Ｓ２０２）。ＰｉｃＷは、観測行列の列の数（第一の列数）に相当し、ＰｉｃＨは、観測行列の行の数（第一の行数）に相当する。
【０１２３】
さらに、制御部１１２は、入力部１１１を介してユーザから、圧縮率（以後ρともいう）の入力を受け、ρを設定する（Ｓ２０３）。制御部１１２は、ＰｉｃＨ、ＰｉｃＷ及びρから観測回数（以後Ｍともいう）を算出し設定する（Ｓ２０４）。
【０１２４】
そして、通信部１１７は、ＰｉｃＨ、ＰｉｃＷ、及び、Ｍを信号復元装置１２０へ送信する（Ｓ２０５）。制御部１１２は、観測回数をカウントするための観測回数カウンタｉの値を０で初期化する（Ｓ２０６）。
【０１２５】
このように初期設定が行われた後、観測行列生成部１１４は、第一の観測行列を生成する（Ｓ１０２）。この第一の観測行列は、ＰｉｃＷｘＰｉｃＨの行数及び列数を有する行列であって、行列要素の全ての値が０である行列である。
【０１２６】
制御部１１２は、生成された第一の観測行列の行列要素と、観測行列表示部１１５を構成する素子との対応付けを行う。対応付けに関しては、図３、図４Ａ〜図４Ｅを用いて説明をしたとおりである。なお、第一の観測行列の行列要素との対応付けがなされなかった観測行列表示部１１５の素子がある場合、制御部１１２は、この素子の状態を予め定めた状態に制御する。予め定めた状態とは制限状態（非透過状態）である。
【０１２７】
制御部１１２は、第一の観測行列に基づいて、観測行列表示部１１５の各画素を制御する（Ｓ１０３）。
【０１２８】
観測行列表示部１１５の各画素を制御するとは、第一の観測行列に基づいて観測行列表示部１１５の各素子を制御することで、各素子に対応する画素の透過率を制御することである。
【０１２９】
このステップＳ１０３における第一の観測行列に基づく観測行列表示部１１５の各素子の制御を第一の制御という。例えば、第一の観測行列を構成する各行列要素の値が一様に０となる場合、制御部１１２は、観測行列表示部１１５の各画素に対応する素子の状態を制限状態に制御する。観測部１１６は、ステップＳ１０３における第一の制御によって観測行列表示部１１５が制御されている間に観測行列表示部１１５を介して伝送された対象信号を観測する（Ｓ１０４）。
【０１３０】
続いて、制御部１１２は、第二の観測行列に基づいて、観測行列表示部１１５の各素子を制御する。具体的には、制御部１１２は、観測回数Ｍと観測回数カウンタｉとを比較する（Ｓ１０５）。ｉがＭ以上である場合（Ｓ１０５のＮ）、信号観測装置１１０は信号観測処理を終了する。
【０１３１】
ｉがＭより小さい場合（Ｓ１０５のＹ）、乱数シード生成部１１３は、乱数シードを生成する（Ｓ１０６）。次に、観測行列生成部１１４は、ＰｉｃＷｘＰｉｃＨの行数及び列数を有する行列であって、この行列に含まれる各行列要素の値を生成された乱数シードに基づき算出した乱数値とする第二の観測行列を生成する（Ｓ１０７）。制御部１１２は、生成された第二の観測行列の行列要素と、観測行列表示部１１５を構成する素子との対応付けを行う。対応付けに関しては、図３、図４Ａ〜図４Ｅを用いて説明をしたとおりである。制御部１１２は、生成された第二の観測行列の各行列要素の値に基づき、観測行列表示部１１５の各素子を制御する。なお、第二の観測行列の行列要素との対応付けがなされなかった観測行列表示部１１５の素子がある場合、制御部１１２は、この素子を制限状態（非透過状態）に制御する（Ｓ１０８）。このステップＳ１０８における第二の観測行列に基づく観測行列表示部１１５の各素子の制御を第二の制御という。観測部１１６は、ステップＳ１０８における第二の制御によって観測行列表示部１１５の各素子の状態が制御されている間に、観測行列表示部１１５を介して伝送された対象信号を観測する（Ｓ１０９）。制御部１１２は、ステップＳ１０９で観測された対象信号の観測値から、ステップＳ１０４で観測された対象信号の観測値を減算することにより、差分値を算出する（Ｓ１１０）。通信部１１７は、ステップＳ１０６で生成された乱数シードと、ステップＳ１１０で算出された差分値とをセットにしたデータセットを信号復元装置１２０へ送信する（Ｓ１１１）。そして、観測回数カウンタｉの値がインクリメントされ（Ｓ１１２）、ステップＳ１０５へ戻る。
【０１３２】
なお、観測回数カウンタｉが観測回数Ｍ以上になるまで、ステップＳ１０５〜ステップＳ１１２が繰り返し行われることにより、第二の観測行列はＭ個生成されることになる。そして、Ｍ個のデータセットが信号復元装置１２０へ送信されることになる。
【０１３３】
このときステップＳ１０６を実行するたびに、異なる乱数シードが生成されるのであれば、生成されるＭ個の第二の観測行列は、互いに異なることが期待できる。
【０１３４】
以上が図５に示した信号観測装置１１０の信号観測処理の手順である。本実施の形態では、図５に示すように、複数回の観測のうち最初の観測で第一の制御が適用される。
【０１３５】
図７は、観測行列生成部１１４による観測行列生成処理の一例を示すフローチャートである。具体的には、図７は、図５のステップＳ１０７の詳細を示す。
【０１３６】
図７に示すように、観測行列生成部１１４は、入力部１１１から入力された対象信号の縦の画素数ＰｉｃＨ及び横の画素数ＰｉｃＷを用いて、第二の観測行列の行列要素の総数Ｎを算出する（Ｓ３０１）。次に、観測行列生成部１１４は、乱数シード生成部１１３が生成した乱数シードを用いて、ベルヌーイ分布に従う、０又は１の値を有するＮ個の乱数を生成する（Ｓ３０２）。そして、観測行列生成部１１４は、生成したＮ個の乱数をラスター順に並べて、ＰｉｃＷ行ＰｉｃＨ列の第二の観測行列を生成する（Ｓ３０３）。
【０１３７】
なお、ここでは、観測行列生成部１１４は、ベルヌーイ分布に従う乱数を生成したが、正規分布あるいはその他の確率分布に従う乱数を生成してもよい。また、乱数がとる値は、０又は１であったが、その他の値でもよい。
【０１３８】
以上が図７に示した、観測行列生成部１１４の観測行列生成処理の手順である。
【０１３９】
図８は、観測行列表示部１１５の制御処理の一例を示すフローチャートである。具体的には、図８は、図５のステップＳ１０３及びステップＳ１０８の詳細を示す。
【０１４０】
図８に示すように、制御部１１２は、入力部１１１から入力された観測行列表示部１１５の横の画素数ＤｅｖＷ及び対象信号の横の画素数ＰｉｃＷを用いて、観測行列の各行列要素の横方向の複製数（以後ＮｕｍＷともいう）を算出する。同様に、制御部１１２は、観測行列表示部１１５の縦の画素数ＤｅｖＨ及び対象信号の縦の画素数ＰｉｃＨを用いて、観測行列の各要素の縦方向の複製数（以後ＮｕｍＨともいう）を算出する（Ｓ４０１）。
【０１４１】
次に、制御部１１２は、処理対象の観測行列の要素の行番号を示すカウンタｊを０で初期化する（Ｓ４０２）。次に、制御部１１２は、カウンタｊが、対象信号の縦の画素数ＰｉｃＨより小さいか否かを判定する（Ｓ４０３）。
【０１４２】
ここで、カウンタｊが縦の画素数ＰｉｃＨ以上である場合（Ｓ４０３のＮ）、観測行列表示部１１５の制御処理を終了する。一方、カウンタｊが縦の画素数ＰｉｃＨより小さい場合は（Ｓ４０３のＹ）、処理対象の観測行列の要素の列番号を示すカウンタｉｉを０で初期化する（Ｓ４０４）。次に、カウンタｉｉが、対象信号の横の画素数ＰｉｃＷより小さいか否かを判定する（Ｓ４０５）。
【０１４３】
ここで、カウンタｉｉが横の画素数ＰｉｃＷ以上である場合（Ｓ４０５のＮ）、カウンタｊは、インクリメントされ（Ｓ４１２）、ステップＳ４０３に戻る。一方、カウンタｉｉが横の画素数ＰｉｃＷより小さい場合は（Ｓ４０５のＹ）、制御部１１２は、観測行列表示部１１５の画素の縦方向の複製数のためのカウンタｍを０で初期化する（Ｓ４０６）。次に、制御部１１２は、カウンタｍがＮｕｍＨよりも小さいか否かを判定する（Ｓ４０７）。
【０１４４】
ここで、カウンタｍがＮｕｍＨ以上である場合（Ｓ４０７のＮ）、カウンタｉｉがインクリメントされ（Ｓ４１３）、ステップＳ４０５に戻る。一方、カウンタｍがＮｕｍＨより小さい場合は（Ｓ４０７のＹ）、制御部１１２は、観測行列表示部１１５の画素の横方向の複製数のためのカウンタｎを０で初期化する（Ｓ４０８）。次に、制御部１１２は、カウンタｎがＮｕｍＷよりも小さいか否かを判定する（Ｓ４０９）。
【０１４５】
ここで、カウンタｎがＮｕｍＷ以上である場合（Ｓ４０９のＮ）、カウンタｍはインクリメントされ（Ｓ４１４）、ステップＳ４０７に戻る。一方、カウンタｎがＮｕｍＷより小さい場合は（Ｓ４０９のＹ）、観測行列（第一の観測行列又は第二の観測行列）のｊ行ｉｉ列目の要素の値を用いて、観測行列表示部１１５の(ｊ×ＮｕｍＨ＋ｍ)行(ｉｉ×ＮｕｍＷ＋ｎ)列目の画素（または素子の状態）を制御する（Ｓ４１０）。次に、カウンタｎがインクリメントされ（Ｓ４１１）、ステップＳ４０９に戻る。
【０１４６】
以上が図８に示した、観測行列表示部１１５の制御処理の手順である。以上の処理により、観測行列の第一の行数及び第一の列数と、観測行列表示部１１５の画素の第二の行数及び第二の列数との関係に応じて、図４Ａ〜図４Ｅを用いて説明をした方法を適用することで観測行列表示部１１５の各素子が制御される。
【０１４７】
図９Ａは、第二の観測行列に基づく観測行列表示部１１５の制御結果の一例を示す図である。図９Ｂは、第一の観測行列に基づく観測行列表示部１１５の制御結果の一例を示す図である。
【０１４８】
図９Ａ及び図９Ｂでは、観測行列（第一の観測行列及び第二の観測行列）は、６４行６４列の行列要素を有する。また、観測行列表示部１１５は、６４行６４列の画素を有する。つまり、観測行列表示部１１５を構成する複数の素子が６４行６４列の行列状に配列されている。
【０１４９】
したがって、観測行列の行列要素と、観測行列表示部１１５の画素とは、一対一で対応している。図９Ａ及び図９Ｂでは、行列要素の値が「０」であれば、観測行列表示部１１５の対応する画素は、制限状態であり、可視光を透過しないので黒で表されている。また、行列要素の値が「１」であれば、観測行列表示部１１５の対応する画素は、非制限状態であり、可視光を透過するので白色で表されている。
【０１５０】
なお、行列要素の値と画素の状態との対応関係は、これに限らない。例えば、「０」が非制限状態に対応し、「１」が制限状態に対応してもよい。
【０１５１】
図９Ｃは、第二の観測行列に基づく観測行列表示部１１５の制御結果の他の一例を示す図である。図９Ｃでは、観測行列（第一の観測行列及び第二の観測行列）は、図９Ａと同様に６４行６４列を有する。一方、観測行列表示部１１５は、７００行７００列の画素を有する。ここでは、観測行列表示部１１５の中心部の６４０ｘ６４０画素が、６４行６４列の観測行列の行列要素に対応付けられている。
【０１５２】
［信号復元処理］
次に、信号復元装置１２０における信号復元処理について、図１０〜図１２を参照して具体的に説明する。
【０１５３】
図１０は、信号復元装置１２０による信号復元処理の一例を示すフローチャートである。
【０１５４】
図１０に示すように、通信部１２３は、信号観測装置１１０から設定値を受信する（Ｓ５０１）。具体的には、通信部１２３は、対象信号の縦の画素数ＰｉｃＨ及び横の画素数ＰｉｃＷ、並びに、観測回数Ｍを信号観測装置１１０から受信する。次に、制御部１２２は、ＰｉｃＨ及びＰｉｃＷの設定を行う（Ｓ５０２）。さらに、制御部１２２は、信号観測装置１１０から受信した観測回数Ｍの設定を行う（Ｓ５０３）。
【０１５５】
そして、制御部１２２は、観測信号がスパース性を示す空間への射影のための基底行列φを入力部１２１から取得する（Ｓ５０４）。なお、基底行列φは、必ずしも入力部１２１から取得されなくてもよく、予め制御部１２２に保持された基底行列φを取得してもよい。基底行列φは、例えば、離散コサイン変換基底、フリーエ変換基底などである。
【０１５６】
次に、通信部１２３は、信号観測装置１１０からデータセットを受信する（Ｓ５０５）。次に、復元行列生成部１２４は、受信したデータセットに基づいて復元行列を生成する（Ｓ５０６）。そして、復元部１２５は、ステップＳ５０４で受信したデータセットとステップＳ５０６で生成した復元行列を用いて復元処理を実施する（Ｓ５０７）。
【０１５７】
具体的には、復元部１２５は、ｓ＝Ｆｅに基づいて対象信号ｅを復元する。ここで、ｓは、受信された複数のデータセットに含まれる複数の観測値からなる観測ベクトル（つまり、圧縮信号）である。また、Ｆは、復元行列である。
【０１５８】
ただし、対象信号ｅは必ずしもスパースな信号とは限らない。対象信号ｅがスパースではない場合は、対象信号ｅがスパース性を示すような空間に基底行列φを用いて射影することでスパース性を持たせる必要がある。具体的には、ｅ’＝φｅという計算で対象信号ｅをスパース性をもつｅ’に変換する。つまり、ｓ＝Ｆｅ＝（Ｆφ-1）ｅ’とすることで、スパース性を持たない対象信号ｅでも、基底行列φを用いることでスパース性を持つ信号ｅ’に変換し、スパース性を信号の復元に利用することが可能になる。ここで、φ-1は、基底行列φの逆行列である。このようにして、復元部１２５は、基底変換後の対象信号ｅ’のスパース性を利用して、ｓ＝Ｆｅ＝（Ｆφ-1）ｅ’を解くことにより対象信号ｅを復元する。
【０１５９】
以上が、図１０に示した信号復元装置１２０の信号復元処理の手順である。
【０１６０】
図１１は、通信部１２３によるデータセット受信処理の一例を示すフローチャートである。具体的には、図１１は、図１０のステップＳ５０５の詳細を示す。
【０１６１】
図１１に示すように、通信部１２３は、まず、信号観測装置１１０から受信したデータセットの数をカウントするためのカウンタｉを０で初期化する（Ｓ６０１）。次に、通信部１２３は、データセットを受信する（Ｓ６０２）。データセットは、信号観測装置１１０の観測部１１６が第二の観測行列を生成するために利用した乱数シード、及び、その第二の観測行列を用いて観測した観測値と第一の観測行列を用いて観測した観測値との差分値を含んでいる。
【０１６２】
そして、通信部１２３は、カウンタｉをインクリメントする（Ｓ６０３）。次に、通信部１２３は、前のデータセットを受信してから一定時間が経過したかどうかを判定する（Ｓ６０４）。ここでの一定時間は、データセット受信待ちのための閾値時間を表す。一定時間は、入力部１２１がユーザから入力を受けてもよいし、制御部１２２によって予め設定されていてもよい。
【０１６３】
一定時間が経過していなければ（Ｓ６０４のＮ）、受信すべきデータセットの個数と一致する観測回数Ｍが、カウンタｉよりも大きいか否かが判定される（Ｓ６０５）。ここで、Ｍがｉ以下である場合（Ｓ６０５のＮ）、又は、一定時間が経過している場合（Ｓ６０４のＹ）、通信部１２３は、実際に受信したデータセット数としてカウンタｉを設定する（Ｓ６０６）。Ｍがｉよりも大きい場合（Ｓ６０５のＹ）、ステップＳ６０２に戻る。
【０１６４】
以上が、通信部１２３のデータセット受信処理の手順である。
【０１６５】
図１２は、復元行列生成部１２４による復元行列生成処理の一例を示すフローチャートである。具体的には、図１２は、図１０のステップＳ５０６の詳細を示す。
【０１６６】
図１２に示すように、復元行列生成部１２４は、通信部１２３が受信した対象信号の縦の画素数ＰｉｃＨ及び横の画素数ＰｉｃＷから復元行列の列数Ｒを算出する（Ｓ７０１）。この例では、復元行列の列数ＲをＰｉｃＨ×ＰｉｃＷとする。そして、復元行列生成部１２４は、図１１のステップＳ６０５で設定されたデータセット数を、復元行列の行数Ｃとして設定する（Ｓ７０２）。
【０１６７】
次に、復元行列生成部１２４は、復元行列の行数をカウントするためのカウンタｉを０で初期化する（Ｓ７０３）。次に、復元行列生成部１２４は、通信部１２３がｉ番目に受信したデータセットに含まれる乱数シードを用いて、ベルヌーイ分布に従う、０又は１の値を有するＲ個の乱数を生成する（Ｓ７０４）。復元行列生成部１２４は、生成したＲ個の乱数を復元行列のｉ行目に採用する（Ｓ７０５）。そして、復元行列生成部１２４は、カウンタｉをインクリメントする（Ｓ７０６）。次に、復元行列生成部１２４は、カウンタｉが復元行列の行数Ｃよりも小さいかどうかを判定する（Ｓ７０７）。カウンタｉが復元行列の行数Ｃ以上であれば（Ｓ７０７のＮ）、復元行列生成処理を終了する。カウンタｉが復元行列の行数Ｃよりも小さければ（Ｓ７０７のＹ）、ステップＳ７０４に戻る。
【０１６８】
以上が、図１２に示した復元行列生成部１２４の復元行列生成処理の手順である。
【０１６９】
以上のような信号観測システム１００による実験結果を図１３Ａ及び図１３Ｂを参照して説明する。
【０１７０】
図１３Ａは、実験で用いた被写体９９を示す図である。図１３Ｂは、実験で復元された画像を示す図である。実験では、図１３Ａに示すように、被写体９９に日本語のかな文字「あ」が表されている。この被写体９９からの光を、図９Ｃのように制御された観測行列表示部１１５を介して観測部１１６が複数回観測した。複数の第二の観測行列を用いた複数回の観測によって得られた複数の観測値の各々と第一の観測行列を用いた観測によって得られた観測値との差分値を用いて、被写体９９を示す画像が図１３Ｂに示すように復元された。
【０１７１】
［効果］
以上のように、本実施の形態に係る信号観測システム１００によれば、予め定められたパターンの第一の観測行列を用いた第一の制御と、乱数シードに基づく第二の観測行列を用いた第二の制御とを選択的に用いて、複数回の観測を行うことができる。したがって、第一の制御によって得られた対象信号の観測値と、第二の制御によって得られた対象信号の観測値との差分値を求めることでき、その差分値を用いて対象信号を復元することが可能となる。その結果、複数回の観測における観測行列による観測値の変化を強調して対象信号の復元を行うことができ、ノイズによる信号の劣化を抑制することができる。
【０１７２】
また、本実施の形態に係る信号観測システム１００によれば、最初の観測で第一の制御を適用することができる。したがって、後続する第二の制御を用いた観測では、既に第一の制御を用いた観測値が得られているので、差分値をすぐに求めることができ、観測値のためのバッファを減少させることができる。
【０１７３】
また、本実施の形態に係る信号観測システム１００によれば、信号観測装置１１０において差分値を算出することができる。したがって、信号復元装置１２０において差分値を算出する必要がなく、対象信号の復元を容易にすることができる。さらに、信号観測装置１１０から信号復元装置１２０への通信量を削減することもできる。
【０１７４】
また、本実施の形態に係る信号観測システム１００によれば、入力部から得られた第一の行数及び第一の列数から決定されるサイズの複数の観測行列に基づいて複数の素子の各々の状態を変更することできる。つまり、観測行列表示部１１５の行数及び列数に固定されずに、観測行列のサイズを任意に設定することができる。対象信号の画素数は、観測行列のサイズに依存する。したがって、観測行列表示部１１５の画素数に依存せずに、対象信号の画素数を自由に設定することが可能となり、観測の自由度を向上させることができる。
【０１７５】
また、本実施の形態に係る信号観測システム１００によれば、観測行列表示部１１５の複数の素子をグループ化し、グループ単位で複数の素子の状態を変更することにより、観測行列の行列要素の数と観測行列表示部１１５の素子（画素）の数とが異なる場合であっても、観測行列に基づいて観測行列表示部１１５の各素子を制御することができる。その結果、任意の画素数で対象信号を復元することができ、観測の自由度を向上させることができる。
【０１７６】
また、本実施の形態に係る信号観測システム１００によれば、乱数シードと、当該乱数シードに基づいて生成された第２の観測行列に基づいて観測行列表示部１１５の各素子を制御している間に観測部１１６で観測された対象信号の量から第１の観測行列に基づいて観測行列表示部１１５の各素子を制御している間に観測部１１６で観測された対象信号の量を減算することで得られる差分値（圧縮信号）とをセットにしたデータセットを信号観測装置１１０から信号復元装置１２０へ送信することができる。また異なる乱数シードに基づいて第二の観測行列を生成することで行列要素の値の異なる第二の観測行列を複数生成することができる。乱数シードと差分値とをセットにしたデータセットを信号観測装置１１０から信号復元装置１２０へ複数、送信することができる。したがって、通信途中で一部のデータセットが失われた場合でも、受信できたデータセットのみを用いて復元行列を生成することができるので、エラー耐性を向上させることができる。
【０１７７】
また、本実施の形態に係る信号観測システム１００によれば、第二の観測行列をそのまま送信せずに、第二の観測行列を生成するために利用した乱数シードを送信することができる。したがって、観測行列のための通信量及び記憶容量を削減することができる。
【０１７８】
（実施の形態２）
次に、実施の形態２について説明する。本実施の形態では、予め定められたパターンの第一の観測行列に基づく第一の制御が複数回の観測のうちの真ん中の観測で適用される点が、実施の形態１と異なる。以下、本実施の形態について、実施の形態１と異なる点を中心に説明する。
【０１７９】
本実施の形態に係る信号観測システムの構成は、図２に示す実施の形態１に係る信号観測システムの構成と実質的に同一であるので、図示及び説明を省略する。
【０１８０】
［信号観測システムの動作］
本実施の形態に係る圧縮信号観測システムの動作を、図１４を参照して具体的に説明する。
【０１８１】
図１４は、実施の形態２における信号観測装置１１０による信号観測処理の一例を示すフローチャートである。図１４において、図５と同一又は類似の処理を行うステップには、同一の符号を付し、適宜説明を省略する。
【０１８２】
初期設定が行われた後（Ｓ１０１）、制御部１１２は、観測回数Ｍと観測回数カウンタｉとを比較する（Ｓ１０５）。ｉがＭより小さい場合（Ｓ１０５のＹ）、制御部１１２は、観測回数カウンタｉを観測回数Ｍの半分で除したときの余りが０であり、かつ、ｉが０より大きいか否かを判定する（Ｓ８０１）。つまり、制御部１１２は、ｉ回目の観測が、Ｍ回の観測のうちの真ん中の観測であるか否かを判定する。
【０１８３】
ここで、観測回数カウンタｉを観測回数Ｍの半分で除したときの余りが０であり、かつ、ｉが０より大きい場合（Ｓ８０１のＹ）、第一の観測行列を用いた観測のための処理が実行される（Ｓ１０２〜Ｓ１０４）。一方、観測回数カウンタｉを観測回数Ｍの半分で除したときの余りが０でない、又は、ｉが０である場合（Ｓ８０１のＮ）、第一の観測行列を用いた観測のための処理をスキップして、第二の観測行列を用いた観測のための処理が実行される（Ｓ１０６〜Ｓ１０９）。つまり、Ｍが偶数の場合には、Ｍ回の観測のうちＭ／２回目の観測の際に第一の観測行列を用いた観測が行われる。すなわち、複数回の観測のうち真ん中の観測で、第一の制御が適用される。その後、観測回数カウンタｉの値がインクリメントされ（Ｓ１１２）、ステップＳ１０５へ戻る。
【０１８４】
ｉがＭ以上である場合（Ｓ１０５のＮ）、制御部１１２は、ステップＳ１０９でＭ回観測された対象信号の観測値の各々から、ステップＳ１０４で観測された対象信号の観測値を減算することにより、Ｍ個の差分値を算出する（Ｓ１１０）。通信部１１７は、ステップＳ１０６で生成された乱数シードと、ステップＳ１１０で算出された差分値とをセットにしたＭ個のデータセットを信号復元装置１２０へ送信する（Ｓ１１１）。
【０１８５】
以上が図１４に示した信号観測装置１１０の信号観測処理の手順である。
【０１８６】
［効果］
以上、本実施の形態に係る信号観測装置１１０によれば、真ん中の観測で第一の制御を適用することができる。したがって、複数回の観測において、実施の形態１よりも時間的に近い観測値同士で差分値を求めることができ、経時変化に伴うノイズの影響を抑制することができる。
【０１８７】
（実施の形態３）
次に、実施の形態３について説明する。本実施の形態では、予め定められたパターンの第一の観測行列に基づく第一の制御が複数回の観測のうちの所定回毎の観測で適用される点が、実施の形態１と異なる。以下、本実施の形態について、実施の形態１と異なる点を中心に説明する。
【０１８８】
本実施の形態に係る信号観測システムの構成は、図２に示す実施の形態１に係る信号観測システムの構成と実質的に同一であるので、図示及び説明を省略する。
【０１８９】
［信号観測システムの動作］
本実施の形態に係る信号観測システムの動作を、図１５を参照して具体的に説明する。
【０１９０】
図１５は、実施の形態３に係る信号観測装置１１０の信号観測処理の一例を示すフローチャートである。図１５において、図５と同一又は類似の処理を行うステップには、同一の符号を付し、適宜説明を省略する。
【０１９１】
初期設定が行われた後（Ｓ１０１）、制御部１１２は、観測回数Ｍと観測回数カウンタｉとを比較する（Ｓ１０５）。ここで、ｉがＭ以上である場合（Ｓ１０５のＮ）、信号観測装置１１０は信号観測処理を終了する。一方、ｉがＭより小さい場合（Ｓ１０５のＹ）、制御部１１２は、観測回数カウンタｉを１０で除したときの余りが０であるか否かを判定する（Ｓ９０１）。つまり、制御部１１２は、ｉ回目の観測が、Ｍ回の観測のうちの所定回数（ここでは１０）周期の観測であるか否かを判定する。なお、本実施の形態では、所定回数を１０として説明するがこれに限定をされるものではない。また、所定回数は、Ｍよりも小さい。
【０１９２】
ここで、観測回数カウンタｉを１０で除したときの余りが０である場合（Ｓ９０１のＹ）、制御部１１２は、第一の観測行列を用いた観測のための処理を実行する（Ｓ１０２〜Ｓ１０４）。一方、観測回数カウンタｉを１０で除したときの余りが０でない場合（Ｓ９０１のＮ）、第一の観測行列を用いた観測のための処理をスキップして、第二の観測行列を用いた観測のための処理が実行される（Ｓ１０６〜Ｓ１０９）。つまり、複数回の観測のうち、所定回数（ここでは１０回）毎の観測で、第一の制御が適用される。
【０１９３】
その後、制御部１１２は、ステップＳ１０９で観測された対象信号の観測値から、ステップＳ１０４で観測された対象信号の観測値を減算することにより差分値を算出する（Ｓ１１０）。通信部１１７は、ステップＳ１０６で生成された乱数シードと、ステップＳ１１０で算出された差分値とをセットにしたデータセットを信号復元装置１２０へ送信する（Ｓ１１１）。その後、観測回数カウンタｉの値がインクリメントされ（Ｓ１１２）、ステップＳ１０５へ戻る。
【０１９４】
以上が図１５に示した信号観測装置１１０の信号観測処理の手順である。
【０１９５】
［効果］
以上、本実施の形態に係る信号観測装置１１０によれば、所定回数毎の観測で第一の制御を適用することができる。したがって、複数回の観測において、比較的時間的に近い観測値同士で差分値を求めることができ、経時変化に伴うノイズの影響を抑制することが可能となる。
【０１９６】
（実施の形態４）
次に、実施の形態４について説明する。本実施の形態では、複数回の観測において、第一の制御と第二の制御とが交互に適用される点が、実施の形態１と異なる。以下、本実施の形態について、実施の形態１と異なる点を中心に説明する。
【０１９７】
（構成）
本実施の形態に係る信号観測システムの構成は、図２に示す実施の形態１に係る信号観測システムの構成と実質的に同一であるので、図示及び説明を省略する。
【０１９８】
［信号観測システムの動作］
本実施の形態に係る信号観測システムの動作を、図１６を参照して具体的に説明する。
【０１９９】
図１６は、実施の形態４に係る信号観測装置１１０の信号観測処理の一例を示すフローチャートである。図１６において、図５と同一又は類似の処理を行うステップには、同一の符号を付し、適宜説明を省略する。
【０２００】
初期設定が行われた後（Ｓ１０１）、制御部１１２は、観測回数Ｍと観測回数カウンタｉとを比較する（Ｓ１０５）。ここで、ｉがＭ以上である場合（Ｓ１０５のＮ）、信号観測装置１１０は信号観測処理を終了する。一方、ｉがＭより小さい場合（Ｓ１０５のＹ）、制御部１１２は、第一の観測行列を用いた観測のための処理を実行する（Ｓ１０２〜Ｓ１０４）。その後、第二の観測行列を用いた観測のための処理が実行される（Ｓ１０６〜Ｓ１０９）。つまり、複数回の観測において、第一の制御及び第二の制御が交互に適用される。
【０２０１】
その後、制御部１１２は、ステップＳ１０９で観測された対象信号の観測値から、ステップＳ１０４で観測された対象信号の観測値を減算することにより差分値を算出する（Ｓ１１０）。通信部１１７は、ステップＳ１０６で生成された乱数シードと、ステップＳ１１０で算出された差分値とをセットにしたデータセットを信号復元装置１２０へ送信する（Ｓ１１１）。その後、観測回数カウンタｉの値がインクリメントされ（Ｓ１１２）、ステップＳ１０５へ戻る。
【０２０２】
以上が図１６に示した信号観測装置１１０の信号観測処理の手順である。
【０２０３】
［効果］
以上、本実施の形態に係る信号観測装置１１０によれば、第一の制御及び第二の制御を交互に適用することができる。したがって、複数回の観測において、時間的に隣り合う観測値同士で差分値を求めることができ、経時変化に伴うノイズの影響を抑制することが可能となる。
【０２０４】
（実施の形態５）
次に、実施の形態５について説明する。本実施の形態では、乱数シードの代わりに乱数の識別番号が信号観測装置から信号復元装置に送信される点が、実施の形態１と異なる。以下、本実施の形態について、実施の形態１と異なる点を中心に説明する。
【０２０５】
［信号観測システムの構成］
図１７は、実施の形態５における信号観測システム２００の構成の一例を示すブロック図である。信号観測システム２００は、信号観測装置２１０、及び、信号復元装置２２０を備える。図１７において、図１と実質的に同一の構成要素については同じ符号を付し、説明を適宜省略する。
【０２０６】
信号観測装置２１０は、圧縮センシングを用いて観測の対象となる対象信号を複数回観測し、観測結果として圧縮信号を出力する。
【０２０７】
信号観測装置２１０は、入力部１１１、制御部１１２、観測行列記憶部２１１、観測行列生成部２１２、観測行列表示部１１５、観測部１１６、及び、通信部１１７を備える。
【０２０８】
観測行列記憶部２１１は、ハードディスクドライブなどの記憶装置である。観測行列記憶部２１１は、乱数列を記憶している。また、観測行列記憶部２１１は、識別番号と、乱数列の読み出し開始位置を示すアドレスとの対応関係を示す情報を記憶している。
【０２０９】
図１８は、識別番号とアドレスとの対応関係を示す情報の一例を示す図である。図１８では、識別番号と一対一で観測行列記憶部２１１のアドレスが対応しており、それぞれのアドレスに対応した乱数が表されている。１つのアドレスは、１ビット単位の位置を示す。観測行列記憶部２１１には、例えば、ビッグエンディアンでデータが格納されている。例えば図１８では、識別番号「０」で１８個の乱数が読み出される場合は、アドレス「０」の最上位ビットから順に１８個の乱数が読み出される。つまり、「０１００１０００１１０１１０１０００」が読み出される。なお、ここでは、識別番号と対応するアドレスは８ビット単位であるが、１ビット単位でもその他の単位でもよい。
【０２１０】
観測行列生成部２１２は、入力部１１１へ入力された対象信号の画素数（行数及び列数）、及び、識別番号に基づき、観測行列記憶部２１１から複数の乱数を読み出し、第二の観測行列を生成する。また、観測行列生成部２１２は、実施の形態１の観測行列生成部１１４と同様に、予め定められたパターンの第一の観測行列を生成する。
【０２１１】
信号復元装置２２０は、信号観測装置２１０が出力した圧縮信号及び識別番号から対象信号を復元する。
【０２１２】
信号復元装置２２０は、入力部１２１、制御部１２２、通信部１２３、観測行列記憶部２２１、復元行列生成部２２２、復元部１２５、及び、出力部１２６を備える。
【０２１３】
観測行列記憶部２２１は、ハードディスクドライブなどの記憶装置である。観測行列記憶部２２１は、信号観測装置２１０の観測行列記憶部２１１と同じ乱数列を記憶している。また、観測行列記憶部２２１は、観測行列記憶部２１１と同様に、識別番号に対応付けて、乱数列の読み出し開始位置を示すアドレスを記憶している。つまり、同一の識別番号が用いられれば、観測行列記憶部２１１及び観測行列記憶部２２１から同一の乱数を読み出すことができる。
【０２１４】
例えば、信号観測装置２１０から受信した行列識別子が０であり、かつ、復元行列列数が１８列である場合は、図１８を参照して、アドレス０の最上位ビットから順に１８個の乱数が読み出される。つまり、「０１００１０００１１０１１０１０００」が読み出される。
【０２１５】
復元行列生成部２２２は、信号観測装置２１０から受信した対象信号の画素数、観測回数、及び、識別番号に基づいて、観測行列記憶部２２１から複数の乱数を読み出し、第二の観測行列を生成する。そして、復元行列生成部２２２は、実施の形態１の復元行列生成部１２４と同様に、複数の第二の観測行列を用いて、対象信号を復元するための復元行列を生成する。
【０２１６】
［圧縮信号観測システムの動作］
次に、以上のように構成された本実施の形態に係る圧縮信号観測システムの動作を図１９〜図２１を参照して具体的に説明する。
【０２１７】
図１９は、実施の形態５における信号観測装置２１０による信号観測処理の一例を示すフローチャートである。図１９において、図５と同一又は類似の処理を行うステップには、同一の符号を付し、適宜説明を省略する。
【０２１８】
図１９に示すように、観測行列生成部２１２は、図５における乱数シードの生成に代えて、識別番号を決定する（Ｓ１００１）。例えば、観測行列生成部２１２は、観測回数カウンタｉの値を識別番号と決定する。観測行列生成部２１２は識別番号に基づいて第二の観測行列を生成する（Ｓ１００２）。
【０２１９】
制御部１１２は、生成された第二の観測行列に基づき、観測行列表示部１１５の各素子を制御する（Ｓ１０８）。つまり、制御部１１２は、第二の制御を実行する。観測部１１６は、ステップＳ１０８における第二の制御が実行されている間に、観測行列表示部１１５を介して伝送される対象信号の量を観測する（Ｓ１０９）。制御部１１２は、ステップＳ１０９で観測された対象信号の観測値から、ステップＳ１０４で観測された対象信号の観測値を減算することにより、差分値を算出する（Ｓ１１０）。通信部１１７は、ステップＳ１００１で決定された識別番号と、ステップＳ１１０で算出された差分値とをセットにしたデータセットを信号復元装置２２０へ送信する（Ｓ１００３）。そして、観測回数カウンタｉの値がインクリメントされ（Ｓ１１２）、ステップＳ１０５へ戻る。
【０２２０】
以上が図１９に示した信号観測装置２１０の信号観測処理の手順である。
【０２２１】
図２０は、観測行列生成部２１２による観測行列生成処理の一例を示すフローチャートである。具体的には、図２０は、図１９のステップＳ１００２の詳細を示す。
【０２２２】
図２０に示すように、観測行列生成部２１２は、入力部１１１から入力された観測信号の縦の画素数ＰｉｃＨ及び横の画素数ＰｉｃＷを用いて、第二の観測行列の要素数Ｎを算出する（Ｓ１１０１）。次に、観測行列生成部２１２は、識別番号を用いて、観測行列記憶部２１１からＮ個の乱数を読み出す（Ｓ１１０２）。具体的には、観測行列生成部２１２は、識別番号に対応するアドレスから乱数の読み出しを開始し、Ｎ個の乱数を観測行列記憶部２１１から読み出す。そして、観測行列生成部２１２は、読み出されたＮ個の乱数をラスター順に並べて、ＰｉｃＷ行ＰｉｃＨ列の第二の観測行列を生成する（Ｓ１１０３）。
【０２２３】
以上が図２０に示した、観測行列生成部２１２の観測行列生成処理の手順である。
【０２２４】
図２１は、復元行列生成部２２２による復元行列生成処理の一例を示すフローチャートである。図２１は、図１０のステップＳ５０６の詳細を示す。図２１において、図１２と同一又は類似の処理を行うステップには、同一の符号を付し、適宜説明を省略する。
【０２２５】
図２１に示すように、カウンタｉが０で初期化された後（Ｓ７０３）、復元行列生成部２２２は、通信部１２３がｉ番目に受信した識別番号に対応する観測行列記憶部２２１のアドレスを導出する（Ｓ１２０１）。復元行列生成部２２２は、導出したアドレスから連続するアドレスに格納されたＲ個の乱数を読み出す（Ｓ１２０２）。復元行列生成部２２２は、読み出したＲ個の乱数を復元行列のｉ行目に採用する（Ｓ１２０３）。
【０２２６】
以上が、図２１に示した復元行列生成部２２２の復元行列生成処理の手順である。
【０２２７】
［効果］
本実施の形態に係る信号観測システム２００によれば、観測行列をそのまま送信せずに、観測行列を生成するための識別番号を送信することができる。したがって、観測行列のための通信量及び記憶容量を削減することができる。
【０２２８】
（他の実施の形態）
以上、１つ又は複数の態様に係る信号観測システムについて、実施の形態に基づいて説明したが、本開示は、この実施の形態に限定されるものではない。本開示の趣旨を逸脱しない限り、当業者が思いつく各種変形を本実施の形態に施したものも、１つ又は複数の態様の範囲内に含まれてもよい。
【０２２９】
例えば、上記各実施の形態において、第一の観測行列に基づく対象信号の観測値と、第二の観測行列に基づく対象信号の観測値との差分値は、信号観測装置において算出されていたが、これに限定されない。差分値の算出は、信号復元装置において行われてもよい。
【０２３０】
また、上記各実施の形態において、観測行列表示部１１５は、透過型液晶ディスプレイであったが、これに限定されない。例えば、観測行列表示部１１５は、デジタルミラーデバイス又はＬＣＯＳ（Liquid Crystal On Silicon）と、レンズ、ミラー又はプリズムと
の組合せであってもよい。また、観測行列表示部１１５は、シャッターアレイであってもよい。例えば、観測行列表示部１１５がデジタルミラーデバイスあるいはＬＣＯＳによって実現される場合は、制御部１１２は、観測行列の行列要素の値が０であれば、観測部１１６へ対象信号が到達しないように対応する画素を制御し、観測行列の行列要素の値が１であれば、観測部１１６へ対象信号が到達するように対応する画素を制御すればよい。
【０２３１】
また、上記各実施の形態において、予め定められたパターンの第一の観測行列として、全ての行列要素の値が０である行列が用いられていたが、これに限定されない。例えば、第一の観測行列として、全ての行列要素の値が１である行列が用いられてもよい。また、第一の観測行列は、値が一様な行列要素を有する行列に限定されない。例えば、第一の観測行列は、０及び１が交互に並ぶ行列であってもよい。
【０２３２】
また、上記各実施の形態において、観測行列の各行列要素は、０又は１の値を有していたが、これに限定されない。例えば、観測行列の各行列要素の値は、整数でなくてもよく、小数であってもよい。この場合、各行列要素の値に対応する透過率で、観測行列表示部１１５の各画素が制御されればよい。
【０２３３】
また、上記各実施の形態において、観測部１１６は、単一の光センサで構成されていたが、これに限定されない。例えば、対象信号がテラヘルツ波信号であれば、観測部１１６は、テラヘルツ波を検知するセンサで構成されてもよい。また、対象信号が赤外線信号であれば、観測部１１６は、赤外線センサで構成されてもよい。また、対象信号が紫外線信号であれば、観測部１１６は、紫外線センサで構成されてもよい。つまり、観測部１１６は、対象信号の種類に適したセンサであればどのようなセンサで構成されてもよい。また、観測部１１６は、複数の光センサで構成されてもよい。また、観測部１１６は、複数種類のセンサ（例えば、光センサ、赤外線センサ、テラヘルツ波センサ、紫外線センサなど）の任意の組み合わせで構成されてもよい。
【０２３４】
また、本開示は、さらに、以下のような場合も含まれる。
【０２３５】
（１）上記の装置は、具体的には、マイクロプロセッサ、ＲＯＭ、ＲＡＭ、ハードディスクユニット、ディスプレイユニット、キーボード、マウスなどから構成されるコンピュータシステムである。前記ＲＡＭまたはハードディスクユニットには、コンピュータプログラムが記憶されている。前記マイクロプロセッサが、前記コンピュータプログラムにしたがって動作することにより、各装置は、その機能を達成する。ここでコンピュータプログラムは、所定の機能を達成するために、コンピュータに対する指令を示す命令コードが複数個組み合わされて構成されたものである。
【０２３６】
（２）上記の装置を構成する構成要素の一部または全部は、１個のシステムＬＳＩ（Large Scale Integration：大規模集積回路）から構成されているとしてもよい。システムＬＳＩは、複数の構成部を１個のチップ上に集積して製造された超多機能ＬＳＩであり、具体的には、マイクロプロセッサ、ＲＯＭ、ＲＡＭなどを含んで構成されるコンピュータシステムである。前記ＲＡＭには、コンピュータプログラムが記憶されている。前記マイクロプロセッサが、前記コンピュータプログラムにしたがって動作することにより、システムＬＳＩは、その機能を達成する。
【０２３７】
（３）上記の装置を構成する構成要素の一部または全部は、各装置に脱着可能なＩＣカードまたは単体のモジュールから構成されているとしてもよい。前記ＩＣカードまたは前記モジュールは、マイクロプロセッサ、ＲＯＭ、ＲＡＭなどから構成されるコンピュータシステムである。前記ＩＣカードまたは前記モジュールは、上記の超多機能ＬＳＩを含むとしてもよい。マイクロプロセッサが、コンピュータプログラムにしたがって動作することにより、前記ＩＣカードまたは前記モジュールは、その機能を達成する。このＩＣカードまたはこのモジュールは、耐タンパ性を有するとしてもよい。
【０２３８】
（４）また、本開示は、上記に示す方法であるとしてもよい。また、これらの方法をコンピュータにより実現するコンピュータプログラムであるとしてもよいし、前記コンピュータプログラムからなるデジタル信号であるとしてもよい。
【０２３９】
（５）また、本開示は、前記コンピュータプログラムまたは前記デジタル信号をコンピュータで読み取り可能な記録媒体、例えば、フレキシブルディスク、ハードディスク、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤ、ＤＶＤ−ＲＯＭ、ＤＶＤ−ＲＡＭ、ＢＤ（Ｂｌｕ−ｒａｙ（登録商標）  Ｄｉｓｃ）、半導体メモリなどに記録したものとしてもよい。また、これらの記録媒体に記録されている前記デジタル信号であるとしてもよい。
【０２４０】
また、本開示は、前記コンピュータプログラムまたは前記デジタル信号を、電気通信回線、無線または有線通信回線、インターネットを代表とするネットワーク、データ放送等を経由して伝送するものとしてもよい。
【０２４１】
また、本開示は、マイクロプロセッサとメモリを備えたコンピュータシステムであって、前記メモリは、上記コンピュータプログラムを記憶しており、前記マイクロプロセッサは、前記コンピュータプログラムにしたがって動作するとしてもよい。
【０２４２】
また、前記プログラムまたは前記デジタル信号を前記記録媒体に記録して移送することにより、または前記プログラムまたは前記デジタル信号を、前記ネットワーク等を経由して移送することにより、独立した他のコンピュータシステムにより実施するとしてもよい。
【産業上の利用可能性】
【０２４３】
本発明は、圧縮センシングを用いた信号観測装置として利用することができる。
【符号の説明】
【０２４４】
１００、２００  信号観測システム
１１０、２１０  信号観測装置
１１１、１２１  入力部
１１２、１２２  制御部
１１３  乱数シード生成部
１１４、２１２  観測行列生成部
１１５  観測行列表示部
１１６  観測部
１１７、１２３  通信部
１２０、２２０  信号復元装置
１２４、２２２  復元行列生成部
１２５  復元部
１２６  出力部
２１１  観測行列記憶部
２２１  観測行列記憶部
【図１】
0006726578.tif 000002
【図２】
0006726578.tif 000003
【図３】
0006726578.tif 000004
【図４Ａ】
0006726578.tif 000005
【図４Ｂ】
0006726578.tif 000006
【図４Ｃ】
0006726578.tif 000007
【図４Ｄ】
0006726578.tif 000008
【図４Ｅ】
0006726578.tif 000009
【図５】
0006726578.tif 000010
【図６】
0006726578.tif 000011
【図７】
0006726578.tif 000012
【図８】
0006726578.tif 000013
【図９Ａ】
0006726578.tif 000014
【図９Ｂ】
0006726578.tif 000015
【図９Ｃ】
0006726578.tif 000016
【図１０】
0006726578.tif 000017
【図１１】
0006726578.tif 000018
【図１２】
0006726578.tif 000019
【図１３Ａ】
0006726578.tif 000020
【図１３Ｂ】
0006726578.tif 000021
【図１４】
0006726578.tif 000022
【図１５】
0006726578.tif 000023
【図１６】
0006726578.tif 000024
【図１７】
0006726578.tif 000025
【図１８】
0006726578.tif 000026
【図１９】
0006726578.tif 000027
【図２０】
0006726578.tif 000028
【図２１】
0006726578.tif 000029