緯度経度から都道府県コードを得る

PostGIS

id:yellow_73:20060406 で ESRI Japan の市区町村界データをPostGISに叩き込みましたが、これを使って、ポイントから都道府県コードを得る方法。

SELECT jcode/1000 AS pcode
FROM city
WHERE WITHIN('POINT(<経度> <緯度>)'::geometry, the_geom);" );

これだけ。
意外とさくさく動いてくれるので、あまりアクセスが無いなら十分使えます。