Определение параметров равноугольной конической проекции Ламберта для растра карты

[lukon]

Доброе время суток, уважаемые геодезисты. У меня возникла проблема: есть электронное изображение карты в равноугольной конической проекции Ламберта (размер карты от Казани до Питера), нужно определить параметры этой проекции (широта двух опорных параллель, центральный меридиан и пр.). Известно, что значения искомых широт/долгот кратны 10 минутам, датум ПЗ-90.02 или WGS-84. Кто может подсказать, в какой доступной ГИС/программе это можно посчитать?

Не совсем по теме: чтоб не было лишних вопросов, хочу преобразовать проекцию карты в псевдомеркатор (EPSG:3857). Можно конечно указать несколько точек и выполнить привязку растра, но тогда наверняка появятся ошибки, а точность мне нужна хотя бы 100 метров

[gamm]

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

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

[lukon]

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

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

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

[gamm]

обычно калькуляторы считают параметры поперечного меркатора (UTM/GK)

1) если карты в одной проекции, то хватит точек на одну 5-9, конвертом (4 по углам примерно, одна в центре, еще 4 на полпути от центра к углам все примерно)
2) можно программу набросать на R, если знаете программирование, то остальное допилите.
3) нужна искомая проекция в формате +proj, с указанием свободных (подбираемых) параметров.

[Boris]

1. для начала стоит перебрать основные советские проекции, только на новой СК, в русской Маинфо есть набор, в который входят проекции :

"--- Прямые конические проекции России и бывш. СССР ---"
"Равноугольная Ламберта-Гаусса (СССР)", 3, 1001, 7, 100, 54, 50, 70, 0, 0
"Равнопромежуточная Каврайского (СССР)", 6, 1001, 7, 100, 54, 47, 62, 0, 0
"Равнопромежуточная Каврайского (Европ. ч. СССР)", 6, 1001, 7, 40, 54, 47, 62, 0, 0
"Равнопромежуточная (РСФСР)", 6, 1001, 7, 100, 54, 49, 68.5, 0, 0
"Lambert Conformal Conic (S-42) with boubds for Euro-RF\p3034", 3, 1001, 7, 40, 54, 47, 62, 7000000, 2800000 Bounds (500000, 1100000) (13500000, 8100000)

Где-то на форуме была давняя статья или ветка про проекции СССР/РСФСР/России, там могла быть и другие варианты.
2. создать список с параметрами PRJ можно в любом минимальном программаторе, хоть в Excel, и скормить эти PRj в gdal_warp указав исходной перебираемые проекции, а результирующей какую-то из тех, что позволяет посмотреть ваша ГИС-программа.
3. вы сказали про точность пересчета, но не сказали про тестовые данные и размер исходного изображения. 100 точности - это минимум 1 пиксель = 100 м (хотя в реальности и близко не так). От Питера до Казани минимум 1000 км по ширине картинки. В карте есть 10 тысяч пикселей в исходнике?
4. при наличии большого числа точек, хотя бы по точке на каждый градус, пересчет через сплайн трансформацию даст приемлемую точность в 2-3 пикселя исходной карты. Главное, чтобы было чем перепроверить.

[tikhpetr]

А можно посмотреть?

[lukon]

Не могу создать вложение, ссылка на одну из карт: https://disk.yandex.ru/i/ZVfgEDIBAcvSeg. Размер 8740 x 5256, при сохранении карты в .png пишет 300 точек на дюйм, но в пдф думаю больше. Масштаб 1 см = 20 км. По поводу ответа Boris, всего подобных карт, которой я скинул: 17 штук (видно по картинке на карте). На каждой разные опорные параллели и наверняка другие параметры (потому что есть еще похожие карты на эти, их 15 штук и там опорные параллели подписаны, почти у всех 15 карт эти параллели уникальны). Итого ~15 карт с разными параметрами проекций (допустим, 2 раза одинаковые параметры встретились). Найдется ли столько проекций LCC в какой-либо из программ? 2-3 пикселя думаю хорошая точность

[Boris]

Масштаб 1 см = 20 км. 2-3 пикселя думаю хорошая точность

Канон российской картографии гласит:"Особо-ответственные объекты карты имеет точность отображения положения равную 0.1 мм в масштабе карты. Основные объекты - 1мм в масштабе карты." Как было неоднократно отмечено специалистами по векторизации карт - "сетка координат точно не является объектом, точность которого нанесена хотя бы с точностью 1 мм". Так что, если заявленный масштаб имеет место быть, точность привязки даже в известной проекции не будет выше 1 мм = 2 км.

[Boris]

Я не нашел в присланном образце указания на проекцию. Она откуда так точно известна? Если это карта ИКАО, то в ней должен в основе быть эллипсоид WGS84 или какой-то из его современных аналогов. Лет 15 назад я пересекался со специалистами местных аэродромов, они жаловались на то что очередные документы на эти аэродромы требовали указания координат в ВПП в WGS84, часто там, где они были в документах указаны в довоенной СК.
В остальном, если известен эллипсоид, пересчитать растр под любую заранее выбранную проекцию, не является проблемой - точек с координатами на листе хватает. Главное понять, в таком ответственном деле, с каким эталоном сравнить.

[lukon]

Всем спасибо, проблема решена путем обращения к составителю карт.Boris, спасибо за интересную информацию

[tikhpetr]

Хорошо, что получилось. Один вопрос. Совершенно не разбираюсь в летных картах. Так и должно быть, что долгота увеличивается на запад? Или я неправильно смотрю?

[lukon]

Север там слева

[tikhpetr]

Принято. Спасибо!