Как добавить на карту линию полярного круга?

[imforlorn]

Начинаю осваивать QGIS. Уже который день бьюсь над вроде бы простой задачей: нужно на карту региона (shape в Quantum GIS 2.4.0, взял отсюда) нанести линию полярного круга. Нашел саму линию в .shp здесь.
Понял, что они физически в разных СК, но неужели никак нельзя преобразовать? Пробовал ogr2ogr, не получилось. Также пробовал менять СК в QGIS, но это вроде вообще не то. Регион в одной стороне, эти линии горизонтальные в другой.
Подскажите, может есть решение с полярным кругом.

[trir]

QGIS вообще то должен спокойно совмещать данные в разных СК
Какие там .prj файлы?

[imforlorn]

у линий:

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

GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.017453292519943295]]

у региона:

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

PROJCS["unnamed",GEOGCS["unnamed ellipse",DATUM["unknown",SPHEROID["unnamed",6378137,0],EXTENSION["PROJ4_GRIDS","@null"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Mercator_2SP"],PARAMETER["standard_parallel_1",0],PARAMETER["central_meridian",0],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1],EXTENSION["PROJ4"," +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs"]]

[trir]

Вообще то у ВСЕГЕИ должен бвыть такой .prj

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

PROJCS["Equidistant_Conic",GEOGCS["GCS_Krassovsky, 1942",DATUM["D_unknown",SPHEROID["krass",6378245,298.3]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Equidistant_Conic"],PARAMETER["standard_parallel_1",46.4],PARAMETER["standard_parallel_2",71.8],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",100],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]

viewtopic.php?t=12624&p=80001

Цифровые географические основы подготовлены с использованием программного обеспечения ArcView и ArcGIS в конической равнопромежуточной проекции на эллипсоиде Красовского с центральным меридианом 100° в.д. и главными параллелями 46,4° и 71,8° с.ш.

P.S. А то, что у вас - это сферический конь в вакууме

[imforlorn]

Да, возможно я перепутал. Уже по-всякому пробовал. Но это не меняет сути вопроса.

[trir]

Это и есть - суть вопроса, пропишите мои параметры в .prj ВСЕГЕИ и попробуйте открыть

[imforlorn]

Попробовал. Все равно линии, горизонтальные, где-то в одном месте. Регион в другом.

[trir]

Попробуйте в новом проекте - вставляя в него .shp'ы с правильными .prj
Установите СК проекта - WGS84LL и посмотрите координаты данных (скриншот тоже не помишает)

[imforlorn]

Вот:

ss1.JPG (124.16 КБ) 13160 просмотров

ss2.JPG (139.04 КБ) 13160 просмотров

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

Создаём пользовательскую проекцию ВСЕГЕИ с параметрами

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

+proj=eqdc +lat_1=46.4 +lat_2=71.8 +lon_0=100 +x_0=0 +y_0=0 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs

Открываем слои, указывая для них эту проекцию. Добавляем географические линии.

capture_20140923_101334.png (114.03 КБ) 13148 просмотров

Сохранив копию одного из слоёв, получим общие для всех файлы CPG, PRJ, QPJ.

[imforlorn]

Не понял, как так у вас получилось, а у меня нет.
Может что не так делаю: QGIS - Установки - Ввод системы координат - Добавить - В поле имя пишу любое - В поле параметры вставляю ваш код - ОК - Создаю проект - Добавляю слои (заранее удалив их .prj файлы) - указываю созданную систему координат. И в итоге у меня те линии так и остаются горизонтальные, вообще не меняются.

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

Перед добавлением слоя ne_10m_geographic_lines.prj в свойствах координатной системы карты поднимите флажок "Enable 'on the fly' CRS transformation" (включить преобразование координат на лету).

capture_20140923_114736.png (60.61 КБ) 13112 просмотров

А у меня в настройках CRS выбрано "Automatically enable 'on the fly' reprojection if layers have different CRS".

capture_20140923_113935.png (51.11 КБ) 13113 просмотров

[imforlorn]

Это всё стоит.

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

Попробуйте открыть файл проекта "project1.qgs". Только сначала исправьте пути к слоям в строках 61 и 165 или поместите слои "polta" и "ne_10m_geographic_lines" в одну папку с файлом проекта.

[imforlorn]

Открылось с вашего проекта!
Что же это значит? Почему в моем проекте не получается..