Систему координат нужно подобрать таким образом, чтобы пиксел изображения меньше всего искажался. Для выбора системы координат, исходное изображение было привязано по углам с использованием координат из файла 3v. После этого файл теперь уже в формате Geo. TIFF был загружен в QGIS и на нем были испробованы различные системы координат с преобразованием. При этом, обращалось внимание, чтобы пиксел изображения оставался. В результате экспериментов было принято решение об использовании системы координат WGS8. UTM соответствующей зоны. Для получения файла рельефа необходимо. На странице SRTM Data Selection Options выписываем номера рядов и столбцов с территорией съемки. Скачиваем с http gis lab. Распаковываем команда для операционной системы UNIXfor i in srtmzip do yes. Получаем несколько в данном случае 4 файлов с данными DEM с именами вида ASTGTM2. Склеиваем все файлы в один Geo. TIFF с помощью gdal. В случае если команды выполнилась с ошибкой, следует прописать систему координат и экстент исходного изображения. Попробовать извлечь из метаданных. Минус варианта выходная система координат будет WGS8. Прописать систему координат и GCP углов изображения. Для этого в файле scene. СК WGS8. 4, пересчитываем их в WGS8. UTM 3. 5N и добавляем их к файлу. Для устранения этого расхождения необходимо определить превышения геоида над эллипсоидом WGS8. Для определения превышения, возьмем координаты средней точки снимка 2. В результате получим 2. EGM9. 6. Если при этом gdalwarp отказывается трансформировать что либо, скачайте соседние тайлы GDEM скорее всего,. Выполним в этом программном продукте аналогичные операции и сравним результаты. Приведу изображение одного и того же участка местности с наложенными на него GPS треками. Результат выполнения ортокоррекции с учетом геоида полностью идентичен результатам работы ПО ENVI EX. Причем, ортокорреция, выполненная в программе wx. GIS, тоже основанной на библиотеке GDAL, не приводит к появлению артефактов. Данное необъяснимое поведение библиотеки версия GDAL 1. Для устранения артефактов в командую строку необходимо добавить ключ et 0. Пример команды. gdalwarp dstnodata 0 srcnodata 0 overwrite t.