Мой вопрос касается методов raster::projectRaster.Согласно документации пакета допускаются две опции: ngb (ближайший сосед) и bilinear (билинейная интерполяция).
В моем текущем приложении я хотел бы преобразовать глобальный растровый слой из 5 минут дугидлинноплановая проекция (+proj=longlat +datum=WGS84) на проекцию Молвейда на равную площадь (+proj=moll +ellps=WGS84):
r_mollweide <- projectRaster(r_longlat, crs="+proj=moll +ellps=WGS84")
Результирующая ячейка сетки в r_mollweide примерно равна (в км2)самые большие клетки в r_longlat.Это означает, что чем дальше мы движемся от экватора, тем больше ячеек r_longlat объединяются в одну ячейку r_mollweide.Что касается меня, метод по умолчанию объединяет эти меньшие ячейки с помощью билинейной интерполяции.Это отлично подходит для температурных сетей.Чем больше r_mollweide, тем больше средняя температура меньших r_longlat пикселей.Однако в некоторых случаях я хочу, чтобы projectRaster вместо этого строил суммы, например, в слоях, которые измеряют общее население, а не температуру.Какой самый эффективный способ сделать это?Можно ли манипулировать projectRaster для использования других функций, кроме двух упомянутых опций?