создание карт для Magellan eXplorist 310, 510, 610, 710,

Дима, ну ты сам себе ответил. Одно из двух - или резать на полоски, или мириться с некоторым ухудшением картинки вверху и внизу.
В принципе, масштаб "3,2 км" больше подходит для обзорной, чем для основной карты, по которой ты непосредственно ориентируешься. А километровки, пятисотки и более подробные карты вряд ли ты будешь делать такими огромными кусками, чтобы между верхом и низом была значительная разница в разрешении.


Ну как сказать. Бывает так, как ты пишешь, а бывает и по-другому. Новые Гармины (Орегоны 600-650, Етрексы 20х-30х и 25-35 Тач, Монтерра) поддерживают те же форматы карт (IMG, KMZ, JNX), что и их предшественники. А бесплатного вектора для Гарминов достаточно много, и вполне можно подобрать себе то, что "влезет на флешку". Так что тут как раз ухудшения от появления новых моделей я не вижу.

Бывают и другие ситуации, когда обновления приносят не совсем те "покращення", на которые рассчитывали. Тот же Ньюман несколько раз менял формат карт для мобильного Озика (ozf2, два вида ozfx3, ozf4), и поздние Озики могут не понимать карты, которые делались для самых ранних Озиков. Или возьмём прошивки для Эксплориста 310. Постепенно прибили Экспорт, Импорт и Усреднение позиции. Из-за усреднения я сижу на 1.19. После майских попробую поставить 1.63 - вдруг вернули. Если нет - опять откачусь на 1.19.

В чём-то ты прав. Я вот на Lowrance Endura облизываюсь. Он мне нравится даже чисто внешне. Да и само сочетание сенсорного экрана и удобных кнопок выглядит очень привлекательно. Но к тому моменту, когда я накоплю на него деньги, он исчезнет из продажи, WinCE вымрет вместе со всеми своими прогами, и останутся только смартфоны с кучей своих приложений и огромными сенсорными экранами, но без кнопочек и без возможности использовать пальчиковые батарейки.

_________________
Карти для GPS-навігаторів: де їх брати? https://encyclopedia-mandriv.blogspot.com/2013/06/karty-dlya-gps-navigatorov.html

Подивився мапу, визначив що з z14 (по SAS). Взяв і переробив, але
1. Границі тайлів SAS мічені
2. Трьома кусками, z12 (1контейнер) z13 (3 контейнера), z14(9).

Бо якщо робити все з z14 то дуже великий діапазон карти в м/пікс: На 85 -1,65м/пікс а на 60- 9,55 м/пікс
3. dв САС поставив більшу якість JPEG
Гляньте, може сподобається. По >>Жаль, не знаю, какое там разрешение. називаєте номер тайлу (він то привязаний до рівня z і широти) і я вам видаю розділову здатність, потім оформлю в екселівський файл.

На вході рівень SAS Z та нумерація тайлу X, Y, на виході все що видно
скрін з екселя та перевірки розміру в SAS додано, координати вернього лівого пікселя тайлу та дозвільну здатнісь у відповідних координатах перевіряйте самі по SAS

В такому розрізі справа лише за підрахунком пікселів на скріні.

Ну что ж... поделюсь результатами своих изысканий.

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

1. Методика измерений.

Если у нас есть в навигаторе путевые точки, то переносим в безопасное место файл \Waypoints\Magellan_New_Waypoints.gpx - чтобы удобнее было работать он должен быть пустым.
Отключаем все карты.
Устанавливаем ориентацию карты "всегда на север".
Выключаем GPS, а текущее положение выставляем вручную на нужной нам широте (хоть вводом координат, хоть маркером на карте).

Дальше устанавливаем желаемый масштабный уровень и ставим четыре точки вблизи углов карты. Делаем скриншот:



Повторяем процедуру для нескольких масштабных уровней.

Сбрасываем скриншоты и файл \Waypoints\Magellan_New_Waypoints.gpx на комп для последующей обработки.

Переходим на другую широту и повторяем все снова.

2. Обработка

В итоге вышеописанной процедуры мы получили набор скриншотов и набор координат точек на них.



