Получение координат прямоугольной сетки Гаусса-Крюгера в градусах на эллипсоиде WGS84#

Для извлечения координат узлов прямоугольной сетки в проекции Гаусса-Крюгера и преобразование их в градусы на эллипсоиде WGS84 необходимо выполнить ряд шагов:

  1. Подготовить сетку точек в узлах прямоугольной сетки Гаусса-Крюгера соответствующей зоны.

  2. Перепроецировать точки в систему координат WGS84 с указанием параметров перехода между эллипсоидами.

  3. Извлечь из координат точек значения в таблицу атрибутов в колонки X и Y.

Подготовка сетки точек#

Для подготовки сетки точек следует использовать ГИС NextGIS QGIS. В NextGIS QGIS устанавливаем систему координат карты (иконка глобуса в левом нижнем углу) и выставляем EPSG:28410 (система координат Гаусса-Крюгера 6-ти градусная зона 10). Далее, в панели «Инструменты анализа» в разделе «Вектор - Создание» выбираем Регулярные точки или вызываем меню Вектор ‣ Анализ ‣ Регулярные точки. В открывшемся диалоге (см. Рис. 125.) вводим следующие данные:

../../_images/regular_points_ru.png

Рис. 125. Диалог «Регулярные точки»#

В поле «Исходный охват» вводим координаты углов прямоугольника, в которые необходимо вписать сетку, и систему координат:

10436000.0000, 10452000.0000, 5664000.0000, 5682000.0000 [EPSG:28410]

Шаг сетки выставляем в 1000 м. Указываем путь и имя файла для сохранения. Далее нажимаем кнопку Выполнить. После формирования файла он будет добавлен в виде слоя в окно карты.

Для проверки корректности формирования точек добавим онлайн подложку. Для этого нужно открыть Интернет ‣ QuickMapServices ‣ Поиск в QMS. Найти OSM Mapnik и нажать Добавить.

Далее необходимо выставить СК карты EPSG:3857 Google merkator (см. howto_select_srs_pic).

../../_images/select_srs_ru.png

Рис. 126. Диалог «Установка системы координат»#

Так как эллипсоиды у подложки и слоя точек отличаются, должен отобразиться диалог выбора параметров перехода, в котором нужно выбрать первый пункт (см. Рис. 127.).

../../_images/parameters_srs_ru.png

Рис. 127. Диалог «Выбор параметров перехода»#

Пример получившейся карты представлен на Рис. 128..

../../_images/sample_map_view_ru.png

Рис. 128. Фрагмент окна карты NextGIS QGIS#

Перепроецирование точек в систему координат WGS84#

Для перепроецирования точек в систему координат WGS84 с указанием параметров перехода необходимо в дереве слоев вызвать контекстное меню и выбрать Сохранить как…. В открывшемся диалоге (см. Рис. 129.) необходимо выбрать формат Shape file ESRI и указать систему координат EPSG:4326 (WGS84).

../../_images/save_as_ru.png

Рис. 129. Диалог сохранения слоя#

Так как мы указали параметры перехода между эллипсоидами, то дополнительно указывать параметры перехода не требуется. В случае, если диалог выбора параметров будет показан, то необходимо выбрать первый пункт, как показано на Рис. 127..

Извлечение координат точек в колонки X и Y таблицы атрибутов#

Для извлечения координат точек в таблицу атрибутов слоя необходимо открыть таблицу атрибутов. В таблице атрибутов активировать калькулятор поля (ctrl+i) - см. Рис. 130..

../../_images/field_calc_add_coord_ru.png

Рис. 130. Диалог «Калькулятор поля»#

Для добавления координаты Х в калькуляторе поля выбрать:

  1. Создать новое поле.

  2. Имя поля - Х.

  3. Тип поля - десятичное целое (real).

  4. Размер - 16, точность - 8 (если этот параметр недоступен, пересохраните слой в формате ESRI Shapefile).

  5. Выражение - $x .

Далее нажимаем ОК.

Для добавления координаты Y в калькуляторе поля выбрать:

  1. Создать новое поле.

  2. Имя поля - Y.

  3. Тип поля - десятичное число (real).

  4. Размер - 16, точность - 8.

  5. Выражение - $y .

Далее нажимаем ОК.