Proj.4 4.6.1 RC1とGoogleメルカトール

　Proj.4 4.6.1 RC1が出たので，以前書いたEPSG:3785(funny epsgish code for the google mercator)に対応しているか調べてみたのですが，epsgファイルには追加されていないようです。ただしesri.extraファイルにはID=900913としてパラメータが追加(4.6.0から?)されているのでprojコマンド等を使用する場合には

    proj +init=esri.extra:900913

とすれば使えます。

　esri.extraの900913のパラメータをepsgファイルにコピーしてIDを900913から3785に変えれば当面は問題なさそうです。

　ところで900913の設定には

    +nadgrids=@null

という見慣れないパラメータがありますが，FAQ(日本語訳)によると球体を使うこと(楕円体ではないこと)による間違ったシフト計算を回避するために，シフトファイル名にnullを渡す裏技のようです。