Искажение при склеивании растра (Raster- Miscellaneous-Merge)

Вопросы по свободной ГИС QGIS. Сообщения об ошибках, предложения по улучшению, локализация.
Ответить
Аватара пользователя
Natalia Novoselova
Гуру
Сообщения: 3020
Зарегистрирован: 15 янв 2013, 20:14
Репутация: 69
Ваше звание: Лиса
Откуда: **
Контактная информация:

Искажение при склеивании растра (Raster- Miscellaneous-Merge)

Сообщение Natalia Novoselova » 27 окт 2017, 00:12

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

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

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

Аватара пользователя
Natalia Novoselova
Гуру
Сообщения: 3020
Зарегистрирован: 15 янв 2013, 20:14
Репутация: 69
Ваше звание: Лиса
Откуда: **
Контактная информация:

Re: Искажение при склеивании растра (Raster- Miscellaneous-Merge)

Сообщение Natalia Novoselova » 27 окт 2017, 03:21

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

Аватара пользователя
Prince660
Активный участник
Сообщения: 164
Зарегистрирован: 25 мар 2011, 20:12
Репутация: 15
Ваше звание: Питоновод
Откуда: г. Симферорполь
Контактная информация:

Re: Искажение при склеивании растра (Raster- Miscellaneous-Merge)

Сообщение Prince660 » 03 ноя 2017, 12:37

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

Ответить

Вернуться в «QGIS»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 20 гостей