Cambodia Population Map Metadata Report

Prediction Weighting Layer Used in Population Redistribution

The data presented below represent the predicted number of people per ~100 m pixel as estimated using the random forest (RF) model as described in Stevens, et al. (In Press). The following pages contain a description of the RF model and its covariates, their sources and any metadata collected for each covariate. The prediction weighting layer is used to dasymetrically redistribute the census counts and project counts to match estimated populations based on UN estimates for the final population maps provided by AfriPop, AsiaPop and AmeriPop.

plot of chunk predict_density

Cambodia Census Data and Observed Population Density

These data are the population density values used to estimate the RF model used to create the prediction weighting layer you see above. Values represent population density as measured by people per hectare and calculated from population counts within each census unit. These values are used as the dependent variable during model estimation.

Cambodia Census, 2008, Adm-level 3

Folder: Census
File Name: khm_popdata_2008.shp
Source: National Institute of Statistics, Cambodia, acquired by Gaughan, et al. for use in AsiaPop data products
Description: These census data were acquired for AsiaPop (Gaughan, et al., 2013). Required fields for map production are ADMINID and ADMINPOP.
Class: polygon
Derived Covariates:
area, buff, zones,

class       : SpatialPolygonsDataFrame 
nfeatures   : 1621 
extent      : 211837, 785034, 1144269, 1625281  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
nvariables  : 153

plot of chunk census_data

Random Forest Model and Diagnostics

These output and figures outline the estimated RF model that is used to predict the population density weighting layer. The model is fitted to the population density values for the preceding census data using covariates aggregatedfrom the ancillary data sources summarized following the model diagnostics.

 randomForest(x = x_data, y = y_data, ntree = popfit$ntree, mtry = popfit$mtry,      nodesize = length(y_data)/1000, importance = TRUE, proximity = TRUE) 
               Type of random forest: regression
                     Number of trees: 500
No. of variables tried at each split: 16

          Mean of squared residuals: 0.18
                    % Var explained: 93

plot of chunk random_forestplot of chunk random_forestplot of chunk random_forest

Covariate Metadata

Remotely-sensed, Classified Landcover

Folder: Landcover
File Name: khm_lc_nibb2_built_wat_rurb_maj_8bit.tif
Source: MDA GeoCover Landcover Product, 30m
Description: Landcover from the Landsat-derived MDA GeoCover product, reclassified to match AfriPop coding and eventually broken down into binary classifications by aggregated land cover type (see Linard, et al., 2010 and Gaughan, et al. 2013 for category information).
Class: raster
Derived Covariates:
prp011, cls011, dst011, prp040, cls040, dst040, prp130, cls130, dst130, prp140, cls140, dst140, prp150, cls150, dst150, prp160, cls160, dst160, prp190, cls190, dst190, prp200, cls200, dst200, prp210, cls210, dst210, prp230, cls230, dst230, prp240, cls240, dst240, prp250, cls250, dst250, prpBLT, clsBLT, dstBLT,

class       : RasterBrick 
dimensions  : 5011, 5933, 29730263, 1  (nrow, ncol, ncell, nlayers)
resolution  : 100, 100  (x, y)
extent      : 201798, 795098, 1134265, 1635365  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=48 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : C:\Users\aegaug02\Research\Population\Data\RF\data\KHM\Landcover\Derived\landcover.tif 
names       : landcover 
min values  :         0 
max values  :       240 

plot of chunk covariate_reports

MODIS 17A3 2010 Estimated Net Primary Productivity, 1km

Folder: NPP
File Name: DEFAULT: MODIS 17A3 2010
Source: United States Geological Survey (USGS)
Description: MODIS 17A3 version-55 derived estimates of net primary productivity for the year 2010, estimated for 1km pixel sizes and subset and resampled to match the available land cover and final population map output requirements.
Class: raster
Derived Covariates:

class       : RasterBrick 
dimensions  : 5011, 5933, 29730263, 1  (nrow, ncol, ncell, nlayers)
resolution  : 100, 100  (x, y)
extent      : 201798, 795098, 1134265, 1635365  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=48 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : C:\Users\aegaug02\Research\Population\Data\RF\data\KHM\NPP\Derived\npp.tif 
names       :   npp 
min values  :     0 
max values  : 13574