land_cover.rmd wrong classification

[hanzlan]

Hi,

after trying this tutorial I noticed that IGBP_raster <- projectRaster(IGBP_raster, crs = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs") creates arithmetic mean values with decimal places. This leads to new land cover classes being created which shouldn't be present in Zimbabwe such as class number 15 "Snow and Ice" resulting from the arithmetic mean of classes 13 "Urban and built-up lands" and 17 "Water bodies", when a pixel is situated right between a city and a river for example. IGBP_raster <- projectRaster(IGBP_raster, crs = "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs", method = "ngb") should help here, but maybe more things need to be changed as less classes will be present, then.