Теперь мы можем для каждого скриншота сопоставить координаты точки (меряем точку у основания древка флажка) в пикселах и в градусах. Для каждого скриншота можно получить 4 разности по широте (в нашем случае это точки 1-2, 1-3, 2-4 и 3-4) и 4 по долготе (1-3, 1-4, 2-3 и 2-4) и посчитать по 4 значения разрешений по широте и долготе, которые затем усреднить.

Разрешения считаются в угловых секундах на пиксель.

Таким образом было промеряно 6 точек для значений широты:

0,248221
19,289529
41,327905
56,827437
71,934742
86,949039

3. Результаты.

Если нанести полученные значения на график, то получается интересная картинка:



Во-первых, разрешения в угловых секундах (ось Y) линейно зависят от масштабного уровня (ось X) для каждого фиксированного значения широты.

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

В-третьих, все линии проходят через 0.

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



где коэффициент K является функцией от широты. Естественно, что эти коэффициенты будут разными по горизонтали и вертикали.

Если использовать МНК и посчитать эти коэффициенты для имеющегося набора значений широты, то получается вот такая картинка:



По оси X значение широты (в радианах), по оси Y - значение коэффициента зависимости масштабный уровень - разрешение.

Синяя линия - горизонтальное разрешение, красная - вертикальное.

Как выяснилось, зависимость линеаризуется подстановкой

(в радианах) или
(в градусах):



Точнее, почти линеаризуется. Наблюдается некоторая остаточная нелинейность, но она может быть связана и с ошибками измерений.

Тем не менее, параметр "R-квадрат" для обоих составляющих достаточно близок к единице чтобы можно было говорить о линейном характере зависимости.

СКО для линейной зависимости составляет 0,002060 по горизонтали и 0,001319 по вертикали.

Если учесть остаточную нелинейность и аппроксимировать эту зависимость квадратичной функцией, то СКО будут меньше и равны, соответственно, 0,000521 и 0,000446

Итак, в конечном итоге получаем формулу зависимости разрешения в угловых секундах на пиксель от широты и масштабного уровня:



линейная зависмость, значения коэффициентов для горизонтального разрешения:

a = 0,679008
b = -0,157750

для вертикального:

a = -0,623140
b = 1,144693

широта в данном случае в радианах.

Более точная формула:



с коэффициентами:

a = 0,084945
b = 0,445941
c = -0,004917

по горизонтали, и

a = 0,052884
b = -0,768238
c = 1,239842

по вертикали.

Этого уже достаточно для расчета парамеров ресамлинга карты в GM под нужный масштабный уровень. При желании несложно пересчитать разрешения из угловых секунд в метры.

_________________
//Cheers, Victor
//Magellan eXplorist 310, Transystem TripMate 852, Columbus V-990, Transystem PersonalMate 860Lite, Garmin Nuvi 50

Вот, быстренько слепил калькулятор. Проверяйте.

_________________
//Cheers, Victor
//Magellan eXplorist 310, Transystem TripMate 852, Columbus V-990, Transystem PersonalMate 860Lite, Garmin Nuvi 50

Смотрел ф-лу
по ходу М у вас в километрах.
То есть если нужно посчитать для 60м, ми вводим значение 0,06?
Делаю тестовые RMP на широте 52,359 под 60м, и 240.

Да. Старался чтобы диапазон значений по x и y при вычислениях был примерно одного порядка. Ну насколько это возможно. Посему километры и угловые секунды. Ну это уже тонкости компьютерной математики.

Кстати, широту можно не переводить в радианы. Просто делить на 90 и от отношения (оно получается безразмерным) сразу брать косинус.

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

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

Kx(lat) = K0 + K(lat)
Ky(lat) = K0 - K(lat)

Где K0 - коэффициент для экватора, а K(lat) универсальная поправка на широту.

Я такой вариант пробовал, но у меня СКО для измеренных данных масштаб-разрешение получились больше чем в "несимметричном" варианте который используется сейчас.

_________________
//Cheers, Victor
//Magellan eXplorist 310, Transystem TripMate 852, Columbus V-990, Transystem PersonalMate 860Lite, Garmin Nuvi 50

