リンク: Paul Ramsey: PostgreSQL 8.3.5 - Important Upgrade.
リンク: PostgreSQL: Documentation: Manuals: PostgreSQL 8.3: Release 8.3.5.
PostgreSQL8.3を使っている人は8.3.5へアップグレードした方がよいようです。
GiSTインデックスの処理にバグがあり，データを削除したときに関係のないデータを"dead"にしてしまい，検索に引っかからなくなる可能性があるとのこと。関係のないデータが削除されてしまうことはないようです。
PostGISはGiSTインデックスを使用するので，重要な更新になりそうです。
| 固定リンク | コメント (0) | トラックバック (0)
リンク: [postgis-users] hexagonal grid.
PostGISでヘキサゴナル・グリッドを生成する例。メーリングリストに挙がっていたのが上のリンクのやり方です。コピペで試してみるとすんなり作成されます。
何に使うの?...と聞かれても特に使い道は無いのですが，少し感動したのでここにメモ。
| 固定リンク | コメント (2) | トラックバック (0)
リンク: [postgis-users] GeoJson export function.
PostGISでGeoJson出力ができるパッチがメーリングリストに投稿されているので試してみました。
上記リンクの内容では最新のSVN用のようですが，PostgreSQL(8.3.1)+PostGIS(1.3.3)のソースでも特に問題はなさそう。PostGISをconifgureする前に，ソースを展開したディレクトリ(POSTGRES_SRC/contrib/postgis-1.3.3/)でパッチをあて，あとは通常通りconfigure, make, make installするだけです。
とするとSELECT ST_AsGeoJson(ST_GeomFromText('POINT(123.4 234.5)'));
{"type":"Point","coordinates":[123.400000000000006,234.500000000000000]}
と出力されます。
パッチを当てた後のlwpostgis.sql(4390行目付近)では
CREATE OR REPLACE FUNCTION ST_AsGeoJson(geometry)
RETURNS TEXT
AS 'SELECT _ST_AsGeoJson(1, $1, 15, 0)'
LANGUAGE 'SQL' IMMUTABLE STRICT; -- WITH (isstrict,iscachable);
となっていて，小数以下はデフォルトで15桁になります。
以下lwpostgis.sqlで定義されている関数と引数のデフォルト値です。
_ST_AsGeoJson(version, geom, precision, options)
ST_AsGeoJson(geom, precision)
version=1
options=0ST_AsGeoJson(geom)
precision=15
version=1
options=0ST_AsGeoJson(version, geom)
precision=15
options=0ST_AsGeoJson(version, geom, precision)
options=0ST_AsGeoJson(geom, precision, options)
version=1ST_AsGeoJson(version, geom, precision,options)
_ST_AsGeoJsonと同じ
| 固定リンク | コメント (0) | トラックバック (0)
別館のPostGISの項目を整理しつつlwpostgis.sqlを見ていて気付いたのですが，AsUKML()という関数があります。確認したバージョンはPostGIS-1.2.1。
定義は
CREATE OR REPLACE FUNCTION AsUKML(geometry, int4, int4)
RETURNS TEXT
AS '$libdir/liblwgeom.so.1.2','LWGEOM_asKML'
LANGUAGE 'C' IMMUTABLE STRICT; -- WITH (isstrict,iscachable);
...
となっていて，これが自分の理解して(思い込んで)いたAsKML()の定義。
AsKML()の定義はどうなっているかというと
-- AsKML(geom, precision, version)
CREATE OR REPLACE FUNCTION AsKML(geometry, int4, int4)
RETURNS TEXT
AS 'SELECT AsUKML(transform($1,4326),$2,$3)'
LANGUAGE 'SQL' IMMUTABLE STRICT; -- WITH (isstrict,iscachable);
...
となっていて，AsUKML()関数の第一引数にtransform($1,4326)を渡しており，WGS84(SRID=4326)へのtransform()を内部で行うようになっています。
処理する元データがWGS84決めうちであれば，AsUKML()を使う方がよいのかもしれません。
もうひとつ。上の定義には引数が三つありますが，最後の引数にはKMLのバージョンを渡せるようになっています。省略した場合には
...
-- AsUKML(geom, precision) / version=2
CREATE OR REPLACE FUNCTION AsUKML(geometry, int4)
RETURNS TEXT
AS '$libdir/liblwgeom.so.1.2','LWGEOM_asKML'
LANGUAGE 'C' IMMUTABLE STRICT; -- WITH (isstrict,iscachable);
...
となっていてversion 2がデフォルト。型がintなので2.1は渡せませんし，ソースコード(lwgeom_kml.c)のLWGEOM_asKML()関数を見ると
if ( version != 2 )
{
elog(ERROR, "Only KML 2 is supported");
PG_RETURN_NULL();
}
となっていて現状では特に意味はもちませんが，今後のKMLのバージョンアップ(3以降)を想定する場合には考慮しておくのも良いかもしれません。
とりあえず別館のページも以上の内容を加筆しておきました。
| 固定リンク | コメント (0) | トラックバック (0)
リンク: Bubble://ちずろぐ/別館/ - B-Wiki - PostGIS.
とある道路関係のWebGISの機能拡張に関係して，PostGISのLRS(Linear Reference System)関連を調べていたり，以前のエントリにコメントをつけたりしつつ気になったのが，別館のPostGISに関する項目が完全に放置状態なこと。
本館・別館のネタではラスタ処理中心でka-MapやGRASSを主に触っていて，仕事でPostGISを使う場合でも，データがシンプルなもの(とはいって点数は億近いレーザープロファイラデータだったりしますが)ばかりだったので，突っ込んでPostGISをいじり倒すことから遠ざかっていたのですが，ちょうど良い機会ということで，久しぶりに以下のPostGIS関数を追加しました。
全てLRSがらみの関数です。分かり易い文章を書くスキルが無いので図で誤魔化してはいますが...。
元々"PostGIS関数リファレンス"という大層なタイトルをつけていたのですが，どう考えても本家のマニュアルには追いつけないので"PostGIS関数サンプル集"と名前も変えて再スタートになります。
"GRASSリファレンスマニュアル"もタイトルが変わりそうな予感がします...。
| 固定リンク | コメント (0) | トラックバック (1)
リンク: [postgis-users] Fun uses for the asKML() function in PostGIS !!.
PostGIS公式サイトの"Case Studies"にある"UC Davis Soil Resource Laboratory"の，PostGISのasKML()を使用したka-MapとGoogleEarthの連動機能について紹介する内容。
実際のka-Mapのサイトはこちら。ページ最下部の"Give it a try!"から入れます。以下操作の概要を紹介。
↑の画像のページで表示する位置を選択し1:35,000以上のスケールにすると，ka-MapのツールバーにあるGoogleEarthのアイコンがクリックできるようになります。
GoogleEarthアイコンをクリックすると，ka-Map上で表示されている位置のkmlデータをダウンロードし，GoogleEarth上に表示。
GoolgeEarth上のPOINTをクリックすると，地質と勾配の概要が吹き出しで表示されます。
吹き出上のリンクをクリックすると，ブラウザに"Map Unit"の情報が表示されます。
"Map Unit"とは何かという説明はここにあります。難しそうなので読んでませんが...。
地質関係でいつも思うのが，GoogleEarthの地表面の透明度を変更できるようにしてほしいということ。
前にも書いたかもしれませんが，数年前に仕事で三次元GISアプリケーション上に柱状図(簡易的な棒グラフ状の3Dモデル)や地質の想定断面，河川断面等を取り込み，地表の半透明化機能を使って表示させるということをやったのですが，データの蓄積が進むと面的(というか三次元的な)構成が分かり易くなりそうです。高さ方向のスケール拡大機能を併用しないと変化が乏しくて見にくいですが...。 究極はこれですかね（そろそろ削除されるようですが)。
| 固定リンク | コメント (0) | トラックバック (0)
リンク: PostGIS 2.0 Wishlist.
PostGISメーリングリストで少し盛り上がっている。かいつまんでみると
We would like to use SQL/MM as our new design road map. Previously we used OGC Simple Features for SQL, but...省略
Perhaps the most common question on the user mailing list is "Why is distance(A,B) returning 2.34? I know A and B are 874 kilometers apart!"...中略...Far better, to just return the right answer, a metric (or imperial) distance.
It would be nice to have GEOMETRY_COLUMNS as a view, or an automatically updated table.
| 固定リンク | コメント (0) | トラックバック (0)
皆さんご存知の話なのかもしれませんが...。
何の気なしにPostGISのクエリで，述部のgeometry columnsにTransform()をかけて，EXPLAINを実行してみました。
使用したテーブルの概要は
という感じ。
このテーブルに対してTransform()無しで普通に約50m×50の矩形範囲に含まれる点を取り出すクエリプランを見てみると...
-- SQL文
EXPLAIN SELECT the_geom
FROM geo_table
WHERE the_geom &&
SetSRID('BOX3D(...,...)'::box3d, 4612);
-- クエリプラン
Index Scan using geo_table_sidx on geo_table (cost=0.00..30.83 rows=1 width=33)
Index Cond: (the_geom && '...'::geometry)
Filter: (the_geom && '...'::geometry)
という感じで，まぁ普通です。
で，WHERE the_geom の部分をあえてWHERE Transform(the_geom, 4612)にかえてみると，
クエリプランがシーケンシャルに変わり，ものすごいコストをはじき出します。
-- SQL文
EXPLAIN SELECT the_geom
FROM geo_table
WHERE Transform(the_geom, 4612) &&
SetSRID('BOX3D(...,...)'::box3d, 4612);
-- クエリプラン
Seq Scan on geo_table (cost=0.00..2025529.48 rows=425 width=33)
Filter: (transform(the_geom, 4612) && '...'::geometry)
元々POINTはSRID=4612なので，わざわざTransform()する必要は無いわけですが，仮にSRID=2444等の公共座標で
入っていて，かつ結果をSRID=4612で取り出したい場合には，POINTをBOX3DのSRIDに合わせるために安易にthe_geom側に
Transform()をかけると，GiSTインデックスは使用されないということですかね...。
GiSTは元々の座標系(この例では2444)で作成されてはずなので，極自然な結果ではあるわけですが。
取り込まれているデータの座標系と，取り出したい座標系が異なる場合は，BOX3Dをthe_geomに合わせ，取得結果にTransform()をかける)という方法が最適な方法なのでしょうね。
-- こんな感じ?
-- BOX3DをはじめからSRID=2444にしても良いわけですが...
SELECT Transform(the_geom, 4612)
FROM geo_table WHERE the_geom &&
Transform(SetSRID('BOX3D(...,...)'::box3d, 4612), 2444);
少なくともクエリプランナはそこまで気を利かせることはしないようです。
まぁ実行した環境の実メモリ量(PostgreSQLの共有メモリ量)等でもプランが変わりそうなので一概には言えないのでしょうけど。
| 固定リンク | コメント (0) | トラックバック (1)
試しただけですが。
-- 点データをKMLで出力する例
--
select AsKML(the_geom, 4) from points limit 3;
kml
------------------------------------------------------------
<Point><coordinates>130.5,33.24,4.19</coordinates></Point>
<Point><coordinates>130.5,33.24,3.71</coordinates></Point>
<Point><coordinates>130.5,33.24,4.61</coordinates></Point>
(3 rows)
| 固定リンク | コメント (7) | トラックバック (1)
AsKML()から始まったPostGIS1.2.1の話。結局今日PostGIS1.2.1のリリースがアナウンスされました。
Bringing all the bug fixes of the last month to you, is the 1.2.1
release!
<http://postgis.refractions.net/download>
The 1.2.1 release of PostGIS is now available. This is a minor bug
fix release.* Fix for bug in Within() caused by point-in-polygon performance
shortcut.
* Fix for bug in indexes with null on PostgreSQL 8.2.
* Fix for JTS handling of multi-dimensional data.
* Fix for GCJ Java support.
* Fix for JDBC compatibility in PostgreSQL 8.2.
* New AsKML() function.
* Better Transform() performance when no transform is actually required!AsKML and Transform courtesy of Eduin Carrillo. JDBC work courtesy
of Markus Schaber, per usual, with an assist from Thomas Marti. Mark
Cave-Ayland tracked down the 8.2 null problem. Mark Leslie cleaned
up the p-i-p problem in within().
ということで，メインはバグフィックスですが，AsKML()の追加とTransform()の改善がされているようです。
| 固定リンク | コメント (3) | トラックバック (0)
昨日PostGIS1.2.0のマニュアルにAsKML()がのっているのだが...という話を書いたのだが...。
その後スレッドは別の方向に流れて，Javaで作った「SimpleKMLExporter」があるよ，という話に。これ，WKTをKMLに変換するというもの。
Have a look at
http://www.zubi.li/wiki/doku.php?id=prog:java
(text german)I written a SimpleKMLExporter
in Java that generate a KML file from a
WKT string.
The javadoc is unfortunately written
in my bad english.
;-)
ソースファイルのサイトはここ。ドイツ語ですか...。
Thanks, Thomas. I'll have a look.
Luckily,I can understand your bad
German as well ;-)
彼らの漫才はとりあえずほっといて，このソースで使われている以下のクラス
は，ここから入手可能。
Javaは何年も触ったことがないし，コンパイルして試す気も無いし，と思いつつソースをざっと読んでみたのだが内部でWGS84へ座標変換しているCH1903Converterクラスって，スイス連邦国土地理局の測地系(CH1903)ですか...。スイスと縁のない私には使えませんね。
| 固定リンク | コメント (0) | トラックバック (0)
さっきメーリングリストに以下のような投稿が...
I just saw that an asKML() function was added to the on-line postgis
manual (1.2.0).
In the release notes, there is no trace to when this was added.
Apparently not in 1.1.5. Does it require an upgrade to 1.2.0 (preferably
soft, but how?).
で調べてみたのだが，確かにVer1.2.0のマニュアルにはAsKML(geometry, [precision])の記述が載っているが，psqlインターフェース上で試してみてもそんな関数は無いと怒られる。
ついでに nm liblwgeom.so.1.2 してもシンボル名は無し。ソースやlwpostgis.sqlをgrepしても見つからず。1.2.1で追加されるのかな?
-- 追記 --
念のため，普段購読していないpostgis-develの12月のアーカイブを覗いてみると，色々やり取りが行われていたみたい。ざっと目を通しただけだが...(以下雰囲気だけの超手抜きダイジェスト)
のような議論があっているようだ。
とりあえず，開発版のソースには取り込まれてるっぽい。
| 固定リンク | コメント (0) | トラックバック (0)
PGdijkstraとはPostGISの最短経路探索用pl/sql関数のこと。
dijkstra(ダイクストラ)法とは，最短経路探索の代表的なアルゴリズム。
Chris氏のページにPGdijkstraを使用した処理速度に関するレポートが乗っている。
私がPGdijkstraを知ったのは，同氏のブログの一連の記事から。
最近はなかなかGIS関連に割り振る時間が無いので，少々GIS界隈の進歩に乗り遅れ気味だが，同記事によると
SELECT astext(the_geom) FROM shortest_path_as_geometry(‘roadstable’, startnode, endnode)
といった感じでSQLを発行すると，探索してくれるようだ。
興味のある方は御一読を。
PGdijkstraに関する日本語のページは
http://www.saruga-tondara.net/GIS/?pgdijkstra
(株)オークニーもPostLBSという独自?に開発PGdijkstraをベースとしたPostGIS関連製品を販売している。
http://www.orkney.co.jp/product/postlbs.html
| 固定リンク | コメント (2) | トラックバック (0)
一つ前のブログも含めて，FDOのコンパイルで気になったのが，AWKTという文字がよく出てくること。
どうもAutoDesk拡張のWKT形式らしい。知らなかった。
grepで色々調べてみると，ソースファイル(テスト用?)内に
GeomFromText('POINT XYZ (10 11 12)')
GeomFromText('POINT XYM (10 11 1.2)')
GeomFromText('POINT XYZM (10 11 12 1.2)')
というような記述がある。
POINTと(座標の並び)の間に，XYZ，XYM，XYZMという記述があるのが拡張だろうか?。これってPostGISのEWKTとも違う。
確かEWKTではPOINTMを導入して構文自体の変更はしていないはず。
オリジナルのWKTにEWKTにAWKTか...ややこしい...。
| 固定リンク | コメント (0) | トラックバック (0)
| 日 | 月 | 火 | 水 | 木 | 金 | 土 |
|---|---|---|---|---|---|---|
| 1 | 2 | 3 | 4 | 5 | 6 | |
| 7 | 8 | 9 | 10 | 11 | 12 | 13 |
| 14 | 15 | 16 | 17 | 18 | 19 | 20 |
| 21 | 22 | 23 | 24 | 25 | 26 | 27 |
| 28 |
Markus Neteler, Helena Mitasova: Open Source Gis: A Grass Gis Approach
GRASSの殆んど全てをこの本で勉強した。自分が持っているのは1st版。これは3rd版。英語版です。
Schuyler Erle: Mapping Hacks: Tips & Tools for Electronic Cartography (HACKS)
地図システム屋さんなら，ぱらぱらめくるだけで楽しくなる本。英語版です。
最近のコメント