1. 課題 　フィールドワークでの測量に問題が生じた。2020年測量のdrone３Dマップと，2018年の光波測距儀を使って得られたネットワークが対応しない。droneマップは光波測距儀で得られたgcpとの対応関係が取られている。そのgcpを使ってdroneマップをコントロールしているので問題はない。この場所には自衛隊のレーダーサイトがあって，コンパスが正常に作動していない危惧があり実際に西偏角度が想定されるものと異なっている。 　光波測距儀の設置点は14カ所で個々の設置点で幾つか測点があって，前視と後視で，位置関係を一応把握したのであるが，手抜きしてGPS計測値を光波測距儀の地球座標値として使用した。これが今抱える問題の元凶だ。そこで，2018年の測量時の測点の現場写真からdrone３Dマップと対応できる測点を利用して，ls_1〜_14の光波測距儀の設置点を求めたいと思う。 2. ls_2の設置点を求める 　ls_2の放射星形区で３Dマップ上での２カ所を現場写真から特定できた。その一つは，図１のls2_20であり，もう一つは図2のls2_11である。後者は光波ネットワークとも一致している。それを図3に示す。 　ls2_20の方は，図２に示すようにかなりズレている。図１の現場写真ではa-rotの下に見えるが，図2では3mほど内陸にずれ込んでいる。 　求めたいのは図１の，赤字で「かわぐち」としたls_2の設置点の３Dマップ上の座標値である。 3. ls_2の3Dマップ上の座標値を求める 3.1 計算法 　当初，3D空間での三角形平面を想定して計算すべきと考えた。しかしながら，熟考するとその必要性が無いことが明らかになった。光波測距儀の構造は台座を水平にして，光波測距儀が装着されている。コンパスで磁北をゼロ度として，台座上を時計回りに回転して，光波測距儀をターゲットに向けて，直線距離を求める。その測距儀の水平角と垂直角を記録する。これだけだ。今，問題なのは３Dマップでのls_2設置点を決定することであり，2D空間に過ぎない。 　連立方程式は簡単にできる。文字表現の簡略化したい。現在の設置点をO’点として，3Dマップに載せうる設置点をO点（Xo, Yo）とする。ls_2放射星形区の２点のうちの１点について，現在の位置をA’点として，3Dマップに載せうる設置点をA点（Xa, Ya）とする。他の１点については，現在の位置をB’点として，3Dマップに載せうる設置点をB点（Xb, Yb）とする。 　光波測量の結果から，次式が成り立つ。赤字表現しているのが未知数である。他の数値はすでに得られた計測値または3Dマップから読み取りできる。|OA|=|OA’|=√((Xa-Xo)^2+(Ya-Yo)^2)　⋯⋯⋯⋯⋯式(1)|OB|=|OB’|=√((Xb-Xo)^2+(Yb-Yo)^2)　⋯⋯⋯⋯⋯式(2)　上記の計算式が成り立つ理由は次のようである。光波測量の結果に基づいて位置ベクトルを３Dマップにプロットするのであるが，ls_2の設置点が決まらないだけであって，その測定結果のうち，水平距離，垂直距離，さらには視準線の間の挟角は変わらない。 　式(1), (2)の連立方程式から，(Xo, Yo)の解が２個現れる。A点を中心として半径|OA|の円と，B点を中心として半径|OB|の円との交点は２点形成されるからである。3Dマップを見ている作者からするといずれが正解かは自明ではある。 　式(1)，(2) は何れも円の方程式である。両円は交わって，その二つの交点を通る直線が成立する。その観点で整理したいと思う。式(2)と(3)はそれぞれ，次のように整理できる。(Xo– Xa)^2+(Yo– Ya)^2 = |OA|^2　⋯⋯⋯⋯⋯式(3)(Xo– Xb)^2+(Yo– Yb)^2 = |OB|^2　⋯⋯⋯⋯⋯式(4) 　二つの交点を通る直線の式は，式(3) – 式(4)，で次のように，得ることができる。2(Xa-Xb)Xo + 2(Ya-Yb)Yo +〔(|OA|^2 – |OB|^2) – (Xa^2 -Xb^2) – (Ya^2 – Yb^2)〕= 0⋯⋯式(5) 　両円の交点の計算は，式(3)と式(4)の何れかの円と式(5)の直線との交点を求めれば良い。ここでは，式(3)の円と式(5)の直線の交点を求めることにする。　式(3)の，中心座標（Xa, Ya），半径 |OA|の円，と，式(5)の直線の交点をM, Nとして，交点を表現したいのであるが，かなり複雑になるので，式(5)を簡潔化して次に表現する。 e=2(Xa-Xb), f=2(Ya-Yb), g=〔(|OA|^2 […]

1. このページ作成の意味 　もともとの関心はGISソフトでの３Dクラウドとgcpの重ね表示の可能性を探ることであった。ここに至る経緯は次のページに示している。 　このページの作成過程で，gisソフトでは実現しないことがわかった。MeshLabに期待したが，これも重い3Dファイルを受け付けない点で問題があった。結局，CloudCompareに頼らざるを得なかった。CloudCompareではMeshLabでどうにもならないものでも何とか作業ができるのである。 　で，問題はレイヤー構造で使用出来るかであるが，Danielの情報で，使えることがわかったのである。ただ，図１のように，Plugins: Cloud layers，は灰色になっている。 　Help > About Plugins，で開いた画面が図２である。Cloud layersを選択するが，これ以上進めない。Referencesとして，サードパーティのWiggins Tech website　が見えるが，ここをクリックするとこの注記が消えてしまう。 2. Plugins: Cloud layersのインストールに向けて 　図３を参照して欲しい。なにか不可解なのだが，CloudCompareのサイトではなくて，3ds-scan.de で公開されている。ぼくがインストールしているのは，図１の v 2.12.4 (Kyiv)なので，2.12.0はより古いバージョンだ。 　キーウで現在活動する企業（サードパーティ）のようだ。 　このサードパーティのリンクをCloudCompareを使用しているWindowsで立ち上げた。ドイツ語が表示される。言語を換える選択肢があり，日本語に替えた。その画面が図３のものである。このサイトは2016年11月から始まっている。図３の「Cloud Layers」プラグインをクリックするとmp4のムーヴィーが始まるがCCもなく，音声も出ない。 　図３の表示のすぐ下には，ダウンロードリンクがある。CloudCompareのミラーサイトの認証を受けているようだが，ヴァージョンが2.12.0のままだ。更新されていない。CloudCompareはフリーだけどサポートが要るというメッセージが色つきになっている。 　次の情報がある。ただし，当方の　Editメニューでは，’Edit > Cloud > Create single point cloud‘しか有効でない。Edit > Scalar fields についても下に示された選択肢はすべて灰色になっているのである。 Here are the new menu entries: 　このWebサイトに入ると，ぼくはすごく怪しい印象を持った。どこにもプラグインのダウンロードのリンクはなく，ただ，コンタクトせよという。しかも種々のページを見ていると，bluetoothに繋げようとする。 　Danielさんが，フォーラムで，Maybe you can use the new ‘qCloudLayers’ plugin?Another way, […]

1. 期待していること 　GrassGISをかつては使ってきたがご無沙汰している。Qgisは１回だけ使って精度面で問題を多く感じてアプリも削除した。Qgisは大衆化が著しく，悪貨が良貨を駆逐する，勢いである。日本の行政はGrassGISを飛び越して，専らQGIS派である。行政のコンテンツがQgisを意識したものになっている。 　さて，3Dメッシュまたは点群ファイルへの対応はどうだろうか。もちろん，GrassGISではこのアプリ公開以来，早くから米軍やNASAのLiDARデータの処理が実行されてきたのであるが，ぼくはこのコンテンツには触れなかった。 　いま，求めるのは，3Dモデルの3D表示能力はもちろんである。CloudCompareはその点では問題はない。しかし，CloudCompareでは，3Dモデルに，他のベクトルデータをオーバーレイすることができない。これはかなり不便なことである。それで，GISソフトに期待してのこのページ作成である。 　繰り返すと，3Dモデル表示が自在で，他のベクトルデータをオーバーレイできるか。この観点で，両アプリを調べたいと思う。 　 2. GrassGISではどうか 2.1 ケース１ 　まずは次の回答を紹介する。 Import Point-Cloud to Grass-GisData Processing – Discussion and Q&A MichaelL Feb ’20 Got it. I know Lidartools (which LAStools is dependent on) was not ported for QGIS 3 so that what I was getting at. 10M points isn’t all that big so if […]