Описание системы координат

Question

Существует местная плановая система координат, описанная и зарегистрированная в ArcGIS (то есть все параметры данной СК известны). Существует производная от нее СК, полученная методом XY-сдвига и поворота на N градусов относительно некоей базовой точки, не являющейся началом системы координат. Необходимо описать производную СК в ArcGIS, то есть получить ее параметры.

Возможно кто-то сталкивался с подобной задачей? Заранее благодарен за ответы.

Answer 1

Итак, если я все правильно понял, то процедура преобразования координат происходит в две итерации. Вначале мы сдвигаем СК, а потом поворачиваем вокруг некой точки. Соответственно для решения мы должны знать сдвиг, координаты центра поворота и угол поворота.

 

Если ничего не напутал, то координаты в новой СК будут равны:

 

X_k=(X_i+ΔX)COSφ-(Y_i+ΔY)SINφ+X_{0<?:NAMESPACE PREFIX = O />}

Y_k=(X_i+ΔX)SINφ+(Y_i+ΔY)COSφ+Y₀  , где

 

X_i, Y_i – координаты в старой системе координат;

ΔX, ΔY – сдвиг координат, значения могут быть как положительными, так и отрицательными;

X₀, Y₀ – координаты центра вращения (поворота) в старой системе координат;

φ – угол поворота;

<?:namespace prefix = st1 ns = "urn:schemas-microsoft-com:office:smarttags" /><?:NAMESPACE PREFIX = ST1 />X_k, Y_k – координаты в новой системе координат

Answer 2

Спасибо, математические формулы пересчета координат мне известны. Но это не выход. Необходимо именно определить параметры проекции. Дабы описав ее однажды уже не возвращаться к проблеме пересчета. То есть "на лету" работать с данными и в исходной и в производной СК. Насколько я понимаю это будет некое преобразование датума. Вот только какое именно - с этим и пытаюсь разобраться.

Answer 3

Необходимо именно определить параметры проекции. 

Но про проекцию не было ни слова.

Преобразование датума это и есть преобразование СК. То что вы спрашивали делается вышеописанными преобразованиями.

Answer 4

Хорошо. А где указанные коэффициенты могут быть заданы при создании собственной СК? В каких полях диалога свойств СК? Какой метод трансформации должен быть при этом задан?

Answer 5

аналогичная проблема )
есть ли варианты решения?
СК на основе проекции Гаусса-Крюгера + смещение по двум координатам и поворот относительно некоторой точки
необходимо описание проекции (*.prj), с помощью которого возможно проецирование "на лету" и совмещение данных в местной системе координат с данными в проекции Г-К

Answer 6

параметры сдвига задать несложно, можно даже напрямую в файле prj. Вот как выглядит сдвиг на юг на нужную величину. Проверено, работает.

"PROJCS["Pulkovo_moved",GEOGCS["GCS_Pulkovo_1942",DATUM["D_Pulkovo_1942",SPHEROID["Krasovsky_1940",6378245.0,298.3]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Gauss_Kruger"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",4785000.0],PARAMETER["Central_Meridian",39.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]"

 

Поворот как параметр отсутствует в большинстве проекций, поэтому так запросто как сдвиг, вам его задать не удастся. Попробуйте для тренировки укротить какую-нибудь косую проекцию (Conic Oblique, например). Однако я не уверен, что есть аналог для Transverse Mercator - а вас именно он небось интересует.

Answer 7

со сдвигом проблем нет ) False_Easting, False_Northing можно и напрямую в файле, и в свойствах проекции в ArcGis...
меня именно поворот и интересует.. я думал, может кто нашел проекцию для Меркатора (Вы правы, именно она мне и нужна)
Может кто знает, что за проекция Local? В ней в качестве параметра присутствует Azimuth. Может как-нибудь с извращениями можно подогнать ее под поворот?
ЗЫ: Извините, если где неправильно что-нибудь сказал. Образование у меня не географическое, а техническое, и в проекциях я пока разбираюсь слабо, но вот, жизнь заставила )

Answer 8

Про проекциии стоит прочитать. Главная особенность - это модель поверхности, на которую идет отображение с геоида. Это разумеется, упрощенное представление (помните школу, наверно - цилиндр, конус, плоскость). Эта модель и задается теми самыми формулами. "Поверхность" жестко вшита в понятие проекции, и параметры этого конуса вам регулировать не дадут. Разве что его положение относительно геоида, т.е. линии касания, сечения и... в некоторых случаях можно отклонить поверхность от стандартного положения на некоторый угол. Следовательно, задача - подобрать именно цилиндрическую проекцию, в которой можно менять угол или как-то еще отклонять ее. Так что "ваша" проекция - Mercator Oblique. Про нее например, тут рассуждают: http://www.navgeocom.ru/forum/index.php?showtopic=67 .

Answer 9

В окне, где задаются проекции, есть кнопочка – Transformations… Там есть и повороты. Только кто бы объяснил, как этим пользоваться.

Answer 10

Вручную подбирать надо. Изменил prj - обновил вид.
На самом деле недолго.