Форумы GIS-Lab.info

gamm писал(а):Исходно SRTM - высоты над геоидом EGM96, нужно их пересчитать в высоты над эллипсоидом.

значение пикселя является высотой над уровнем моря в данной точке

Я в сметении)

Модель геоида примерно соответствует уровню моря там, где море есть. С внутренними морями, конечно, всё сложнее, но не думаю, что в вашей задаче это существенно.

От артефактов избавился, косяк был в последовательности соединения вершин. А вот с искаженными очертаниями вопрос до сих пор открыт.
Может по-диагонали отсечь треугольник и перенести его вправо? получится параллелограмм. Или это я бред несу?)

Тоже больше похоже на ошибку в коде при обработке растра. Не случайно же он у вас наискосок разрезан. Где-нибудь переменную цикла перепутали или переполнение словили.
Получившийся из исходного квадрата фрагмент должен напоминать трапецию.

Selot писал(а):Или это я бред несу?)

для начала выплюньте не стрипы, а облако точек (центров SRTM, для начала лучше заменить высоту на константу), и посмотрите, что получится - должна быть вытянутая трапеция. Похоже на баг ...

Глупый косяк и странный. Файл сртм представлен матрицей 1201х1201, я читал только 1200х1200. Стоило прочитать весь файл и все встало на свои места. Не совсем понимаю почему оно так. Да и трапеция не вышла, а обычный прямоугольник

В средних широтах трапеция и не должна сильно отличаться от прямоугольника, так что, я думаю, всё правильно.
Если вы читали только 1200 значений из строки матрицы и не пропускали 1201-е, а сразу начинали читать следующую строку - понятно, что это значение оказывалось в начале следующей строки. Это всё объясняет.

Рассмотрим участок поверхности Земли размером 1000 × 1000 км. Расстояние от середины стороны до её концов 500 км, что в угловой мере равно ρ = 500 / R = 500 км / 6371 км = 0,078 рад = 4,5°. Вычислим стрелку прогиба середины по отношению к краям: h = R (1 - cos ρ) ≈ R ρ² / 2 = 19,6 км. В отношении к длине это 2%. Такая кривизна видна на трёхмерном изображении:Смотрите на ближнюю и дальнюю стороны, это дуги меридианов. Правая и левая стороны образованы параллелями, они в плане кривые.

Другое дело градусные тессеракты SRTM, — они в девять раз меньше в меридианальном направлении и почти в тринадцать раз в субширотном на широте Крыма. Относительный прогиб дуги меридиана 0,2%, так что кривизна не видна глазом:

Всем спасибо за советы и помощь)