送電線下横断面図生成方法および装置
【課題】平面地図上で指示された方向の断面線に沿った地盤上の樹高線を容易にしかも高精度に生成する。
【解決手段】送電線が張られている径間方向とは垂直な方向に指示された断面線を長さ方向に、指定された内挿間隔で区切った点を注目点とし、各注目点を中心とする検索半径の円からなる検索範囲内を設定し、（ａ）ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点についての３次元座標を有する三角メッシュデータから、該検索範囲内にある複数の頂点を決定し、（ｂ）該複数の頂点が有する高さ値をｚｑ１，ｚｑ２，…，ｚｑｎとし、該注目点から該検索範囲内の各頂点までのそれぞれの距離に反比例した重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、該注目点における高さ方向の値ｚｐをｚｐ＝（ｗ１・ｚｑ１＋ｗ２・ｚｑ２＋…＋ｗｎ・ｚｑｎ）で示される内挿計算により決定し、（ｃ）決定された各注目点の高さを用いて、樹高線を生成する。
【解決手段】送電線が張られている径間方向とは垂直な方向に指示された断面線を長さ方向に、指定された内挿間隔で区切った点を注目点とし、各注目点を中心とする検索半径の円からなる検索範囲内を設定し、（ａ）ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点についての３次元座標を有する三角メッシュデータから、該検索範囲内にある複数の頂点を決定し、（ｂ）該複数の頂点が有する高さ値をｚｑ１，ｚｑ２，…，ｚｑｎとし、該注目点から該検索範囲内の各頂点までのそれぞれの距離に反比例した重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、該注目点における高さ方向の値ｚｐをｚｐ＝（ｗ１・ｚｑ１＋ｗ２・ｚｑ２＋…＋ｗｎ・ｚｑｎ）で示される内挿計算により決定し、（ｃ）決定された各注目点の高さを用いて、樹高線を生成する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば航空写真から取得した三次元地図データを処理する技術に係り、特に三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成方法およびシステムに関する。
【背景技術】
【０００２】
図化機で取得された３次元地図データから、ＣＡＤ技術を利用して目的とする断面図を作成する方法が特許文献１に開示されている。
【０００３】
特許文献１では、（ａ）内挿間隔と検索半径及び重みを含む内挿パラメータを入力し、（ｂ）入力された三次元地図データを用いて平面地図を表示し、（ｃ）該平面地図の表示画面上で指示された断面位置の両側に各々一つ以上の平行な第２の断面位置を生成し、（ｄ）該三次元地図データから該第１及び第２の断面位置の各断面と交わる各点の地形データを抽出し、（ｅ）抽出された地形データのうち、該第１の断面位置を長さ方向に該内挿間隔で区切った注目点を中心とする該検索半径の円からなる検索範囲内の各点における地形データの高さ方向の値をＺｑ１，Ｚｑ２，…，Ｚｑｎとし、該注目点から該検索範囲内の各点までのそれぞれの距離に反比例した該重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、該注目点における高さ方向の値ＺｐをＺｐ＝（ｗ１・Ｚｑ１＋ｗ２・Ｚｑ２＋…＋ｗｎ・Ｚｑｎ）で示される内挿計算により求め、（ｆ）該Ｚｐを用いて該第１の断面位置の断面線を出力する。
【０００４】
特許文献１では、指示された断面位置の地盤高の変化を示す断面図を高い精度で生成・出力することができる。しかし、当該断面位置に沿った地盤上にある樹木や建物などの地物の高さの変化を表示できない。
【０００５】
従って、例えば、送電線下付近にある樹木や建物が送電線に触れることなく、適切な距離に保たれているかどうかを確認する場合、作業員が実際に現地に赴いて、地盤上の樹木や建物の高さを確認する必要があった。
【０００６】
また、特許文献２には、支障樹木を自動的に割り出して知らせる送電線下支障樹木表示装置が開示されている。特許文献２には、
・３次元データから、送電線が張られている経路の径間が指定されると、当該経間の樹木情報や送電線の鉄塔位置をデータベースから読み出して、鉄塔間を結ぶ送電線と樹木位置との離隔距離を求めること、
・ｘ−ｙ軸面（平面）、ｘ−ｚ軸面（縦面）、ｙ−ｚ軸面（横面）のうちのいずれかの面上における離隔距離が、数週間以内または数ヶ月以内に送電線に所定以上近接することを表示すること、
・樹木の成長高を年月日と成長率とから計算して、この成長高をデータベースの樹木情報のｚ値に加算することにより、樹木情報を更新すること、
が開示されている。
【０００７】
特許文献２の図６（ｂ）には、送電線が張られている径間の断面に、樹木の高さを示す線（樹高線）が示されている。しかし、特許文献２には、樹高線の生成手法は記載されていない。しかも送電線が張られている径間方向と垂直な方向の地盤高線や樹木の高さの変化を示す樹高線を含む断面図を生成する手法も何ら開示されていない。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特許第３２６７５８９号公報
【特許文献２】特許第３７１４０１７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
引用文献１および２では、送電線が張られている径間方向と垂直な方向に指示された断面線、または平面地図上の任意の方向に指示された断面線に沿った地盤上の樹高線を容易にしかも高精度に生成することができないという問題点があった。
【００１０】
そこで本願発明は、平面地図上で指示された任意の方向の断面線に沿った地盤上の樹高線を容易にしかも高精度に生成することができる送電線下横断面図生成方法および装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
上記課題を解決するために、本発明に係る、送電線下横断面図生成方法は、入力された三次元地図データを用いてｘ−ｙ平面地図を表示し、前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成し、前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成する。
【００１２】
好ましくは、上記樹高線を生成するステップは、（ａ）前記断面線を長さ方向に、指定された内挿間隔で区切った点を注目点とし、各注目点を中心とする検索半径の円からなる検索範囲内を設定し、（ｂ）前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点についての３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記検索範囲内にある複数の頂点を決定し、（ｃ）前記複数の頂点が有する高さ値をｚｑ１，ｚｑ２，…，ｚｑｎとし、前記注目点から前記検索範囲内の各頂点までのそれぞれの距離に反比例した重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、前記注目点における高さ方向の値ｚｐをｚｐ＝（ｗ１・ｚｑ１＋ｗ２・ｚｑ２＋…＋ｗｎ・ｚｑｎ）で示される内挿計算により決定し、（ｄ）決定された各注目点の高さを用いて、前記樹高線を生成する。
【００１３】
あるいは、前記樹高線を生成するステップは、（ａ）前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点についての３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記断面線と前記三角メッシュとが交わる複数の交点を決定し、（ｂ）前記交点が存在する前記三角メッシュ上の三角形の辺の両端の頂点における高さの差を前記両端の頂点間のｘｙ平面上の距離と、前記両端の頂点のうちの一方と前記交点との間のｘｙ平面上の距離との比に応じて比例配分することにより、前記交点における高さを決定し、（ｃ）決定された前記複数の交点の高さを用いて、前記樹高線を生成する。
【００１４】
また、前記地盤高線を生成するステップは、（ａ）前記断面線と等高線とが交わる複数の第１の交点のｘ、ｙ座標と、前記断面線と地物とが交わる第２の交点のｘ、ｙ座標を前記三次元地図データから抽出し、（ｂ）前記第１の交点における高さを前記等高線と等しい値に決定し、（ｄ）前記第２の交点における高さを、当該第２の交点に隣接する２つの第１の交点間の高さの差を前記２つの第１の交点間のｘｙ平面上の距離と、前記２つの第１の交点のうちの一方と前記第２の交点との間のｘｙ平面上の距離との比に応じて比例配分することにより決定し、（ｅ）決定された前記第１の交点の高さおよび前記複数の第２の交点の高さを用いて、前記断面線に沿った地盤高の変化を表す地盤高線を生成する。
【００１５】
本発明によれば、三次元地図データを用いて表示された平面地図の表示画面上でオペレータが計測したい断面線を指示することによって、自動的に三次元地図データから必要な地形データ、三角メッシュデータが抽出され、これに基づいて所望の断面線に沿った地盤高線の他に、高精度の樹高線を生成および出力できる。さらにはこれに関連する付随情報を簡単な手順で出力することができる。
【００１６】
また、断面線上の地盤高線は、断面線と交わる等高線の値を用いた比例配分計算により、地盤高を計算し、当該地盤高線上の樹高線は、三角メッシュ状に配置された各点のもつｘ、ｙ、ｚ座標値を用いた再送計算または比例配分計算により、容易かつ高精度に計算することが可能となる。
【００１７】
さらに、三次元地図データを基にして、断面線と交わる建物の位置や高さを含めた断面図を得ることができる。
【発明の効果】
【００１８】
本発明によれば、平面地図上で指示された断面線に沿った地盤上の樹高線を容易にしかも高精度に生成することができる。
【図面の簡単な説明】
【００１９】
【図１】本発明の一実施形態に係るシステムの構成を示すブロック図。
【図２】本実施形態に係る送電線下横断面図生成手順を示すフローチャート。
【図３】平面地図の一例を示す図。
【図４】平面地図上に指示された断面線と、伸長された断面線の例とを示す図。
【図５】地盤高線を生成するために決定さる断面線と等高線や地物との交点の例を示す図。
【図６】交点の高さ計算について説明するための図。
【図７】平面地図上に配置された三角メッシュを示す図。
【図８】図７に示す三角メッシュを各点の高さを用いて三次元的に示した図。
【図９】三角メッシュと指示された断面線との交点の例を示す図。
【図１０】三角メッシュと断面線との交点Ｊ1〜Ｊ5を模式的に示した図。
【図１１】地盤高線と樹高線とを含む断面図の表示例を示す図。
【図１２】建物が散在する平面地図上に指示された断面線の例を示す図。
【図１３】地盤高線と建物の高さの変化を示す線とを含む断面図の表示例を示す図。
【図１４】他の送電線下横断面図生成手順を示すフローチャート。
【図１５】内挿計算を用いた樹高線の生成方法を説明するための図。
【発明を実施するための形態】
【００２０】
以下、本発明の実施の形態を説明するが、まず最初に本発明の前提となる技術の一例について述べる。航空写真から地図データを取得するために、図化機（実体図化機）と呼ばれる装置が開発されている。図化機は、例えば６０％ずつ重複して撮影された一対の航空写真フィルムを写真架台に載せ、それを接眼鏡を通して左右の目で個別に観察した立体像を作製する装置である。具体的には、まず写真架台上の航空写真フィルムに対して、航空写真フィルムに対する撮影時のカメラの相対位置関係と相似の相対位置関係をもって投影機を設置し、この投影機により空間に地表と同様の立体像を再現する。この操作を相互標定と呼ぶ。
【００２１】
この相互標定後に、絶対標定という作業に移る。絶対標定では、図化機内のデータ取得システムに座標値ｘ，ｙ、および高さ（標高値）ｚの基準点のデータを入力して、航空写真フィルムの立体像を表示するディスプレイの表示画面上に展開させ、この表示画面上の基準点と立体像の中の基準点を一致させる。これらの標定作業は、図化機で定められた観測用測標（メスマーク）を用いて行われる。
【００２２】
図化に際しては、図化するモデル内にｘ，ｙ，ｚが既知である三つ以上の基準点が必要である。このために、航空写真上に地上の基準点が写るように対空標識を設置し、これを基に各モデルに基準点を増設する作業を行う。この作業を空中三角測量といい、コンピュータを使用して解析計算を行う。
【００２３】
道路、建物などの航空写真の図化は、図化しようとする対象に測標をタッチさせ、角部などの変化点毎のポイントデータを取得することにより行われる。等高線については、測標を航空写真上の地上部分にタッチさせて移動させることにより、任意の指定した間隔で自動的にデータが取得される。
【００２４】
このようにして、三次元座標を有する三次元地図データの取得、いわゆるディジタルマッピングを行い、図面化する。なお、ディジタルマッピングとは、既製地図のアナログ情報及び各種の付加情報（地形、土地利用種別等）をディジタル化して電子記憶媒体に記録したもの及び電子記憶媒体から地形図図形として再生出力したもの（建設省・国土地理院による）、あるいは航空写真測量等により、地形・地物等に関わる地図情報をディジタル形式で測定し、コンピュータ技術により体系的に整理されたディジタル地形図を新たに構築する作業をいい、地形図等の原図の作製を含む（建設省公共測量作業規程のうちのディジタルマッピング作業要領による）。
【００２５】
このディジタルマッピングにより得られた三次元座標（ｘ、ｙ、ｚ）を有する三次元地図データを三次元ＣＡＤシステムに取り込み、要求仕様に応じた図面として編集する。建設省で定められた線や記号等はＣＡＤシステムの基本性能だけでは対応できないため、カスタマイズによりＣＡＤシステムに新たな機能を持たせることで対応する。
【００２６】
このカスタマイズは大別して、編集を目的としたものとアプリケーションを目的としたものの二種類があるが、本発明は特に後者のアプリケーションについてのカスタマイズに関する。以下、上記の事項を踏まえて本発明の一実施形態を説明する。
【００２７】
図１は、本発明の一実施形態に係るシステムの構成を示している。このシステムは大きくわけて、図化機１とＣＡＤシステム２及び自動製図機９からなり、これらがバス１０により接続されている。なお、図では各種のインタフェースは省略している。
【００２８】
図化機１は前述した通りのものであり、航空写真フィルム１１を写真架台上に載せ、これを投影機で立体像として投影し（相互標定）、さらに絶対標定を行った後、三次元地図データを取得する。
【００２９】
三次元地図データは、航空写真中の地形や地物をｘ、ｙ座標値および高さ値ｚを有する。例えば、建物などの地物については、その当該地物の最も高い点（標高点）のｘ、ｙ座標と高さ値（標高値）を有する。一方、平面地図（２次元地図）では、地物はｘ、ｙ座標値をもつが、高さ値（標高値）ｚを持たない。
【００３０】
本実施形態の単一断面計測手順では、指示された断面線上の地盤の高さの変化を示す地盤高線とともに、断面線上の樹高や建物の高さの変化を示す線を含む断面図を生成し出力する。
【００３１】
樹高や建物の高さは、例えば、図化機１を用いて、航空写真から測定される。ｘ−ｙ平面上に、相互に接続させた三角形の面のセットとして表される三角メッシュ（図７参照）を配置する。三角形メッシュ中の各三角形の頂点について、ｘ、ｙ座標値と、高さｚとを図化機１を用いて測定する。従って、三角メッシュには、等間隔または任意の間隔（例えば５〜１０ｍ）で、高さｚの測定された点が設定されている。ここでは、ｘ−ｙ平面地図上に三角メッシュ状に設定された各点（各三角形の各頂点）についての座標値（ｘ、ｙ）および高さ値（ｚ）を含むデータを三角メッシュデータと呼ぶ。三角メッシュデータは、３次元地図データのもつデータの一部として作成され、記憶される。
【００３２】
三角メッシュのｘ、ｙ座標と平面地図のｘ、ｙ座標とは、同じ航空写真から測定されたものであるから、両者は一致する。従って、図７に示すように、ｘ−ｙ平面地図上に三角メッシュを配置することができ、平面地図上の三角メッシュの各点に対応するｘ、ｙ座標値に対応する高さ（標高）が得られる。図８は、図７に示す三角メッシュを三角メッシュの各点の高さを用いて三次元的に示したものである。
【００３３】
図化機１で作成された３次元地図データは、ＣＡＤシステム２に入力され、図１の主記憶装置４およびまたは外部記憶装置５に記憶される。
【００３４】
ＣＡＤシステム２は、汎用コンピュータを用いて構成され、具体的にはＣＰＵ３、ＲＯＭやＲＡＭ等の主記憶装置４、ＨＤＤ（ハードディスクドライブ）などの外部記憶装置５、キーボードやポインティングデバイス（例えばマウス）を含む入力装置６、ディスプレイ７及びスキャナ８等からなり、図化機１から入力される３次元地図データを要求仕様に応じた地図データとして編集する。ＣＡＤシステム２で作製された地図データは、自動製図機９により紙ベースの地図として出力される。
【００３５】
ここで、ＣＡＤシステム２は本発明に基づくアプリケーションプログラムの一つとして設定された送電線下横断面図生成プログラムに従って、指示された断面位置の地盤の高さ、樹高や建物などの高さを計算し、ディスプレイ７上に単一断面図として表示したり、自動製図機９を介して紙ベースの単一断面図１２として出力する処理を行う。この送電線下横断面図生成プログラムは外部記憶装置５に予めインストールされており、ここから主記憶装置４にダウンロードされて使用されるものとする。
【００３６】
以下、この送電線下横断面図生成プログラムに基づく送電線下横断面図生成手順について、図２に示すフローチャートを用いて説明する。まず、送電線下横断面図生成に関わる各種パラメータを入力する（ステップＳ１１）。パラメータとしては、例えば、断面図の縦軸、横軸の縮尺、基準標高値Ｈ、後述する樹高線などを生成するために必要なパラメータなどがある。
【００３７】
次に、図化機１から三次元地図データを取り込む（ステップＳ１２）。取り込まれた３次元地図データは、主記憶装置４または外部記憶装置５に記憶される。次に、取り込んだ三次元地図データのうちのｘ，ｙデータを用いて図３に示すようなｘ−ｙ平面地図２０をディスプレイ７で表示する（ステップＳ１３）。
【００３８】
そして、この平面地図２０の表示画面上で横断面図を生成したい断面位置２１をオペレータが指示する（ステップＳ１４）。具体的には、入力装置６に含まれるマウスなどのポインティングデバイスを用いて、平面地図２０上の表示画面上で断面位置２１となる線（断面線）の両端点ａ，ｂを指示することにより、送電線がはられている経路の径間に垂直な方向に断面線２１が指示される。指示された断面線２１は、平面地図２０の表示画面上に表示され、これによりオペレータは指示した断面線（断面位置）２１を確認できる。なお、断面線は、ここでは、径間に垂直な方向に指定されるが、これに限定するものではなく、任意の方向に断面線が指定されても、本実施形態に係る手法が適用可能である。
【００３９】
次に、上記のように指示された断面線２１と交わる地形データを三次元地図データから抽出し、これを用いて、断面線２１に沿った地盤の高さの変化を示す線（地盤高線）を生成する。また、この断面線上の樹木の高さ（樹高）や建物の高さの変化を示す線を生成する。
【００４０】
まず、断面線２１上の地盤高線の生成について説明する。
【００４１】
指示された断面線２１の両端ａ、ｂを拡張して断面線を長く延ばす。図４に示すように伸長後の新たな断面線の両端をａ´、ｂ´とする。以下、両端ａ´、ｂ´を有する断面線を断面線（ａ´、ｂ´）と表す。また、両端ａ、ｂを有する断面線を断面線（ａ、ｂ）と表す。
【００４２】
断面線（ａ´、ｂ´）の位置（ｘ、ｙ）と、等高線や、道路、河川、鉄道などの地物の位置（ｘ、ｙ）とから、断面線（ａ´、ｂ´）と交差する等高線や地物との交点を決定する（ステップＳ１６）。例えば、図５に示すように、断面線（ａ´、ｂ´）と等高線や地物との交点Ｉ1〜Ｉ5が決定される。
【００４３】
そして、ステップＳ１７では、各交点（例えばＩ1〜Ｉ5）の座標（ｘn、ｙn、ｚn）を決定する（ｎは、各交点を識別するインデックスであり、ここでは、例えばｎ＝１，２，…、５）。各交点についてのｘ、ｙ座標値は、３次元地図データの地形データから抽出する。断面線（ａ´、ｂ´）と等高線との交点の高さｚとして、３次元地図データから当該等高線の値を抽出する。３次元地図データから抽出できない道路、河川、鉄道などの地物との交点の高さｚは、等高線の値を用いて計算する。
【００４４】
図５において、断面線と等高線との交点（例えば、Ｉ1、Ｉ3、など）の高さは、その等高線の値とする。また、３次元地図データから抽出できない地物と断面線の交点（例えば、Ｉ2）の高さは、図６に示すように、断面線上で当該交点Ｉ2と隣接する（高さ決定済みの）交点Ｉ1、Ｉ3間の傾斜を考慮して、次のようにして計算する。
【００４５】
図６において、交点Ｉ1の座標を（ｘ1、ｙ1、ｚ1）、交点Ｉ3の座標を（ｘ3、ｙ3、ｚ3）、とする。ここで交点Ｉ1、Ｉ3のｘ、ｙ、ｚ座標値、Ｉ2のｘ、ｙ座標値は上記のように決定されるから、交点Ｉ2のｚ座標は、式（１）に示すように、Ｉ1とＩ3との間のｚ座標値の差分を、Ｉ1とＩ3との間のｘｙ平面上の距離と、Ｉ1とＩ2との間のｘｙ平面上の距離との比に応じて比例配分することで求めることができる。
【数１】
【００４６】
断面線の端点ａ、ｂのｚ座標も、断面線上で隣接する（高さ値が決定済みの）２つの交点のｘ、ｙ、ｚ座標を用いて、上記同様に計算する。
【００４７】
以上のようにして、各交点Ｉ1〜Ｉ5と、断面線の両端点ａ，ｂのｘ、ｙ、ｚ座標値を決定した後、断面線上の中点または始点ａを基準点（０，０）として、各交点Ｉ1〜Ｉ5について、基準点からの距離を計算する。交点Ｉ1〜Ｉ5と、断面線の両端点ａ，ｂを距離の小さいものから順に並べる。並べられた各点の高さ（ｚ）を用いて、断面線２１に沿った地盤高の変化を示す地盤高線１０１を生成する。生成された地盤高線１０１は、図１１に示すように、ディスプレイ７により表示され、さらには、自動製図機９により紙で出力される（ステップＳ１８）。
【００４８】
次に、上記指示された断面線に沿った地盤上に存在する樹木や建物などの地物の高さ（樹高等）を計算し、樹高線を生成する。ここでは、まず、樹高線を生成するための第１の方法について説明する。
【００４９】
ステップＳ２１では、ステップＳ１２で主記憶装置４または外部記憶装置５に記憶された３次元地図データから、現在表示されている平面地図に対応するｘ−ｙ領域の三角メッシュデータを読み出し、当該三角メッシュとステップＳ１４で指示された断面線２１の位置とから、三角メッシュと断面線２１との交点を決定する（図９参照）。
【００５０】
そして、決定された各交点の（ｘ、ｙ、ｚ）を、地盤高の場合と同様にして計算する（ステップＳ２２）。図１０は、三角メッシュと断面線２１との交点Ｊ1〜Ｊ6を模式的に示したものである。前述したように、三角メッシュデータには、三角形の各頂点のｘ、ｙ座標の他に高さｚをもつ。従って、交点が三角メッシュの三角形の頂点と一致する場合には、当該交点の（ｘ、ｙ、ｚ）は、三角メッシュデータに含まれる当該頂点の（ｘ、ｙ、ｚ）値に決定する。
【００５１】
断面線２１と三角メッシュの三角形の辺との交点（例えば、図１０の交点Ｊ1〜Ｊ6）のｘ、ｙ座標値は、３次元地図データから抽出する。また、これらの交点の高さは、当該交点が位置する辺の２つの端点（頂点）のもつ（ｘ、ｙ、ｚ）値を用いて計算する。
【００５２】
例えば、交点Ｊ2は、三角メッシュの頂点Ｏ2、Ｏ3を端点とする辺上にあるから、交点Ｊ2の高さは、この２つの頂点Ｏ2、Ｏ3のｘ、ｙ、ｚを用いて、次のように計算することができる。
【００５３】
頂点Ｏ2の座標を（ｘo2、ｙo2、ｚo2）、頂点Ｏ3の座標を（ｘo3、ｙo3、ｚo3）、交点Ｊ2の座標を（ｘJ2、ｙJ2、ｚJ2）と表し、交点Ｊ2のｚ座標ｚJ2は、次式（２）に示すように、Ｏ2とＯ3との間のｚ座標値の差分を、Ｏ2とＯ3との間のｘｙ平面上の距離と、Ｏ2とＯJ2との間のｘｙ平面上の距離との比に応じて比例配分することで求めることができる。
【数２】
【００５４】
他の交点Ｊ1、Ｊ3、Ｊ4、Ｊ5、Ｊ6についても交点Ｊ2と同様にして、高さを計算できる。すなわち、例えば、交点Ｊ3の場合には、頂点Ｏ3、Ｏ4のｘ、ｙ、ｚ座標値を用いる。上式（２）のＯ2、Ｏ3のｘ、ｙ、ｚ座標値を頂点Ｏ3、Ｏ4のｘ、ｙ、ｚ座標値に置き換えれば、交点Ｊ3の高さが計算できる。また、交点Ｊ4の場合には、頂点Ｏ4、Ｏ5のｘ、ｙ、ｚ座標値を用いる。上式（２）のＯ2、Ｏ3のｘ、ｙ、ｚ座標値を頂点Ｏ4、Ｏ5のｘ、ｙ、ｚ座標値に置き換え、交点Ｊ2の代わりに交点Ｊ4の座標を（ｘJ4、ｙJ4、ｚJ4）と表すと、交点Ｊ4の高さが計算できる。
【００５５】
以上のようにして、各交点のｘ、ｙ、ｚ座標値を決定した後、断面線上の中点または始点ａを基準点（０，０）として、各交点について、基準点からの距離を計算する。そして、交点を当該距離の小さいものから順に並べる。並べられた各点の高さを用いて、断面線２１に沿った樹木や建物などの高さの変化を示す線１０２を生成する。生成された樹木や建物などの高さの変化を示す線１０２は、図１１に示すように、ステップＳ１８で出力される地盤高線１０１に重ねて出力される（ステップＳ２３）。
【００５６】
以上説明したような樹木や建物の高さの変化を表す線を生成・出力する手順は、例えば、ほぼ連続的に樹木や建物などが連なって存在する地域において有効であるが、例えば平坦な（高さ変化の少ない）場所に建物などが散在するような地域においても適用することも可能である。後者（すなわち、例えば建物が散在するような地域）の場合、次に示すような手順で、高さ変化を表す線を生成・出力することもできる。
【００５７】
例えば、図１２に示すような建物が散在するような地域の平面地図上に、断面線２１が指示されたとする。前述したように、ステップＳ１５〜ステップＳ１８で、地盤高線２０１が生成・出力される。
【００５８】
次に、この断面線２１の位置（ｘ、ｙ）と、各建物の位置（ｘ、ｙ）とから、断面線２１と交差する建物（例えば、図１２では、建物４０，４１，４２）と、当該建物と断面線２１との交点を決定する。
【００５９】
そして、図１２において、建物４０の図形の頂点数（ここでは、８）と、各頂点ｃ1〜ｃ8の座標（ｘ、ｙ、ｚ）を３次元地図データから抽出する。建物４０と断面線２１との交点は５１，５２である。各交点の（ｘ、ｙ）は、建物４０の頂点の座標（ｘ、ｙ、ｚ）から計算する。交点５１は辺［ｃ１、ｃ２］上にあるから、頂点ｃ１のｘ、ｙ座標と、頂点ｃ２のｘ、ｙ座標とから交点５１の交点のｘ、ｙ座標を計算する。また、交点５２は辺［ｃ５、ｃ６］上にあるから、頂点ｃ５のｘ、ｙ座標と、頂点ｃ６のｘ、ｙ座標とから交点５２の交点のｘ、ｙ座標を計算する。
【００６０】
３次元地図データには、１つの建物について、その建物の最も高い点（標高点）のｘ、ｙ、ｚ値が含まれている。例えば、図１２において標高点はＨ0〜Ｈ3である。どの標高点が注目の建物（例えばここでは建物４０）の内部の点であるかを決定するために次のような処理を行う。
【００６１】
図１２において、建物４０の各頂点の座標値から、その近傍にある標高点Ｈ0、Ｈ3を検出する。そして、標高点が建物４０の各辺の進行方向に対し、右側にあるか左側にあるかを決定し、当該標高点が右側にある場合と左側にある場合とをそれぞれカウントする。ただし、標高点から各辺に垂線をおろす場合に垂線との交点がないような辺はカウントしない。
【００６２】
例えば、建物４０の場合、辺［ｃ1,ｃ2］の進行方向は、図１２に示されているようにｃ1からｃ２の方向であり、辺［ｃ5,ｃ6］の進行方向は、図１２に示されているようにｃ5からｃ6の方向である。従って、標高点Ｈ0は、辺［ｃ1,ｃ2］、辺［ｃ6,ｃ7］、辺［ｃ7,ｃ8］、辺［ｃ8,ｃ1］の進行方向に対し右側にあるが、標高点Ｈ0が左側にあるような辺はない。同様に、標高点Ｈ3は、辺［ｃ1,ｃ2］の進行方向に対し右側にあるが、辺［ｃ7,ｃ8］の進行方向に対しては左側にある。
【００６３】
すなわち、
標高点Ｈ0の場合、「右側」のカウント数＝４、「左側」のカウント数＝０
標高点Ｈ3の場合、「右側」のカウント数＝１、「左側」のカウント数＝１
である。
【００６４】
「右」のカウント数＞０かつ「左」のカウント数＝０、または、「右」のカウント数＝０かつ「左」のカウント数＞０である場合は、当該標高点は建物内にあり、そうでない場合は、当該標高点は建物の外部にあると決定できる。従って、標高点Ｈ0は建物４０の内部にあり、標高点Ｈ3は建物４０の外部にあると決定できる。
【００６５】
以上のようにして、建物４０の高さは標高点Ｈ0のｚ値（標高値）に決定する。断面線２１と交差する他の建物４１，４２の高さｚについても上記建物４０の場合と同様にして、標高点Ｈ1、Ｈ2のｚ値（標高値）に決定される。
【００６６】
以上のようにして、各建物について、断面線２１との交点の（ｘ、ｙ、ｚ）を決定した後、断面線上の中点または始点ａを基準点（０，０）として、各建物について、基準点からの距離を計算する。そして、図１３に示すように、建物を表す矩形（建物４０に対応する２０２ａ、建物４１に対応する２０２ｂ、建物４２に対応する２０２ｃ）を距離の小さいものから順に、地盤高線２０１上に並べる。建物を表す矩形の幅は、当該建物を通過する断面線２１の長さ（すなわち、当該建物と断面線との２つの交点間の距離）に対応し、矩形の高さは、上記のように決定された当該建物の高さに対応する。
【００６７】
生成された建物の高さ２０２ａ〜２０２ｃの追加された断面図は、図１３に示すように、ディスプレイ７により表示され、さらには、自動製図機９により紙で出力される。
【００６８】
送電線路沿いでは、送電線と樹木や建物との離隔距離が適切に保たれている必要がある。本実施形態によれば、指示された断面線２１が、例えば、図３に示すように、送電線路を横切る場合には、この送電線の位置（高さ）や風等によって送電線が揺れる範囲も断面図に追加して表示する（図１１および１３参照）。このような表示により、送電線と樹木や建物との間に適切な離隔距離が保たれているかどうかを目視で容易に確認できる。
【００６９】
送電線路に関する情報（例えば、鉄塔番号、鉄塔座標および高さ（ｘ、ｙ、ｚ）、送電線の種類、送電電力、最下電線の横振れの中心点の座標等）は、３次元地図データに含まれている。また、送電線の横振は、鉄塔の断面の最下電線の振れの中心となる点を中心とした半円で表示され、この半円の半径（離隔）は、送電線の種類などにより予め定められた値であり、この値も予め主記憶装置４および外部記憶装置５に記憶されていてもよい。
【００７０】
図１１および図１３に示すように、地盤高線や樹高線とともに送電線の横振れ断面を表示する場合には、ユーザは、例えば、図２のステップＳ１１において、送電電圧や、当該送電電圧に対応する送電線との離隔、鉄塔番号や、送電線の線巾などのパラメータを入力する。これにより、上述のように、地盤高線や樹高線（または建物に高さ線）が生成・表示されるとともに、指示された断面線２１と交差する送電線の横振断面１０３、２０３を主記憶装置４および外部記憶装置５に記憶されている情報を基に生成し、図１２や図１３のように表示する（ステップＳ３１）。
【００７１】
さらに、送電線と樹木や建物との離隔距離を図１２や図１３に示すような断面図に、離隔線（１０４、２０４）とともに表示する（ステップＳ３２）。離隔距離は、３次元地図データに含まれる送電線や樹木・建物のｘ、ｙ座標や、上記のように高さｚを決定することにより、計算することができる。
【００７２】
次に、指示された断面線に沿った地盤上の樹高線を生成するための第２の方法を説明する。
【００７３】
上述の第１の方法に従って生成された樹高線は、三角メッシュ上の単位面積当たりの頂点の数が多いほど精度が高くなる。例えば、ステップＳ１４で指示された断面線に沿った地形に、崖、道路、河川その他が存在することにより、地盤の高さ（Ｚ値）に大きな変化があれば、それに伴い樹高も大きく変化する。三角メッシュ上の単位面積当たりの頂点の数が多いほど、この地盤高の大きな変化に伴う樹高の急激な変化も樹高線に反映させることができる。しかし、単位面積当たりの頂点の数が多いほど、各頂点について高さを測定するための手間がかかる。
【００７４】
第１の方法では、断面線と三角メッシュとの交点の高さｚを、当該交点に隣接する２つの頂点の高さを用いて計算していた。一方、第２の方法は、断面線上の（高さを求めるべき）複数の点（例えば、指定された断面線と三角メッシュとの交点）の高さｚを、当該注目点の周辺の複数の点（三角メッシュの頂点または標高点）の高さを用いて内挿計算により求めることで、高精度に樹高線を生成することができる。この第２の方法によれば、三角メッシュの単位面積当たりの頂点の数（あるいは基準とする高さ値が既知の点の数）がそれほど多くなくとも上述のような地盤高の急激な変化に伴う樹高の変化を高精度に表した樹高線を生成できる。
【００７５】
樹高線を生成するための第２の方法を用いた場合の送電線下横断面図生成手順を図１４に示す。図２のステップＳ２１およびＳ２２が、図１４ではステップＳ２０に置き換わっている。以下、図１４において、図２と異なる部分について説明する。
【００７６】
まず、ステップＳ１１では、内挿間隔ｄ、検索半径ｒ、および重みｗを含む内挿パラメータが入力される。以下、ステップＳ２０の処理について図１５を参照して説明する。
【００７７】
図１５に示すように、断面線２１をその長さ方向に内挿間隔ｄで区切った点をそれぞれ注目点とし、各注目点における高さ方向の値ｚを内挿計算により計算する。この場合、まず、主記憶装置４または外部記憶装置５に記憶された３次元地図データから、現在表示されている平面地図に対応するｘ−ｙ領域の三角メッシュデータを読出す。
【００７８】
図１５において、断面線２１上の注目点ｐの高さを計算する場合、この点ｐを中心として半径ｒの検索範囲１４１を設定し、三角メッシュデータから、この検索範囲１４１内にある頂点ｑを決定する。検索範囲１４１内にあるｎ個の三角メッシュの頂点ｑ(例えば、ここでは４個の頂点Ｏ4、Ｏ5、Ｏ6、Ｏ7)のｚ値から、交点ｐのz値を内挿により計算する。この内挿計算に際し、頂点のｚ値に対しては交点ｐからの距離に反比例した重みｗを与える。
【００７９】
この内挿計算を数式で表すと、次式となる。
【００８０】
ｚp＝（ｗ１・ｚq1＋ｗ２・Ｚq2＋…＋ｗｎ・ｚqn） 式（３）
ここで、zpは交点ｐの求めようとするｚ値、ｚｑ１，ｚｑ２，…，ｚｑｎは検索範囲内のｎ個の頂点ｑのｚ値、そしてｗ１，ｗ２，…，ｗｎは点ｐからｎ個の点ｑまでのそれぞれの距離に反比例した重み（重み係数）であり、ｗ１＋ｗ２＋…＋ｗｎ＝１である。すなわち、点ｑのうち点ｐに対してより近い点のｚ値に、より大きな重みを与える。
【００８１】
以上のようにして、各注目点のｘ、ｙ、ｚ座標値を決定した後、ステップＳ２３において、断面線上の中点または始点ａを基準点（０，０）として、各注目点について、基準点からの距離を計算する。そして、注目点を当該距離の小さいものから順に並べる。並べられた各注目点の高さを用いて、断面線２１に沿った樹木や建物などの高さの変化を示す線１０２を生成する。生成された樹高線１０２は、図１１に示すように、ステップＳ１８で出力される地盤高線１０１に重ねて出力される（ステップＳ２３）。
【００８２】
以上の説明では、断面線２１をその長さ方向に内挿間隔ｄで区切った点をそれぞれ注目点として定めて内挿計算により当該注目点の高さ方向の値を計算したが、高さを求めるべき注目点の定め方はこの場合に限らない。例えば、断面線と三角メッシュとの各交点を注目点として定め、内挿計算により当該注目点の高さ方向の値を計算してもよい。この場合、少なくとも、検索半径ｒ、および重みｗを含む内挿パラメータが、図２のステップＳ１１において入力される。そして、ステップＳ２０では、まず、地盤高の場合と同様にして、現在表示されている平面地図に対応するｘ−ｙ領域の三角メッシュとステップＳ１４で指示された断面線２１の位置とから、三角メッシュと断面線２１との交点を決定する（図９参照）。決定された各交点（例えば、図１０の交点Ｊ1〜Ｊ6）のｘ、ｙ座標値は、３次元地図データから抽出する。これらの交点の高さは、当該交点を中心とする半径ｒの円（検索範囲）１４１内の三角メッシュの頂点の高さ値を用いて内挿計算により求める。
【００８３】
例えば、図１５において、断面線２１上の交点ｐ（例えば交点Ｊ5）の高さを計算する場合、この点ｐを中心として半径ｒの検索範囲１４１を設定し、三角メッシュから、この検索範囲１４１内にある頂点ｑを決定する。検索範囲１４１内にあるｎ個の三角メッシュの頂点ｑ(例えば、ここでは４個の頂点Ｏ4、Ｏ5、Ｏ6、Ｏ7)のｚ値から、交点ｐのz値を、式（３）に従って計算する。
【００８４】
以上のようにして、各交点（例えば、図１０の交点Ｊ1〜Ｊ6）のｘ、ｙ、ｚ座標値を決定した後、ステップＳ２３において、断面線上の中点または始点ａを基準点（０，０）として、各交点について、基準点からの距離を計算する。そして、交点を当該距離の小さいものから順に並べる。並べられた各点の高さを用いて、断面線２１に沿った樹木や建物などの高さの変化を示す線１０２を生成する。生成された樹高線１０２は、図１１に示すように、ステップＳ１８で出力される地盤高線１０１に重ねて出力される（ステップＳ２３）。
【００８５】
原理的には、同じ精度の樹高線を生成しようとする場合、三角メッシュの単位面積当たりの頂点の数が多いほど（密であるほど）、検索範囲の大きさｒは小さくてすむ。しかし、三角メッシュの単位面積当たりの頂点の数が相対的に少ない（疎）場合も、検索範囲の大きさｒを大きくしたり、上記のような間隔ｄを小さくすることで、頂点数が密の場合と同様な精度の樹高線を生成することができる。
【技術分野】
【０００１】
本発明は、例えば航空写真から取得した三次元地図データを処理する技術に係り、特に三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成方法およびシステムに関する。
【背景技術】
【０００２】
図化機で取得された３次元地図データから、ＣＡＤ技術を利用して目的とする断面図を作成する方法が特許文献１に開示されている。
【０００３】
特許文献１では、（ａ）内挿間隔と検索半径及び重みを含む内挿パラメータを入力し、（ｂ）入力された三次元地図データを用いて平面地図を表示し、（ｃ）該平面地図の表示画面上で指示された断面位置の両側に各々一つ以上の平行な第２の断面位置を生成し、（ｄ）該三次元地図データから該第１及び第２の断面位置の各断面と交わる各点の地形データを抽出し、（ｅ）抽出された地形データのうち、該第１の断面位置を長さ方向に該内挿間隔で区切った注目点を中心とする該検索半径の円からなる検索範囲内の各点における地形データの高さ方向の値をＺｑ１，Ｚｑ２，…，Ｚｑｎとし、該注目点から該検索範囲内の各点までのそれぞれの距離に反比例した該重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、該注目点における高さ方向の値ＺｐをＺｐ＝（ｗ１・Ｚｑ１＋ｗ２・Ｚｑ２＋…＋ｗｎ・Ｚｑｎ）で示される内挿計算により求め、（ｆ）該Ｚｐを用いて該第１の断面位置の断面線を出力する。
【０００４】
特許文献１では、指示された断面位置の地盤高の変化を示す断面図を高い精度で生成・出力することができる。しかし、当該断面位置に沿った地盤上にある樹木や建物などの地物の高さの変化を表示できない。
【０００５】
従って、例えば、送電線下付近にある樹木や建物が送電線に触れることなく、適切な距離に保たれているかどうかを確認する場合、作業員が実際に現地に赴いて、地盤上の樹木や建物の高さを確認する必要があった。
【０００６】
また、特許文献２には、支障樹木を自動的に割り出して知らせる送電線下支障樹木表示装置が開示されている。特許文献２には、
・３次元データから、送電線が張られている経路の径間が指定されると、当該経間の樹木情報や送電線の鉄塔位置をデータベースから読み出して、鉄塔間を結ぶ送電線と樹木位置との離隔距離を求めること、
・ｘ−ｙ軸面（平面）、ｘ−ｚ軸面（縦面）、ｙ−ｚ軸面（横面）のうちのいずれかの面上における離隔距離が、数週間以内または数ヶ月以内に送電線に所定以上近接することを表示すること、
・樹木の成長高を年月日と成長率とから計算して、この成長高をデータベースの樹木情報のｚ値に加算することにより、樹木情報を更新すること、
が開示されている。
【０００７】
特許文献２の図６（ｂ）には、送電線が張られている径間の断面に、樹木の高さを示す線（樹高線）が示されている。しかし、特許文献２には、樹高線の生成手法は記載されていない。しかも送電線が張られている径間方向と垂直な方向の地盤高線や樹木の高さの変化を示す樹高線を含む断面図を生成する手法も何ら開示されていない。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特許第３２６７５８９号公報
【特許文献２】特許第３７１４０１７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
引用文献１および２では、送電線が張られている径間方向と垂直な方向に指示された断面線、または平面地図上の任意の方向に指示された断面線に沿った地盤上の樹高線を容易にしかも高精度に生成することができないという問題点があった。
【００１０】
そこで本願発明は、平面地図上で指示された任意の方向の断面線に沿った地盤上の樹高線を容易にしかも高精度に生成することができる送電線下横断面図生成方法および装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
上記課題を解決するために、本発明に係る、送電線下横断面図生成方法は、入力された三次元地図データを用いてｘ−ｙ平面地図を表示し、前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成し、前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成する。
【００１２】
好ましくは、上記樹高線を生成するステップは、（ａ）前記断面線を長さ方向に、指定された内挿間隔で区切った点を注目点とし、各注目点を中心とする検索半径の円からなる検索範囲内を設定し、（ｂ）前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点についての３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記検索範囲内にある複数の頂点を決定し、（ｃ）前記複数の頂点が有する高さ値をｚｑ１，ｚｑ２，…，ｚｑｎとし、前記注目点から前記検索範囲内の各頂点までのそれぞれの距離に反比例した重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、前記注目点における高さ方向の値ｚｐをｚｐ＝（ｗ１・ｚｑ１＋ｗ２・ｚｑ２＋…＋ｗｎ・ｚｑｎ）で示される内挿計算により決定し、（ｄ）決定された各注目点の高さを用いて、前記樹高線を生成する。
【００１３】
あるいは、前記樹高線を生成するステップは、（ａ）前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点についての３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記断面線と前記三角メッシュとが交わる複数の交点を決定し、（ｂ）前記交点が存在する前記三角メッシュ上の三角形の辺の両端の頂点における高さの差を前記両端の頂点間のｘｙ平面上の距離と、前記両端の頂点のうちの一方と前記交点との間のｘｙ平面上の距離との比に応じて比例配分することにより、前記交点における高さを決定し、（ｃ）決定された前記複数の交点の高さを用いて、前記樹高線を生成する。
【００１４】
また、前記地盤高線を生成するステップは、（ａ）前記断面線と等高線とが交わる複数の第１の交点のｘ、ｙ座標と、前記断面線と地物とが交わる第２の交点のｘ、ｙ座標を前記三次元地図データから抽出し、（ｂ）前記第１の交点における高さを前記等高線と等しい値に決定し、（ｄ）前記第２の交点における高さを、当該第２の交点に隣接する２つの第１の交点間の高さの差を前記２つの第１の交点間のｘｙ平面上の距離と、前記２つの第１の交点のうちの一方と前記第２の交点との間のｘｙ平面上の距離との比に応じて比例配分することにより決定し、（ｅ）決定された前記第１の交点の高さおよび前記複数の第２の交点の高さを用いて、前記断面線に沿った地盤高の変化を表す地盤高線を生成する。
【００１５】
本発明によれば、三次元地図データを用いて表示された平面地図の表示画面上でオペレータが計測したい断面線を指示することによって、自動的に三次元地図データから必要な地形データ、三角メッシュデータが抽出され、これに基づいて所望の断面線に沿った地盤高線の他に、高精度の樹高線を生成および出力できる。さらにはこれに関連する付随情報を簡単な手順で出力することができる。
【００１６】
また、断面線上の地盤高線は、断面線と交わる等高線の値を用いた比例配分計算により、地盤高を計算し、当該地盤高線上の樹高線は、三角メッシュ状に配置された各点のもつｘ、ｙ、ｚ座標値を用いた再送計算または比例配分計算により、容易かつ高精度に計算することが可能となる。
【００１７】
さらに、三次元地図データを基にして、断面線と交わる建物の位置や高さを含めた断面図を得ることができる。
【発明の効果】
【００１８】
本発明によれば、平面地図上で指示された断面線に沿った地盤上の樹高線を容易にしかも高精度に生成することができる。
【図面の簡単な説明】
【００１９】
【図１】本発明の一実施形態に係るシステムの構成を示すブロック図。
【図２】本実施形態に係る送電線下横断面図生成手順を示すフローチャート。
【図３】平面地図の一例を示す図。
【図４】平面地図上に指示された断面線と、伸長された断面線の例とを示す図。
【図５】地盤高線を生成するために決定さる断面線と等高線や地物との交点の例を示す図。
【図６】交点の高さ計算について説明するための図。
【図７】平面地図上に配置された三角メッシュを示す図。
【図８】図７に示す三角メッシュを各点の高さを用いて三次元的に示した図。
【図９】三角メッシュと指示された断面線との交点の例を示す図。
【図１０】三角メッシュと断面線との交点Ｊ1〜Ｊ5を模式的に示した図。
【図１１】地盤高線と樹高線とを含む断面図の表示例を示す図。
【図１２】建物が散在する平面地図上に指示された断面線の例を示す図。
【図１３】地盤高線と建物の高さの変化を示す線とを含む断面図の表示例を示す図。
【図１４】他の送電線下横断面図生成手順を示すフローチャート。
【図１５】内挿計算を用いた樹高線の生成方法を説明するための図。
【発明を実施するための形態】
【００２０】
以下、本発明の実施の形態を説明するが、まず最初に本発明の前提となる技術の一例について述べる。航空写真から地図データを取得するために、図化機（実体図化機）と呼ばれる装置が開発されている。図化機は、例えば６０％ずつ重複して撮影された一対の航空写真フィルムを写真架台に載せ、それを接眼鏡を通して左右の目で個別に観察した立体像を作製する装置である。具体的には、まず写真架台上の航空写真フィルムに対して、航空写真フィルムに対する撮影時のカメラの相対位置関係と相似の相対位置関係をもって投影機を設置し、この投影機により空間に地表と同様の立体像を再現する。この操作を相互標定と呼ぶ。
【００２１】
この相互標定後に、絶対標定という作業に移る。絶対標定では、図化機内のデータ取得システムに座標値ｘ，ｙ、および高さ（標高値）ｚの基準点のデータを入力して、航空写真フィルムの立体像を表示するディスプレイの表示画面上に展開させ、この表示画面上の基準点と立体像の中の基準点を一致させる。これらの標定作業は、図化機で定められた観測用測標（メスマーク）を用いて行われる。
【００２２】
図化に際しては、図化するモデル内にｘ，ｙ，ｚが既知である三つ以上の基準点が必要である。このために、航空写真上に地上の基準点が写るように対空標識を設置し、これを基に各モデルに基準点を増設する作業を行う。この作業を空中三角測量といい、コンピュータを使用して解析計算を行う。
【００２３】
道路、建物などの航空写真の図化は、図化しようとする対象に測標をタッチさせ、角部などの変化点毎のポイントデータを取得することにより行われる。等高線については、測標を航空写真上の地上部分にタッチさせて移動させることにより、任意の指定した間隔で自動的にデータが取得される。
【００２４】
このようにして、三次元座標を有する三次元地図データの取得、いわゆるディジタルマッピングを行い、図面化する。なお、ディジタルマッピングとは、既製地図のアナログ情報及び各種の付加情報（地形、土地利用種別等）をディジタル化して電子記憶媒体に記録したもの及び電子記憶媒体から地形図図形として再生出力したもの（建設省・国土地理院による）、あるいは航空写真測量等により、地形・地物等に関わる地図情報をディジタル形式で測定し、コンピュータ技術により体系的に整理されたディジタル地形図を新たに構築する作業をいい、地形図等の原図の作製を含む（建設省公共測量作業規程のうちのディジタルマッピング作業要領による）。
【００２５】
このディジタルマッピングにより得られた三次元座標（ｘ、ｙ、ｚ）を有する三次元地図データを三次元ＣＡＤシステムに取り込み、要求仕様に応じた図面として編集する。建設省で定められた線や記号等はＣＡＤシステムの基本性能だけでは対応できないため、カスタマイズによりＣＡＤシステムに新たな機能を持たせることで対応する。
【００２６】
このカスタマイズは大別して、編集を目的としたものとアプリケーションを目的としたものの二種類があるが、本発明は特に後者のアプリケーションについてのカスタマイズに関する。以下、上記の事項を踏まえて本発明の一実施形態を説明する。
【００２７】
図１は、本発明の一実施形態に係るシステムの構成を示している。このシステムは大きくわけて、図化機１とＣＡＤシステム２及び自動製図機９からなり、これらがバス１０により接続されている。なお、図では各種のインタフェースは省略している。
【００２８】
図化機１は前述した通りのものであり、航空写真フィルム１１を写真架台上に載せ、これを投影機で立体像として投影し（相互標定）、さらに絶対標定を行った後、三次元地図データを取得する。
【００２９】
三次元地図データは、航空写真中の地形や地物をｘ、ｙ座標値および高さ値ｚを有する。例えば、建物などの地物については、その当該地物の最も高い点（標高点）のｘ、ｙ座標と高さ値（標高値）を有する。一方、平面地図（２次元地図）では、地物はｘ、ｙ座標値をもつが、高さ値（標高値）ｚを持たない。
【００３０】
本実施形態の単一断面計測手順では、指示された断面線上の地盤の高さの変化を示す地盤高線とともに、断面線上の樹高や建物の高さの変化を示す線を含む断面図を生成し出力する。
【００３１】
樹高や建物の高さは、例えば、図化機１を用いて、航空写真から測定される。ｘ−ｙ平面上に、相互に接続させた三角形の面のセットとして表される三角メッシュ（図７参照）を配置する。三角形メッシュ中の各三角形の頂点について、ｘ、ｙ座標値と、高さｚとを図化機１を用いて測定する。従って、三角メッシュには、等間隔または任意の間隔（例えば５〜１０ｍ）で、高さｚの測定された点が設定されている。ここでは、ｘ−ｙ平面地図上に三角メッシュ状に設定された各点（各三角形の各頂点）についての座標値（ｘ、ｙ）および高さ値（ｚ）を含むデータを三角メッシュデータと呼ぶ。三角メッシュデータは、３次元地図データのもつデータの一部として作成され、記憶される。
【００３２】
三角メッシュのｘ、ｙ座標と平面地図のｘ、ｙ座標とは、同じ航空写真から測定されたものであるから、両者は一致する。従って、図７に示すように、ｘ−ｙ平面地図上に三角メッシュを配置することができ、平面地図上の三角メッシュの各点に対応するｘ、ｙ座標値に対応する高さ（標高）が得られる。図８は、図７に示す三角メッシュを三角メッシュの各点の高さを用いて三次元的に示したものである。
【００３３】
図化機１で作成された３次元地図データは、ＣＡＤシステム２に入力され、図１の主記憶装置４およびまたは外部記憶装置５に記憶される。
【００３４】
ＣＡＤシステム２は、汎用コンピュータを用いて構成され、具体的にはＣＰＵ３、ＲＯＭやＲＡＭ等の主記憶装置４、ＨＤＤ（ハードディスクドライブ）などの外部記憶装置５、キーボードやポインティングデバイス（例えばマウス）を含む入力装置６、ディスプレイ７及びスキャナ８等からなり、図化機１から入力される３次元地図データを要求仕様に応じた地図データとして編集する。ＣＡＤシステム２で作製された地図データは、自動製図機９により紙ベースの地図として出力される。
【００３５】
ここで、ＣＡＤシステム２は本発明に基づくアプリケーションプログラムの一つとして設定された送電線下横断面図生成プログラムに従って、指示された断面位置の地盤の高さ、樹高や建物などの高さを計算し、ディスプレイ７上に単一断面図として表示したり、自動製図機９を介して紙ベースの単一断面図１２として出力する処理を行う。この送電線下横断面図生成プログラムは外部記憶装置５に予めインストールされており、ここから主記憶装置４にダウンロードされて使用されるものとする。
【００３６】
以下、この送電線下横断面図生成プログラムに基づく送電線下横断面図生成手順について、図２に示すフローチャートを用いて説明する。まず、送電線下横断面図生成に関わる各種パラメータを入力する（ステップＳ１１）。パラメータとしては、例えば、断面図の縦軸、横軸の縮尺、基準標高値Ｈ、後述する樹高線などを生成するために必要なパラメータなどがある。
【００３７】
次に、図化機１から三次元地図データを取り込む（ステップＳ１２）。取り込まれた３次元地図データは、主記憶装置４または外部記憶装置５に記憶される。次に、取り込んだ三次元地図データのうちのｘ，ｙデータを用いて図３に示すようなｘ−ｙ平面地図２０をディスプレイ７で表示する（ステップＳ１３）。
【００３８】
そして、この平面地図２０の表示画面上で横断面図を生成したい断面位置２１をオペレータが指示する（ステップＳ１４）。具体的には、入力装置６に含まれるマウスなどのポインティングデバイスを用いて、平面地図２０上の表示画面上で断面位置２１となる線（断面線）の両端点ａ，ｂを指示することにより、送電線がはられている経路の径間に垂直な方向に断面線２１が指示される。指示された断面線２１は、平面地図２０の表示画面上に表示され、これによりオペレータは指示した断面線（断面位置）２１を確認できる。なお、断面線は、ここでは、径間に垂直な方向に指定されるが、これに限定するものではなく、任意の方向に断面線が指定されても、本実施形態に係る手法が適用可能である。
【００３９】
次に、上記のように指示された断面線２１と交わる地形データを三次元地図データから抽出し、これを用いて、断面線２１に沿った地盤の高さの変化を示す線（地盤高線）を生成する。また、この断面線上の樹木の高さ（樹高）や建物の高さの変化を示す線を生成する。
【００４０】
まず、断面線２１上の地盤高線の生成について説明する。
【００４１】
指示された断面線２１の両端ａ、ｂを拡張して断面線を長く延ばす。図４に示すように伸長後の新たな断面線の両端をａ´、ｂ´とする。以下、両端ａ´、ｂ´を有する断面線を断面線（ａ´、ｂ´）と表す。また、両端ａ、ｂを有する断面線を断面線（ａ、ｂ）と表す。
【００４２】
断面線（ａ´、ｂ´）の位置（ｘ、ｙ）と、等高線や、道路、河川、鉄道などの地物の位置（ｘ、ｙ）とから、断面線（ａ´、ｂ´）と交差する等高線や地物との交点を決定する（ステップＳ１６）。例えば、図５に示すように、断面線（ａ´、ｂ´）と等高線や地物との交点Ｉ1〜Ｉ5が決定される。
【００４３】
そして、ステップＳ１７では、各交点（例えばＩ1〜Ｉ5）の座標（ｘn、ｙn、ｚn）を決定する（ｎは、各交点を識別するインデックスであり、ここでは、例えばｎ＝１，２，…、５）。各交点についてのｘ、ｙ座標値は、３次元地図データの地形データから抽出する。断面線（ａ´、ｂ´）と等高線との交点の高さｚとして、３次元地図データから当該等高線の値を抽出する。３次元地図データから抽出できない道路、河川、鉄道などの地物との交点の高さｚは、等高線の値を用いて計算する。
【００４４】
図５において、断面線と等高線との交点（例えば、Ｉ1、Ｉ3、など）の高さは、その等高線の値とする。また、３次元地図データから抽出できない地物と断面線の交点（例えば、Ｉ2）の高さは、図６に示すように、断面線上で当該交点Ｉ2と隣接する（高さ決定済みの）交点Ｉ1、Ｉ3間の傾斜を考慮して、次のようにして計算する。
【００４５】
図６において、交点Ｉ1の座標を（ｘ1、ｙ1、ｚ1）、交点Ｉ3の座標を（ｘ3、ｙ3、ｚ3）、とする。ここで交点Ｉ1、Ｉ3のｘ、ｙ、ｚ座標値、Ｉ2のｘ、ｙ座標値は上記のように決定されるから、交点Ｉ2のｚ座標は、式（１）に示すように、Ｉ1とＩ3との間のｚ座標値の差分を、Ｉ1とＩ3との間のｘｙ平面上の距離と、Ｉ1とＩ2との間のｘｙ平面上の距離との比に応じて比例配分することで求めることができる。
【数１】
【００４６】
断面線の端点ａ、ｂのｚ座標も、断面線上で隣接する（高さ値が決定済みの）２つの交点のｘ、ｙ、ｚ座標を用いて、上記同様に計算する。
【００４７】
以上のようにして、各交点Ｉ1〜Ｉ5と、断面線の両端点ａ，ｂのｘ、ｙ、ｚ座標値を決定した後、断面線上の中点または始点ａを基準点（０，０）として、各交点Ｉ1〜Ｉ5について、基準点からの距離を計算する。交点Ｉ1〜Ｉ5と、断面線の両端点ａ，ｂを距離の小さいものから順に並べる。並べられた各点の高さ（ｚ）を用いて、断面線２１に沿った地盤高の変化を示す地盤高線１０１を生成する。生成された地盤高線１０１は、図１１に示すように、ディスプレイ７により表示され、さらには、自動製図機９により紙で出力される（ステップＳ１８）。
【００４８】
次に、上記指示された断面線に沿った地盤上に存在する樹木や建物などの地物の高さ（樹高等）を計算し、樹高線を生成する。ここでは、まず、樹高線を生成するための第１の方法について説明する。
【００４９】
ステップＳ２１では、ステップＳ１２で主記憶装置４または外部記憶装置５に記憶された３次元地図データから、現在表示されている平面地図に対応するｘ−ｙ領域の三角メッシュデータを読み出し、当該三角メッシュとステップＳ１４で指示された断面線２１の位置とから、三角メッシュと断面線２１との交点を決定する（図９参照）。
【００５０】
そして、決定された各交点の（ｘ、ｙ、ｚ）を、地盤高の場合と同様にして計算する（ステップＳ２２）。図１０は、三角メッシュと断面線２１との交点Ｊ1〜Ｊ6を模式的に示したものである。前述したように、三角メッシュデータには、三角形の各頂点のｘ、ｙ座標の他に高さｚをもつ。従って、交点が三角メッシュの三角形の頂点と一致する場合には、当該交点の（ｘ、ｙ、ｚ）は、三角メッシュデータに含まれる当該頂点の（ｘ、ｙ、ｚ）値に決定する。
【００５１】
断面線２１と三角メッシュの三角形の辺との交点（例えば、図１０の交点Ｊ1〜Ｊ6）のｘ、ｙ座標値は、３次元地図データから抽出する。また、これらの交点の高さは、当該交点が位置する辺の２つの端点（頂点）のもつ（ｘ、ｙ、ｚ）値を用いて計算する。
【００５２】
例えば、交点Ｊ2は、三角メッシュの頂点Ｏ2、Ｏ3を端点とする辺上にあるから、交点Ｊ2の高さは、この２つの頂点Ｏ2、Ｏ3のｘ、ｙ、ｚを用いて、次のように計算することができる。
【００５３】
頂点Ｏ2の座標を（ｘo2、ｙo2、ｚo2）、頂点Ｏ3の座標を（ｘo3、ｙo3、ｚo3）、交点Ｊ2の座標を（ｘJ2、ｙJ2、ｚJ2）と表し、交点Ｊ2のｚ座標ｚJ2は、次式（２）に示すように、Ｏ2とＯ3との間のｚ座標値の差分を、Ｏ2とＯ3との間のｘｙ平面上の距離と、Ｏ2とＯJ2との間のｘｙ平面上の距離との比に応じて比例配分することで求めることができる。
【数２】
【００５４】
他の交点Ｊ1、Ｊ3、Ｊ4、Ｊ5、Ｊ6についても交点Ｊ2と同様にして、高さを計算できる。すなわち、例えば、交点Ｊ3の場合には、頂点Ｏ3、Ｏ4のｘ、ｙ、ｚ座標値を用いる。上式（２）のＯ2、Ｏ3のｘ、ｙ、ｚ座標値を頂点Ｏ3、Ｏ4のｘ、ｙ、ｚ座標値に置き換えれば、交点Ｊ3の高さが計算できる。また、交点Ｊ4の場合には、頂点Ｏ4、Ｏ5のｘ、ｙ、ｚ座標値を用いる。上式（２）のＯ2、Ｏ3のｘ、ｙ、ｚ座標値を頂点Ｏ4、Ｏ5のｘ、ｙ、ｚ座標値に置き換え、交点Ｊ2の代わりに交点Ｊ4の座標を（ｘJ4、ｙJ4、ｚJ4）と表すと、交点Ｊ4の高さが計算できる。
【００５５】
以上のようにして、各交点のｘ、ｙ、ｚ座標値を決定した後、断面線上の中点または始点ａを基準点（０，０）として、各交点について、基準点からの距離を計算する。そして、交点を当該距離の小さいものから順に並べる。並べられた各点の高さを用いて、断面線２１に沿った樹木や建物などの高さの変化を示す線１０２を生成する。生成された樹木や建物などの高さの変化を示す線１０２は、図１１に示すように、ステップＳ１８で出力される地盤高線１０１に重ねて出力される（ステップＳ２３）。
【００５６】
以上説明したような樹木や建物の高さの変化を表す線を生成・出力する手順は、例えば、ほぼ連続的に樹木や建物などが連なって存在する地域において有効であるが、例えば平坦な（高さ変化の少ない）場所に建物などが散在するような地域においても適用することも可能である。後者（すなわち、例えば建物が散在するような地域）の場合、次に示すような手順で、高さ変化を表す線を生成・出力することもできる。
【００５７】
例えば、図１２に示すような建物が散在するような地域の平面地図上に、断面線２１が指示されたとする。前述したように、ステップＳ１５〜ステップＳ１８で、地盤高線２０１が生成・出力される。
【００５８】
次に、この断面線２１の位置（ｘ、ｙ）と、各建物の位置（ｘ、ｙ）とから、断面線２１と交差する建物（例えば、図１２では、建物４０，４１，４２）と、当該建物と断面線２１との交点を決定する。
【００５９】
そして、図１２において、建物４０の図形の頂点数（ここでは、８）と、各頂点ｃ1〜ｃ8の座標（ｘ、ｙ、ｚ）を３次元地図データから抽出する。建物４０と断面線２１との交点は５１，５２である。各交点の（ｘ、ｙ）は、建物４０の頂点の座標（ｘ、ｙ、ｚ）から計算する。交点５１は辺［ｃ１、ｃ２］上にあるから、頂点ｃ１のｘ、ｙ座標と、頂点ｃ２のｘ、ｙ座標とから交点５１の交点のｘ、ｙ座標を計算する。また、交点５２は辺［ｃ５、ｃ６］上にあるから、頂点ｃ５のｘ、ｙ座標と、頂点ｃ６のｘ、ｙ座標とから交点５２の交点のｘ、ｙ座標を計算する。
【００６０】
３次元地図データには、１つの建物について、その建物の最も高い点（標高点）のｘ、ｙ、ｚ値が含まれている。例えば、図１２において標高点はＨ0〜Ｈ3である。どの標高点が注目の建物（例えばここでは建物４０）の内部の点であるかを決定するために次のような処理を行う。
【００６１】
図１２において、建物４０の各頂点の座標値から、その近傍にある標高点Ｈ0、Ｈ3を検出する。そして、標高点が建物４０の各辺の進行方向に対し、右側にあるか左側にあるかを決定し、当該標高点が右側にある場合と左側にある場合とをそれぞれカウントする。ただし、標高点から各辺に垂線をおろす場合に垂線との交点がないような辺はカウントしない。
【００６２】
例えば、建物４０の場合、辺［ｃ1,ｃ2］の進行方向は、図１２に示されているようにｃ1からｃ２の方向であり、辺［ｃ5,ｃ6］の進行方向は、図１２に示されているようにｃ5からｃ6の方向である。従って、標高点Ｈ0は、辺［ｃ1,ｃ2］、辺［ｃ6,ｃ7］、辺［ｃ7,ｃ8］、辺［ｃ8,ｃ1］の進行方向に対し右側にあるが、標高点Ｈ0が左側にあるような辺はない。同様に、標高点Ｈ3は、辺［ｃ1,ｃ2］の進行方向に対し右側にあるが、辺［ｃ7,ｃ8］の進行方向に対しては左側にある。
【００６３】
すなわち、
標高点Ｈ0の場合、「右側」のカウント数＝４、「左側」のカウント数＝０
標高点Ｈ3の場合、「右側」のカウント数＝１、「左側」のカウント数＝１
である。
【００６４】
「右」のカウント数＞０かつ「左」のカウント数＝０、または、「右」のカウント数＝０かつ「左」のカウント数＞０である場合は、当該標高点は建物内にあり、そうでない場合は、当該標高点は建物の外部にあると決定できる。従って、標高点Ｈ0は建物４０の内部にあり、標高点Ｈ3は建物４０の外部にあると決定できる。
【００６５】
以上のようにして、建物４０の高さは標高点Ｈ0のｚ値（標高値）に決定する。断面線２１と交差する他の建物４１，４２の高さｚについても上記建物４０の場合と同様にして、標高点Ｈ1、Ｈ2のｚ値（標高値）に決定される。
【００６６】
以上のようにして、各建物について、断面線２１との交点の（ｘ、ｙ、ｚ）を決定した後、断面線上の中点または始点ａを基準点（０，０）として、各建物について、基準点からの距離を計算する。そして、図１３に示すように、建物を表す矩形（建物４０に対応する２０２ａ、建物４１に対応する２０２ｂ、建物４２に対応する２０２ｃ）を距離の小さいものから順に、地盤高線２０１上に並べる。建物を表す矩形の幅は、当該建物を通過する断面線２１の長さ（すなわち、当該建物と断面線との２つの交点間の距離）に対応し、矩形の高さは、上記のように決定された当該建物の高さに対応する。
【００６７】
生成された建物の高さ２０２ａ〜２０２ｃの追加された断面図は、図１３に示すように、ディスプレイ７により表示され、さらには、自動製図機９により紙で出力される。
【００６８】
送電線路沿いでは、送電線と樹木や建物との離隔距離が適切に保たれている必要がある。本実施形態によれば、指示された断面線２１が、例えば、図３に示すように、送電線路を横切る場合には、この送電線の位置（高さ）や風等によって送電線が揺れる範囲も断面図に追加して表示する（図１１および１３参照）。このような表示により、送電線と樹木や建物との間に適切な離隔距離が保たれているかどうかを目視で容易に確認できる。
【００６９】
送電線路に関する情報（例えば、鉄塔番号、鉄塔座標および高さ（ｘ、ｙ、ｚ）、送電線の種類、送電電力、最下電線の横振れの中心点の座標等）は、３次元地図データに含まれている。また、送電線の横振は、鉄塔の断面の最下電線の振れの中心となる点を中心とした半円で表示され、この半円の半径（離隔）は、送電線の種類などにより予め定められた値であり、この値も予め主記憶装置４および外部記憶装置５に記憶されていてもよい。
【００７０】
図１１および図１３に示すように、地盤高線や樹高線とともに送電線の横振れ断面を表示する場合には、ユーザは、例えば、図２のステップＳ１１において、送電電圧や、当該送電電圧に対応する送電線との離隔、鉄塔番号や、送電線の線巾などのパラメータを入力する。これにより、上述のように、地盤高線や樹高線（または建物に高さ線）が生成・表示されるとともに、指示された断面線２１と交差する送電線の横振断面１０３、２０３を主記憶装置４および外部記憶装置５に記憶されている情報を基に生成し、図１２や図１３のように表示する（ステップＳ３１）。
【００７１】
さらに、送電線と樹木や建物との離隔距離を図１２や図１３に示すような断面図に、離隔線（１０４、２０４）とともに表示する（ステップＳ３２）。離隔距離は、３次元地図データに含まれる送電線や樹木・建物のｘ、ｙ座標や、上記のように高さｚを決定することにより、計算することができる。
【００７２】
次に、指示された断面線に沿った地盤上の樹高線を生成するための第２の方法を説明する。
【００７３】
上述の第１の方法に従って生成された樹高線は、三角メッシュ上の単位面積当たりの頂点の数が多いほど精度が高くなる。例えば、ステップＳ１４で指示された断面線に沿った地形に、崖、道路、河川その他が存在することにより、地盤の高さ（Ｚ値）に大きな変化があれば、それに伴い樹高も大きく変化する。三角メッシュ上の単位面積当たりの頂点の数が多いほど、この地盤高の大きな変化に伴う樹高の急激な変化も樹高線に反映させることができる。しかし、単位面積当たりの頂点の数が多いほど、各頂点について高さを測定するための手間がかかる。
【００７４】
第１の方法では、断面線と三角メッシュとの交点の高さｚを、当該交点に隣接する２つの頂点の高さを用いて計算していた。一方、第２の方法は、断面線上の（高さを求めるべき）複数の点（例えば、指定された断面線と三角メッシュとの交点）の高さｚを、当該注目点の周辺の複数の点（三角メッシュの頂点または標高点）の高さを用いて内挿計算により求めることで、高精度に樹高線を生成することができる。この第２の方法によれば、三角メッシュの単位面積当たりの頂点の数（あるいは基準とする高さ値が既知の点の数）がそれほど多くなくとも上述のような地盤高の急激な変化に伴う樹高の変化を高精度に表した樹高線を生成できる。
【００７５】
樹高線を生成するための第２の方法を用いた場合の送電線下横断面図生成手順を図１４に示す。図２のステップＳ２１およびＳ２２が、図１４ではステップＳ２０に置き換わっている。以下、図１４において、図２と異なる部分について説明する。
【００７６】
まず、ステップＳ１１では、内挿間隔ｄ、検索半径ｒ、および重みｗを含む内挿パラメータが入力される。以下、ステップＳ２０の処理について図１５を参照して説明する。
【００７７】
図１５に示すように、断面線２１をその長さ方向に内挿間隔ｄで区切った点をそれぞれ注目点とし、各注目点における高さ方向の値ｚを内挿計算により計算する。この場合、まず、主記憶装置４または外部記憶装置５に記憶された３次元地図データから、現在表示されている平面地図に対応するｘ−ｙ領域の三角メッシュデータを読出す。
【００７８】
図１５において、断面線２１上の注目点ｐの高さを計算する場合、この点ｐを中心として半径ｒの検索範囲１４１を設定し、三角メッシュデータから、この検索範囲１４１内にある頂点ｑを決定する。検索範囲１４１内にあるｎ個の三角メッシュの頂点ｑ(例えば、ここでは４個の頂点Ｏ4、Ｏ5、Ｏ6、Ｏ7)のｚ値から、交点ｐのz値を内挿により計算する。この内挿計算に際し、頂点のｚ値に対しては交点ｐからの距離に反比例した重みｗを与える。
【００７９】
この内挿計算を数式で表すと、次式となる。
【００８０】
ｚp＝（ｗ１・ｚq1＋ｗ２・Ｚq2＋…＋ｗｎ・ｚqn） 式（３）
ここで、zpは交点ｐの求めようとするｚ値、ｚｑ１，ｚｑ２，…，ｚｑｎは検索範囲内のｎ個の頂点ｑのｚ値、そしてｗ１，ｗ２，…，ｗｎは点ｐからｎ個の点ｑまでのそれぞれの距離に反比例した重み（重み係数）であり、ｗ１＋ｗ２＋…＋ｗｎ＝１である。すなわち、点ｑのうち点ｐに対してより近い点のｚ値に、より大きな重みを与える。
【００８１】
以上のようにして、各注目点のｘ、ｙ、ｚ座標値を決定した後、ステップＳ２３において、断面線上の中点または始点ａを基準点（０，０）として、各注目点について、基準点からの距離を計算する。そして、注目点を当該距離の小さいものから順に並べる。並べられた各注目点の高さを用いて、断面線２１に沿った樹木や建物などの高さの変化を示す線１０２を生成する。生成された樹高線１０２は、図１１に示すように、ステップＳ１８で出力される地盤高線１０１に重ねて出力される（ステップＳ２３）。
【００８２】
以上の説明では、断面線２１をその長さ方向に内挿間隔ｄで区切った点をそれぞれ注目点として定めて内挿計算により当該注目点の高さ方向の値を計算したが、高さを求めるべき注目点の定め方はこの場合に限らない。例えば、断面線と三角メッシュとの各交点を注目点として定め、内挿計算により当該注目点の高さ方向の値を計算してもよい。この場合、少なくとも、検索半径ｒ、および重みｗを含む内挿パラメータが、図２のステップＳ１１において入力される。そして、ステップＳ２０では、まず、地盤高の場合と同様にして、現在表示されている平面地図に対応するｘ−ｙ領域の三角メッシュとステップＳ１４で指示された断面線２１の位置とから、三角メッシュと断面線２１との交点を決定する（図９参照）。決定された各交点（例えば、図１０の交点Ｊ1〜Ｊ6）のｘ、ｙ座標値は、３次元地図データから抽出する。これらの交点の高さは、当該交点を中心とする半径ｒの円（検索範囲）１４１内の三角メッシュの頂点の高さ値を用いて内挿計算により求める。
【００８３】
例えば、図１５において、断面線２１上の交点ｐ（例えば交点Ｊ5）の高さを計算する場合、この点ｐを中心として半径ｒの検索範囲１４１を設定し、三角メッシュから、この検索範囲１４１内にある頂点ｑを決定する。検索範囲１４１内にあるｎ個の三角メッシュの頂点ｑ(例えば、ここでは４個の頂点Ｏ4、Ｏ5、Ｏ6、Ｏ7)のｚ値から、交点ｐのz値を、式（３）に従って計算する。
【００８４】
以上のようにして、各交点（例えば、図１０の交点Ｊ1〜Ｊ6）のｘ、ｙ、ｚ座標値を決定した後、ステップＳ２３において、断面線上の中点または始点ａを基準点（０，０）として、各交点について、基準点からの距離を計算する。そして、交点を当該距離の小さいものから順に並べる。並べられた各点の高さを用いて、断面線２１に沿った樹木や建物などの高さの変化を示す線１０２を生成する。生成された樹高線１０２は、図１１に示すように、ステップＳ１８で出力される地盤高線１０１に重ねて出力される（ステップＳ２３）。
【００８５】
原理的には、同じ精度の樹高線を生成しようとする場合、三角メッシュの単位面積当たりの頂点の数が多いほど（密であるほど）、検索範囲の大きさｒは小さくてすむ。しかし、三角メッシュの単位面積当たりの頂点の数が相対的に少ない（疎）場合も、検索範囲の大きさｒを大きくしたり、上記のような間隔ｄを小さくすることで、頂点数が密の場合と同様な精度の樹高線を生成することができる。
【特許請求の範囲】
【請求項１】
三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成方法であって、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示するステップと、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成するステップと、
を含み、
前記樹高線を生成するステップは、
前記断面線を長さ方向に、指定された内挿間隔で区切った点を注目点とし、各注目点を中心とする検索半径の円からなる検索範囲内を設定するステップと、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記検索範囲内にある複数の頂点を決定するステップと、
前記複数の頂点が有する高さ値をｚｑ１，ｚｑ２，…，ｚｑｎとし、前記注目点から前記検索範囲内の各頂点までのそれぞれの距離に反比例した重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、前記注目点における高さ方向の値ｚｐを
ｚｐ＝（ｗ１・ｚｑ１＋ｗ２・ｚｑ２＋…＋ｗｎ・ｚｑｎ）で示される内挿計算により決定するステップと、
決定された各注目点の高さを用いて、前記樹高線を生成するステップと、
を含む送電線下横断面図生成方法。
【請求項２】
三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成方法であって、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示するステップと、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成するステップと、
を含み、
前記樹高線を生成するステップは、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記断面線と前記三角メッシュとが交わる複数の交点を決定するステップと、
前記交点が存在する前記三角メッシュ上の三角形の辺の両端の頂点における高さの差を前記両端の頂点間のｘｙ平面上の距離と、前記両端の頂点のうちの一方と前記交点との間のｘｙ平面上の距離との比に応じて比例配分することにより、前記交点における高さを決定するステップと、
決定された前記複数の交点の高さを用いて、前記樹高線を生成するステップと、
を含む送電線下横断面図生成方法。
【請求項３】
前記地盤高線を生成するステップは、
前記断面線と等高線とが交わる複数の第１の交点のｘ、ｙ座標と、前記断面線と地物とが交わる第２の交点のｘ、ｙ座標を前記三次元地図データから抽出するステップと、
前記第１の交点における高さを前記等高線と等しい値に決定するステップと、
前記第２の交点における高さを、当該第２の交点に隣接する２つの第１の交点間の高さの差を前記２つの第１の交点間のｘｙ平面上の距離と、前記２つの第１の交点のうちの一方と前記第２の交点との間のｘｙ平面上の距離との比に応じて比例配分することにより決定するステップと、
決定された前記第１の交点の高さおよび前記複数の第２の交点の高さを用いて、前記断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
を含む請求項１または２記載の送電線下横断面図生成方法。
【請求項４】
三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成装置であって、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示する手段と、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成する手段と、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成する手段と、
を備え、
前記樹高線を生成する手段は、
前記断面線を長さ方向に、指定された内挿間隔で区切った点を注目点とし、各注目点を中心とする検索半径の円からなる検索範囲内を設定し、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記検索範囲内にある複数の頂点を決定し、
前記複数の頂点が有する高さ値をｚｑ１，ｚｑ２，…，ｚｑｎとし、前記注目点から前記検索範囲内の各頂点までのそれぞれの距離に反比例した重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、前記注目点における高さ方向の値ｚｐを
ｚｐ＝（ｗ１・ｚｑ１＋ｗ２・ｚｑ２＋…＋ｗｎ・ｚｑｎ）で示される内挿計算により決定し、
決定された各注目点の高さを用いて、前記樹高線を生成する、
送電線下横断面図生成装置。
【請求項５】
三次元地図データから送電線下の横断面を生成する送電線下横断面図生成装置であって、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示する手段と、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成する手段と、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成する手段と、
を備え、
前記樹高線を生成する手段は、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記断面線と前記三角メッシュとが交わる複数の交点を決定し、
前記交点が存在する前記三角メッシュ上の三角形の辺の両端の頂点における高さの差を前記両端の頂点間のｘｙ平面上の距離と、前記両端の頂点のうちの一方と前記交点との間のｘｙ平面上の距離との比に応じて比例配分することにより、前記交点における高さを決定し、
決定された前記複数の交点の高さを用いて、前記樹高線を生成する、
送電線下横断面図生成装置。
【請求項６】
前記地盤高線を生成する手段は、
前記断面線と等高線とが交わる複数の第１の交点のｘ、ｙ座標と、前記断面線と地物とが交わる第２の交点のｘ、ｙ座標を前記三次元地図データから抽出し、
前記第１の交点における高さを前記等高線と等しい値に決定し、
前記第２の交点における高さを、当該第２の交点に隣接する２つの第１の交点間の高さの差を前記２つの第１の交点間のｘｙ平面上の距離と、前記２つの第１の交点のうちの一方と前記第２の交点との間のｘｙ平面上の距離との比に応じて比例配分することにより決定し、
決定された前記第１の交点の高さおよび前記複数の第２の交点の高さを用いて、前記地盤高線を生成する、
請求項４または５記載の送電線下横断面図生成装置。
【請求項７】
三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成プログラムを記憶するコンピュータ読取可能な記憶媒体であって、
前記送電線下横断面図生成プログラムは、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示するステップと、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成するステップと、
を含み、
前記樹高線を生成するステップは、
前記断面線を長さ方向に、指定された内挿間隔で区切った点を注目点とし、各注目点を中心とする検索半径の円からなる検索範囲内を設定するステップと、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記検索範囲内にある複数の頂点を決定するステップと、
前記複数の頂点が有する高さ値をｚｑ１，ｚｑ２，…，ｚｑｎとし、前記注目点から前記検索範囲内の各頂点までのそれぞれの距離に反比例した重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、前記注目点における高さ方向の値ｚｐを
ｚｐ＝（ｗ１・ｚｑ１＋ｗ２・ｚｑ２＋…＋ｗｎ・ｚｑｎ）で示される内挿計算により決定するステップと、
決定された各注目点の高さを用いて、前記樹高線を生成するステップと、
を含むコンピュータ読取可能な記憶媒体。
【請求項８】
三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成プログラムを記憶したコンピュータ読取可能な記憶媒体であって、
前記送電線下横断面図生成プログラムは、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示するステップと、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成するステップと、
を含み、
前記樹高線を生成するステップは、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記断面線と前記三角メッシュとが交わる複数の交点を決定するステップと、
前記交点が存在する前記三角メッシュ上の三角形の辺の両端の頂点における高さの差を前記両端の頂点間のｘｙ平面上の距離と、前記両端の頂点のうちの一方と前記交点との間のｘｙ平面上の距離との比に応じて比例配分することにより、前記交点における高さを決定するステップと、
決定された前記複数の交点の高さを用いて、前記樹高線を生成するステップと、
を含む、コンピュータ読取可能な記憶媒体。
【請求項９】
前記地盤高線を生成するステップは、
前記断面線と等高線とが交わる複数の第１の交点のｘ、ｙ座標と、前記断面線と地物とが交わる第２の交点のｘ、ｙ座標を前記三次元地図データから抽出するステップと、
前記第１の交点における高さを前記等高線と等しい値に決定するステップと、
前記第２の交点における高さを、当該第２の交点に隣接する２つの第１の交点間の高さの差を前記２つの第１の交点間のｘｙ平面上の距離と、前記２つの第１の交点のうちの一方と前記第２の交点との間のｘｙ平面上の距離との比に応じて比例配分することにより決定するステップと、
決定された前記第１の交点の高さおよび前記複数の第２の交点の高さを用いて、前記断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
を含む請求項７または８記載のコンピュータ読取可能な記憶媒体。
【請求項１】
三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成方法であって、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示するステップと、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成するステップと、
を含み、
前記樹高線を生成するステップは、
前記断面線を長さ方向に、指定された内挿間隔で区切った点を注目点とし、各注目点を中心とする検索半径の円からなる検索範囲内を設定するステップと、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記検索範囲内にある複数の頂点を決定するステップと、
前記複数の頂点が有する高さ値をｚｑ１，ｚｑ２，…，ｚｑｎとし、前記注目点から前記検索範囲内の各頂点までのそれぞれの距離に反比例した重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、前記注目点における高さ方向の値ｚｐを
ｚｐ＝（ｗ１・ｚｑ１＋ｗ２・ｚｑ２＋…＋ｗｎ・ｚｑｎ）で示される内挿計算により決定するステップと、
決定された各注目点の高さを用いて、前記樹高線を生成するステップと、
を含む送電線下横断面図生成方法。
【請求項２】
三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成方法であって、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示するステップと、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成するステップと、
を含み、
前記樹高線を生成するステップは、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記断面線と前記三角メッシュとが交わる複数の交点を決定するステップと、
前記交点が存在する前記三角メッシュ上の三角形の辺の両端の頂点における高さの差を前記両端の頂点間のｘｙ平面上の距離と、前記両端の頂点のうちの一方と前記交点との間のｘｙ平面上の距離との比に応じて比例配分することにより、前記交点における高さを決定するステップと、
決定された前記複数の交点の高さを用いて、前記樹高線を生成するステップと、
を含む送電線下横断面図生成方法。
【請求項３】
前記地盤高線を生成するステップは、
前記断面線と等高線とが交わる複数の第１の交点のｘ、ｙ座標と、前記断面線と地物とが交わる第２の交点のｘ、ｙ座標を前記三次元地図データから抽出するステップと、
前記第１の交点における高さを前記等高線と等しい値に決定するステップと、
前記第２の交点における高さを、当該第２の交点に隣接する２つの第１の交点間の高さの差を前記２つの第１の交点間のｘｙ平面上の距離と、前記２つの第１の交点のうちの一方と前記第２の交点との間のｘｙ平面上の距離との比に応じて比例配分することにより決定するステップと、
決定された前記第１の交点の高さおよび前記複数の第２の交点の高さを用いて、前記断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
を含む請求項１または２記載の送電線下横断面図生成方法。
【請求項４】
三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成装置であって、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示する手段と、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成する手段と、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成する手段と、
を備え、
前記樹高線を生成する手段は、
前記断面線を長さ方向に、指定された内挿間隔で区切った点を注目点とし、各注目点を中心とする検索半径の円からなる検索範囲内を設定し、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記検索範囲内にある複数の頂点を決定し、
前記複数の頂点が有する高さ値をｚｑ１，ｚｑ２，…，ｚｑｎとし、前記注目点から前記検索範囲内の各頂点までのそれぞれの距離に反比例した重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、前記注目点における高さ方向の値ｚｐを
ｚｐ＝（ｗ１・ｚｑ１＋ｗ２・ｚｑ２＋…＋ｗｎ・ｚｑｎ）で示される内挿計算により決定し、
決定された各注目点の高さを用いて、前記樹高線を生成する、
送電線下横断面図生成装置。
【請求項５】
三次元地図データから送電線下の横断面を生成する送電線下横断面図生成装置であって、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示する手段と、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成する手段と、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成する手段と、
を備え、
前記樹高線を生成する手段は、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記断面線と前記三角メッシュとが交わる複数の交点を決定し、
前記交点が存在する前記三角メッシュ上の三角形の辺の両端の頂点における高さの差を前記両端の頂点間のｘｙ平面上の距離と、前記両端の頂点のうちの一方と前記交点との間のｘｙ平面上の距離との比に応じて比例配分することにより、前記交点における高さを決定し、
決定された前記複数の交点の高さを用いて、前記樹高線を生成する、
送電線下横断面図生成装置。
【請求項６】
前記地盤高線を生成する手段は、
前記断面線と等高線とが交わる複数の第１の交点のｘ、ｙ座標と、前記断面線と地物とが交わる第２の交点のｘ、ｙ座標を前記三次元地図データから抽出し、
前記第１の交点における高さを前記等高線と等しい値に決定し、
前記第２の交点における高さを、当該第２の交点に隣接する２つの第１の交点間の高さの差を前記２つの第１の交点間のｘｙ平面上の距離と、前記２つの第１の交点のうちの一方と前記第２の交点との間のｘｙ平面上の距離との比に応じて比例配分することにより決定し、
決定された前記第１の交点の高さおよび前記複数の第２の交点の高さを用いて、前記地盤高線を生成する、
請求項４または５記載の送電線下横断面図生成装置。
【請求項７】
三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成プログラムを記憶するコンピュータ読取可能な記憶媒体であって、
前記送電線下横断面図生成プログラムは、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示するステップと、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成するステップと、
を含み、
前記樹高線を生成するステップは、
前記断面線を長さ方向に、指定された内挿間隔で区切った点を注目点とし、各注目点を中心とする検索半径の円からなる検索範囲内を設定するステップと、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記検索範囲内にある複数の頂点を決定するステップと、
前記複数の頂点が有する高さ値をｚｑ１，ｚｑ２，…，ｚｑｎとし、前記注目点から前記検索範囲内の各頂点までのそれぞれの距離に反比例した重みをｗ１，ｗ２，…，ｗｎ（ｗ１＋ｗ２＋…＋ｗｎ＝１）として、前記注目点における高さ方向の値ｚｐを
ｚｐ＝（ｗ１・ｚｑ１＋ｗ２・ｚｑ２＋…＋ｗｎ・ｚｑｎ）で示される内挿計算により決定するステップと、
決定された各注目点の高さを用いて、前記樹高線を生成するステップと、
を含むコンピュータ読取可能な記憶媒体。
【請求項８】
三次元地図データから送電線下の横断面を生成するための送電線下横断面図生成プログラムを記憶したコンピュータ読取可能な記憶媒体であって、
前記送電線下横断面図生成プログラムは、
入力された三次元地図データを用いてｘ−ｙ平面地図を表示するステップと、
前記平面地図の表示画面上で、送電線が張られている径間方向とは垂直な方向に指示された断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
前記断面線に沿った地盤上の樹木の高さの変化を表す樹高線を生成するステップと、
を含み、
前記樹高線を生成するステップは、
前記ｘ―ｙ平面地図上に設定された三角形を相互に接続させた三角メッシュ中の各三角形の頂点について３次元座標（ｘ、ｙ、ｚ）を有する三角メッシュデータから、前記断面線と前記三角メッシュとが交わる複数の交点を決定するステップと、
前記交点が存在する前記三角メッシュ上の三角形の辺の両端の頂点における高さの差を前記両端の頂点間のｘｙ平面上の距離と、前記両端の頂点のうちの一方と前記交点との間のｘｙ平面上の距離との比に応じて比例配分することにより、前記交点における高さを決定するステップと、
決定された前記複数の交点の高さを用いて、前記樹高線を生成するステップと、
を含む、コンピュータ読取可能な記憶媒体。
【請求項９】
前記地盤高線を生成するステップは、
前記断面線と等高線とが交わる複数の第１の交点のｘ、ｙ座標と、前記断面線と地物とが交わる第２の交点のｘ、ｙ座標を前記三次元地図データから抽出するステップと、
前記第１の交点における高さを前記等高線と等しい値に決定するステップと、
前記第２の交点における高さを、当該第２の交点に隣接する２つの第１の交点間の高さの差を前記２つの第１の交点間のｘｙ平面上の距離と、前記２つの第１の交点のうちの一方と前記第２の交点との間のｘｙ平面上の距離との比に応じて比例配分することにより決定するステップと、
決定された前記第１の交点の高さおよび前記複数の第２の交点の高さを用いて、前記断面線に沿った地盤高の変化を表す地盤高線を生成するステップと、
を含む請求項７または８記載のコンピュータ読取可能な記憶媒体。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【公開番号】特開２０１２−１３８０２２（Ｐ２０１２−１３８０２２Ａ）
【公開日】平成２４年７月１９日（２０１２．７．１９）
【国際特許分類】
【出願番号】特願２０１０−２９１１８６（Ｐ２０１０−２９１１８６）
【出願日】平成２２年１２月２７日（２０１０．１２．２７）
【出願人】（５０７４０１５１３）株式会社 ケイソク (3)
【Ｆターム（参考）】
【公開日】平成２４年７月１９日（２０１２．７．１９）
【国際特許分類】
【出願日】平成２２年１２月２７日（２０１０．１２．２７）
【出願人】（５０７４０１５１３）株式会社 ケイソク (3)
【Ｆターム（参考）】
[ Back to top ]