Вычисление левого верхнего угла растра

Question

Здравствуйте, помогите пожалуйста ! Имеется рисунок, размером 500x375 известны координаты центра в десятичных градусах и размер пикселя в метрах.
Нужно вычислить:
•     координаты левого верхнего угла;
•     размер пикселя в десятичных градусах.

Answer 1

Тебе быстро или точно? :)

Если точно, жди, кто-нить из геодезистов проснется отпишет. Было не так давно, в праздники лень вспоминать, поверь. А быстро 1 градус=111 км по меридиану, уже одна цифра есть :) Долготу считаем из длины дуги - простая тригонометрия. Это если быстро. Нннооо...

Точность зависит от размеров шарика, какой взять - то есть по сути от выбора проекций. Так что надо сразу определяться, КАКИЕ у тебя градусы, красовские или WGS например - разница будет метров 100. Мое IMHO: лучше в таких задачах схитрять и использовать встроенный механизм проекций - всетки мы в ArcMap работаем, 21-й век на дворе. Например:

- выбираем географическую проекцию, строим точку-центр по широте-долготе.

- Меняем проекцию на метровую, например ГК. Привязываем растр углом по центру (а какая разница?), задавая пиксель в метрах. Если надо, сдвигаем растр на нужное расстояние, целясь центром на точку-центр

- перепроецируем обратно в нужную географию, замеряем размер пикселя в градусах, если еще не отпало :)

 

Answer 2

Точность, прежде всего. Проекция GCS_WGS_1984. Использование ArcMap метод хорош, но для большого количества рисунков занимает много времени, поэтому меня и интересует именно вычисление.

Answer 3

Координаты левого верхнего угла (в метрах) можно вычислить только если есть проекция...

А так как координаты центра в десятичных градусах и размер пикселя в метрах, то ерунда получается... и каким образом удалось получить размер пикселя ещё тот вопрос

я так понял необходимо создать файл привязки растра для ArcGIS там как раз и есть размер пикселя, левый угол и поворот...

ширина= 500*размер пикселя (по горизонтали, если они различаются)

высота= 375*размер пикселя (по вертикали)

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

Answer 4

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

Answer 5

Характеристики не проблема, описания эллипсоидов есть в составе описания проекции. Но тут вы попадаете в область достаточно сложной математики, где придется на все вопросы отвечать самому, например:

 - как ориентированы полуоси эллипсоида и можно ли их использовать для вычисления длины окружности;

- как вычислить длину полуокружности эллипса и годится ли она для вычисления длины дуги;

- надо ли и как учесть датум при переходе от российских топокарт к мериканским;

Так что лучше, наверно, тогда начать не с форума, а с учебника по геодезии, там это все есть. Вряд ли стоит уповать на его упрощенное изложение тут в теме. А так вобще задача благородная, и наверняка решается "скупыми средствами".

 

Answer 6

Вот нашел кое-что "как это делается". Есть скрипт ESRI для пересчетов, оказывается, можно использовать как пример.

http://giscraft.ru/programs/calcarea.shtml