光波測量結果とグーグルアース画像とを平面座標系で補正する試み aim for suitability correction btw electro-optical distance measurements and Google Earth’s images using the coordinate plane
1. Google Earthとのつきあい方 光波測距儀設置点をGEPro情報から捉える:後方交会法 では,後方交会法の限界を思い知った結果になった。ドローン測量の際に,公共基準点を捉えられなかったことが今なお,響いている。現地に再度行った方が,よほど人生を生かすことになるはずだが,だらだらと不完全なドローン測量結果を生かすべく,闘っている? ドローン測量と併せて光波測量を実施したが,光波測距儀の設置位置が数メートル規模で確定できない。gpsはこの精度では信頼できない。それではと,Google Earthを使おうとのめり込んだ。この限界については,グーグルアースプロの限界 にぼくが得た情報をまとめている。その結果,Google Earth上で何らかの作業をするのは,狭い範囲の測量情報を取り扱うのは,時間をドブに捨てるのと同じだと,遅ればせながら,わかった。GEが提供するプログラムがあるが,これは児童のお遊びツールとしてのみ,何とか価値があるかも知れない。 2. 平面上での作業を選択 現在の悩みは,前述のように,もともとのドローン測量の周辺的な問題もあるが,Google Earthの画像と光波測距儀のgcpのズレを気にしてしまったことにある。Google Earthの画像が当方の使用目的には不適当であることがわかっても,ドローン測量から得られたortho.tifをGoogleアースに落とした際にほぼGoogle Earthの画像と一致していて,違和感を感じさせないということに意味があると,グーグルアースの画像がモザイク写真であっても,思っている。 このページ作成の前に,前述の後方交会法のページで最後に掲載したGoogle Earthの画像から後方交会法で得たプッシュピンを削除し画像を図1に示す。赤いプッシュピン(電柱byKなどの個別名称)はGoogle Earth上で,現地調査時の現場写真やドローン測量で得られた空中写真に基づいて,手作業で落としたものである。これをここでは仮にGoogle Earth系A〜Fとする。 茶色のプッシュピン(a-GE〜g-GE,-GEを付加しているのは光波測距儀の設置位置を現場写真なろからGE上に落としたことに由来する)は光波測距儀で得られたものである。これをここでは仮に光波系a〜fとする。 今後の作業過程を次に。 ステップ1 光波系a〜f(計7点)とbase-GE(GE画像にプロットした光波測距儀の想定設置位置)(以上,計8点。平面直角座標系1の形)。Google Earth系A〜F(計7点)についてはkml出力して,経緯度(十進角)を平面直角座標系1に一括変換する。 Wikipedia由来か。次の情報を得て,今後は,十進角で表現したい。「十進角(じゅうしんかく、英:Decimal degrees (DD) )とは、角度の表現法の一つである。緯度と経度の地理座標を度分秒法が北緯a度b分c.d秒、東経e度f分g.h秒と表現するのに対して、十進角ではa.bcd度、経度e.fgh度といったように分秒の部分を小数、十進法で表す。」 ステップ2 この二つの系を一つの散布図にプロットして画像出力して,それをイラストレーターに配置する。この一連の作業の目的はGoogle Earthのモザイク写真との整合性を取ることなので,光波系a〜f(計7点)のみ,トレースしてベクトル化して,手作業?で回転および平行移動をする。 なお,この回転と平行移動は現実的でなければならない。平行移動は最大2mほどだろうし,回転は数度以内であろう。光波測距儀では回転軸はコンパスで実行されるが,大山の自衛隊基地のレーダーの影響が考えられるので,数度の回転を要する可能性はある。 ステップ3 ステップ2の結果に基づいて,幾何学的な数値計算を実施する。 ステップ4 その結果に基づいて,光波測量の結果を回転および平行移動を実施する。 ステップ5 ステップ4で得られた結果を十進角に変換してKMLファイルを作成してGoogle Earthに取り込む。 ステップ6 Google Earth上で,ノータッチのGoogle Earth系A〜F(7点)と変換後の光波系a〜fとベース(8点)の整合性を一応,確認する。この作業を通じて,当然,みかけ状の改善は見られるはずである。 3. 計算作業 3.1 ステップ1 Google Earthから出力したKMLをみると,図2のように, 経緯度値が2セット見られる。<LookAt></LookAt>内に,128.5300162650217 27.35264129535775<Point></Point>内に,128.5300076128268,27.35263507066106KML リファレンス には,<LookAt></LookAt>について,次の説明がある。 Feature から派生する要素に関連付ける仮想カメラを定義します。LookAt 要素は、表示対象のオブジェクトに対する「カメラ」の位置を設定します。Google Earth では、ユーザーが [場所] パネルのアイテムや 3D ビューア内のアイコンをダブルクリックすると、「空を飛ぶように」視点が移動し、LookAt で指定した視点のビューに切り替わります。 というわけで,建物南西角の場所は,<Point></Point>内の経緯度を採用する必要がある。 […]
