Нахождение параметров преобразования между разными СК

Системы координат, проекции, преобразования, привязка
Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Сообщение Максим Дубинин » 31 янв 2008, 00:15

администрация поддерживает и сама упорно пишет :)
если хотите поучаствовать как авторы - милости прошу

по форме: алгоритмы с описанием + формулы, библиотеки - по желанию.
пристегивайтесь, турбулентность прямо по курсу

Andrey1
Активный участник
Сообщения: 209
Зарегистрирован: 13 июл 2007, 18:58
Репутация: 1

Сообщение Andrey1 » 31 янв 2008, 07:01

EXE в архиве

Аватара пользователя
kitos
Интересующийся
Сообщения: 49
Зарегистрирован: 09 апр 2007, 09:47
Репутация: 0
Откуда: Екатеринбург
Контактная информация:

Сообщение kitos » 01 фев 2008, 14:53

sim писал(а): по форме: алгоритмы с описанием + формулы, библиотеки - по желанию.
я правильно понял, что алгоритм - схема расчета, а не код на С++?
Andrey1 писал(а): EXE в архиве
- самораспаковывающийся архив чтоли?

дошел до формул преобразования прямоугольных координат в географические. на руках у меня было 2 верси ГОСТа 51794-2001: сделаный из скаченых с вашего сайта djvu и в ворде. как показали мои расчеты, ближе к истине первый вариант(в вордовском варианте в формуле 35 вместо sin(2B0) стоит (sin(B0))^2). также в обоих вариантах в 33-й формуле не понятно, к числителю или знаменателю относится последний сомножитель (cosB0). к знаменятелю!(с) если нужно, могу прислать этот вордовские вариант, исправленный. (когда пройду все что мне нужно и исправлю еще ошибки которые найду). жаль, не помню откуда скачал(

Andrey1
Активный участник
Сообщения: 209
Зарегистрирован: 13 июл 2007, 18:58
Репутация: 1

Сообщение Andrey1 » 01 фев 2008, 16:41

Я имею в виду скомпилированные (Delphi) файлы *.exe (т.е. программа-калькулятор), упакованные винраром.

Идёт под любую винду, начиная с Win95, (русскую, со стандартными размерами элементов интерфейса и шрифтов), никаких дополнительных библиотек или софта не требуется. Так как эта программа - прототип, созданный для какого-никакого решения моих собственных задач, то возможны некоторые упущения, ошибки и недоработки.

Оптимальное разрешение экрана - 1024*768 или 1280*1024.

Программ на самом деле две - пересчёт B,L <=> x,y и B1,L1,H1 <=> B2,L2,H2. Для прототипа такая разбивка по функциям сойдёт.

Итак, если интересно и выложишь на сайт - почему бы и нет?

P.S. Секретных алгоритмов, ключей, коэффициентов и т.п. нет - все данные и формулы взяты из открытых и весьма распространённых книг типа "Справочник геодезиста".

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Сообщение Максим Дубинин » 02 фев 2008, 04:29

да, алгоритм - схема расчета, однако, если он сопровождается еще и кодом - еще лучше. Если обратите внимание, обычно я так и делаю - расчет в чем-нибудь общеупотребимом, Excel, например + код, тут уж в чем удобнее.

Andrey1 - калькулятор тоже хорошо, но важно все-таки показать, что именно как и зачем он считает (см. первую часть).

Общий принцип простой - решение должно быть как можно более полным и экономить время нам самим, если мы вдруг решаем вернуться к этой задаче. Это требует хорошей документированности, ссылок на источники и т.д.
пристегивайтесь, турбулентность прямо по курсу

Andrey1
Активный участник
Сообщения: 209
Зарегистрирован: 13 июл 2007, 18:58
Репутация: 1

Сообщение Andrey1 » 02 фев 2008, 11:21

Обычный WinEXE, совместимый с любой виндой, а формулы - из книг Бугаевского, справочника геодезиста и прочих общедоступных источников, также были использованы источники в интернете. Так мне сбрасывать калькулятор тебе на мыло?

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Сообщение Максим Дубинин » 02 фев 2008, 18:31

да
пристегивайтесь, турбулентность прямо по курсу

Аватара пользователя
kitos
Интересующийся
Сообщения: 49
Зарегистрирован: 09 апр 2007, 09:47
Репутация: 0
Откуда: Екатеринбург
Контактная информация:

Сообщение kitos » 04 фев 2008, 11:28

перешел к пересчету СК42 <-> WGS 84.
может кто-то написать точные значения(координаты одной точки в обеих СК)? или параметры(если они не те, что в ГОСТе)?

Аватара пользователя
kitos
Интересующийся
Сообщения: 49
Зарегистрирован: 09 апр 2007, 09:47
Репутация: 0
Откуда: Екатеринбург
Контактная информация:

Сообщение kitos » 04 фев 2008, 16:07

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

Andrey1
Активный участник
Сообщения: 209
Зарегистрирован: 13 июл 2007, 18:58
Репутация: 1

Сообщение Andrey1 » 04 фев 2008, 17:10

Какая-то страшная у тебя задача...

А нахождение параметров преобразования между разными СК для меня тоже интересно (исходные данные для решения такой задачи - координаты некоторого множества точек в обеих СК).

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Сообщение Максим Дубинин » 05 фев 2008, 02:52

про точку: ну вот например +towgs84=23.92,-141.27,-80.9,0,0,0,0
50 50
50.0001958957 49.9984778591
пристегивайтесь, турбулентность прямо по курсу

Аватара пользователя
kitos
Интересующийся
Сообщения: 49
Зарегистрирован: 09 апр 2007, 09:47
Репутация: 0
Откуда: Екатеринбург
Контактная информация:

Сообщение kitos » 05 фев 2008, 09:28

Какая-то страшная у тебя задача...
дааа... :twisted:
50 50
50.0001958957 49.9984778591
я правильно понял, что это трехпараметрическое преобразование из СК-42 в WGS84? и его точность 3м? (мне бы точнее на самом деле :oops: ) в обратную сторону если проверять, эти значения должны сработать?

Andrey1
Активный участник
Сообщения: 209
Зарегистрирован: 13 июл 2007, 18:58
Репутация: 1

Сообщение Andrey1 » 05 фев 2008, 10:02

Sim, наверное, имеет в виду параметры 7-параметрического преобразования?

Если параметры преобразования (3,7) правильные, то его точность порядка 0,1 мм. Но 3-параметрическое преобразование адекватно только в некоторых окрестностях данной точки, размер области зависит от требуемой точности.

Аватара пользователя
kitos
Интересующийся
Сообщения: 49
Зарегистрирован: 09 апр 2007, 09:47
Репутация: 0
Откуда: Екатеринбург
Контактная информация:

Сообщение kitos » 05 фев 2008, 10:11

не понял, что за окрестность?!? откуда она взялась, если параметры - для преобразования СК вцелом? про 3-параметрическое преобразование я предположил потому что в конце строки с пар-ми стоят нули.
еще, если 3 параметра дают точность в 0.1 ММ, то зачем существует и приводится в литературе 7-параметрическое? по идее оно точнее должно быть... куда уж точнее? сдвиги плит замерять?))

Andrey1
Активный участник
Сообщения: 209
Зарегистрирован: 13 июл 2007, 18:58
Репутация: 1

Сообщение Andrey1 » 05 фев 2008, 10:54

Kitos, возьми и посчитай сам - пересчитай координаты точек с помощью 3-параметрического преобразования и с помощью 7-параметрического преобразования. Навскидку мне сложно сказать, на каком расстоянии расхождение станет неприличным, дельта зависит от твоих потребностей. Но, тем не менее, 3-параметрическое преобразование применяется только для относительно небольших районов, обычно метровая ошибка набегает в круге радиусом километров этак сто, может, и немного ошибаюсь.

Считай, разбирайся.

Ответить

Вернуться в «Координаты и привязка»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 0 гостей