потоковый Change Detection

[Boris]

В чем бы быстро без раздумий выполнить потоковое сравнение каналов одинаковой размерности, с двумя дополнительными условиями:
а) игнорировать пиксели с NoData
б) сравнение проводить с приведением гистограмм без учета NoData
дополнительное пожелание:
в) с возможностью использования векторной маски, для задания зон интереса (ROI)

На выходе желательно получить диапазон отрицательных и положительных значений.

[gamm]

если снимки помещаются в памяти, то проще всего - в R

[SS_Rebelious]

если снимки помещаются в памяти, то проще всего - в R

код в студию, пожалуйста

[gamm]

если снимки помещаются в памяти, то проще всего - в R

код в студию, пожалуйста

Не вопрос - cran.r-project.org/

поставьте все 3000 пакетов, и запустите поиск по ключевым словам ...

[SS_Rebelious]

поставьте все 3000 пакетов, и запустите поиск по ключевым словам ...

Ну это вы меня совсем далеко послали Я думал у вас есть готовый пример соответствующего кода...

[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)        

[Voltron]

Если снимки не помещаются в памяти (а R даже жалкую сцену Landsat не осиливает) — OTB.

[Александр Мурый]

Если снимки не помещаются в памяти (а R даже жалкую сцену Landsat не осиливает) — OTB.

Если снимки не помещаются в памяти --- можно попробовать пакет raster, а если это еще и Landsat --- landsat

[Boris]

Я извиняюсь, что я опять о своем, но я же написал "быстро без раздумий".
А

А дальше - обычная статистика, для исходного поста нужно совмещение гистограмм (квантильная регрессия)

особенно,

(квантильная регрессия)

- это за пределами начальной постановки задачи (для меня).

PS
Если трудозатраты на создание скрипта и написание инструкции к нему - "тыкай сюда", "на ввод 'подавай! - результат получай!" имеют указанный порядок, то я готов их компенсировать любым удобным путем.
Поскольку это действительно Landsat'ы. Порезанные специальным образом и приправленные контурами, того, что априори неинтересно или известно. Суть работы в анализе изменений, а большая часть трудозатрат уходит на их выявление и отделение от "шума".

[Максим Дубинин]

http://gis-lab.info/qa/rgdal.html

[gamm]

Если снимки не помещаются в памяти (а R даже жалкую сцену Landsat не осиливает) — OTB.

если на 1Гб, то таки да. На 4Гб уже осиливает. И при нынешних ценах на память нужно постаить 16Гб (на х64, конечно) - тогда все работает.