Добавление местной координатной системы в GIS
Введение
Под местной системой координат (МСК) будет подразумеваться так называемая «городская» система, построенная независимо от государственной системы (ГСК) и включенная в неё заданием ключей перехода к СК-42 или СК-63. МСК крупных территорий, сравнимых с размерами субъектов Федерации, не являются предметом данной статьи, поскольку относятся к классическим картографическим проекциям.
Постановка задачи
Некоторые программы позволяют реализовать работу в МСК непосредственно. Так, в MapInfo любая проекция может быть дополнена аффинным преобразованием. ArcGIS в качестве МСК предлагает локальную проекцию: аналог ортометрической проекции на эллипсоиде, дополненный разворотом и масштабированием.
Другие программы, включая QGIS, работают только с «обыкновенными» проекциями, не допуская дополнительных геометрических преобразований. Задача статьи — показать, как сконструировать проекцию, позволяющую работать в МСК в таких средах, как QGIS или, скажем, бортовой софт приёмников GARMIN.
В качестве рабочей среды будем использовать командную строку UNIX. Это идеальный инструмент для экспериментирования, позволяющий непринуждённо сочетать PROJ.4, findkey и утилиты обработки текстовых потоков awk, pr.[1]
Подготовка данных
| i | X | Y | x | y | p |
|---|---|---|---|---|---|
| 1 | 7383477.64 | 6087377.60 | 1334.71 | 285.94 | 1.0 |
| 2 | 7382557.14 | 6081916.51 | 563.67 | −5197.34 | 1.0 |
| 3 | 7386610.19 | 6088160.39 | 4444.27 | 1153.79 | 1.0 |
| 4 | 7381962.05 | 6090016.34 | −252.07 | 2881.90 | 1.0 |
1 7383477.64 6087377.60 1334.71 285.94 1.0 2 7382557.14 6081916.51 563.67 -5197.34 1.0 3 7386610.19 6088160.39 4444.27 1153.79 1.0 4 7381962.05 6090016.34 -252.07 2881.90 1.0
- ↑ Команды воспроизводятся в среде MinGW, что входит в состав QGIS под MS Windows. Две особенности:
- нужно заменить команду awk на gawk, или лучше определить awk через gawk;
- после некоторых команд (proj, pr) придётся добавить в пайп команду удаления лишних символов CR:
| tr -d '\r'