Для зональной статистики (чтобы не разбивать анализ на много кусков) склеиваю несколько растров (продукт MODIS - температура в кельвинах) при помощи Raster- Miscellaneous-Merge (QGIS 2.18). Несмотря на то, что тайлов продукта было 8, он успешно и быстро склеил все в один файл.
Но при точной проверке вижу, что значения пикселей немного отличаются (н-р 301.14 и 300.06 К) Также произошло смещение (что заметно на краях, когда NoData сместилась).
Это нормальное явление при этой операции для склеивания растров? Мне казаось, что оператор Merge просто должен склеивать без искажений. Хоть вроде оно и небольшое, но почему вообще оно произошло? Ведь я не задавала никаких доп. операций. Может быть, настройки были не те?
Сам этот продукт сам поделен тайлами очень точно. На первый взгляд нет наложений.
Вообще, мне важно понять, насколько корректно после такого склеивания растров делать зональную статистику. Или лучше ее делать по исходным продуктам, а склеивание дает слишком большое искажение..
Вопрос точности очень важен. От этого я цепляюсь к деталям.
Искажение при склеивании растра (Raster- Miscellaneous-Merge)
- Natalia Novoselova
- Гуру
- Сообщения: 3020
- Зарегистрирован: 15 янв 2013, 20:14
- Репутация: 69
- Ваше звание: Лиса
- Откуда: **
- Контактная информация:
- Natalia Novoselova
- Гуру
- Сообщения: 3020
- Зарегистрирован: 15 янв 2013, 20:14
- Репутация: 69
- Ваше звание: Лиса
- Откуда: **
- Контактная информация:
Re: Искажение при склеивании растра (Raster- Miscellaneous-Merge)
Также в операции Merge он превращает NoData в "0", что вносит искажение в зональную статистику (NoData ей просто не видится, а 0 считается за число). В настройке операции Merge не нашла способ оставить NoData . Сравнение по оригиналу и по растру merged показало значимые различия в зон. статистике. Вероятно, они связаны как раз с превращением NoData в "0".
Но при этом медиана и максимум обоих растров - почти одинаковы. Что говорит все же о том, что склеенный растр можно использовать для зональной статистики. Но нужно брать не среднее, а медиану... (так уйдет ошибка от превращения NoData в "0").
Но при этом медиана и максимум обоих растров - почти одинаковы. Что говорит все же о том, что склеенный растр можно использовать для зональной статистики. Но нужно брать не среднее, а медиану... (так уйдет ошибка от превращения NoData в "0").
- Вложения
-
- statistic.jpg (129.86 КБ) 2018 просмотров
- Prince660
- Активный участник
- Сообщения: 164
- Зарегистрирован: 25 мар 2011, 20:12
- Репутация: 15
- Ваше звание: Питоновод
- Откуда: г. Симферорполь
- Контактная информация:
Re: Искажение при склеивании растра (Raster- Miscellaneous-Merge)
В свойствах исходных тайлов посмотрите значение "NoData" и когда будуте сшивать в меню укажите его, а не "0".Natalia Novoselova писал(а): ↑27 окт 2017, 03:21Также в операции Merge он превращает NoData в "0", что вносит искажение в зональную статистику (NoData ей просто не видится, а 0 считается за число). В настройке операции Merge не нашла способ оставить NoData . Сравнение по оригиналу и по растру merged показало значимые различия в зон. статистике. Вероятно, они связаны как раз с превращением NoData в "0".
Но при этом медиана и максимум обоих растров - почти одинаковы. Что говорит все же о том, что склеенный растр можно использовать для зональной статистики. Но нужно брать не среднее, а медиану... (так уйдет ошибка от превращения NoData в "0").
В итоге получите сшитый растр со значениями "NoData" как в исходниках.
Кто сейчас на конференции
Сейчас этот форум просматривают: Bing [Bot] и 11 гостей