同僚、こんにちは！
解決すべき問題-座標を通常のカートグラフィック度（WGS84）に変換するための100％信頼できるパラメーターを取得します。同僚は私が何を意味するのかをすでに理解していますが、好奇心旺盛な人に説明します。たとえば、USRNの抽出からGPS受信機の座標への座標を再計算するためのパラメータを使用してインターネット上を歩くアプリケーションとアルゴリズムは、圧倒的に「地図上の野郎」になるという事実です。家の大きさのオブジェクトを検索する場合、これは問題にはなりませんが、エンジニアリングベンチャーの場合、精度はすでに弱いです。たとえば、地形の境界をセンチメートルの精度でマークしたり、地下のパイプやケーブルを見つけたり、フラット座標のマップを使用してドローンを起動したり、インターネットのオンラインカートグラフィックベースを使用してフラット座標で図面を描画したりするなど、サブメートルの精度が必要です。
正確な座標が不正確になる理由
私たちの不動産に関する情報やプロジェクトや図面から私たちによく知られているフラットメーターの座標は、ローカルでは非常に正確ですが、数学だけでは地球に結び付けるのに十分ではありません。事実、文書からの「フラット、メーター」座標系の数学モデルは、当時の技術の精度で、最初に測地点の形で地上に実装されました（ロシア連邦では、ほとんどの座標系はソビエト時代に開発され、運用されています今日まで）。そして、その時だけ、これらのファーストクラスの地理的ポイントから、敵を混乱させるために、それらすべてから、座標グリッドの回転と歪みを伴うSK63などの秘密の座標系を導き出した他のものが作成されました。そのような変換のたびに、歪みは許容され、重要ではありませんが、変換の数に関して直線的に増加することはありません。しかし、はるかに進歩的です。その結果、ほとんどの座標グリッドは、わずかにしわくちゃになって片方の端から引っ張られたシートのように見えます。そのため、99％の地理計算機では、座標グリッドの「しわくちゃのシート」からあなたを救うことはできません。有料の座標を再計算するためのいくつかの測地サービスがあります。そこにいる人々は、座標系の理論的パラメータに従って数えられませんが、「しわくちゃのシート」のすべてのパラメータを持っていると推測できます。ほとんどのロシア連邦では、小さな領域の座標系のパラメータを計算する必要があります。これらの領域の半径は、多くの場合15kmを超えません。このような小さな領域では、座標グリッドの歪みが1センチメートルを超えないことが多く、座標系は非常に正確に地球上に「存在」します。関心が20〜30 kmのスペースを超えて忍び寄る場合は、より小さな領域の移行のいくつかのローカルパラメータに依存する必要があります。座標系をより小さなサブゾーンに分割します。
自転車を再発明しますか？
座標系の真のパラメータを計算するために、それはソフトウェアの「車」であることがわかります。確かに、これらは通常、ジオデシーとカートグラフィ用の複雑で高価なプロフェッショナルソフトウェアシステムのコンポーネントであり、すべてのプロフェッショナルが必要とするわけではありません。一部のソフトウェアは、理解できない、使いにくいファイル形式でWindows上で実行されます。より多くのコマンドラインオプションがあります。プロファイルフォーラムの記事に触発されて、ブラックジャックとWebを使用して独自のソフトウェアを作成することにしました。
それがどのように機能するかを簡単に説明します
パラメータを計算するには、WGSと同じローカル座標系の2つの座標系で、パラメータを知りたい3つのポイント（できるだけ多く、より正確になります）が必要です。カタログ、テリトリーのカダストラルプラン、リファレンスカードなどからフラットメーター座標を取得します。度座標または同じソースから、またはGNSS基地局の座標を使用することが多いか、衛星受信機を使用して、平面座標で知られているこれらの点を単純に測定します。
Webフォームは、座標系のパラメーターを計算し、それをGISシステムの99％に適用できる2つの一般的な形式（proj stringとWKT）で表示します。
ここに、まさにパラメータと少しの用語についてのいくつかの物語があります
多くの理解できない手紙
, . x,y,z lat long h. h , , , ( ).
— «». , . , , .
() — , . , , ..
, . PROJ4 (MapInfp, ArcGIS . . , ): +proj=omerc +lat_0=59.8338730825 +lonc=33 +alpha=-0.0001 +gamma=-1.771957267229058 +k=0.9996584453038837 +x_0=2365031.423134961 +y_0=426397.2888527482 +ellps=krass
(+ellps=krass)— . : , , . . . , 42 .
(+proj=omerc) — . « », , .
. , , , . . . , 42 . (omerc).
(+lat_0=59.8338730825 +lonc=33) — , « » ( ), . . .
(+alpha=-0.0001 +gamma=-1.771957267229058) .
(+k=0.9996584453038837), . , . , .
(+x_0=2365031.423134961 +y_0=426397.2888527482), .
— «». , . , , .
() — , . , , ..
, . PROJ4 (MapInfp, ArcGIS . . , ): +proj=omerc +lat_0=59.8338730825 +lonc=33 +alpha=-0.0001 +gamma=-1.771957267229058 +k=0.9996584453038837 +x_0=2365031.423134961 +y_0=426397.2888527482 +ellps=krass
(+ellps=krass)— . : , , . . . , 42 .
(+proj=omerc) — . « », , .
. , , , . . . , 42 . (omerc).
(+lat_0=59.8338730825 +lonc=33) — , « » ( ), . . .
(+alpha=-0.0001 +gamma=-1.771957267229058) .
(+k=0.9996584453038837), . , . , .
(+x_0=2365031.423134961 +y_0=426397.2888527482), .
オープンソースパッケージに基づく-
- 地理的変換のためのproj4
- 地図上に情報を表示するためのリーフレット
- センチメートルの精度でパラメーターの効果の領域を計算するためのgeophp（この記事の執筆時点では、実装されていません）
Webフォームのソースコードは、オープンリポジトリのAGPLライセンスの下で入手できます。
tutuのウェブフォームの議論。