Форумы GIS-Lab.info

да, этот модуль нуждается в целом ряде исправлений и улучшений. Что касается его авторов (нашей группы) мы будем что-то исправлять только когда появятся ресурсы или необходимость, это касается всех наших проектов. Тем временем, код открыт, править его могут любой другой человек. Это конечно не означает, что не нужно сообщать об ошибках и предлагать улучшения. Когда и если у нас или у кого-то еще дойдут руки это исправить, ему будет полезна эта тема.

Про сплющенность вам стоит почитать вот тут: http://gis-lab.info/qa/proj-sk-faq.html#16

спасибо большое за статью и за модуль. Все получилось сделать прям с первого раза

Куда сохраняется файл после трансформации? В указанном месте его нет. И как его отобразить в окне программы?

Здравствуйте!
Делала по инструкции привязку - получилось нечто непонятное. Видимо, буду первым человеком с 2013 года, кто задаёт подобного рода вопросы )
1. почему в статье есть фраза "Если координаты находятся в формате отличном от dd.dddddd (десятичные градусы), их нужно будет к нему привести", но при этом в вылезающем окне указан иной формат: [я не являюсь носителем мега-знаний о форматах координат, поэтому могу что-то не допонять - если поясните, в чём тут дело - буду благодарна]

2. всё же я в свой растр внесла данные в формате dd.dddddd (тыкала в карты гугла и переносила так точки на свой растр).
Что я сделала не так?? Вот так было в процессе: После привязки растров вдруг стало два и их вид меня совсем не впечатляет (((
[т.к. ограничение по вложениям - 3 шт, то объединила два изображения] Я что-то опять по незнанию с системами координат накосячила или что-то ещё, неизвестное мне, натворила?
[система координат растра WGS 84 / UTM zone 37N, версия QGIS 2.14.0]

Похоже, что вы перепутали Х с Y.
Еще не очень понятно, что вы привязывете и к чему.

Во-первых, Вы пытаетесь привязать снимок Landsat-8 от USGS. Зачем? Ведь он уже с геопривязкой. Кстати, это видно по значениям "X источника", "Y источника" в Вашей таблице привязки.

Я скачал эту же сцену и просто открыл: Границы административных районов вместо банана.

Во-вторых, указав проекцию растра UTM zone 37N, следует вводить координаты в единицах именно этой проекции, которые метры, а не градусы долготы/широты.

В-третьих, darsvid прав, координаты всегда вводятся в порядке "восточная—северная", т. е. долгота—широта.

Несколько вопросов по статье:

1) Если выбрать тип трансформации Линейная, то становится доступной опция создания world-файла. А в чём собственно отличие этого типа трансформации от Полинома 1 степени. Судя по формулам - это одно и то же. Или нет?
2) В каких случаях лучше применять трансформацию Гельмерта, а не полиномиальные преобразования?

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

Denis Rykov писал(а):1) Если выбрать тип трансформации Линейная, то становится доступной опция создания world-файла. А в чём собственно отличие этого типа трансформации от Полинома 1 степени. Судя по формулам - это одно и то же. Или нет?

Не одно. Полином 1 степени - аффинное преобразование:
x = A x + B y + C
y = D x + E y + F
"Линейная трансформация" - полином 1. степени с ограничениями:
x = A x + (0 y) + C
y = (0 x) + E y + F

Denis Rykov писал(а):2) В каких случаях лучше применять трансформацию Гельмерта, а не полиномиальные преобразования?

В случаях, когда преобразование должно быть конформным. Преобразование Гельмерта - это полином 1. степени с ограничениями:
x = A x + B y + C
y = -B x + A y + F

Denis Rykov писал(а):в статье приведён пример привязки карты, когда в соответствие пиксельным координатам ставятся градусы. Так ведь делать неправильно?

Если пользователь считает, что проекция карты "Широта/долгота", он должен использовать соответствующие единицы - градусы.

То есть QGIS при создании world-файла - обеспечивает только часть его возможностей, так как часть коэффициентов заведомо равны 0? Коэффициенты world-файла в общем случае - это и есть параметры полинома 1 степени?

Denis Rykov писал(а):Коэффициенты world-файла в общем случае - это и есть параметры полинома 1 степени?

Точно так.

Denis Rykov писал(а):То есть QGIS при создании world-файла - обеспечивает только часть его возможностей, так как часть коэффициентов заведомо равны 0?

Цель трансформирования растра - получить грид: колонки за X, строки за Y. Поэтому после трансформирования коэффициенты B и D непременно станут нулями.

В этом плане добавленная "линейная трансформация" очень полезна. Она позволяет нивелировать ошибки (нашей) привязки для заведомо трансформированных (кем-то ранее) растров.

Насчёт других возможностей (открыть растр с ненулевыми коэффициентами в world-файле) не экспериментировал.

При открытии растра в окне привязчика он спрашивает систему координат. На что влияет выбор системы координат на данном шаге? Выходная система координат задаётся в настройках трансформации, а эта для чего - не понятно.

Denis Rykov писал(а):Выходная система координат задаётся в настройках трансформации, а эта для чего - не понятно.

На выходе можно сразу получить растр в проекции, отличной от исходной: "gdalwarp -t_srs"…

Denis Rykov писал(а):При открытии растра в окне привязчика он спрашивает систему координат. На что влияет выбор системы координат на данном шаге?

…но для этого надо знать исходную проекцию: "gdal_translate -a_srs". Да и при сохранении трансформированных растров в исходной проекции GDAL должен что-то писать в их тэги SRS.

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

Так и есть, если привязчик используется только для создания world-файла. Если же производится трансформация, GDAL пишет в тэги выходного растра GeoTIFF проекцию, которую пользователь должен указать.

Так я её и указываю в настройках привязчика ("Target SRS") и с этой проекцией всё понятно. Непонятно что за проекцию с меня спрашивают при открытии растра в привязчике. Я ведь не трансформирую растр из СК в СК, для чего мне нужно было бы знать исходную систему координат, а просто трансформирую растр по точкам и назначаю ему выходную проекцию.