三次元空間における流体の流れを表す流線を、任意の視点位置から見た二次元画像として描画する描画装置において、所定の開始位置から流線を生成する流線生成部と、前記視点位置を前記三次元空間内に設定するビュー設定部と、複数の多角形面を連結した帯状面を、前記設定された視点位置に基づき、前記流線に沿うように形成する帯状面形成部と、前記三次元空間の所定位置に配置された光源に対する前記帯状面の反射光を表すテクスチャを、前記複数の多角形面の各々にマッピングするテクスチャマッピング部とを有することを特徴とする描画装置。前記テクスチャマッピング部は、前記三次元空間の所定位置に配置された光源への法線を表す法線情報を用いることにより、前記三次元空間の所定位置に配置された光源に対する前記帯状面の反射光を表すテクスチャを、前記複数の多角形面の各々にマッピングすることを特徴とする請求項1記載の描画装置。前記帯状面形成部は、前記設定された視点位置から見た前記複数の多角形面の各々の面積が最大となるよう、又は、前記複数の多角形面の各々の上に存在する特定の点から延ばした法線に対して、前記設定された視点位置から下ろした垂線の長さが最小となるよう、又は、前記複数の多角形面の各々の各頂点と前記設定された視点位置間の距離の総和が最小となるように、前記複数の多角形面の各々の傾きを決定することを特徴とする請求項1記載の描画装置。前記帯状面形成部は、前記複数の多角形面において互いに隣接する第1の多角形面と第2の多角形面とが共有する辺が、前記流線と交叉するように前記帯状面を形成することを特徴とする請求項1記載の描画装置。前記テクスチャマッピング部は、テクスチャをマッピングする前記複数の多角形面の各々に対応する前記流線の位置における前記流体の物理値に基づいて、テクスチャをマッピングすることを特徴とする請求項1記載の描画装置。前記帯状面形成部は、前記複数の多角形面の各々に対応する前記流線の位置における前記流体の物理値に基づいて、前記複数の多角形面の各々の形状を決定することを特徴とする請求項1記載の描画装置。三次元空間における流体の流れを表す流線を、任意の視点位置から見た二次元画像として描画する描画方法において、コンピュータが、所定の開始位置から流線を生成し、前記視点位置を前記三次元空間内に設定し、複数の多角形面を連結した帯状面を、前記設定された視点位置に基づき、前記流線に沿うように形成し、前記三次元空間の所定位置に配置された光源に対する前記帯状面の反射光を表すテクスチャを、前記複数の多角形面の各々にマッピングすることを特徴とする描画方法。三次元空間における流体の流れを表す流線を、任意の視点位置から見た二次元画像として描画する描画プログラムにおいて、コンピュータに、所定の開始位置から流線を生成させ、前記視点位置を前記三次元空間内に設定させ、複数の多角形面を連結した帯状面を、前記設定された視点位置に基づき、前記流線に沿うように形成させ、前記三次元空間の所定位置に配置された光源に対する前記帯状面の反射光を表すテクスチャを、前記複数の多角形面の各々にマッピングさせることを特徴とする描画プログラム。1つの側面では、本発明は、数値シミュレーション結果を表示する際の処理負荷の軽減を図ることを目的とする。本装置の一観点によれば、三次元空間における流体の流れを表す流線を、任意の視点位置から見た二次元画像として描画する描画装置において、所定の開始位置から流線を生成する流線生成部と、前記視点位置を前記三次元空間内に設定するビュー設定部と、複数の多角形面を連結した帯状面を、前記設定された視点位置に基づき、前記流線に沿うように形成する帯状面形成部と、前記三次元空間の所定位置に配置された光源に対する前記帯状面の反射光を表すテクスチャを、前記複数の多角形面の各々にマッピングするテクスチャマッピング部を有することを特徴とする描画装置が提供される。数値シミュレーション結果を表示する際の処理負荷の軽減を図ることができる。図3は、三次元空間上に満たされた流体不図示の流線L300を示している。一般に、ある瞬間における各点の流体分子の速度ベクトルに沿ってできる曲線を流線という。また、同一の開始点から次々と流れる複数の流体分子をある瞬間に結んでできる曲線を流脈線という。そして、一つの流体分子の移動経路を継続追跡して描かれる曲線を流跡線という。定常な流れでは、これら三者は同一となるが、非定常な流れでは異なる線となる。本明細書においては、上記の流線、流脈線及び流跡線を、便宜上総称して流線の語を用いて説明する。本明細書に記載された複数の実施例は、これら三者のいずれにも適用できることは言うまでもない。図3において、流線L300は、各時間ステップのベクトル場を入力データとして、本出願の出願日時点で当業者に知られているRungeKutta法、あるいは、DormandPrince法により、流線L300を計算することができる。流線L300には、曲率の大きな特徴部分に点節点P301乃至P306を配置してもよい。図4A乃至図4Dは、流線に沿って四辺形を連結した帯状面を形成する処理を示している。なお、本実施形態では、四辺形矩形を用いているが、四辺形以外の多角形面を用いてもよい。図4Aにおいて、三次元空間上に視点I及び視線ベクトルVを設定する。この視点I及び視線ベクトルVは、ユーザによって指定されてもよい。ユーザからの指定がない場合は、予め定められた視点Iの位置及び視点ベクトルVの値を用いてもよい。また、様々な角度から、流線を観察するために、これらの値は、ユーザにより適宜変更されてもよい。指定された視点I及び視線ベクトルVの情報は、記憶部に記憶され、本実施例のシステムによって利用される。開始位置P401から終端位置Pendまで、流線L400が示されている。まず、流線L400の開始位置P401と、隣の節点P402を含むように、四辺形M410を生成する。四辺形M410の辺H401上に開始位置P401が存在し、辺H402上に節点P402が存在するようにしてもよい。節点を通る辺H401、H402は、予め定められた一定の長さとしてもよい。あるいは、節点における物理値例えば、流体の速度、温度に応じて、辺の長さを決定してもよい。四辺形M410は、4つの頂点P401a、P401b、P402a、P402bを持つ。四辺形M410の面の傾きは、視点Iから四辺形M410が見える面積が最大となるように決定される。同様に、節点P402と節点P403を用いて、四辺形M420を生成する。四辺形M420は、節点P402を通る辺H402と、節点P403を通る辺H403を持ち、4つの頂点P402a、P402b、P403a、P403bを有する。そして、四辺形M410と四辺形M420は、辺H402を共有する。なお、各々の四辺形は、必ずしも長方形又は正方形である必要はなく、平行四辺形又は台形となってもよい。あるいは四辺形以外の三角形やその他の多角形でもよい。以上の操作を繰り返し、終端位置Pendまで、四辺形を生成する。この操作によって、複数の四辺形により接続された帯状面Q400が流線L400に沿って形成される。なお、各四辺形の面の向きをどのように決定するかのアルゴリズムの詳細は後述する。図4Bは、着目する流線と平行になるように新たな流線をコピーする操作を示している。この操作は、図4Cで後述するように、四辺形の数を増やし、帯状面Q400を滑らかにするための前処理として位置づけられる。四辺形を適切な数だけ生成する処理は、ここで述べる方法以外にもあり得る点に留意すべきである。例えば、図4Aにおいて、予め節点の数を増やしておき、四辺形を生成してもよい。流線L400を挟むように、流線L400に近似した近似流線L400a及びL400bを描く。L400aは、四辺形の頂点P401a、P402a、P403a、、Pendaを通るように、流線L400を平行移動して生成した曲線としてもよい。同様に、四辺形の頂点P401b、P402b、P403b、、Pendbを通る曲線L400bを生成する。なお、本実施形態では、着目する流線を挟むように2つの新たな曲線を生成したが、新たな流線を1つだけ生成してもよい。図4Cは、辺を追加することにより、四辺形の数を増加させる操作を示している。例えば、辺H401と辺H402の間に、辺H401aが追加されている。同様に、辺H402と辺H403の間に、辺H402aが追加されている。辺の追加は、例えば以下のように行えばよい。1各流線において隣り合う2つの節点を結ぶ線分の勾配勾配ベクトルを求める。22つの節点を結ぶ線分勾配の差が、予め定められた閾値thより大きい場合、例えば、その2つの節点の中間点に新たな節点を生成する。3各流線において生成した新たな節点を結ぶ新たな辺を生成する。4新たな辺と、隣接する既存の辺とで構成される四辺形を生成する。以上の操作により、例えば、辺H401と、新たな辺H401aとを持つ新たな四辺形を形成することができる。同様に辺H401aと辺H402とで、新たな四辺形を形成する。この操作を繰り返すことにより、四辺形の数を増加させることができる。なお、上述のように、四辺形の増加の方法は、この実施形態に限定されない。なお、上記2の処理において、2つの節点の中間点以外に新たな節点を生成しても良い。図4Dは、三角形のポリゴンを形成する操作を示している。ポリゴンとは、立体形状を三次元的に表現する場合に用いられる当該立体形状の表面を分割する多角形である。計算のしやすさを考慮した場合、ポリゴンとして三角形が用いられる。本実施形態では、三角形のポリゴンを形成している。例えば、四辺形に一本の対角線m401を引くことにより、三角形のポリゴンPg401とPg402が形成される。対角線m401乃至m404を引くことにより、ポリゴンPg401乃至Pg408が形成される。三角形のポリゴンを形成することにより、汎用的な三次元グラフィクスソフト及びレンダリングハードウエアを用いて処理を行うことが容易となる。なお、この段階で、四辺形の辺のうち、流線と交叉している辺の長さを、流線のその位置での物理値に基づいて伸縮させてもよい。辺の長さを伸縮させることにより、チューブの太さを変化させることができる。辺の長さの伸縮は、適切な他の処理段階で行ってもよいことは言うまでもない。なお、最終的なチューブ形状に見える二次元画像をレンダリングするために、ポリゴンを生成することは必須ではない。従来技術では、流線1本あたりのポリゴン数がより多い。このため、処理を行う計算機の画像処理性能によっては、流線の可視化により多大な時間が必要となる。このため、入力の段階において単純にデータ量を減らす方法が知られている。シミュレーションのデータは既に時間空間で離散化されているが、再度サンプリングして、より疎なデータとする。この方法では、観察者が見たい現象を正確に再現できないというリスクがある。すなわち、この方法は、処理前に情報を削減してしまうため、流線の構造が、データの削減前のものと同じにならない可能性が高い。これに比して、本実施例では、データ量を削減せずに、流線の1本あたりのポリゴン数を大幅に削減できる。図5Aは、法線マップの模式図を示している。四辺形500は面であるため、このままの状態でレンダリングを行うと、単純な面として描画されてしまう。この四辺形をチューブ形状に見えるようにするために、四辺形の各領域c1乃至c4に、それぞれ法線ベクトルu1乃至u4をマッピングする。この法線ベクトルu1乃至u4は、チューブ形状の法線をシミュレートするよう設定される。図5Bは、法線マップ情報をメモリ上に記憶する際のテーブルの例を示している。それぞれの位置に対応づけて法線ベクトルの値が格納される。なお、この実施形態では、4つの位置に法線ベクトルを対応付けたが、更に多くのベクトルを持つ法線マップを用いて細かい対応付けを行ってもよい。また、法線マップの代わりに、チューブ形状の高さを示す高さマップを保持してもよい。なお、法線マップは、長方形500に対応して保持しておけば、平行四辺形、台形等の四辺形矩形に対する法線マップを計算することができる。たとえば、四辺形500の四隅の頂点を台形の四隅の頂点と対応付けて、モーフィングと同様の手法で法線マップを台形にマッピングすればよい。同様に、三角形や四辺形以外の多角形に対しても、長方形500の法線マップを利用することができる。なお、三角形やその他の多角形用のための法線マップ情報を、ライブラリとして別途記憶しておいてもよいことは言うまでもない。なお、上述の例では、ポリゴンとして三角形を用いているが、ポリゴンは、三次元形状のモデリングの手法として用いるものであり、法線マップのデータ構造とは別個に扱うことができる。図6Aは、反射光のテクスチャのマッピングを示している。上述の法線マップを適用して光源650からの光の反射を考慮し、視点Iから視線ベクトルV方向に見た四辺形600をレンダリングした結果が示されている。例として、適用したテクスチャは単色である。図6Bは、それぞれの流線の位置における物理値例えば、速度を色にマッピングしたテクスチャを用いてレンダリングした結果を示している。図6Bは、白黒の階調表示で描かれているが、実際には、色を変化させたテクスチャをマッピングしている。このように、流線をチューブ形状によって可視化すると共に、流線の個々の位置における物理値例えば、温度を色によって表現することができる。更に、別の物理値例えば、速度に応じてチューブ形状の太さを変えることにより、チューブ形状の色と太さにより、二つの物理値例えば、温度と速度を同時に表現することができる図9参照。図7は、物理値とテクスチャの色とのマッピングを示している。例えば物理値として流体の温度を四辺形に対して色のテクスチャとしてマッピングするための関係を示している。温度が低い場合には、青にマッピングされ、温度が高い場合には赤にマッピングされる。温度が中程度であれば黄色にマッピングされる。図8は、物理値とチューブの太さとのマッピングを示している。例えば、流体の速度流速が小さければチューブの太さを太くするすなわち四辺形の辺のうち、流線と交叉している辺の長さを長くする。逆に、流体の速度が大きければチューブの太さを捕捉するすなわち四辺形の辺のうち、流線と交叉している辺の長さを短くする。図9は、図7及び図8の関係を用いて、帯状面を形成し、テクスチャマッピングを施し、レンダリングを行った結果の例を示している。図9において、951乃至954は、それぞれ、四辺形を表している。描かれたチューブは、901及び903の領域で細くなっており、902の領域で太くなっている。このことから、902の領域では、他の部分に比べて流速が小さいことが分かる図8参照。また、図9においては、白黒となっているが、901及び903の領域では、赤い色を示しており、902の領域では黄色い色を示している。したがって、902の領域では、他の部分に比べて温度が低いことが分かる。図10A及び図10Bは、視点を視点Iaから視点Ibに移動させた場合に、四辺形M100の面の向きを決定する処理を示している。まず、図10Aにおいて、視点は、視点Iaに位置している。四辺形M100において、2つの対向する辺は、流線L100の節点P101とP102において接している。この場合、四辺形M100は、方向D120又は方向D130に、面の傾きを変化させることができる。四辺形M100の面の傾きは、視点が視点Iaの位置にある場合には、例えば以下のように決定することができる。[視点が視点Iaの位置にある場合の四辺形M100の向きの決定]1四辺形M100上に単位面積T150を設定する。そして、視点Iaから見てT150の面積が最大に見えるように、四辺形M100の傾きを決定する。又は、2四辺形M100上の任意の点S160から四辺形M100に対する法線ベクトルVsを延ばす。そして、視点Iaから法線ベクトルVsに垂線を下ろし、その垂線の長さLaが最小となるように、四辺形M100の傾きを決定する。3四辺形M100の四隅の点から視点までの距離の総和が最小となるように、四辺形M100の傾きを決定する。以上が、四辺形M100の面の傾きを決定する処理の例である。上記の決定処理以外にも、当業者であれば、M100の面の傾きを決定する手法は見いだせるであろう。図10Bは、視点が視点Iaから視点Ibに移動した場合の、四辺形M100の面の傾きを決定する処理を示している。上述のように、視点Ibから見て、四辺形M100上の単位面積T150が最大に見えるように、四辺形M100の面の傾きを決定すればよい。又は、四辺形M100上の任意の点S160から四辺形M100に対する法線ベクトルVsを延ばす。そして、視点Ibから法線ベクトルVsに垂線を下ろし、その垂線の長さLbが最小となるように、四辺形M100の傾きを決定すればよい。図10Cは、視点が視点Iaから視点Ibに移動した場合の、四辺形M100の面の傾きを決定する別の処理を示している。視点Ibと、四辺形M100の各頂点K181、K182、K183、及びK184とを結んだ直線を、それぞれN181、N182、N183、及びN184とする。そして、直線N181、N182、N183、及びN184の長さの総和が最も小さくなるように、四辺形M100の傾きを決定すればよい。以上の操作を他の四辺形の面の傾きの決定にも行えば、帯状面の全体の傾きが変更される。その後、テクスチャマッピングを施し、レンダリングを行えば、新たな視点Ibから見たチューブ状の流線の二次元画像が得られる。このように、帯状面を生成し、テクスチャマッピングを利用して、帯状面をチューブ状に見えるようレンダリングする。三次元空間上にチューブを実際に生成するよりも、帯状面を利用する方が、ポリゴン数を従来よりも減少させることができる。したがって、流線を可視化する処理を従来よりも短縮することが可能である。また、視点を変えた場合にも、上述のように四辺形の面の向きを変更し、テクスチャマッピングを施し、レンダリングを行えばよい。上述のように、少ないポリゴン数のため、視線を変更しても従来よりも高速に二次元画像が得られる。なお、四辺形の面の傾きをより高速に決定するには、更に視線ベクトル及び視野を考慮してもよい図示せず。すなわち、視線ベクトルと視野によって、最終的にレンダリングされるべき三次元モデルの部分が特定される。この特定された三次元モデル部分だけに着目して、四辺形の傾きの変更、テクスチャマッピング、レンダリング等の処理を限定して行えばよい。図11Aは、一実施形態の処理フローを示している。以下に示す処理フローは、矛盾のない限り、順番を入れ替えてもよい。そして、また、プログラムに係る請求項及び、方法に係る請求項における、各ステップについても、矛盾のない限り順番を入れ替えてもよいことは言うまでもない。そして、請求項内のステップの順番を入れ替えた場合でも、その請求項の技術的範囲に包含される。また、既に記憶部にデータが記憶されており、改めてユーザからの指示を受け付ける必要がない場合がある点に留意すべきである。加えて、ユーザからの指示がない場合には、デフォルト値を用いて処理を行ってもよい。ステップ1102で、まず、三次元空間上の視点位置、視線方向、流線開始位置の指示をユーザから受け取る。なお、これらのデータは、予め記憶装置に記憶しておいてもよい。ステップ1104で、シミュレーションなどによって得られたベクトル場のデータを読み込む。ステップ1106で、記憶部に記憶された流線の開始位置を読み込む。ステップ1108で、流線を生成する。ステップ1110で、記憶部に記憶された視点位置を読み込む。ステップ1112で、流線に沿って、帯状面を生成する。ここでは、まず、複数の四辺形を生成する。次に、流線をコピーして新たな流線を生成する。そして、四辺形の両端の辺において、流線に沿った勾配ベクトルの変化量が、予め定められた閾値thよりも大きい場合には、例えば中間点において、四辺形を分割する。なお、既に述べたように、流線をコピーせずに、流線における節点を適切に増加させ、隣接する節点間によって四辺形を形成させてもよい。なお、例えば、節点における物理値例えば、速度、温度等に基づいて、節点を含む四辺形の辺の長さを、予め定めた関数を用いて決定してもよい。ステップ1114で、四辺形に対角線を一本引いて、2つの三角形を生成する。この処理を、全ての四辺形に行う。この処理によって、帯状面は、三角形のポリゴンで表現されることとなる。なお、ポリゴンは、必ずしも三角形である必要はないため、四辺形のままでポリゴンとして扱ってもよい。なお、三角形のポリゴンを用いる場合には、三次元モデルの処理が容易化される。ステップ1116で、法線マップのデータを読み込む。法線マップは、長方形に対する法線マップのデータを記憶させておいてもよい。実際の四辺形がひし形、台形である場合には、長方形の法線マップをひし形、台形に適用することができる。例えば、長方形の各頂点を、台形の各頂点と対応させ、変形量を比例配分して台形の法線マップを計算すればよい。ステップ1118で、法線マップを利用して、反射を表現するテクスチャを生成する。また、流線の各位置における物理値例えば、速度、温度等に対応付けて色を加えたテクスチャを生成してもよい。ステップ1120において、未処理の流線があるか否かがチェックされる。未処理の流線がある場合には、ステップ1106に戻る。なお、流線の数は、流線の開始位置の個数と一致する。そして、未処理の流線がない場合には、全ての流線に対する帯状面の作成及びテクスチャマッピングが終了したこととなる。この場合には、ステップ1130に進む。ステップ1130で、視点位置から、視線方向に見たチューブ状の流線が、レンダリングされる。レンダリングされたデータは、表示装置に表示される。そして、処理はAを経由してステップ1142に進む。ステップ1142において、視点を変更するユーザからの指示があるか否かがチェックされる。ユーザは、視点位置を変更して、流線の状況を観察したい場合があるからである。このチェックで、視点位置を変更しない場合には、ステップ1162に進む。ステップ1132で、視線方向を変更するユーザからの指示があるか否かがチェックされる。このチェックで視線方向の変更が指示されていることが分かった場合には、ステップ1164に進む。ステップ1164では、ユーザから視線方向の指示を受け取る。この場合には、視点位置は、変更されていないため、視界の変更の指示がなされたことになる。そして、この場合には、ステップ1130において指示された視点及び視線方向に対する二次元画像をレンダリングすればよい。ステップ1142のチェックで、YESである場合には、ステップ1144に進む。ステップ1144では、視点位置、視線方向の変更指示をユーザから受け取り、記憶部に記憶する。この場合、視線方向については、変更の指示がなく、視点位置のみが変更される場合もある。ステップ1146で、受け付けた視点位置を読み込む。ステップ1148で、視点位置に基づいて、帯状面の各々の四辺形の傾きを修正する。ステップ1150で、ポリゴンの生成が行われる。この処理は、ステップ1114と同様である。ステップ1152で、法線マップデータが読み込まれる。この処理は、ステップ1116と同様である。ステップ1154で、法線マップを利用して、反射を表現するテクスチャを生成する。また、流線の各位置における物理値例えば、速度、温度等に対応付けて色を加えたテクスチャを生成してもよい。この処理は、ステップ1118と同様である。ステップ1156で、未処理の流線があるか否かがチェックされる。未処理の流線がある場合には、ステップ1146に戻る。そして、未処理の流線がない場合には、全ての流線に対する帯状面の各四辺形の傾きの修正及びテクスチャマッピングが終了したこととなる。この場合には、ステップ1130に進み、レンダリングが行われる。ステップ1162において、視線方向を変更しない場合NOには、ユーザは、視点位置及び視線方向の変更を望んでいないため、更なる処理は不要であるため、一連の処理が終了する。図11Bは、一実施形態の処理装置のブロック図を示す。まず、データ入力部1172から流体の流れのシミュレーションデータ等のベクトルデータを読み込み、記憶部1192に記憶する。ビュー設定部1174では、ユーザからの指示を受け付ける。例えば、視点位置、視線方向、焦点距離、並びに、三次元空間上での光源の位置等を受付け、記憶部に記憶する。開始位置設定部1176では、流体の流れを示す流線の開始位置に関するユーザの指定を受け付ける。このデータは、記憶部に記憶される。色、チューブ形状設定部1178では、チューブ形状の色をどのように設定するかに関するユーザからの指定を受け付ける。例えば、流線上の各位置での流体の温度をチューブの色にマッピングする場合には、例えば図7に示したようなマッピング関係を設定する。また、チューブの太さを流体の物理値と関連付けることができる。この場合には、図8に示したように、例えば、流体の流速をチューブの太さに関連付けることができる。なお、チューブの太さは、帯状面を形成する各四辺形の辺の幅を伸縮させることによりコントロールできる。流線生成部1180は、開始位置設定部1176から受け付けられて記憶部に記憶された流線の開始位置と、シミュレーションデータとを用いて、流線を生成する。複数の開始位置が設定されていれば、複数の流線が生成される。帯状面形成部1182は、帯状面を形成する。面形状決定手段1184において、流線に沿って複数の四辺形が決定される。面傾き決定部1183において、この複数の四辺形の傾きが決定される。この処理を経て、帯状面が形成される。また、帯状面を三角形のポリゴンに分割する。なお、面傾き決定部1183は、視点の位置等が変更された場合にも、面の傾きを変更するために利用される。テクスチャマッピング部1185は、複数の四辺形の各々にテクスチャをマッピングする。色適用手段1186は上述のように、物理値と色とのマッピングから、テクスチャに色を適用する。法線マップ適用手段1187は、法線マップ1193を記憶部1192から読み出す。そして、光源の位置と法線マップの法線ベクトルを四辺形の傾きに適用し、反射をシミュレートしたテクスチャを四辺形にマッピングする。レンダリング部1189は、二次元平面に流線を可視化したチューブを描画する。表示部1190は、レンダリングされた画像を表示する。図12は、一実施形態のハードウエア構成1200を示している。ハードウエアには、CPU1202、メモリ1204、動的メモリ1206、媒体1209を読み書きするドライブ1208、通信インターフェース1210、表示装置1213等に画像信号を送るビデオインターフェース1212、I/Oインターフェース1214、I/Oとしてのキーボード1216、マウス1218、プリンタ1220が含まれる。CPU1202は、メモリ1204や動的メモリ1206に格納されている上記説明した種々の機能を実現するプログラムを実行する。また、プログラムに係る本発明は、機械読み取り可能な記録媒体に格納されてもよい。機械読み取り可能な記録媒体には、磁気記録媒体、光ディスク、光磁気記録媒体、半導体メモリなどがある。磁気記録媒体には、HDD、フレキシブルディスクFD、磁気テープMTなどがある。光ディスクには、DVDDigitalVersatileDisc、DVD-RAM、CD-ROMCompactDisc-ReadOnlyMemory、CD-RRecordable/RWReWritableなどがある。また、光磁気記録媒体には、MOMagneto-Opticaldiskなどがある。図13は、四辺形の傾きを変更した場合の四辺形の修正を示す図である。図13aは、流線L1300に対し、四辺形M1310及びM1320が、節点P1320を通る辺H1350を共有している。視点図示せずの位置を変更した場合、これらの四辺形の傾きを変更する処理が行われ、図13bのように、二つの四辺形がM1311及びM1321のように傾きを変えたと仮定する。この場合、二つの四辺形M1311及びM1321は、辺を共有しなくなってしまう場合がある。そのため、節点P1302を通り四辺形M1311及びM1321のそれぞれの面が交わる直線α1380を求める。図13cに示すように、先に求めた直線α1380を共有する辺H1352として、四辺形M1312及びM2322に示すように修正すればよい。この処理を他の四辺形と共有する辺に対しても行えばよい。図14は、四辺形の傾きを変更した場合の四辺形の別の修正を示す図である。図14aは、四辺形M1410、M1420、及びM1430が、流線L1400に沿って帯状面を形成している。この3つの四辺形は、それぞれ隣り合う四辺形と辺を共有している。図14bは、視点図示せずの位置が変更されたために、1つ置きに四辺形の傾きを変更させた状態を示している。すなわち、四辺形M1431及びM1421が四辺形の傾きを変えている。そして、この状態で、両四辺形の頂点K1491とK1494とを結ぶ。同様に両四辺形の頂点K1492とK1493とを結ぶ。図14cは、その結果、新たな四辺形M1411が形成された様子を示している。この処理により、簡便に四辺形の傾きを変えることができる。なお、新たに作成された四辺形M1411は、厳密に平面にならず、捻れた面となる場合があるが、平面に近似して例えば、対角線で2分して、2つの三角形ポリゴンに置き換えて後のレンダリング処理等を行ってもよい。