UTM座標上での半径検索続き

前の記事にレスポンスをいただきました。

2点間距離(1点からの距離)による検索方法3：MugeSoの日記

これはあるゾーンでの第２象限における等距離円とUTMで得た外接矩形(赤)、経緯度座標での外接矩形(青)、expand_by_meterで得られるポリゴンのバウンディングボックスの関係を示したものです。このように、青で示した外接矩形よりも大きくはなりますが、ポリゴンのバウンディングボックスは等距離円を内包しています。

本当だ...スゲエ。
私のほうが頭の中だけで考えて、メルカトル脳的に結論だしちゃってました。
やっぱりちゃんと図に描かないと判らないですね。
感動、感服いたしました。

中央経線を跨ぐ場合にはバウンディングボックスの北側の辺が等距離円に食い込みますが

これはUTM上の緯線の、もっとも凸なところを跨いじゃうから、ということですね。
で、解決法が、

これは前のエントリの方法で回避できます。

＞解決方法としてはUTM上で中央経線と正方形の交点を求め、元の座標系でこの点も含めたバウンディングボックスを得ればＯＫなはずです。

と。
なるほど、その通りじゃないかなと思います。

ただ、そこまでして経緯度座標系上でのバウンディングボックスを求める理由が、

私がなぜこのような方法に固執するのかというと、検索対象レコードが等距離円を持っていて条件として与えた点をその等距離園内に含むすべてのレコードを得たいからです。等距離円は中心を示すgeometry型列と半径を示すdouble型列を用いて表現しています。このとき、中心点の列で空間インデックスを張っても半径を得るために逐次検索となってしまうのです。

考えてみるとUTM上で外接矩形を用いたインデックスを張る手もありますが、国内で動かす場合6つのゾーンすべてについてインデックスを張らなければいけないのが気になります。検索対象が頻繁に移動することがあるため避けたいところなのです。

だとすれば、UTMではなくメルカトル図法を用いる方法を検討されてはどうでしょうか。

元ネタはここ↓

Geographic/Photogrammetric/Cartographic - 1点からの半径検索は、Google 投影に変換してから検索するのがベストプラクティス？

第1回: Google 投影にふれる
第2回: Mercator 投影の縮尺係数

メルカトル図法なので、検索中心の緯度が変わるごとに検索範囲に補正係数をかけてやる必要がありますが、関数インデックスは単一で済みます。

1エントリ：1アイヌリンク

CHIKAR STUDIO

アサンカラ（旭川）アイヌの血を引く、カラフト・アイヌの伝統弦楽器「トンコリ」の奏者、OKIさんのページ。
DUB AINU BANDを率いて活躍中とのこと。

初めて聞きましたけど、幻想的な感じですごくいいです。
Fuji Rock Festとかにも出てるんですね。
恥ずかしながら知りませんでした。

• DUB AINU BANDの楽天ダウンロード検索結果