Я пытаюсь преобразовать таблицу данных о лесных пожарах в растровые данные. Каждая строка таблицы wildfire является событием пожара и содержит информацию о долготе, широте, году и fire_size.
Я хочу, чтобы растр имел разрешение 0,5 и был способен преобразовать обратно в таблицу, подобную этой: Таблица
#dat3 <- read.csv("USFires3.csv")
coordinates(dat3) <- ~LONGITUDE+LATITUDE
rast <- raster(res=0.5)
projection(rast) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
projection(dat3) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
dat3
dat3ras <- rasterize(dat3, res=0.5)
dat3ras
dat3raspts <- rasterToPoints(dat3ras)
hasValues(dat3ras)
extent(rast) <- extent(dat3)
rast
rast2 <- rasterize(dat3, rast, field=c('year', 'FIPS_CODE'), fun='count')
rast2pts <- rasterToPoints(rast2)
Я использовал приведенный выше код, но Финальная таблица 'rast2pts' не выглядит так, как я хочу.
Я очень плохо знаком с R и пространственными данными или кодированием, поэтому я был бы признателен, если бы вы могли сделать ваше объяснение простым и легким для понимания.
Кроме того, long & lat в таблице огня имеет NAD83, но я хочу перенести его в WGS84.