５月から公開されるらしい。 http://www.kaerucafe.co.jp/ippenshonin/ 時宗含め特定の宗派は信仰してないけど、先日の善光寺のエントリでも書いたように一遍上人は好きだから、公開されたら地雷原に特攻一応見に行こう。いまんとこ近畿で公開が決まってるの…

1960年没なんで著作権切れてんのね。私が持ってる岩波文庫版の「古寺巡礼」では、一部の古地名に勘違い？誤植？*1と思われる箇所があって、でもずいぶんと版を重ねてる本だしそういう史料があったかもしれず確信がもてなくてもんもんとしてたんだが、青空文…

本題とは外れますが、高安の先で太子道と交差するR25も古代道で、業平道（横大路）といいます。業平道は東西に一直線で名阪国道の天理IC付近までつづいてます。天理IC横には天理スタミナラーメン（天スタ）の本店がありますのでついでにラーメンくってきます…

大和川渡ったあたりで太子道の痕跡がわからなくなります。ここはマジで水没してて道ははなかったんじゃないだろか。とりあえず富雄川沿いの道をR25と交差してるところまでいくと、太子道も終わりです。 太子道の終点付近は高安という集落です。伊勢物語や謡…

太子道の大和側につきあたる付近の建物がない部分は、明治の地図だとびっしりと地図記号が書き込まれているんですが、ほとんどが通常の「田」の地図記号になっています。 「沼田」でもなければ、「草地」や「荒地」でもありません。当然白洲正子さんが歩いた…

黒田を過ぎると三宅町に入ります。万葉集に「打久津 三宅の原ゆ 直土に 足踏み貫き 夏草を 〜」から始まる歌がありますが、足場が悪い夏草の原をずぶずぶ踏みわけて行ったということですから、万葉集の頃は水田や森林じゃなかったことは確かでしょう。このあ…

R24を西に入って田原本駅を過ぎて少し行くと、右手に不安なぐらい細い道が現れます。これが聖徳太子が通ったという太子道です。細い路地の両側には古い民家が建ち並んでいます。 短い距離のところどころに聖徳太子が休んだという石とかがあります。どんだけ…

GPSもって一番下の矢印、保津環濠集落のあたりから太子道に入って、大和川まで抜けてみます。かなーり寂れた細い道、かつ四輪進入禁止の区間などもあるのでもし真似する人は充分注意してください。 今回の GPS ログは Android 4.0端末の 104SH + MyTracks で…

いきなりお外に飛び出す前に、どっからどこまで太子道を辿れそうか地図で確認しておきましょう。 例によって陸軍陸地測量部の二万分の一地形図（明治４１年）と、現在の国土地理院の二万五千分の一地形図の、太子道付近を拡大して並べてみました。デジカメ撮…

標高データとかを地図上でごにょごにょすると、実際にその場にいるだけではわからないことに気付けて楽しいです。でも実際に対象の場所を見ないでデータだけみてるのもつまんねえ*1から、実際の太子道がどうなってんのかちょっとだけ紹介しておきます*2。 *1…

今日の奈良はいい天気でした。こんな日に自宅警備とかしてるのもどうかと思うので、外に出てスポーツGPSログ取りましょう。

昨夏は紀伊半島南部の水害で大きな被害が出ました。今でもR168/R169周辺はちょこちょこ新たな崩落がおきて通せんぼされることもしばしば。仕方ないから諦めてはるか山の彼方の道を通らされる（「広域迂回路」というステキな名前でご案内）ことになるわけです…

こんな見た目になりました。 各自なめるように見てください。マイマップへのリンクはこれでいいのかな？。 お寺や神社、古墳群がきれいに岸辺に並んでて嬉しくなります。太子道も歩けるし。 属性ごとにわけたポリゴンが必要ならgdal_cont… geomval 型の戻り…

立派なアーリーアダプターさんです。PostGIS 2.0が正式リリースされた後ぐらいにこの日記をぐぐって来た人であればたぶん上記のやりかたで問題ないはずですが、この日記を書いてる時点のGDAL/OGR(1.9.0)では上記コマンドがエラーになるはず。これはPostGIS 2…

