2017年9月27日水曜日

地盤沈下観測等における衛星活用マニュアル

「地盤沈下観測等における衛星活用マニュアル」 平成29年3月 環境省水・大気環境局

SAR を使った広域観測に関するマニュアルです。実務に取り入れられるのはまだ先かと思っていましたが、すでにマニュアル化されていました。私のような初心者にもちょうど良いレベルになっています。ありがたい。
意外と内容が面白く、一気に読んでしまいました。以下、備忘録です。
p27
ALOS-2 は、非常に高い精度で軌道の調整を行っているため基線長は 500m以下に抑えられており十分に短く、どの組合せでも高い干渉性が得やすいため、ほとんど全 ての観測された衛星データで干渉 SAR解析を行うことができる
p28
ALOS-2 の場合、1ペアの干渉 SAR 解析の2つの 衛星データの間隔は、3 年未満とすることが望ましい
p28
解析に用いる衛星データの選定にあたっては、データ観測時間帯の気象条件、特 に降雨データ(X バンド降雨データ)などを確認し、降雨が確認されれば解析対象デ ータから除外することが望ましい。また、用いる場合は解析精度が悪化しやすいこと に注意する必要がある。
p37
ALOS-2 では、処理レベル 1.1 データ(SLC:Single Look Complex)を入手
p65
ALOS-2 と Sentinel-1 を地盤高の観測に用いた場合の長所・短所

雲の有無に左右されないと聞いていましたので、降雨にも大きく影響されないのではと漠然と考えていましたが、そうではないようです。マニュアルに書かれていますので、やはり経験上よくなかったのでしょう。(先日、データを購入したのですが、確認すると1シーンのみ数時間前に降雨がありました。)

Sentinel-1 に関しても書かれています。
実施例では C バンドにもかかわらず、意外と変位量を把握できているように見えます。沈下は都市域の問題だからでしょうか?
esa のサイトでシーンを検索てみました。が、ALOS2 と同様に2014年以降しか引っ掛かりません。古いデータが安価で手に入ると、過去の変動を多量に処理して把握できるかと思ったのですが。残念。ま、近年のデータがあれば、ALOS2 との比較はできます。山中や広いグラウンドなど、都市域以外でどこまで把握できるか知りたいですね。

あと、InSAR は初歩で、そのスタッキング処理、2.5次元解析までの実施がマニュアル化されています。これ、重要ですね。
試行したことがあるのは 最初のステップのみ。スタッキング処理、2.5次元解析までできて当たり前、といった感覚なのでしょうか?早急にスタッキング処理だけでもクリアーしたいですね。

InSAR では多くの実施例を見かけるのですが、すべて広域です。七五三掛の例でも1km弱。解像度は3mと十分に高いので、色のつながりを見るにしても、もっと狭い範囲の変動を把握できないかと期待しています。
今後の動向としては、3次元解析はもちろん、過去にさかのぼったり、上記のように狭い箇所に着目したりする、などでしょうか?
今後の動向に着目しながら、早く現状に追いつきましょう。


2017年9月24日日曜日

ジオリファレンサーで位置合わせできない QGIS2.14

QGIS で踏査結果をまとめようとして、また躓きました。

踏査で使った地図をラスターで貼り付けたいのですが、できません。ジオリファレンサープラグインでラスターを読み込み、位置を「マップキャンパスより追加」するのですが、まったく違う位置に張り付きます。

違う PC でもダメ。
手法を変えてもダメ。

晩ご飯を食べた後、もしやと思い画像を反転させてからポイントを追加してみると、、、
できました。
移動や拡縮は容易にできるようですが、回転は得意でないのでしょうか?今回、ほぼ180度回さないといけなかったので、うまく機能しなかったようです。北が上になるように画像を調整してからジオリファレンサーで位置合わせ、といった手順になります。今まで使ったソフトに比べてクセがありました。

使っていると、まだまだ出てくるのでしょう。ひとつづつ、クセを覚えていきましょう。

2017年9月23日土曜日

air lift method

簡易索道を張ってもらった若い助手の方にコンプレッサーを借りていました。

それを使って、エアリフト法による孔内洗浄をテストしました。結果は上々。
15mくらいの孔内水位でも汲み出すことができました。

その結果を助手の方に報告していたのですが、横からオペさんが「昔はよくやっていた」と入ってこられました。

2人で「え?!」と驚き。
「エアリフト」という言葉は知っていましたが、その知識と映像、具体的手順が結びついたのがつい最近。しかもUSGSのサイト。「海外では豪快な洗い方をするなあ」と思っていた程度です。
https://www.usgs.gov/media/videos/usgs-groundwater-monitoring-well-redevelopment-using-air-lift-method
https://water.usgs.gov/ogw/video/gwpd.html

オペさん曰く、日本でも「さく井では昔から実施している」とのこと。昔は会社でも井戸掘りの仕事を多く請け負っていたそうで、大型のコンプレッサーを持ち込んで実施されていたそうです。圧をかけすぎるとVPが割れるので注意が必要とか、VP100を入れるべきところをVUを入れてしまい、割って掘りなおした先輩がいるとか、教えていただきました。
この応用で、「エアリフトポンプ」といったものがあるということも御存知でした。井戸掘りを多く請け負っていた支店では複数台所有していたとのこと。もっと早く知りたかったですね。

