計測領域の背後に背景画像を、前方側にデジタルカメラを配置し、密度分布のない状態で計測領域を撮影した参照画像と、密度分布が発生した状態での計測画像の2つの画像から画素毎の差分と密度勾配を演算することによって得た結果画像から密度勾配を可視化する方法。結果画像は次式によって得られるものである請求項1に記載の密度勾配を可視化する方法。Rresult=Rrun-Rref×∠RrunOffsetここでRはi、j番地の画素の輝度値を表し、添え字refは参照画像、runは計測画像、resultは結果の画像を表している。上式中の∠Rrunは、密度変化の方向を与える項であり、画素点周りの1次微分値を用いる。密度変化の方向を与える項∠Rrunは、次式によって得られるものである請求項2に記載の密度勾配を可視化する方法。∠Rrun=Rrun-Rrun2Rref-Rref22得られた結果画像に対して移動平均処理でスムージングを施す請求項3に記載の密度勾配を可視化する方法。計測領域の背後に配置する背景画像と、計測領域の前方側に配置するデジタルカメラと、密度分布のない状態で計測領域を撮影した参照画像と、密度分布が発生した状態での計測画像の2つの画像から画素毎の差分と密度勾配を演算する手段とからなる密度勾配を可視化する装置。背景画像には連続的な輝度分布の縞模様を用いたものである請求項5に記載の密度勾配を可視化する装置。前記演算手段はデジタルカメラのボディ内に設置したものである請求項5又は6に記載の密度勾配を可視化する装置。計測領域は風洞内であって、該領域は背景画像側の後方と前方カメラ側は透明な窓部とされている請求項5乃至7のいずれかに記載の密度勾配を可視化する装置。本発明の課題は、凹面鏡や凸レンズを必要せず、平易な画像演算によって気流の密度勾配の可視化を可能にする手法を提供することにある。本発明の密度勾配を可視化する方法は、計測領域の背後に背景画像を、前方側にデジタルカメラを配置し、密度分布のない状態で計測領域を撮影した参照画像と、密度分布が発生した状態で撮影した計測画像の2つの画像から画素毎の差分と密度勾配を演算することによって得た結果画像を用いるものとした。そして、結果画像を得る前記の演算は次式によって得られるものとした。Rresult=Rrun-Rref×∠RrunOffsetここでRはi、j番地の画素の輝度値を表し、添え字refは参照画像、runは計測画像、resultは結果の画像を表している。上式中の∠Rrunは、密度変化の方向を与える項であり、画素点周りの1次微分値を用いる。更に、例えば今、i方向の密度勾配を可視化するため、縞がi方向にあるとすると、最も平易な演算法として、前記の密度変化の方向を与える項∠Rrunは、次式によって得られるものとした。∠Rrun=Rrun-Rrun2Rref-Rref22更に、得られた結果画像に対して移動平均処理でスムージングを施すようにした。本発明の密度勾配を可視化する装置は、計測領域の背後に配置する背景画像と、計測領域の前方側に配置するデジタルカメラと、密度分布のない状態で計測領域を撮影した参照画像と、密度分布が発生した状態での計測画像の2つの画像から画素毎の差分と密度勾配を演算する手段とからなるものとした。そして、前記背景画像には連続的な輝度分布の縞模様を用いるものとした。また、一つの形態として前記演算手段はデジタルカメラのボディ内に設置する形態を採用した。本発明の密度勾配を可視化する装置は、その計測領域は風洞内であって、該領域は背景画像側の後方と前方カメラ側は透明な窓部とされている構成を採用した。本発明では、従来のシュリーレン法で不可欠な平行光線を作るための凹面鏡等による大がかりな光学系を必要としないシンプルな装置構成で実施できる。その結果、凹面鏡等の光学系の設置に依存する撮影方向の制約がなく、あらゆる方向からの可視化が可能である。このことによって従来、被写体の影になり可視化することができなかった部分の可視化が容易になる。更に、室内外等適用場所を問わない可搬性を持った計測法である。また従来のBOS法で用いる直接相互相関法やフーリエ変換によるクロススペクトル計算といった複雑な演算や、それに伴って統計的に発生するエラーの除去を必要としない。演算方法は式に基づく単純な画像処理であるため、撮像装置内に組み込み、単体で可視化が可能な撮像装置を製品化することもできる。また光学的装置の設置調整も従来技術に比べて簡単に済むという効果を奏する。以下、本発明の実施の形態について、詳細に説明する。可視化したい密度勾配の方向に分布を持つ背景画像を用意する。例えば、モノクロカメラを用いて鉛直方向の密度勾配を可視化しようとする場合、図2に示すような縞状に輝度値分布を持つ画像を印刷したスクリーンを準備する。密度勾配を持つ対象の流れ場1、図1に示す例ではバナーの炎による高温雰囲気領域を挟んで、背景画像スクリーン2に対して反対側に撮影用のデジタルカメラまたはデジタルビデオカメラ3を設置する。本発明では画像情報取得のための装置は以上の構成で足りる。カメラ3と対象と背景2の位置関係は、従来のBOS法と同じ位置関係で最適値を決める。この際、想定される縞の移動量より、背景2の縞の間隔を大きくしておく。流れ場を生成する前に、予め背景画像2を参照画像として取得しておく。次いで流れ場1を生成し、参照画像2の取得と同様に計測画像を取得する。この計測画像では、発生した密度勾配によって光は屈折し、歪んだ縞の画像となってカメラ3に撮像される。次いで、得られた計測画像と参照画像に対して以下の演算を行う。この演算は特殊な装置を必要とせず汎用のパーソナルコンピュータで行うことが可能であり、デジタルカメラのボディ内に内蔵させる形態も可能である。次式1はi、j番地の画素の計測画像と参照画像との差分とその密度変化の方向を示すものである。Rresult=Rrun-Rref×∠RrunOffset............ここでRはi、j番地の画素の輝度値を表し、添え字refは参照画像、runは計測画像、resultは結果の画像を表している。式中の∠Rrunは、密度変化の方向を与える項であり、点周りの1次微分値を用いる。例えば今、i方向の密度勾配を可視化するため、縞がi方向にあるとすると、最も平易な演算法として、計測画像と参照画像の点周りの空間差分の平均は次式2∠Rrun=Rrun-Rrun2Rref-Rref22............を用いることができる。式の作用により、式の第1項は、密度の変化があった場所について、その方向と変化量に従った輝度値を与える。末尾のOffset項は、変化のない点の輝度値を与えるために用いる付与項である。モノクロ8ビットの場合、付与項の値をグレースケールの中間値である128とすると、従来のシュリーレン法のように、変化のない場合は灰、変化の方向に従って白または黒に変化する画像が得ることができる。カラー画像を用いる場合は、RGB値それぞれについて同様の演算を行えばよい。温められた気流1の可視化の実施例を説明する。図1のようにデジタルカメラ3、背景画像2、ガスバーナー4を設置する。前記背景画像2はアクリルフレームに取付けられている。鉛直方向の密度勾配を得るために背景画像は図2に示す連続的な輝度分布の縞模様を用いた。前述したように従来のシュリーレン法ではこの現象を可視化する場合、図8のような複雑な光学装置を必要としていた。しかし、本発明ではこのような複雑な光学系は一切必要がなく、測定領域温められた気流場1の背後に背景画像2を設置するだけで前方所望角からデジタルカメラ3によって該測定領域を撮影すればよい。まず、ガスバーナー4を点火させない時点で図3Aに示す参照画像を取得し、ガスバーナーを点火後、図3Bに示す計測画像を取得する。取得したこの2つの画像情報から式にしたがって画像処理をし、各点についてRresultを演算し、演算結果の合成画像として図4Aが得られる。この演算は図示していないパーソナルコンピュータを用いて実行する。この図4Aにおいては画像上向きに正の密度勾配がある場所は黒方向、負の密度勾配がある場合は白方向で階調化し可視化されている。高い解像度の撮像機材を用いると、図4Bに示すよう補間のために移動平均処理を行った場合、滑らかな自然な画像が十分な解像度でシュリーレン写真として得られることが分かる。次に本発明の手法を超音速風洞試験における密度勾配の可視化に適用した実施例を説明する。図5のように供試体5を風洞6内の測定領域1内の設置し、その背後に背景画像2、前方に風洞にデジタルカメラ3を設置する。風洞6に流す気流はM=30の超音速流れであり、流れの中では円錐の供試体5周りには、衝撃波による圧縮領域と、膨張波による膨張領域が生じる。背景画像は図6Aに示す縞模様を用いた。気流を発生させる前に参照画像を取得し、気流を発生させた後、図6Bに示す計測画像を取得する。先の実施例と同様に各点についてRresultを演算し、演算結果の合成画像として式にしたがって画像処理をし、結果の画像として図6Cを得た。画像では、圧縮領域が黒色に、膨張領域が白く可視化されている。図7に本発明による可視化画像を従来のシュリーレン画像と比較して示す。上側画像が本発明による画像で、下側画像が従来のシュリーレン画像である。図5に見られるように本発明は、設置が容易かつ平演算も平易である上、従来のシュリーレン法と同等の画像を得ることができる。図7では上側の本発明による画像の撮像感度が下側の従来のシュリーレン画像よりよく、その結果、衝撃波による圧縮領域と、膨張波による膨張領域がより鮮明に得られている。