Страница 1 из 1
потоковый Change Detection
Добавлено: 30 авг 2011, 02:53
Boris
В чем бы быстро без раздумий выполнить потоковое сравнение каналов одинаковой размерности, с двумя дополнительными условиями:
а) игнорировать пиксели с NoData
б) сравнение проводить с приведением гистограмм без учета NoData
дополнительное пожелание:
в) с возможностью использования векторной маски, для задания зон интереса (ROI)
На выходе желательно получить диапазон отрицательных и положительных значений.
Re: потоковый Change Detection
Добавлено: 30 авг 2011, 05:31
gamm
если снимки помещаются в памяти, то проще всего - в R
Re: потоковый Change Detection
Добавлено: 30 авг 2011, 09:20
SS_Rebelious
gamm писал(а):если снимки помещаются в памяти, то проще всего - в R
код в студию, пожалуйста
Re: потоковый Change Detection
Добавлено: 30 авг 2011, 10:44
gamm
SS_Rebelious писал(а):gamm писал(а):если снимки помещаются в памяти, то проще всего - в R
код в студию, пожалуйста
Не вопрос -
cran.r-project.org/
поставьте все 3000 пакетов, и запустите поиск по ключевым словам ...
Re: потоковый Change Detection
Добавлено: 30 авг 2011, 10:58
SS_Rebelious
gamm писал(а):поставьте все 3000 пакетов, и запустите поиск по ключевым словам ...
Ну это вы меня совсем далеко послали
Я думал у вас есть готовый пример соответствующего кода...
Re: потоковый Change Detection
Добавлено: 30 авг 2011, 11:27
gamm
SS_Rebelious писал(а):gamm писал(а):поставьте все 3000 пакетов, и запустите поиск по ключевым словам ...
Ну это вы меня совсем далеко послали
Я думал у вас есть готовый пример соответствующего кода...
Код - в хелпе к rgdal, там несколько строк. После прочтения снимков через rgdal они превращаются в обычные таблицы, а point_in_poly() позволяет создавать/использовать векторные маски. А дальше - обычная статистика, для исходного поста нужно совмещение гистограмм (квантильная регрессия). поэтому для написания соотвествующего кода нужно примерно 30 минут
Вот кусок, который читает растр:
Код: Выделить всё
fn.region<-"L7.TIF"
work.region.tif <- new("GDALReadOnlyDataset", fn.region)
p.b<-getRasterTable(work.region.tif)
Re: потоковый Change Detection
Добавлено: 30 авг 2011, 14:11
Voltron
Если снимки не помещаются в памяти (а R даже жалкую сцену Landsat не осиливает) — OTB.
Re: потоковый Change Detection
Добавлено: 30 авг 2011, 14:25
Александр Мурый
Voltron писал(а):Если снимки не помещаются в памяти (а R даже жалкую сцену Landsat не осиливает) — OTB.
Если снимки не помещаются в памяти --- можно попробовать пакет
raster, а если это еще и Landsat ---
landsat
Re: потоковый Change Detection
Добавлено: 30 авг 2011, 19:31
Boris
Я извиняюсь, что я опять о своем, но я же написал "быстро без раздумий".
А
А дальше - обычная статистика, для исходного поста нужно совмещение гистограмм (квантильная регрессия)
особенно,
(квантильная регрессия)
- это за пределами начальной постановки задачи (для меня).
PS
Если трудозатраты на создание скрипта и написание инструкции к нему - "тыкай сюда", "на ввод 'подавай! - результат получай!" имеют указанный порядок, то я готов их компенсировать любым удобным путем.
Поскольку это действительно Landsat'ы. Порезанные специальным образом и приправленные контурами, того, что априори неинтересно или известно. Суть работы в анализе изменений, а большая часть трудозатрат уходит на их выявление и отделение от "шума".
Re: потоковый Change Detection
Добавлено: 30 авг 2011, 23:46
Максим Дубинин
Re: потоковый Change Detection
Добавлено: 31 авг 2011, 17:29
gamm
Voltron писал(а):Если снимки не помещаются в памяти (а R даже жалкую сцену Landsat не осиливает) — OTB.
если на 1Гб, то таки да. На 4Гб уже осиливает. И при нынешних ценах на память нужно постаить 16Гб (на х64, конечно) - тогда все работает.