GDAL, сравнение двух Geotiff

[Rumato]

Добрый день, у меня такой вопрос: есть два геотифа, полученных из одного HDF-файла разными программами (одна программа точно правильно работает, нужно проверить другую, самописную), они должны быть одинаковыми, но нужно сравнить чтобы убедиться.

В GDAL есть какие-либо функции чтобы это проделать? Или чем ещё можно так проверить?


Заранее большое спасибо за помощь!

[Дмитрий Барышников]

В gdal если только сравнить вывод gdalinfo, но сами пикселы так не проверишь.

[Максим Дубинин]

если выхлоп gdalinfo одинаков, то дальше тупо корреляция или вычесть растры друг из друга

[Rumato]

Максим, а в смысле вычесть растры, типа один двумерный массив из другого? Если ноль то тогда равны, так?

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

Сначала сделать diff для вывода gdalinfo, затем перегнать тифы в текстовый формат (например, Arc/Info ASCII Grid) и сделать diff уже по самим данным.

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

Если по-простому сравнить (отличаются или нет), то хватит обычного diff:

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

$ diff -q rast.tif rast1.tif
Файлы rast.tif и rast1.tif различаются

Или можно исп-ть подобные функции, например, на Питоне (модуль filecmp).

[gamm]

Если по-простому сравнить (отличаются или нет), то хватит обычного diff:

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

$ diff -q rast.tif rast1.tif
Файлы rast.tif и rast1.tif различаются

Или можно исп-ть подобные функции, например, на Питоне (модуль filecmp).

лучше сравнивать содержимое (вывалив в формат без метаданных типа Surfer), там могут какие-то метаданные не совпасть

[Rumato]

Большое спасибо за ответы, все ответы пригодились(или ещё пригодятся). Спасибо!

[ericsson]

можно еще сделать наложение difference в ImageMagick и считать статистику по нему