問題タブ [proj]

Alberta 10 TM Projection で緯度/経度座標を東/北座標に変換する必要があります。

10 TM 投影法は UTM に似ていますが、カナダのアルバータ州のカスタム投影法です。私は（多少の努力で）自分でコーディングできると思いますが、すでに行われている場合は車輪を再発明したくありません.

「NAD 1983 StatePlane Michigan South FIPS 2113 Feet」(別名ESRI 102690 ) に x、y 座標の大規模なデータセットがあります。それらを緯度経度ポイントに変換したいと思います。

理論的には、これはprojが処理するように構築されたものですが、ドキュメントでは手がかりが得られませんでした。はるかに複雑なケースが説明されているようです。

次のように、 python インターフェイスを使用してみました。

しかし、それは非常に間違った出力をもたらします。

Javascript ライブラリがありますが、処理するポイントが 50,000 まであるため、適切ではないようです。

---

私のために働いたもの：

ptest次のように、x 座標と y 座標をスペースで区切って、各ペアを 1 行に配置したファイルを作成しました。

次に、そのファイルをコマンドに入力し、結果を出力ファイルにパイプしました。

Java のPROJ.4または JavaScript のProj4jsを使用して、SWEREF99 TM の地図投影座標に変換したい WGS84 の GPS 座標があります。

PROJ.4とその使用方法のドキュメントを見つけるのは難しいです。良いリンクがあれば、コメントとして投稿してください。

SWEREF99 TMの PROJ.4 パラメータは次のとおりです。+proj=utm +zone=33 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

変換にPROJ.4 Javaライブラリを使用しようとLat: 55° 00’ N, Long: 12° 45’ Eしましたが、次のコードで試しました:

これは私に要点を与えますPoint2D.Double[5197915.86288144, 1822635.9083898761]が、私はN: 6097106.672, E: 356083.438 何を間違っているのですか? 代わりにどのメソッドとパラメーターを使用する必要がありますか?

正しい値はLantmäterietから取得されます。

proj.transform(testLatLng, new Point2D.Double());使用する正しい方法であるかどうかはわかりません。

kmlに変換するシェープファイルがたくさんあります。

それらを開くことはできますが、それらの座標系がわかりません (作成者から提供されたものでもありません)。

たとえば、これはアムステルダム市の境界上の点である必要があります: 'x': 106980.0, 'y': 491810.0.

理想的には、ウェブサイト、ライブラリ、スクリプト、オラクル、タコなど、これらのデータがどの投影法に基づいているかを教えてくれるものを見つけたいと思っています。

このようなものは存在しますか？私の現在の必要性は、この特定の投影システムを見つけることですが、可能であれば、未知の形式の座標がいくつかあるが、それらがどこにマップされるべきかを知っているすべてのケースでそれを推測する方法を希望します. この質問の読者もそれから恩恵を受けるでしょうが、私の本当に具体的な質問への答えから恩恵を受けることはありません.

投影 SRID 27700 で PostGIS データベースに保存されたポリゴン データがあります。

OpenStreetMap、つまり SRID 900913 で形状を表示したい (と思いますか?)。

だから、2つの質問：

1. 以下のコードを OpenStreetMap の正しい SRID で出力するように変更するにはどうすればよいですか?
2. 以下の Django コードを変更して、ポリゴンとして表示する準備が整った素敵な json オブジェクトを提供するにはどうすればよいですか?

この質問が意味をなさない場合はお詫び申し上げます - 私は GeoDjango を初めて使用します。

これは古き良きものですが、解決策が見つからないようです。

900913 座標で 4326 システムに st_transform を実行したい場合、y 座標がシフトします。

例：

SELECT AsText( Transform( Transform( GeomFromText( 'POINT( 449760.25168159 6790560.4594059 )', 900913), 4326 ), 900913 ) )

ここで、元の 900913 開始点は st_stransformed で 4326 になり、900913 に戻ります。結果は元の点ではなく、y が異なります。(後で結果を挿入します。ここにはありません)。

どこかで読んだように +nadgrids=@null を追加して、4326 の proj4text を変更してみました

現在、srid 4326 の proj4text は次のとおりです。

"srid=4326 の spatial_ref_sys から proj4text を選択" +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs

現在、srid 900913 の proj4text は次のとおりです。

"srid=900913 の spatial_ref_sys から proj4text を選択" +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@ null +no_defs

900913 から 4326 への別の投影への投影も試みましたが、900913 から 4326 への直接変換とまったく同じポイントが得られます。

誰でもアイデアはありますか？

EJ

座標をWGS84からMGAゾーン53に再投影しようとしています。これは、GDA94データムに基づくUTM投影です。結果として無限大になりますが、これは間違いなく正しくありません。私は次のproj4ようにRのパッケージを使用しています：

ここで何が問題になっていますか？

ライブラリ Proj4js を使用して座標を変換しています。座標を EPSG: 23029 から EPSG: 4326 に変換しています。問題は、座標を変換することによって、処理負荷で使用される精度から多少の偏差があることです。これは、座標系を起点から終点へ、またはその逆に変換することでわかりました。変換に使用するコードは次のとおりです。

両方向に変換すると、原点が同じではありません。

誰でもこの問題を解決する方法を知っていますか?

ありがとう。