Более эффективная реклассификация растров
У меня есть raster, где первый столбец нужно переклассифицировать, оставив все остальные значения как NA. Я написал код для этого (ниже), но он очень неэффективен, так как размер растра значительно превышает этот игрушечный пример.
Может ли кто-нибудь порекомендовать более эффективный способ переписать этот код, чтобы он мог обрабатывать гораздо большие растровые файлы? Я подозреваю, что apply() может быть здесь полезен, но я не уверен, в каком направлении это делать.
r<-raster::raster(nrow=10, ncol=10)
r[]<-NA
for(i in 1:r@nrows){
r[(r@nrows*i)-(r@ncols-1)]<-2
}
Стоит ли изучать PHP в 2023-2024 годах?
Привет всем, сегодня я хочу высказать свои соображения по поводу вопроса, который я уже много раз получал в своем сообществе: "Стоит ли изучать PHP в...
Поведение ключевого слова "this" в стрелочной функции в сравнении с нормальной функцией
В JavaScript одним из самых запутанных понятий является поведение ключевого слова "this" в стрелочной и обычной функциях.
Приемы CSS-макетирования - floats и Flexbox
Здравствуйте, друзья-студенты! Готовы совершенствовать свои навыки веб-дизайна? Сегодня в нашем путешествии мы рассмотрим приемы CSS-верстки - в...
Тестирование функциональных ngrx-эффектов в Angular 16 с помощью Jest
В системе управления состояниями ngrx, совместимой с Angular 16, появились функциональные эффекты. Это здорово и делает код определенно легче для...
Концепция локализации и ее применение в приложениях React ⚡️
Локализация - это процесс адаптации приложения к различным языкам и культурным требованиям. Это позволяет пользователям получить опыт, соответствующий...
Пользовательский скаляр GraphQL
Листовые узлы системы типов GraphQL называются скалярами. Достигнув скалярного типа, невозможно спуститься дальше по иерархии типов. Скалярный тип...
0
0
100
1
Перейти к ответу
Данный вопрос помечен как решенный
Ответы 1
Ответ принят как подходящий
Сначала лучший способ делать то, что вы делаете
library(raster)
r <- raster::raster(nrow=10, ncol=10)
for(i in 1:nrow(r)){
r[i, 1] <- 2
}
Гораздо лучшая альтернатива
r <- raster::raster(nrow=10, ncol=10)
r[,1] <- 2
Но для больших растров лучше всего сделать
r <- raster::raster(nrow=10, ncol=10)
xy <- cbind(xFromCol(r, 1), yFromRow(r, 1:nrow(r)))
r <- rasterize(xy, r, field=2)
И, вероятно, менее эффективен, как это
r <- raster::raster(nrow=10, ncol=10)
r <- init(r, "col")
r <- reclassify(r, rbind(c(1,1,2), c(2,Inf,NA)), right=TRUE, include.lowest=TRUE)
Та же идея, лаконично
r <- subs(init(r, "col"), data.frame(from=1, to=2))
Другие вопросы по теме
Добавьте стрелки к графику уровня
Рассчитать максимальную длину последовательных дней выше определенного порога в растровом стеке
Будущее многопроцессорного параллелизма в Windows
Чтение определенных растровых файлов и создание среднего растра в R
Как создать временной ряд (NDVI) для каждого полигона, содержащегося в наборе функций в движке Google Earth?
Почему rasterToPolygons создает горизонтальные линии?
Векторная память исчерпана с помощью функции projectRaster в R
Добавление растровых значений при построении всего набора растров
Затухающий буфер вокруг точек
Как наложить два растра, используя ifelse (условные операторы) в R?
Похожие вопросы
Как остановить автоотключение блестящей приборной панели
Отфильтровать первые n самых больших групп в data.frame
Rstuido Сбой после изменения графических настроек и не открывается снова даже после переустановки
Ggplot2 position_stack(vjust) с абсолютной единицей измерения, а не пропорциональной
Точный критерий Q Кокрана в R
Удалить глобальную переменную и освободить память от функции
Как извлечь столбцы, которые содержат много значений 0, с помощью R?
Python для возврата определенного столбца с использованием понимания списка или применения()
Корреляции между кадром данных и списком кадров данных в R
Создать блестящее боковое меню приборной панели из фрейма данных
Спаси меня еще раз... спасибо Роберту за все возможные подходы к этому вопросу!