ГОСТ Р 51794-2008: Преобразование геодезических координат в плоские прямоугольные и обратно

Восточнее... аж почти на две зоны...
Александр Юрьевич, наверное, на автопилоте отбросил от игреков первые цифирьки...

 

Что называется, оговорка по Фрейду. Правильно читать - восточнее.

 

Не-е, если бы отбросил, то расстояние не посчитал бы.
Перевел в геодезические, и увидел, что восточнее. - 109°26'14" и 100°32'43"

 

Epic fail
Достану козырь: если в таблице есть еще такие данные, как номер датума (я так понимаю) и смещение:

Код:

project_name map_type Xmin Xmax Ymin Ymax resolution GK Смещение
burjat_1 clutter 18001236 18998736 5494759 6382959 100 18 18500000
burjat_2 clutter 20067479 20535879 5708988 6406188 100 20 20500000
chita_1 clutter 19988298 20927998 5399026 6548426 100 20 20500000
irkutsk_1 clutter 17998477 19000677 5612595 7204995 100 18 18500000
irkutsk_2 clutter 17228563 17538163 5844787 6460887 100 17 17500000
irkutsk_3 clutter 20090971 20674371 6226568 6773368 100 20 20500000
piter clutter 6149400 6698300 6446000 6863300 100 6 6500000

Это что-нибудь изменит?
Я думал они не нужны, достаточно координат.

 

О-о-о!

Это в корне меняет дело! (раскрыть) Это в корне меняет дело! (свернуть)

Однажды у Генри Киссинджера спросили:
— Что такое "челночная дипломатия"?
Киссинджер ответил:
— О! Это универсальный еврейский метод! Поясню на примере:
Вы хотите методом челночной дипломатии выдать дочь Рокфеллера
замуж за простого парня из русской деревни.
— Каким образом?
— Очень просто.

Я еду в русскую деревню, нахожу там простого парня и спрашиваю:
— Хочешь жениться на американской еврейке?
Он мне:
— Нахрена?! У нас и своих девчонок полно.
Я ему:
— Да. Но она — дочка миллиардера!
Он:
— О! Это в корне меняет дело!

Тогда я еду в Швейцарию, на заседание правления банка и спрашиваю:
— Вы хотите иметь президентом сибирского мужика?
— Фу, — говорят мне в банке.
— А если он, при этом, будет зятем Рокфеллера?
— О! Это в корне меняет дело!

И таки–да, я еду домой к Рокфеллеру и спрашиваю:
— Хотите иметь зятем простого русского мужика?
Он мне:
— Что вы такое говорите, у нас в семье все — финансисты!
Я ему:
— А он, как раз, — президент правления Швейцарского банка!
Он:
— О! Это в корне меняет дело!
Сюзи! Пойди сюда.
Мистер Киссинджер нашел тебе
жениха. Это президент Швейцарского банка!

Сюзи:
— Фи… Все эти финансисты — дохляки или педики!
А я ей:
— Да! Но этот — здоровенный сибирский мужик!
Она:
— О–о–о! Это в корне меняет дело!

(правда, это не совсем то, что некоторые, пытаясь выглядеть профи, называют новомодным словом "датум".)
Сейчас попробую подоходчивее сформулировать, как считать.

Оффтоп

Проверил свою программку, тестовые данные брал с Test data for the transverse Mercator projection

Результатом полностью удовлетворён.

 

По мин/макс видно, что ваши карты поболее чем одна зона (первые цифры игреков, по недоразумению названные тут иксами).
GK (это номер не датума, а зоны в стандартной 6-градусной разграфки проекции Гаусса-Крюгера) дублируется в первых цифрах смещения - именно в эти зоны приведены те ваши карты, которые выходят за пределы одной зоны (piter никуда приводить не надо - он и так в одной зоне). Если эти карты необходимо каким-то образом совместить со стандартными листами, то либо стандартные трансформировать, либо ваши резать и тоже трансформировать надо. Резать не по прямоугольным осям, а по меридианам/параллелям - на "трапеции". Соответственно, и анализировать что куда попадает следует не по прямоугольным координатам, а по геодезическим (долгота/широта). Впрочем, я уже повторяюсь...

 

К сути: как определить, попадает ли точка в квадрат, заданный геодезическими координатами?
И не проще ли тогда будет всё-таки переводить координаты точки в плоские, но с учетом номера зоны карты?

 

Не проще. Хотя бы потому, что границы карты не прямоугольник - присмотритесь к листу карты-километровки...
Переведите ваши прямоугольные точки ;-)) в широту-долготу (программа у вас уже есть) и сравнивайте с пределами широтно-долготного "квадрата" простым сложением-вычитанием...

 

Код:

function XYtoBL(X, Y, GK, Shift) {
var z0 = (Y - Shift) / (6378245*Math.cos(B0));
…
var res = new Coordinates((B0 + dB)*180/Math.PI, 6*GK-3 + l*180/Math.PI);
}
 
function BLtoXY(B, L, GK, Shift) {
…
var n = GK;
…
var y = Shift +
    l*Math.cos(B)*…
}

Вроде бы не ошибся.
И, да, 57.29577951 = 180/Math.PI
А по сути ответить сложно. Насколько мне известно, никто ещё не нарезал карты в прямоугольных координатах.
Поэтому и возникают вопросы по происхождению этих кластеров.
Кроме того, та же теорема Пифагора справедлива (до определённой степени) только для точек в одной зоне.
В проекции Гаусса-Крюгера, в отличии от привычной многим Декартовой системы координат, масштаб не постоянен, а увеличивается с удалением от осевого меридиана. И чем дальше от осевого меридиана, тем быстрее изменяется масштаб.
Поэтому и стремятся не работать в зонах шире стандартных 6°. У вас же для для первого района ширина зоны аж 14°.

 

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

 

Оффтоп

ИТЦ "СканЭкс"... ;-))

 

Что вы понимаете под нарезал? К электронным картам это имеет отношение? Почему бы их и не сделать в прямоугольных?

 

Оффтоп

ага...Andrey_NM, Вы имеете отношение к упомянутой компании?... ну объясните же пожалуйста на популярном языке что такое "тайл" применительно к векторной карте и/или зачем вектор растрировать, теряя объекты, и нарезать на квадратики разных масштабов... был объект, например здание 5КН с адресом, стало прямоугольное пятно на растре...

 

а не все ли равно в каких координатах хранить векторную карту? и вообще, зачем их нарезать?...

 

Не имею. Соответственно объяснить ничего не могу :-)

 

нет.
вот широта-долгота (типа проекция в мапинфо):

вот "Пулково-42" (Гаусса-Крюгера):


На первой, разумеется большие искажения будут, если там контуры материков нарисовать... Непривычно, непонятно... Но нам же не смотреть, а измерять надо, верно?... Причем в вашей задачке "куда попадает" именно по карте. И без ее изображения, а лишь по цифрам...
И что проще?....
Собсно, в мапинфо в любой проекции есть пространственные операторы "внутри", "содержит" - там вообще ваша задача решается одним запросом...
Но вы же программу свою пишите... Не... Делайте как хотите, время ваше... Изобретайте свой велосипед, создавайте себе трудности...

 

кстати, на вашей картинке линии 2-3 и 4-1 должны быть дугами...

 

Так вы оппоненты, а я за помощью пришел :)

А тут отбиваться приходится...
Уже на троллинг становится похоже.

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