Добавление местной координатной системы в GIS: различия между версиями

Материал из GIS-Lab
Перейти к навигации Перейти к поиску
Нет описания правки
Нет описания правки
Строка 11: Строка 11:
Другие программы, включая QGIS, работают только с «обыкновенными» проекциями, не допуская дополнительных геометрических преобразований. Задача статьи — показать, как сконструировать проекцию, позволяющую работать в МСК в таких средах, как QGIS или, скажем, бортовой софт приёмников GARMIN.
Другие программы, включая QGIS, работают только с «обыкновенными» проекциями, не допуская дополнительных геометрических преобразований. Задача статьи — показать, как сконструировать проекцию, позволяющую работать в МСК в таких средах, как QGIS или, скажем, бортовой софт приёмников GARMIN.


В качестве рабочей среды будем использовать командную строку UNIX. Это идеальный инструмент для экспериментирования, позволяющий непринуждённо сочетать PROJ.4, [[%D0%9A%D0%BE%D0%BD%D1%84%D0%BE%D1%80%D0%BC%D0%BD%D0%BE%D0%B5_%D0%BF%D1%80%D0%B5%D0%BE%D0%B1%D1%80%D0%B0%D0%B7%D0%BE%D0%B2%D0%B0%D0%BD%D0%B8%D0%B5#.D0.9F.D1.80.D0.B8.D0.BC.D0.B5.D1.80_.D0.BF.D1.80.D0.BE.D0.B3.D1.80.D0.B0.D0.BC.D0.BC.D0.BD.D0.BE.D0.B9_.D1.80.D0.B5.D0.B0.D0.BB.D0.B8.D0.B7.D0.B0.D1.86.D0.B8.D0.B8|'''findkey''']] и утилиты GNU обработки текстовых потоков.<ref>Команды воспроизводятся в среде MinGW, что входит в состав QGIS под MS Windows. Две особенности:
В качестве рабочей среды будем использовать командную строку UNIX. Это идеальный инструмент для экспериментирования, позволяющий непринуждённо сочетать PROJ.4, [[%D0%9A%D0%BE%D0%BD%D1%84%D0%BE%D1%80%D0%BC%D0%BD%D0%BE%D0%B5_%D0%BF%D1%80%D0%B5%D0%BE%D0%B1%D1%80%D0%B0%D0%B7%D0%BE%D0%B2%D0%B0%D0%BD%D0%B8%D0%B5#.D0.9F.D1.80.D0.B8.D0.BC.D0.B5.D1.80_.D0.BF.D1.80.D0.BE.D0.B3.D1.80.D0.B0.D0.BC.D0.BC.D0.BD.D0.BE.D0.B9_.D1.80.D0.B5.D0.B0.D0.BB.D0.B8.D0.B7.D0.B0.D1.86.D0.B8.D0.B8|'''findkey''']] и утилиты обработки текстовых потоков '''awk''', '''pr'''.<ref>Команды воспроизводятся в среде MinGW, что входит в состав QGIS под MS Windows. Две особенности:
* нужно заменить команду '''awk''' на '''gawk''', или лучше определить '''awk''' через '''gawk''';
* нужно заменить команду '''awk''' на '''gawk''', или лучше определить '''awk''' через '''gawk''';
* после некоторых команд ('''proj''', '''pr''') придётся добавить в пайп команду удаления лишних символов CR:<pre>| tr -d '\r'</pre>
* после некоторых команд ('''proj''', '''pr''') придётся добавить в пайп команду удаления лишних символов CR:<pre>| tr -d '\r'</pre>
</ref>
</ref>
== Подготовка данных ==
{| class="wikitable"
|-
! ''i'' !! ''X'' !! ''Y'' !! ''x'' !! ''y'' !! ''p''
|- align="right"
| 1 || 7383477.64 || 6087377.60 || 1334.71  || 285.94 || 1.0
|- align="right"
| 2 || 7382557.14 || 6081916.51  || 563.67 || −5197.34 || 1.0
|- align="right"
| 3 || 7386610.19 || 6088160.39 || 4444.27  || 1153.79 || 1.0
|- align="right"
| 4 || 7381962.05 || 6090016.34 || −252.07  || 2881.90 || 1.0
|}
<pre>
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
</pre>


<references />
<references />

Версия от 13:20, 10 марта 2013

Эта страница является черновиком статьи.


Введение

Под местной системой координат (МСК) будет подразумеваться так называемая «городская» система, построенная независимо от государственной системы (ГСК) и включенная в неё заданием ключей перехода к СК-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
  1. ↑ Команды воспроизводятся в среде MinGW, что входит в состав QGIS под MS Windows. Две особенности:
    • нужно заменить команду awk на gawk, или лучше определить awk через gawk;
    • после некоторых команд (proj, pr) придётся добавить в пайп команду удаления лишних символов CR:
      | tr -d '\r'
Источник — https://wiki.gis-lab.info/index.php?title=Добавление_местной_координатной_системы_в_GIS&oldid=10913