сделал 4 тестовые RMP, комментарии и замечания завтра.
маленькие для 240м тут же

два других растра из z18 SAS для 60м и геотифы из z16 для 240м и из z18 для 60м тут

Основание - это нижняя часть? Насколько я заметил, у Магелланов на точку указывает центральная часть иконки. Поэтому со флажком получается не совсем корректно. Лучше в качестве иконки использовать крестик.

Программка у вас получилась интересная. И результаты, которые она выдаёт, получаются близкими к реальности. Только хотелось, чтобы вы пояснили, в чём отличие между "Linear" и "Square". Непонятно, что лучше использовать. Мне показалось, что "Square" даёт более реальные результаты. А ещё у меня не получается ввести дробную широту. Запятая не вводится, а введение точки выдаёт ошибку и зависание программы.

_________________
Карти для GPS-навігаторів: де їх брати? https://encyclopedia-mandriv.blogspot.com/2013/06/karty-dlya-gps-navigatorov.html

Linear - линейная зависимость коэффициента K от широты, Square - квадратичная. См. Формулы выше - в первом случае два параметра (a и b), во втором - три (a, b и c)

С точкой - поправлю завтра. Опять забыл десятичный разделитель поменять.

Меня больше беспокоит разность разрешения в меирахипо осям. Сомневаюсь что это правильно.

_________________
//Cheers, Victor
//Magellan eXplorist 310, Transystem TripMate 852, Columbus V-990, Transystem PersonalMate 860Lite, Garmin Nuvi 50

Выше вы сами пишете "более точная формула" про формулу с тремя коэффициентами. Так зачем в калькулятор включать Linear?

_________________
Карти для GPS-навігаторів: де їх брати? https://encyclopedia-mandriv.blogspot.com/2013/06/karty-dlya-gps-navigatorov.html

Потому что, как написал ввше, не уверен что остаточная нелинейность не связана с ошибками измерения.

"Более точная" она только с точки зрения СКО измеренных данных. Так что все это еще нуждается в проверке.

_________________
//Cheers, Victor
//Magellan eXplorist 310, Transystem TripMate 852, Columbus V-990, Transystem PersonalMate 860Lite, Garmin Nuvi 50

Ну коль есть наработка и из практики формула и утилитка, то нужно проверить на практике.
По ф-ле Засядьько можна определить некоторые интересные моменты

Один из них , это широта , на которой мерные отрезки прибора не врут, она равна : 31,86273883.
Если ф-ра справедлива и работает, этот факт стоит проверить.
Другой момент, это широты, на которых тайлы из SAS один к одному, отображаются в навигаторе.
Для масштаба отображения 240м, это широта 52,35898812, тайлы уровня Z16 и таже широта но тайлы уровня z18.


Если к тому же справедлива ф-ла Поморцева, то в принципе значения Width of tile in degree: 0,0027466 Height of tiles in degree: 0,0016773 из .tlm для малых (по высоте) rmp из SAS z16 (z18) должны на этой широте совпадать с данными ф-лы Поморцева для 240м (60м).
Поэтому первым делом возле 52,36 пометил тайлы кеша САС z18? z16, получилась сетка 256х256.
Сделал тестовый rmp в SAS 4х4 тайла из z16 и та-же область в z18 16х16 тайлов.
В GM оформил геотифы вбивая данные от утилиты Поморцева для 52,35898812 и соответственно 240м и 60м и отдельно геотифы от САС ,
но чтоб нарезка тайлов в RMP 1:1 совпадала с тайлами в САС (разбирал контейнеры,проверял-100%cовпадение, то есть в z16_berlinSASgeotiff.tif.rmp и z18_berlinSASgeotiff.tif.rmp границы тайлов физически такие как я пометил ).
Сделал 4 тестовых RMP . , они же с геотифами прилагаются выше.
Смотрел геотиф z16 по данным калькулятора Поморцева в фотошопе, имеем несовпадение пикселов сетки, вместо 256х256 имеем 256х253....
Нужна практика,то есть скрины выложенных мной rmp