|
||
|
Работа с растровыми данными в R: rgdal
R - статистический пакет с открытым кодом, обладающий огромным набором функций, которые несомненно могли бы пригодится любым специалистам имеющим дело с анализом данных, в том числе и пространственных. Целью данной статьи является демонстрация того, как можно осуществлять статистические операции с многоканальными растровыми данными, с которыми мы повседневно работаем в GIS, не переводя их для этого в табличный формат. Работа напрямую с растром экономит время, способствует более чистой логике процесса и значительно облегчает программирование. Работу с растрами напрямую делает возможной пакет RGDAL. Перед использованием примера необходимо убедиться что библиотека rgdal установлена, для этого из меню необходимо выбрать Packages\Set CRAN Mirror..., Packages\Install package и загружена Packages\Load package.... Примечание: если выбор этих пунктов из меню приводит к ошибке Error: unexpected input in "packages\", следует, значит вы выбираете не из меню программы, а пытаетесь ввести эти операции как команды в командную строку. Другим способом загрузки и установки пакета является работа из командной строки:
Чтение растровых данных Введем в командной строке R следующие команды и выберем файл для анализа, например вот этот.
Работа с растровыми данными Для начала, посмотрим на распределение данных:
Построим регрессию и изучим ее результаты:
Для того, чтобы получить из изображения фрейм данных (dataframe в терминологии R) необходимый для некоторых статистический операций можно указать все каналы изображения (в данном случае их 15) следующим образом:
Как видно из примеров, все операции работают с растром напрямую, без обычного преобразования в другие форматы, что удобно. Сохранение данных Самый простой вариант - сохранение данных в текстовый файл с разделителями
Если в результате некоторой операции у вас получилась матрица значений, ее можно сохранить как растр. Такая матрица часто является результатом обработки исходных растровых данных. Если результатом некой операции (над импортированным растром или новыми данными) является вектор (как например при применении деревьев классификации), то его легко сконвертировать в матрицу с нужными размерами, вот так:
Где resultvector - результирующий вектор, height - высота (количество рядов) матрицы, width - ширина (количество колонок). rgdal позволяет читать и сохранять растры со следующей разрядностью согласно спецификации GeoTIFF:
rgdal позволяет читать и сохранять данные в следующих растровых форматах:
Для сохранения данных необходимо инициализировать драйвер соответствующего формата и создать растр с разрядностью данных соответствующей данным в матрице:
Если матрица одна то и канал у результирующего растра будет только один:
После этого остается только назвать результирующий файл и сохранить его:
После окончания работы стоит также удалить использованные объекты: файл с данными и драйвер формата:
Таким образом, данный пример иллюстрирует, как применить всю математическую мощность популярного бесплатного пакета для работы с данными для обработки данных в растровых форматах. Язык R легко позволяет работать со скриптами, что позволяет еще больше увеличить производительность. Если у вас есть свои интересные примеры использования R в анализе пространственных данных, пожалуйста, пишите, мы с удовольствием их опубликуем у нас на сайте. Обсудить в форуме (Комментариев - 0)
Последнее обновление: April 29 2008 |