いままで Google Earth で KMZ をみていましたが、「そもそもGoogle Mapsに重ねろよ」という声が聞こえてきそう。でも Google Maps のマイマップは Raster 入りのKMLは対応してないようでそのままだと表示できないんです＞＜。 Google Mapsと重ねるためには…

奈良湖をずっと真っ黒で描いてましたが、あんまりな見た目なんで、ここは葦原をイメージした緑色＋深いところは青、にしてみましょう。 CREATE TABLE [destテーブル] (rid serial PRIMARY KEY, rast raster); INSERT INTO [destテーブル](rast) SELECT ST_Ad…

今回はお絵かきの話です。

データを奈良湖だけに絞ったところで、最大の水深と平均の水深を測ってみます。ピクセル値を取り出して集計...とかしなくても大丈夫です。 SELECT (stat).mean 平均, (stat).max 最大, (stat).stddev 偏差 FROM (SELECT ST_SummaryStats('[テーブル]','rast'…

奈良と大阪の県境付近、大和川が大阪に流れ出す場所である亀の瀬付近は地盤が弱く、昔からずーっと治水工事してます。奈良湖は亀の瀬のがけ崩れで水がせき止められてできたに違いないと適当に断定して、ここをせき止めてダム化してしまいましょう。 ・・・と…

大きいラスタの操作はけっこう時間かかります。ラスタをインポートするとき、普通に入れると元ファイルの図郭が丸々入りますが、インポート後にいろいろ調べ物をするには必要最小限のラスタタイルだけに絞って扱ったほうが何かと便利です。邪魔な生駒山西側…

初回はデータのお掃除です。

ところで上のSQL例はさも自明なことのようにさらりと書いてますが、誰にでも直観的な書き方とは言えないと思います。 ピクセル操作しようとした人がまずマニュアルを見ると、ST_Value() と ST_SetValue() というのが最初に目にはいってきます。名前からして…

今回はとりあえず、古墳時代〜飛鳥時代あたりの水面をターゲットにして水面を調整しようと思います。 前回の「お試し版」奈良湖でとくに気になったところは、島の山古墳が水没してること。島の山古墳周辺はその名前から考えて、かつて湖に浮かぶ"島"だったよ…

湖水面の高さは時代によって上下します。前回の地図で出た水面がいつの時代相当かはともかく、唐古鍵遺跡が水没*1しちゃってるんで、少なくとも弥生時代には湖の面積はこの範囲より狭かったんちゃうかと。その前提だと、石上〜三輪〜初瀬あたりは相当昔から…

タイトルの分母は飾りです。

で、せっかく作ったデータを見たり見せたりしたいわけですが、一般の人がQGISをインストールしてることなんて期待できません。ここは誰でもお手軽にGoogle Earthで表示できるよう、gdal_translate で KMZ にエクスポートします。 gdal_translate -of KMLSUPE…

インポートしたままだとピクセル値が32bit floatで扱いづらい。そんな精度もいらないので 8bit uint（メートル単位）のテーブル用意します。 CREATE TABLE [8BUIテーブル名] (rid serial PRIMARY KEY, rast raster); 簡単すぎて説明することがないｗ。PostGI…

つづきです。※PostGISでぐぐって来てしまった人向け注記：2.0 beta 基準で書いてますので正式版では手順が変更になるかもしれません。また、最適な処理方法にはなってないような。ツッコミは歓迎。

とりあえずQGISで表示します。 メニューの[プラグイン(P)]->[Pythonプラグインを呼び出す...]->"Load PostGIS Raster to QGIS" を選んでラスタ表示ツールを追加。レイヤの表示スタイルをちょちょいといじれば、 見える…見えるぞ！私にも盆地が見える！まあこ…

PostGIS添付のraster2pgsqlコマンドを使います。たとえばこんな感じで。 $ raster2pgsql -s 4326 -t 250x250 -I -M [入力元ファイル] [出力先テーブル名] | psql -d [DB名] 拍子抜けするほどに簡単！