Формулы перевода координат на изображении

ДД!
Нужно решить задачу расчета на растровом рисунке карты положения точки, которая соответствует определенной геогр.координате.
Проекция - коническая равнопромежуточная, все параметры есть. А фот формул x,y = F (L,B) нету :(

помогите, пож, сделать такие формулы.
Можно и на небольшой возмездной основе! :)

Заранее спасибо!

 

спасибо за ссылки на книги, но я сейчас не справлюсь с изучением всего учебника :( нашел раздел про этот тип проекции - яснее ничего не стало.
я уже пробовал реализовать одну мат.модель проекции, но результат не попадает в координатную сетку. Видимо, эта модель подразумевает параллели строго круговые, а на карте это скорее дуги

 

Вам для каких целей это надо?

 

программно нарисовать на растровой карте точки в географически верных местах.

 

Озик пробовали?

 

нет, а что это?

 

http://www.oziexplorer.com/ Только название темы уже пора поменять...

 

"программно" - самому запрограммировать.
мне нужен алгоритм, чтобы засунуть его в свою программу.

 

мож, проекция другая?
выложить бы картиночку-то...
вообще, неплохо бы сформулировать задачу с самого начала - что за карты, что за точки... мож, какой другой путь подскажет коллективный разум?... ;-)

 

задача простая - есть куча геогр.координат в таблице - нужно автоматизировано эти точки поставить на растровой карте.
Карта сделана в профессиональном ГИС, т.ч. подозрений в ее криворукости нет. Зато есть параметры проекции - см.ниже.
Из них мне пока непонятна цифра 0,297 в описании сфероида, и что цифра после UNIT DEGREE
Как можно понять параметры кривизны линий параллелей этой карты?

Europe_Equidistant_Conic
SPHEROID["International_1924",6378388.0,297.0]]
UNIT["Degree",0.0174532925199433]]
PROJECTION["Equidistant_Conic"],
PARAMETER["False_Easting",0.0],
PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",10.0],
PARAMETER["Standard_Parallel_1",43.0],
PARAMETER["Standard_Parallel_2",62.0],
PARAMETER["Latitude_Of_Origin",30.0],
UNIT["Meter",1.0]]

 

Это не 0,297 - это 297.0 (сжатие)
Эллипсоид Хейфорда (1910 г.)

 

этот параметр эллипсоида мне сможет помочь понять кривизну линий параллелей на проецкии?

 

Вряд ли. Это отношение большой и малой полуосей. На таком эллипсоиде возможна проекция с прямыми параллелями. Скорее всего, дело в конической эквидистантной проекции. Эквидистантная - эквивалентные расстояния. При этом направления сильно искажаются. Прямая линия будет дуговой. Это предположнение (не работал с коническими проекциями) Выложите сюда хотя бы скриншот фрагмента вашей карты, чтобы не гадать.

непонятно... может это угол разворота в градусах?.. Не масштаб же...
Может, поможет:

Раскрыть Спойлер Свернуть Спойлер

x= xp + rsinb y =yc + r cosb
r = ( qmin - q )/(90 - q min) yp + xpctga b = a (j - jc )/(j -jmin),x = -x'sinj c + y'cosj c + xc
y = -x' sinjc - y' sinjc sinqc + z' cosqc + yc
x' = cosq cosj CA /s
y' = cosq sinj CAB /s
z' = sinq CB /s
гдеj cиqc, - долгота и широта центральной точки карты, xc и yc - смещение по горизонтали и по вертикали центральной точки карты относительно начала координат s - масштаб, CA, CBи CAB - соответственно минимальный, максимальный и средний радиус земного эллипсоида. Геодезические координаты центральной точки карты определяют единственную точку (j c,qc), в которой нет искажений проекции, и при удалении от которой искажения возрастают.

Нашел на http://pc601s.vigg.ru/atlas/B_3_1_2.htm

 

я кажется понял свою ошибку
неверно предполагал, что в конической проекции полюс (90 широта) - это центральная точка для кругов параллелей. на самом же деле 90я широта тоже круговая, а центр.точка немного выше.

 

LegacyGT, наткнулся на такую штуку - http://trac.osgeo.org/proj/ Возможно, там будут формулы преобразования широта-долгота (градусы) в север-восток (метры). Программка PROJ4 может быть там, ее можно включить как библиотеку или просто нарыть формулы в исходном коде.

 

Нет, конечно.
π/180 = 0.01745329251994329576923690768489.

Библиотека PROJ.4 - довольно старая, проверенная и очень уважаемая в мире ГИС. С сырцами.
Если нужны формулы, я бы посоветовал Coordinate Conversions and Transformations including Formulas

 

А-а-а-а... Теперь понятно, объявление юнита угловых единиц, градус, выраженный в радианах. А то голову ломал.

 

Ура!
с помощью "калькулятора, транспортира и такой-то матери", т.е. экселя, визио и ряда грубых слов, мне удалось подогнать формулы и параметры так, что погрешность на самых дальних меридианах карты составляет несколько пЫкселей.

Пришлось схалтурить и измерить реальный угол наклона меридианов на рисунке, т.к. угол наклона линий меридианов не равен самому меридиану.
Но почему-то в параметрах проекции нет никаких данных о том, какой угол наклона относительно базового меридиана должен соответствовать одному градусу долготы. По факту оказалось, 10гр долготы ~~7,88 градусов на карте.