JP6291189B2 - ハイパースペクトル画像の分解能を制御する方法 - Google Patents

ハイパースペクトル画像の分解能を制御する方法 Download PDF

Info

Publication number
JP6291189B2
JP6291189B2 JP2013172769A JP2013172769A JP6291189B2 JP 6291189 B2 JP6291189 B2 JP 6291189B2 JP 2013172769 A JP2013172769 A JP 2013172769A JP 2013172769 A JP2013172769 A JP 2013172769A JP 6291189 B2 JP6291189 B2 JP 6291189B2
Authority
JP
Japan
Prior art keywords
hyperspectral
pixel
window
image
pixels
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2013172769A
Other languages
English (en)
Other versions
JP2014049126A (ja
JP2014049126A5 (ja
Inventor
エリック・ダニエル・ビューラー
コンラッド・ロバート・クチンスキー
ベンジャミン・トーマス・オチピンティ
Original Assignee
ジーイー・アビエイション・システムズ・エルエルシー
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ジーイー・アビエイション・システムズ・エルエルシー filed Critical ジーイー・アビエイション・システムズ・エルエルシー
Publication of JP2014049126A publication Critical patent/JP2014049126A/ja
Publication of JP2014049126A5 publication Critical patent/JP2014049126A5/ja
Application granted granted Critical
Publication of JP6291189B2 publication Critical patent/JP6291189B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04NPICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/10Circuitry of solid-state image sensors [SSIS]; Control thereof for transforming different wavelengths into image signals
    • H04N25/11Arrangement of colour filter arrays [CFA]; Filter mosaics
    • H04N25/13Arrangement of colour filter arrays [CFA]; Filter mosaics characterised by the spectral characteristics of the filter elements
    • H04N25/135Arrangement of colour filter arrays [CFA]; Filter mosaics characterised by the spectral characteristics of the filter elements based on four or more different wavelength filter elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/2823Imaging spectrometer
    • G01J2003/2826Multispectral imaging, e.g. filter imaging
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10032Satellite or aerial image; Remote sensing
    • G06T2207/10036Multispectral image; Hyperspectral image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20021Dividing image into blocks, subimages or windows
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20212Image combination
    • G06T2207/20216Image averaging

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Image Processing (AREA)
  • Spectrometry And Color Measurement (AREA)

Description

本発明は、ハイパースペクトル画像の分解能を制御する方法に関する。
従来型およびハイパースペクトルの両方のイメージングシステムは、複数の相違するスペクトル帯域で電磁エネルギーを収集かつ記録する。結果として得られるイメージ(imagery)は、スペクトル帯域情報を1つまたは多くのチャネルに組み合わせて画像のグレースケールまたは色表現を形成することにより表示される。ハイパースペクトルイメージ(HSI)デバイスは、画像センサ上で画素と呼ばれる多数の空間画像素子で同時に多くの離散的なスペクトル帯域または色でエネルギーを記録する分光計の部類である。標準的な広帯域撮像素子は、広範囲のスペクトルにわたってすべての検出される入射エネルギーに対して各々の画素で1つの値を記録し、検出器の2次元配列から2つの空間次元での画像を創出する。HSIデバイスは、追加的なスペクトル次元を伴う画像を創出することにより標準的な広帯域撮像素子とは異なる。各々のHSI画素は、10から数百の波長値が記録される場合がある。
米国特許第7979209号明細書
本発明の一態様は、画素、および各々の画素内部で部分画素を規定する少なくとも1つのフィルタを有する画像センサからのハイパースペクトル画像の分解能を制御する方法に関する。方法は、部分画素の行および列の配列を伴う画像センサ上に窓を規定するステップ、ハイパースペクトル画像の1つまたは複数の既定のパラメータに基づいて、既定のパラメータに対する配列に対する重み付き平均に対する値を確定するために、窓内部の部分画素を重み付けするステップ、既定の数の行または列だけ窓をシフトするステップであって、既定の数が配列内のそれぞれの行または列の数未満であるステップ、画像センサ上のすべての可能な窓に対して重み付けするステップおよびシフトするステップを反復するステップ、ならびに、重み付き平均に基づいてハイパースペクトル画像を処理するステップを含む。
ハイパースペクトルデータキューブからのハイパースペクトル画像の典型的な既知のサンプリングを明示する代表的な例図である。 本発明の実施形態によるハイパースペクトルデータキューブから導出されるハイパースペクトル画像をアップサンプリングする方法を明示するフローチャートである。 図2の方法でのハイパースペクトルデータキューブから導出されるハイパースペクトル画像の第1の画素のアップサンプリングを明示する代表的な例図である。 図2の方法でのハイパースペクトルデータキューブから導出されるハイパースペクトル画像の第2の画素のアップサンプリングを明示する代表的な例図である。 図2の方法でのハイパースペクトルデータキューブから導出されるハイパースペクトル画像の第3の画素のアップサンプリングを明示する代表的な例図である。 図2の方法でのハイパースペクトルデータキューブから導出されるハイパースペクトル画像の第4の画素のアップサンプリングを明示する代表的な例図である。 図2の方法でのハイパースペクトルデータキューブから導出されるハイパースペクトル画像の完全なアップサンプリングを明示する代表的な例図である。
背景技術および以下の説明において説明目的で、数多くの特定の詳細を、本明細書で説明する技術の徹底した理解をもたらすために記載する。しかしながらこれらの特定の詳細がなくとも、例示的な実施形態が実践可能であることは当業者には明白であろう。他の事例では、構造およびデバイスを、例示的な実施形態の説明を容易にするために線図の形式で示す。
例示的な実施形態を、図面を参照して説明する。これらの図面は、本明細書で説明するモジュール、方法、またはコンピュータプログラム製品を実装する特定の実施形態の決まった詳細を図示する。しかしながら図面を、図面に存在し得る何らかの制限を加えるものであると解釈すべきではない。方法およびコンピュータプログラム製品を、それらの動作を達成するために任意の機械可読媒体に付して提供することができる。実施形態を、既存のコンピュータプロセッサを使用して、またはこのもしくは別の目的で組み込まれる専用コンピュータプロセッサにより、またはハードワイヤードのシステムにより実装することができる。
上記で指摘したように、本明細書で説明する実施形態は、機械実行可能命令もしくはそこで記憶されたデータ構造体を搬送するまたは有するための、機械可読媒体を備えるコンピュータプログラム製品を含み得る。そのような機械可読媒体は、プロセッサを伴う汎用または専用のコンピュータまたは他の機械によりアクセスされ得る任意の利用可能な媒体であり得る。例としてそのような機械可読媒体は、RAM、ROM、EPROM、EEPROM、CD−ROMもしくは他の光ディスク記憶装置、磁気ディスク記憶装置もしくは他の磁気記憶デバイス、または、所望のプログラムコードを機械実行可能命令もしくはデータ構造体の形式で搬送もしくは記憶するために使用され得る、および、プロセッサを伴う汎用もしくは専用のコンピュータもしくは他の機械によりアクセスされ得る、任意の他の媒体を備え得る。情報をネットワークまたは別の通信接続(ハードワイヤード、ワイヤレス、またはハードワイヤードもしくはワイヤレスの組み合わせのいずれか)を介して機械に伝送または提供するとき、機械は当然のことながらその接続を機械可読媒体とみなす。したがって任意のそのような接続を、当然のことながら機械可読媒体と呼ぶ。上記の組み合わせもまた、機械可読媒体の範囲内に含まれる。機械実行可能命令は例えば、汎用コンピュータ、専用コンピュータ、または専用処理機械に、決まった1つの機能または機能の群を遂行させる命令およびデータを含む。
実施形態を、例えばネットワーク化された環境での機械により実行されるプログラムモジュールの形式でのプログラムコードなどの機械実行可能命令を含むプログラム製品により一実施形態において実装され得る、方法ステップの全体的な脈絡で説明する。一般にプログラムモジュールは、個別のタスクを遂行する技術的効果を有する、または個別の抽象データ型を実装する、ルーチン、プログラム、オブジェクト、コンポーネント、データ構造体等を含む。機械実行可能命令、関連するデータ構造体、およびプログラムモジュールは、本明細書で開示する方法のステップを実行するためのプログラムコードの例を表す。そのような実行可能命令の個別のシーケンスまたは関連するデータ構造体は、そのようなステップで説明する機能を実装するための対応する作用の例を表す。
実施形態を、プロセッサを有する1つまたは複数のリモートコンピュータへの論理接続を使用して、ネットワーク化された環境で実践することができる。論理接続は、ここでは例としてであり限定としてではなく提示される、ローカルエリアネットワーク(LAN)および広域ネットワーク(WAN)を含み得る。そのようなネットワーク化環境は、職場内または企業内のコンピュータネットワーク、イントラネットおよびインターネットではありふれたものであり、多種多様の異なる通信プロトコルを使用することができる。そのようなネットワークコンピューティング環境が典型的には、パーソナルコンピュータ、ハンドヘルドデバイス、マルチプロセッサシステム、マイクロプロセッサベースまたはプログラマブルの民生用電子機器、ネットワークPC、ミニコンピュータ、メインフレームコンピュータ等を含む多くのタイプのコンピュータシステム構成を包含することになることを当業者ならば理解するであろう。
実施形態を、通信ネットワークを介して(ハードワイヤードリンク、ワイヤレスリンクによって、またはハードワイヤードもしくはワイヤレスのリンクの組み合わせによってのいずれかで)リンクされる、ローカルおよびリモートの処理デバイスによりタスクが遂行される分散コンピューティング環境でさらに実践することができる。分散コンピューティング環境ではプログラムモジュールは、ローカルおよびリモートの両方のメモリ記憶デバイスに位置し得る。
例示的な実施形態の全体または一部分を実装するための例示的なシステムは、処理装置、システムメモリ、および、システムメモリを含む様々なシステム構成要素を処理装置に結合するシステムバスを含む、コンピュータの形式での汎用コンピューティングデバイスを含み得る。システムメモリは、リードオンリーメモリ(ROM)およびランダムアクセスメモリ(RAM)を含み得る。コンピュータは、磁気ハードディスクに関する読み出しおよび書き込みを行うための磁気ハードディスクドライブ、リムーバブル磁気ディスクに関する読み出しまたは書き込みを行うための磁気ディスクドライブ、ならびに、CD−ROMまたは他の光学媒体などのリムーバブル光ディスクに関する読み出しまたは書き込みを行うための光ディスクドライブをさらに含み得る。ドライブおよびそれらの関連する機械可読媒体は、機械実行可能命令、データ構造体、プログラムモジュール、およびコンピュータ用の他のデータの不揮発性の記憶をもたらす。
実施形態で開示する方法の技術的効果には、特に物体を検出かつ追跡する方法がセンサと連関して使用されるときに、センサの範囲を向上させることと直接の相関関係がある、既存の情報から画像の空間分解能を向上させることがある。さらにこの方法は、システムにより収集されるデータをより効率的に使用することにより、既存のスペクトル複合画像技法を改善する。この技法を、スペクトルキューブ配列から複合イメージを生成する任意のシステムで使用することができる。スペクトルキューブ配列システムでこの技法を実装することによって、スペクトルキューブの大きさに等しく、ある倍数だけセンサの空間分解能が向上することになる。
図1は、ハイパースペクトルデータキューブ12からのハイパースペクトル画像14のパンクロマチック表現の既知のサンプリングを明示する代表的な例図である。ハイパースペクトルデータキューブ12から形成されるハイパースペクトル画像14のパンクロマチック表現は、HSIデバイスの画素から収集されるスペクトル帯域のすべてからのデータを単一の配列に組み合わせることにより形成される広帯域画像である。ハイパースペクトルデータキューブ12は、HSIデバイスの画像センサ上の画素からの値の配列からなり、ハイパースペクトルデータキューブ12の22などの各々の画素は、16、18、20などの、この実施形態では総計9個となり、番号1〜9と識別される部分画素の配列からなる。画素22内の16、18、20などの各々の部分画素は、特定のスペクトル帯域を表し、単一の画素22の内部に包含される。特定の用途に応じて各々の画素に、より多い、またはより少ない部分画素が存在し得ることが理解されよう。部分画素により表されるスペクトル帯域は、ハイパースペクトルデータキューブ12の9個のスペクトル帯域である場合があり、またはSurface Optics Corporation製のSOC710などのHSIデバイスの16個のスペクトル帯域である場合があり、または周知のベイヤーフィルタの赤、青および緑の色に対応する3つの相違するスペクトル帯域である場合があるが、部分画素により表されるスペクトル帯域が画素内の複数の部分画素に対して同じである構成を含む特定の実装形態の要求に応じて、他のスペクトル帯域フィルタ構成が可能である。本願の好適な実施例において、少なくとも1つのフィルタは、マイクロレンズの配列のレンズレットを備えるハイパースペクトルフィルタである。
ハイパースペクトルデータキューブ12は、所与のハイパースペクトル画素22の、この実施形態では番号1〜9のすべての部分画素の値を単一のパンクロマチック画素28の値に組み合わせる変換方法24により、ハイパースペクトル画像14のパンクロマチック表現に変換される。この実施形態では画素22内の9個の部分画素は、1から9で番号付けされた一意的な波長での強度を表す値を各々包含し、変換方法24によって組み合わされる。各々のハイパースペクトル画素22の部分画素値のすべてを組み合わせるために使用される変換方法24は、ハイパースペクトルデータキューブ12内のすべてのハイパースペクトル画素に対して反復され、この実施形態では各々の画素は、波長5に対応する部分画素(16、26、44、46、48、50、52、54)を中心とし、波長5に対応する部分画素および8個の最近傍の部分画素からなる。
単一のパンクロマチック画素28の値を形成するために使用される1つの変換方法24は、ハイパースペクトル画素22の部分画素のすべての9個の値を平均化することである。所定の変換方法24を使用してハイパースペクトルデータキューブ12内のデータをサンプリングすることにより、ハイパースペクトル画像14は、ハイパースペクトルデータキューブ12と同じ分解能または画素の数を有することになる。波長5に対応する部分画素(16、26、27、44、46、48、50、52、54)を中心とするハイパースペクトルデータキューブ12の各々の画素は、28および30などの単一のパンクロマチック画素を生成するためにサンプリングされることになる。
図2は、本発明の実施形態によるハイパースペクトルデータキューブから導出されるハイパースペクトル画像の分解能を制御するアップサンプリング方法80を明示するフローチャートである。方法80を用いると、ハイパースペクトル画像の各々のパンクロマチック画素を、ハイパースペクトルデータキューブの画素の境界は考慮せずに、ハイパースペクトルデータキューブの部分画素のサンプリングされた部分集合により生成することができる。方法80は、窓を規定するステップ82、所与の窓位置に対して窓内の部分画素値を重み付けするステップ86、重み付き平均を出力ハイパースペクトル画像に割り当てるステップ88、および、第1の窓位置で開始し(84)、すべての窓位置が処理された(92)状態まで、連続する窓位置のはじめから終わりまでインクリメントして(90)、窓位置の集合のはじめから終わりまで繰り返すことにより、窓が設定された部分画素を反復して重み付けかつシフトするステップからなる場合がある。
アップサンプリング方法80の第1のステップは、ハイパースペクトルデータキューブのどの部分画素が、ハイパースペクトル画像の結果として得られる画素に対する値を計算するために使用可能であるかを決定するために、画像処理ではマスクまたはテンプレートとしても既知である窓を規定すること82であり得る。窓は、含まれる部分画素の数および位置の観点から画素と同じサイズであり得る。しかしながら、任意のサイズの窓を実装形態に応じて使用することができる。窓は、画素に包含されるより多いまたは少ない部分画素を含む場合があり、窓は、画素を規定する部分画素の集合の形状とまったく異なる形状である場合がある。例えば部分画素の長方形の配置構成からなる窓を、部分画素の正方形の配置構成を用いる画素からなるハイパースペクトルデータキューブに適用する場合がある。
窓を第1の窓位置に適用し(84)、続いて窓内の部分画素値を重み付けする(86)ことができる。各々の窓が設定された位置で、重み付けを、窓が設定された部分画素の各々の値に適用することができる。数学的にはこのことは、次式のように表すことができる。
ただしzは結果として得られるハイパースペクトル画像の画素値であり、yは窓内のハイパースペクトルデータキューブの部分画素の値であり、wは部分画素に適用されるべき重みであり、Pは窓内の部分画素の総数である。その最も単純な形式では、平均化演算に影響を与える等しい重み付けを適用することができる、すなわち、すべてのwpの重みが1/Pの値に設定される。
重み付けを、結果として得られる画素の値に対する各々の部分画素の相対的寄与を修正するために使用することができ、その場合重みの集合は、既定のパラメータに基づく場合がある。既定のパラメータを、結果として得られる画像の対象の特徴を有利に表示するように選択することができる。本願の好適な実施例において、少なくとも1つの規定のパラメータは、赤、緑および青のスペクトル帯域のうちの少なくとも1つまたは全ての強度値である。
他の重みを適用することができ、窓内のすべての重みの集合は典型的には総計が1になり得るが、本発明はその制約には限定されない。太陽放射の周波数応答、および地球の大気内の吸収帯に起因する選択的な波長のさらなる減衰を補償するように、重みの好ましい集合が導出される。重み付けの体系を、結果として得られる画像の対象の特徴を有利に表示するように設計することができる。
重み付けされた、窓が設定された部分画素の値を、部分画素の単一の重み付き平均に組み合わせることができる。次いで重み付き平均を、現在の窓位置に対応する位置での、結果として得られるハイパースペクトル画像内の画素に割り当てることができる(88)。
窓内の部分画素値を重み付けするステップ86、および、重み付き平均を結果として得られるハイパースペクトル画像内の出力画素に割り当てるステップ88の処理を、すべての位置が処理された(92)状態まで、すべての可能な窓位置に対して反復することができる。順次的にはこの処理は、既定の数の行または列で窓をシフトするステップに相当し得る。本発明の好ましい実施形態では、垂直または水平のいずれかの方向に1つの部分画素だけ窓をシフトすることができるが、結果として得られるハイパースペクトル画像の分解能を制御するために他のシフトを使用することができる。次いで重み付け処理を反復して、新しい窓が設定された近隣の部分画素を使用してハイパースペクトル画像内の次の画素を形成することができる。本発明の別の実施形態では、同時の窓について、それらが垂直または水平のいずれかの方向に1つの部分画素だけシフトされる窓の集合であり、ハイパースペクトル画像のすべての結果として得られる出力画素を同時に生成することが可能であるように、例示する場合がある。
例として図3は、図2の方法によるハイパースペクトルデータキューブ12から導出されるハイパースペクトル画像114の第1の画素128のアップサンプリング100を明示する代表的な例図である。ハイパースペクトルデータキューブ12から導出されるハイパースペクトル画像114の第1の画素128のアップサンプリング100を明示する変換方法124は、ハイパースペクトルデータキューブ12のどの部分画素が、ハイパースペクトル画像114の結果として得られる画素128に対する値を計算するために使用可能であるかを決定するために、窓22を使用することができる。
次いで、既定の数の行または列だけ窓22をシフトすることができる。本発明の一実施形態では、垂直または水平のいずれかの方向に1つの部分画素だけ窓をシフトすることができるが、結果として得られるハイパースペクトル画像114の分解能を制御するために他のシフトを使用することができる。図4〜6は、単一の部分画素のシフトを用いて窓22の反復されるシフトを使用することによる、ハイパースペクトル画像114の画素226、328、430のアップサンプリングを明示する代表的な例図である。
図4は、図3の実施形態によるハイパースペクトルデータキューブ12から導出されるハイパースペクトル画像114の第2の画素226のアップサンプリング200を明示する代表的な例図である。ハイパースペクトルデータキューブ12から導出されるハイパースペクトル画像114の第2の画素226のアップサンプリング200を明示する変換方法224は、ハイパースペクトルデータキューブ12の部分画素216を中心とするどの部分画素が、ハイパースペクトル画像114の結果として得られる第2の画素226に対する値を計算するために使用可能であるかを決定するために、図3での窓位置から右に1列シフトされた窓22を使用することができる。
図5は、図3および4の実施形態によるハイパースペクトルデータキューブ12から導出されるハイパースペクトル画像114の第3の画素328のアップサンプリング300を明示する代表的な例図である。ハイパースペクトルデータキューブ12から導出されるハイパースペクトル画像114の第3の画素328のアップサンプリング300を明示する変換方法324は、ハイパースペクトルデータキューブ12のどの部分画素が、ハイパースペクトル画像114の結果として得られる第3の画素328に対する値を計算するために使用可能であるかを決定するために、図4での窓位置から右に1列シフトされた窓22を使用することができる。
図6は、図3〜5の実施形態によるハイパースペクトルデータキューブ12から導出されるハイパースペクトル画像114の第4の画素430のアップサンプリング400を明示する代表的な例図である。ハイパースペクトルデータキューブ12から導出されるハイパースペクトル画像114の第4の画素430のアップサンプリング400を明示する変換方法424は、ハイパースペクトルデータキューブ12のどの部分画素が、ハイパースペクトル画像114の結果として得られる第4の画素430に対する値を計算するために変換方法424により使用可能であるかを決定するために、図3での窓位置から下に1行シフトされた窓22を使用することができる。
図7は、図2の方法によるハイパースペクトルデータキューブ12から導出されるハイパースペクトル画像114の完全なアップサンプリング500を明示する代表的な例図である。ハイパースペクトル画像114のパンクロマチック表現を、ハイパースペクトルデータキューブ12の部分画素のすべての可能な配列に対して重み付けするステップおよびシフトするステップを反復することにより形成することができる。ハイパースペクトルデータキューブ12から導出されるハイパースペクトル画像114を創出するためのアップサンプリング500を明示する変換方法524は、図3〜6に示したように、ハイパースペクトルデータキューブ12のどの部分画素が、ハイパースペクトル画像114の各々の結果として得られる画素に対する値を計算するために使用可能であるかを決定するために、1度に1つの部分画素だけシフトされる窓を反復して使用することができる。図3〜6に示した重み付けされた、窓が設定された平均からの、結果として得られる画素128、226、328、430は、ハイパースペクトル画像114に対する画素のすべてを構築するために使用される、重み付けされた、窓が設定された平均のすべてに対しての最終的なハイパースペクトル画像114に配置されるように示される。
本発明の分解能制御方法の実施形態を使用してハイパースペクトルデータキューブ12から形成されるハイパースペクトル画像114のパンクロマチック表現は、518、520などの境界のハイパースペクトルデータキューブの部分画素を除外することが可能な状態で、ハイパースペクトルデータキューブの部分画素とハイパースペクトル画像の画素との間の1対1の写像を確定することができる。境界のハイパースペクトルデータキューブの部分画素518、520は、結果として得られるハイパースペクトル画像114内に類似物を有さない場合があり、そのことによって、ハイパースペクトル画像に対する最終的な分解能が、MおよびNがそれぞれ1行および1列内のハイパースペクトルデータキューブの部分画素の数である場合にM−2×N−2個の画素として確定し、スペクトルキューブの大きさに等しく、ある倍数だけセンサの空間分解能が実質的に向上する。
単一の重み付けの体系を使用することで、グレースケールとして表示される画像をもたらすことができる。複数の異なる重み付けの体系を、結果として得られる画像の相違するチャネルを生成するために同時に用いることができる。例えばハイパースペクトル画像114のパンクロマチック表現の表示は、複数の色チャネルを有し得る。典型的にはRGBとして表現される結果として得られるチャネルは、画像の対象が実際に発現する実際の可視色、または虚偽の色の画像を表すことができ、そのことによって、画像の対象の実際の色特性を表さない場合があるが、それでも観視者に視覚的に訴える可視色特性で、ハイパースペクトル画像の重み付けされたパンクロマチック表現が着色される。
本発明の方法とデモザイキングの周知の技法との間の重要な相違は、デモザイキングで一般に使用される技法は、各々の部分画素で本来的に欠損している赤、緑および青の部分画素値を推定するために周囲の画素を補間するということである。今回の発明の方法は、周囲の部分画素の分析により欠損した画素値に対して補間を試みるのではなく、代わりに、表示用に使用されるRGB画素の中心を変化させる。色チャネルの各々にわたって繰り返すことにより、この方法は、欠損した画素値を補間しようとすることなく、ベイヤーフィルタの色フィルタ配列パターン規定を絶えず再規定することによりRGB画素値の補間を遂行することができる。
記述した本説明では、最良の形態を含めて本発明を開示するために、さらには、任意のデバイスまたはシステムを作製かつ使用すること、および任意の組み込んだ方法を遂行することを含めて、本発明を実践することを当業者ならば誰でも可能にするために例を使用する。本発明の特許的な範囲は、特許請求の範囲により定義され、当業者が想到する他の例を含み得る。そのような他の例は、それらの例が、特許請求の範囲の文字通りの文言と異ならない構造要素を有するならば、またはそれらの例が、特許請求の範囲の文字通りの文言と実質的な違いのない等価の構造要素を含むならば、特許請求の範囲の範囲内にあることが意図される。
10 ハイパースペクトル画像のサンプリングの従来技術の方法
12 ハイパースペクトルデータキューブ
14 ハイパースペクトル画像
16、18、20 ハイパースペクトル部分画素
22 ハイパースペクトル画素、窓
24 従来技術の変換方法
26、27 ハイパースペクトル部分画素
28、30 パンクロマチック画素
44、46、48、50、52、54 ハイパースペクトル部分画素
80 アップサンプリング方法
82 窓を規定するステップ
84 第1の窓位置で開始するステップ
86 窓内の部分画素値を重み付けするステップ
88 重み付き平均を出力画像画素に割り当てるステップ
90 窓カウンタをインクリメントするステップ
92 処理が完了したかどうかを決定するステップ
100 第1の画素をアップサンプリングする
114 ハイパースペクトル画像
124 第1の画素に対するアップサンプリング方法
128 ハイパースペクトル画像の第1の画素
200 第2の画素をアップサンプリングする
216 シフトされた窓の中心部分画素
224 第2の画素に対するアップサンプリング方法
226 ハイパースペクトル画像の第2の画素
300 第3の画素をアップサンプリングする
316 シフトされた窓の中心部分画素
324 第3の画素に対するアップサンプリング方法
328 ハイパースペクトル画像の第3の画素
400 第4の画素をアップサンプリングする
416 シフトされた窓の中心部分画素
424 第4の画素に対するアップサンプリング方法
426 ハイパースペクトル部分画素
430 ハイパースペクトル画像の第4の画素
500 ハイパースペクトル画像をアップサンプリングする
518、520 ハイパースペクトル部分画素
524 アップサンプリング方法

Claims (7)

  1. 画素、および各々の画素内部で部分画素を規定する少なくとも1つのフィルタを備える画像センサからのハイパースペクトル画像の分解能を制御する方法であって、
    各々が他の部分画素のスペクトル帯域から独立したスペクトル帯域を表す部分画素の行および列の配列を備える前記画像センサ上に窓を規定するステップと、
    前記ハイパースペクトル画像の少なくとも1つの既定のパラメータに基づいて、前記少なくとも1つの既定のパラメータに対する前記配列に対する重み付き平均に対する値を確定するために、前記窓内部の前記部分画素を重み付けするステップと、
    既定の数の行または列だけ前記窓をシフトするステップであって、前記既定の数が前記画素内の部分画素のそれぞれの行または列の数未満であるステップと、
    前記画像センサ上のすべての可能な配列に対して前記重み付けするステップおよびシフトするステップを反復するステップと、
    前記重み付き平均に基づいて前記ハイパースペクトル画像を処理するステップと、
    を含む、方法。
  2. 前記少なくとも1つの既定のパラメータが、赤、緑および青のスペクトル帯域のうちの少なくとも1つの強度値である、請求項1に記載の方法。
  3. 前記少なくとも1つの既定のパラメータが、赤、緑および青のスペクトル帯域のすべての強度値である、請求項2に記載の方法。
  4. 前記少なくとも1つのフィルタがベイヤーフィルタである、請求項1乃至3のいずれかに記載の方法。
  5. 前記少なくとも1つの既定のパラメータが、隣接するスペクトル帯域の強度値である、請求項1乃至4のいずれかに記載の方法。
  6. 前記強度値が16個のスペクトル帯域から決定される、請求項5に記載の方法。
  7. 前記少なくとも1つのフィルタが、マイクロレンズの配列のレンズレットを備えるハイパースペクトルフィルタである、請求項1乃至6のいずれかに記載の方法。
JP2013172769A 2012-08-29 2013-08-23 ハイパースペクトル画像の分解能を制御する方法 Expired - Fee Related JP6291189B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US13/598,029 2012-08-29
US13/598,029 US8797431B2 (en) 2012-08-29 2012-08-29 Method of controlling the resolution of a hyperspectral image

Publications (3)

Publication Number Publication Date
JP2014049126A JP2014049126A (ja) 2014-03-17
JP2014049126A5 JP2014049126A5 (ja) 2016-09-29
JP6291189B2 true JP6291189B2 (ja) 2018-03-14

Family

ID=49054431

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013172769A Expired - Fee Related JP6291189B2 (ja) 2012-08-29 2013-08-23 ハイパースペクトル画像の分解能を制御する方法

Country Status (6)

Country Link
US (1) US8797431B2 (ja)
EP (1) EP2703792A1 (ja)
JP (1) JP6291189B2 (ja)
CN (1) CN103679758B (ja)
BR (1) BR102013020692A2 (ja)
CA (1) CA2824231A1 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014125804A1 (ja) * 2013-02-13 2014-08-21 パナソニック株式会社 マルチスペクトル撮像装置およびマルチスペクトル撮像方法
US20220404531A1 (en) * 2021-06-17 2022-12-22 Drs Network & Imaging Systems, Llc Method and system for fabrication and use of a spectral basis filter
CN114383668B (zh) * 2022-03-24 2022-05-24 北京航空航天大学 一种基于可变背景的流场测量装置及方法
CN115060367B (zh) * 2022-06-13 2023-04-21 华东师范大学 基于显微高光谱成像平台的全玻片数据立方体采集方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6847737B1 (en) * 1998-03-13 2005-01-25 University Of Houston System Methods for performing DAF data filtering and padding
US6639665B2 (en) * 2001-02-09 2003-10-28 Institute For Technology Development Multispectral imaging system for contaminant detection
US7231084B2 (en) * 2002-09-27 2007-06-12 Motorola, Inc. Color data image acquistion and processing
US7751594B2 (en) * 2003-04-04 2010-07-06 Lumidigm, Inc. White-light spectral biometric sensors
CA2547359C (en) 2003-11-26 2012-11-27 Florida Environmental Research Institute, Inc. Spectral imaging system
US7242478B1 (en) * 2003-12-05 2007-07-10 Surface Optics Corporation Spatially corrected full-cubed hyperspectral imager
US7570810B2 (en) 2005-02-24 2009-08-04 Seiko Epson Corporation Method and apparatus applying digital image filtering to color filter array data
US7979209B2 (en) 2005-04-15 2011-07-12 Mississippi State University Research And Technology Corporation Temporal mapping and analysis
US20070058050A1 (en) * 2005-09-12 2007-03-15 Manuel Innocent Reconstructing a full color image from an image encoded by a bayer pattern
US7826685B2 (en) * 2006-12-01 2010-11-02 Harris Corporation Spatial and spectral calibration of a panchromatic, multispectral image pair
EP2174176A1 (en) * 2007-07-17 2010-04-14 Explay Ltd. Coherent imaging method of laser projection and apparatus thereof
JP5305992B2 (ja) * 2009-03-10 2013-10-02 キヤノン株式会社 画像処理装置及び画像処理方法
US8351031B2 (en) * 2009-06-05 2013-01-08 Spectral Sciences, Inc. Single-shot spectral imager
US8611603B2 (en) * 2012-02-14 2013-12-17 The United States Of America As Represented By The Secretary Of The Army Method and apparatus for object tracking via hyperspectral imagery

Also Published As

Publication number Publication date
CA2824231A1 (en) 2014-02-28
EP2703792A1 (en) 2014-03-05
JP2014049126A (ja) 2014-03-17
US8797431B2 (en) 2014-08-05
BR102013020692A2 (pt) 2016-03-15
CN103679758B (zh) 2018-01-02
US20140063298A1 (en) 2014-03-06
CN103679758A (zh) 2014-03-26

Similar Documents

Publication Publication Date Title
Abraham et al. The Gemini Deep Deep Survey. I. Introduction to the survey, catalogs, and composite spectra
JP2021506168A5 (ja)
CN105096286B (zh) 遥感图像的融合方法及装置
JP6291189B2 (ja) ハイパースペクトル画像の分解能を制御する方法
US20110115812A1 (en) Method for colorization of point cloud data based on radiometric imagery
CA2931741C (en) Analysis of a multispectral image
CN104052979B (zh) 用于图像处理的装置和技术
Viticchié et al. Imaging spectropolarimetry with IBIS. II. On the fine structure of G-band bright features
US9898953B2 (en) Offset method and equipment of RGBW panel subpixel
US20190139189A1 (en) Image remosaicing
US9336570B2 (en) Demosaicking system and method for color array based multi-spectral sensors
US20180366075A1 (en) Sub-pixel rendering method and rendering apparatus
CN104170377A (zh) 图像处理装置、摄像装置及图像处理程序
Cai et al. Feature-driven multilayer visualization for remotely sensed hyperspectral imagery
Zhu et al. HCNNet: A hybrid convolutional neural network for spatiotemporal image fusion
JP2014116706A (ja) 二次元カラーコード表示装置、表示方法、読み取り装置、及び読み取り方法
Hu et al. A spatial constraint and deep learning based hyperspectral image super-resolution method
EP2517172B1 (en) Filter setup learning for binary sensor
CN106470335B (zh) 基于亚像素采样的图像处理方法及图像显示方法
KR101329134B1 (ko) 부화소 렌더링을 이용한 영상 변환 시스템 및 방법
KR101671927B1 (ko) 이미지 처리 장치 및 방법
Humphries et al. A multi-wavelength analysis of small-scale brightenings observed by IRIS
EP3002731A2 (en) Determing color information using a binary sensor
JP6426909B2 (ja) 色情報補完装置およびそのプログラム
CN210578954U (zh) 一种具有多通道窄带滤色片阵列的图像传感器

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160815

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160815

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20170519

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170530

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170807

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20180124

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180209

R150 Certificate of patent or registration of utility model

Ref document number: 6291189

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees