wgs​​84とnad83間の変換

これは、wgs84とnad83の間の変換に関する初心者の質問です。航空写真があり、緯度経度座標をできる限り厳密に計算しようとしています。疑わしいパラメータなど、通常のタイプの修正をトリッキーにするいくつかの要因があります。そのため、試行錯誤が伴います。

したがって、修正を行うときに進捗状況を確認する方法が必要です。最初はgoogle earthをチェックするだけで問題ありませんでしたが、この時点でデータ（nad83）とgoogle earth（wgs84）の2.5フィート程度の差が制限要因になりました。したがって、nad83をwgs84に変換する必要があります。

Rでrgdalを使用し、C＃でdotspatialを使用して投影を行っていますが、nad83とwgs84を切り替えることは「投影」とは見なされません（そうですか？）。「+ proj = NAD83 ellps = WGS84」でプロジェクトを呼び出そうとすると、エラーがスローされます。

私は自分で7パラメータ座標の変更を書いて大丈夫ですが、いずれかを（X、Y、Z）座標に変換する方法がわかりません（これは、彼らが基づいている楕円の表面に緯度経度から表面）。nad83とwgs84に関する既存の質問は、私のような非常に基本的な質問に答えるのではなく、どの方法を使用するかに関するものです。つまり、どのようにそれを行うのでしょうか。

通常、NAD83とWGS84は互いに1メートル以内にあります。メートル未満の2.5フィートの違いに関する懸念は、このデータ変換を実行する必要があることを示しています。

簡単に言うと、この計算には、座標の経時的な動きを説明できるように座標が収集された時期の知識が必要です（主に、基盤となる構造プレートの動きと、「地震に伴う突発的な動き」による）。最初に、生データを1997.0や2002.0などの特定の標準日付を参照する同等の座標に変換します。次に、（b）時間に関するこれらのパラメーターの一次導関数とともに、特定の日付のデータムを変更するために必要な（a）オフセット、回転、再スケーリング（つまり、7パラメーター変換）を説明する14個のパラメーターを検索します。（b）に経過時間を掛けて（a）に加算すると、必要な7パラメーター変換の1次近似が得られ、それを適用します。

手順は数学的には難しくありませんが、データが最初に収集および処理された方法の詳細にかなり関係しています。NOAA論文、Transformations Between NAD83とWGS84は、明確で読みやすい説明を提供します。NOAA Webサイト（http://www.ngs.noaa.gov/TOOLS/Htdp/Htdp.shtml）で入手可能な対話型の「水平時間依存位置」ソフトウェアHTDPを使用することをお勧めします。Fortranのソースコードは、主にプレート構造モデルの係数といくつかのインターフェイスコードで構成されていますが、実際に計算を行うのは約3,000行です。移植するのは大きな苦痛かもしれません。

Google Earthの航空写真の精度が2.5フィートオフセットに下がった場合、私は確信が持てません。それは、正確に直交する写真ではないという問題かもしれません。

必要なのは、写真を写真に合わせるためのアフィン変換だけです。