To model the population density, geospatial data representing natural and built environment has been collected. These variables are land use data from ESA, Global Human Settlement Layer (GHSL), road network from Open Street Map, DEM, and nightlight dataset. After masking water bodies, distance is calculated from the centroids to the 100x100 meters grids to the edge of eight land uses, and road layer. In order to distinguish between built settlement and for GHSL the distances inside the layer has been marked negative and outside the layer marked positive. DEM layer is converted to elevation and slope. Table 1 shows the descriptive statistics of covariates.