伝承が途切れていた技術でした。危ないですね。
ひとまず、ノウハウをお持ちの方が身近にいたというこがわかりました。

2017年9月22日金曜日

簡易索道 その2

簡易索道でボーリング資材を搬入していただいています。

ポイントとの通りに良い立木がなく「難しいかな?」と思っていたのですが、オペさん「楽勝です」とのこと。

実際見てみると、三又などをうまく活用して搬入されていました。確かに、難なく作業されていましたが、私には無理。何度も経験してノウハウを身に着けていないと、現場で「楽勝です」なんて言えないと思います。感謝です。

若い方々は見たこともないでしょう。そうなると、計画もできません。災害対応もあるので現場に引き連れてくる時間はありませんが、今回の搬入方法を写真や動画で残しておき、あとで教育に使うことにしましょう。

とにかく、オペさんに助けていただき、感謝、感謝の1日でした。

2017年9月21日木曜日

ReCap Photo

台風18号による被災写真を Recap (photo) に投げ込みました。

できたモデルを表示しようとすると、以下の情報が表示されました。

写真プロジェクトの作成は処理機能が強化された新しい ReCap Photo ™ に移行します。 詳細は、ここをクリックしてください。

3か月くらい前に使用した際には表示されていなかったように思います。

内容を見てみると、いくつかの変更点がありました。気になったのは以下の2点。
・使用できる写真が1000枚に増えた
・ReCap pro のサブスクリプションにバンドルされる
長い時間・コストをかけて実用に耐えうる機能を実装したので、有料にしますよ、ということでしょうか?

従来の Recap (photo) でも、細々した仕様変更がありました。
e-mail で招待し、モデルをシェアできるのは以前からでしょうか?知りませんでした。

一般的な測量でも UAVが 使われだしています。SfM の精度向上を求め、メーカーさんは努力されていることでしょう。
まだまだ仕様変更は続くと思われます。ついていきましょう。

2017年9月18日月曜日

QGIS で 走向傾斜 その2

QGIS で 走向傾斜のシンボルを DXF へ変換した際、歪んでしまう件の続きです。https://phreeqc.blogspot.jp/2017/08/qgis.html

新たな情報を探しているうちに、以下の情報が引っ掛かりました。プリントコンポーザにて SVG や PDF 出力に対応しているとのこと。
https://docs.qgis.org/2.14/ja/docs/user_manual/print_composer/overview_composer.html

そういえば、Civil3D も 2017 から PDF のジオメトリ読み込みに対応していました。
試してみると、DXF 経由ではないので座標やスケールは飛びますが、読み込み自体はOK。歪まずに Civil3D へ取り込むことができました。(Civil 側にて取り込み時のスケール設定はあるのですが、まだよくわかりません。)現段階ではこのあたりが落としどころでしょう。

GIS を使って少し作業を進めると、ラインの表記が楽になったと気づきました。これ、Arcを使われている先輩が、大昔に言われていました。
遷急線や flow line 等のラインが折れていても、等間隔で破線や進行方向のマーク(マーカーライン)を表示してくれます。小崩壊頭部の ┬┬┬ マークも綺麗に並んでくれます。欲を言えば、┬┬┬  のケバが自動で崩積土の範囲まで伸びてくれたら大きな崩壊にも使えるのですが。ま、これはCADと同様、手作業です。

当面は QGIS でルートマップを作成し、Civil3D + GEORAMAで 3次元分布を考慮した地質図を作成することになるでしょうか?
もう少し触ってみましょう。


2017年9月17日日曜日

縦断図 測点間隔の変更(Civil3D 2017)

弾性波探査で、逆解析に遠隔の情報を取り入れたいと考えました。

地形は広域の3次元データ(LP)になっているため、任意のピッチで地盤高を把握することができます。探査測線上は既に Civil3D で線形を設定し、縦断を切っていますので、あとは遠隔まで線形を伸ばし、5mピッチで高さを表示させるだけです。

計画も現地盤もサーフェスから縦断に取り込んでいますので、レポート機能が使えません。今回は以下の手順としました。

①既存の線形を複写し、端部を 50m 延長。
 この時点で縦断図も複写され、延長されます。
②開始測点を-0+50に設定。
 線形を選択し右クリック - 線形プロパティ - 測点 -0+50
 縦断図で複写前の地盤高と見比べる。ずれていなければOK。
③縦断図の測点表示を5mピッチに変更
 縦断図を選択し右クリック-縦断ビュープロパティ
 「情報」タブ - オブジェクトスタイル - 現在の選択をコピー、リネーム
 必要により表示する標高を「標高」タブで変更
 右クリック - 縦断ビュースタイルを編集
 「水平軸」- 副ティックの詳細 - 間隔:5.00m

サーフェスから作成した縦断の測点毎の地盤高を、テキストに吐き出す方法までは思いつきませんでした。ひとまず、これで先に進みましょう。


**************************************************************
おまけ1
 線形に測点ラベルを付ける場合
 線形を右クリック-線形ラベルを編集
 主測点 追加:100m
 副測点 追加:5m

おまけ2
 測点ピッチにかかわらず、縦断図内の地形や計画ラインの折れ点で高さを抜きたい場合は以下の通り(235の技の125番目に詳細が書かれています)。安定計算などではこちらの方が良いでしょう。
 縦断図内の地形線等を選択し右クリック - 「縦断ジオメトリを編集」
 縦断レイアウトツール - 縦断グリッドビュー