Extract Population Information.R

#SetWD
setwd("C:/Users/misba/Desktop/PROCESS/")

#Load libraires
library(dplyr)
library(raster)
library(rgdal)
library(data.table)

#Read object
raster_object <- raster(tiff_file)

#Extract XY eastings and northing values and population values
population_df <- as.data.frame(raster_object, xy=TRUE)
names(population_df) <- c("X","Y","Population")

#Remove rows with 0 population to reduce size
population_df <- population_df %>% filter(Population!=0)

#Write to dataframe
fwrite(population_df, "XY_POP_DATA_FULL.csv", row.names=FALSE)