Привязка листа из атласа четвертичных отложений

[xen87]

Доброе время суток!
Столкнулся с задачей привязки листа (№7) из атласа четвертичных отложений европейской части СССР 1971 года, М 1:1500000.
Привязывать буду в модуле привязки растров QGIS.
Возникло ряд вопросов :
1. используемая система географических координат для подобного атласа, предполагаю, что - Альберса?
2. если и в правду Альберса, то какие параметры необходимо задавать при создании пользовательской СК? Поискал по сайту нашел здесь: http://gis-lab.info/qa/grids-qgis.html для ХМАО (+proj=aea +lat_1=60 +lat_2=64 +lat_0=0 +lon_0=72 +x_0=18500000 +y_0=0 +ellps=krass +units=m +towgs84=28,-130,-95,0,0,0,0 +no_defs), а для скажем Московской области??
3. Если есть, поделитесь пожалуйста ссылкой - каким образом формируется подобная строка указания СК, чтоб хоть чуть-чуть разбираться в подобной "китайской грамоте" ...

Заранее большое спасибо.

[Александр Мурый]

1. Было бы неплохо посмотреть на саму карту для начала..
2. Может быть, Albers-Europe:

Код: Выделить всё

+proj=aea +lat_1=52 +lat_2=64 +lat_0=0 +lon_0=45 +x_0=8500000 +y_0=0 +ellps=krass +units=m +towgs84=28,-130,-95,0,0,0,0 +no_defs

3. Это строка описания СК/проекции в формате PROJ4; здесь написано про основные параметры + поиск на форуме даст кое-что

[xen87]

Большое спасибо за отклик!!!
Прикрепил к первому посту лист из атласа.

[Донецков]

То что это не Альберса это 100%, посмотрите в выходных данных карты, если не ошибаюсь на последнем листе набора...
Скорее всего это Гаусса-Крюгера, только нужно создавать свои параметры с ЦМ по середине кажого листа...

Вы собираетесь привязать только один лист или весь набор?

Если интересно, посмотрите на сайте Геокнига в разделе карты приведено много тематических карт, в т.ч. есть и 1.5 млн. смежных тематик на эту же территорию...

[Ariki]

Для Гаусса-Крюгера слишком большая протяженность по долготе. Проекция похожа на коническую. Центральный меридиан, безусловно, 36-й.

[Игорь Белов]

Подходит коническая равнопромежуточная Каврайского для европейской части СССР с переносом осевого меридиана на 36°в.д. Соответственно набор параметров выглядит так:

Код: Выделить всё

+proj=eqdc +lat_1=47 +lat_2=62 +lat_0=54 +lon_0=36 +x_0=0 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +no_defs

P. S. И ещё раз повторю: не ищите на советских картах проекцию Альберса. В ГУГК обходились без неё. Если нужна была равноплощадная проекция, использовали равноплощадную азимутальную.

[xen87]

Большое спасибо всем откликнувшимся!!!
Отвечу на вопросы:
1. Мне необходимо привязать только один лист из этого атласа (лист №7).
2. Весь атлас просмотрел, но данных о системе координат и проекции не нашел
3. Большое спасибо за ссылку на ГеоКнигу.
Проблема:
Воспользовался параметрами СК/проекции из поста ErnieBoyd, но тут у меня возникли проблемы:
1. карта деформировалась ... см. вложение.
2. деформированная карта "улетела" относительно точек привязки на 2 градуса на запад и на 10 градусов на юг (видно в теплые страны собралась ...).
Что делать? В свойствах проекта и при привязке модулем QGIS выставил пользовательскую СК/проекцию Каврайского из предыдущего поста.

[Донецков]

Так как это геологическая карта, то и обращаемся к ИСТОКАМ, находим сайт ВСЕГЕИ, скачиваем картографическую основу, правда ближайшая только 1:2500000, и для начала привязываем к ней...
Параметры СК этих основ на сайте ВСЕГЕИ опубликованы...
Также можете посмотреть параметры отечественных СК для тематических карт на сайте ЦГИ ИГ РАН...

[Донецков]

Соответственно набор параметров выглядит так:
Код: Выделить всё+proj=eqdc +lat_1=47 +lat_2=62 +lat_0=54 +lon_0=36 +x_0=0 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +no_defs

почему использованы параметры из ГОСТ 51794-2001 года, а не ГОСТ 51794-2008:

23.57,-140.95,-79.8,0,-0.35,-0.79,-0.22

?

[Игорь Белов]

почему использованы параметры из ГОСТ 51794-2001 года, а не ГОСТ 51794-2008?

• не хочется в QGIS заводить кучу пользовательских CRS вместо существующих "Pulkovo 1942";
• разница для большинства задач несущественна;
• к тому же в данном случае разрешение растра 400-500 метров;

Ну а на пленэре обе редакции ГОСТа одинаково не работают без коррекции к имеющимся пунктам ГГС.

[Игорь Белов]

1. карта деформировалась...
2. деформированная карта "улетела" относительно точек привязки...

Пока не пробовал трансформировать карту в QGIS'е. Сделал афинное преобразование в Global Mapper'е по шести пересечениям градусной сетки. Создал пользовательскую проекцию, открыл в ней растр и слой пересечений.

[Игорь Белов]

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

Код: Выделить всё

+lat_1=49 +lat_2=64

Растр на выходе получается несколько краше, однако тоже не фонтан.
Донецков прав. Если есть возможность добраться до выходных данных карты, это надо сделать.
И по поводу картографической основы ВСЕГЕИ Донецков тоже прав, только она нужна не для предварительной привязки, а для чистовой: большое количество геоточек, равномерно покрывающих весь лист, и трансформирование локальными сплайнами или локально-афинное с триангуляцией Делоне.

[xen87]

Добрый день!!!
Большое спасибо ВСЕМ за помощь!!!
Последовал советам ErnieBoyd:
1. привязывал лист атласа в ГлобалМаппере, данные о пользовательской СК/проекции взял из файла вложения предыдущего поста. Карта получилась чуть-чуть деформированной (очень похожа на карту из предыдущего поста ).
2. затем в том же самом ГлобалМаппере изменил СК/проекцию на ВГС84 на нужную мне 38 зону (интересует только Костромская область).
Результат во вложении.
На рисунке вектор озер (сверху Чухломское, а снизу Галичское) наложен на лист атласа. Видно что лист атласа смещен относительно озер примерно на 10 км на запад (по широте практически совпадают). Нельзя ли каким-то простеньким инструментом сместить подложку (лист атласа) на эти 10 км ( подобный инструмент в Ердасе вроде есть, но я забыл как он называется) ... Конечно понимаю, что это халтура и делать надо как пишет Донецков, но супер точности от меня (пока) не требуют ...
Заранее большое ВАМ спасибо!!!

[Александр Мурый]

На рисунке вектор озер (сверху Чухломское, а снизу Галичское) наложен на лист атласа. Видно что лист атласа смещен относительно озер примерно на 10 км на запад (по широте практически совпадают). Нельзя ли каким-то простеньким инструментом сместить подложку (лист атласа) на эти 10 км

Если вы выложите привязанный растр и слой озёр, могу в GRASS подвинуть растр куда надо (там это делается в одну команду) и выдам в исходном формате.

[xen87]

amuriy большое Вам спасибо за помощь.
Шейп файл озер прикладываю к сообщению, а подложка (лист атласа) лежит по адресу (файл большой 55.7 МБ): http://files.mail.ru/GGEP1R