Форумы GIS-Lab.info

Для зональной статистики (чтобы не разбивать анализ на много кусков) склеиваю несколько растров (продукт MODIS - температура в кельвинах) при помощи Raster- Miscellaneous-Merge (QGIS 2.18). Несмотря на то, что тайлов продукта было 8, он успешно и быстро склеил все в один файл.
Но при точной проверке вижу, что значения пикселей немного отличаются (н-р 301.14 и 300.06 К) Также произошло смещение (что заметно на краях, когда NoData сместилась).

Это нормальное явление при этой операции для склеивания растров? Мне казаось, что оператор Merge просто должен склеивать без искажений. Хоть вроде оно и небольшое, но почему вообще оно произошло? Ведь я не задавала никаких доп. операций. Может быть, настройки были не те?
Сам этот продукт сам поделен тайлами очень точно. На первый взгляд нет наложений.

Вообще, мне важно понять, насколько корректно после такого склеивания растров делать зональную статистику. Или лучше ее делать по исходным продуктам, а склеивание дает слишком большое искажение..
Вопрос точности очень важен. От этого я цепляюсь к деталям.

Также в операции Merge он превращает NoData в "0", что вносит искажение в зональную статистику (NoData ей просто не видится, а 0 считается за число). В настройке операции Merge не нашла способ оставить NoData . Сравнение по оригиналу и по растру merged показало значимые различия в зон. статистике. Вероятно, они связаны как раз с превращением NoData в "0".
Но при этом медиана и максимум обоих растров - почти одинаковы. Что говорит все же о том, что склеенный растр можно использовать для зональной статистики. Но нужно брать не среднее, а медиану... (так уйдет ошибка от превращения NoData в "0").

Natalia Novoselova писал(а): ↑

27 окт 2017, 03:21

Также в операции Merge он превращает NoData в "0", что вносит искажение в зональную статистику (NoData ей просто не видится, а 0 считается за число). В настройке операции Merge не нашла способ оставить NoData . Сравнение по оригиналу и по растру merged показало значимые различия в зон. статистике. Вероятно, они связаны как раз с превращением NoData в "0".
Но при этом медиана и максимум обоих растров - почти одинаковы. Что говорит все же о том, что склеенный растр можно использовать для зональной статистики. Но нужно брать не среднее, а медиану... (так уйдет ошибка от превращения NoData в "0").

В свойствах исходных тайлов посмотрите значение "NoData" и когда будуте сшивать в меню укажите его, а не "0".
В итоге получите сшитый растр со значениями "NoData" как в исходниках.