Using datum (TM rectangle to WGS-84)

Using datum in transformations

Now, let's complicate our case and transform rectangle TM coordinates to geodetic coordinates at ellipsoid WGS-84. As a work file we continue using file at last case (SPb_Demo.dwg). Теперь усложним задачу и переведем плоские прямоугольные координаты в проекции Поперечная Меркатора в геодезические координаты на эллипсоиде WGS-84. При этом в качестве рабочего файла продолжим использовать пример из предыдущего кейса (SPb_Demo.dwg).

Script for getting geodetic coordinates on other ellipsod (WGS-84)

As test point we retry using our previous point (116479.50982048,99975.21929290,16.0) В качестве тестовой точки повторим использование нашей прежней точки (116479.50982048,99975.21929290,16.0).

Cheme of script

Group 0 - input parameters/Начальные параметры

Here we using 3 coordinates (X,Y,Z), because we change elipsoids and Height of point also changings (when we will use datum's calculations). In that script a also using codeblock for initialize values of coordinates:

Здесь используются уже 3 координаты (X,Y,Z), так как при смене эллипсоидов высота должна быть пересчитана (когда мы перейдем к вычислениям с датумом). В этом скрипте значения координат также представлены в формате codeblock:

Input coordintes of point

Group 1 - Info about CS/Информация о системе координат

As opposed to previous case we use only one CS (source CS), because needing result we get as geodedic coordinates at ellipsoid WGS-84.

В отличие от предыдущего случая, тут мы используем только одну СК (исходную СК), поскольку нам требуется получить результат - геодезические координаты на эллипсоиде WGS-84.

Input source CS

Group 2 - Info about ellipsoids/Информация по эллипсоидам

As opposed to previous case we need select source ellipsoid and target ellipsoid (WGS-84). In ellipsoid parameters are missing - look at page above.

В отличие от предыдущего случая, здесь необходимо выбрать исходный эллипсоид и целевой эллипсоид (WGS-84). Если параметры эллипсоидов отсутствуют - см. предыдущую страницу.

Choosing ellipsoid parameters

Group 3 - Info about datum/Информация о датуме

Here we need select one of datum (number's parameters to go from one ellipsoid to another). These parameters also there are at EPSG web-site or in government laws. In our case we need use EPSG:1267. After that we will take that string as input data to node GeneralTerms.GetDatumInfo. If needing datum is missing - loot at pink nodes and write parameters as string in order below:

"ΔX,ΔY,ΔZ,ω_x_rad,ω_y_rad,ω_z_rad,mE-6";

Здесь необходимо выбрать один из датумов (числовые параметры трансформации для перехода от одного эллипсоида к другому). Эти параметры также могут быть найдены на веб-сайте EPSG или опубликованы в государственных нормативных документах. Для нашего случая, это датум EPSG:1267. После этого данную строку необходимо подать на вход ноду GeneralTerms.GetDatumInfo.Если требуемый датум отсутствует, см. на розовые ноды и запиши параметры в строку согласно порядку выше.

Group 4 - Node TM_FromRectangleToGeodetic (was decribed early)

Group 5 - Node Datum_Recalculation

View of node Datum_Recalculations

That node takes as input parameters info about datum (Group 3), about source and tarhet ellipsoids (Group 2) and coordinates of point (results of Group 4 ang Height from input parameters).

Данный нод беретв качестве входных значений информацию по датуму (Шаг 3), информацию по исходному и целевому эллипсоиду (Шаг 2), и координаты точки (пересчитанные широту/долготу в Шаге 4) и высоту из начальных параметров.

Group 6 - Results and comparing accuracy/Результаты и сравнение точности

Results of script

As result we get a recalculated coordinates and Height of point. Node CoordTransform.StrFormatOfGraduses returm more viewing result. В качестве результата мы получаем пересчитанные координаты точки и ее высоту. Для преобразования результата в более удобный формат можно использовать нод CoordTransform.StrFormatOfGraduses.

We can check accuracy in comparing results witj standard system function _geomarkpoint. As we can see, accuracy of calculation = system's results.

Мы можем проверить точность в сравнении результатов со стандартной системной функцией - маркером геоположения _geomarkpoint.Как мы можем видеть на скриншоте ниже - вычисленные значения равны системным.

Comparing results and system info

Last updated

Was this helpful?