Bhutan 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. (2015). 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 WorldPop.

Stevens, F. R., Gaughan, A. E., Linard, C., & Tatem, A. J. (2015). Disaggregating Census Data for Population Mapping Using Random Forests with Remotely-Sensed and Ancillary Data. PLOS ONE, 10(2), e0107042. doi:10.1371/journal.pone.0107042

plot of chunk predict_density

Bhutan 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.

Bhutan Census, 2005

Folder: Census
File Name: BTN_adm1_pop.shp
Source: Bhutan Census, 2005
Description: Census numbers obtained from geohive were combined with a shapefile from GADM. More information found in the license file in the ! raw data/census folder.
Class: polygon
Derived Covariates:
area, buff, zones,

class       : SpatialPolygonsDataFrame 
nfeatures   : 20 
extent      : 375742, 710659, 2953288, 3124656  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
nvariables  : 22

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.

[1] "Random Forest model is a merged RF model using models from:"
NPL,   

Call:
 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: 12

          Mean of squared residuals: 0.13
                    % Var explained: 91

plot of chunk random_forestplot of chunk random_forestplot of chunk random_forest

Covariate Metadata

ICIMOD 30M Bhutan Landcover 2010

Folder: Landcover
File Name: lc_r_n_prj455.img
Source: ICIMOD, 2010
Description: Land cover information is provided by ICIMOD, original landcover coding was converted to GlobCover coding using a conversion recorded in the raw data folder.
Class: raster
Derived Covariates:
cls011, dst011, cls040, dst040, cls130, dst130, cls140, dst140, cls150, dst150, cls160, dst160, cls190, dst190, cls200, dst200, cls210, dst210, cls230, dst230, cls240, dst240, cls250, dst250, clsBLT, dstBLT,

class       : RasterBrick 
dimensions  : 1793, 3394, 6085442, 1  (nrow, ncol, ncell, nlayers)
resolution  : 100, 100  (x, y)
extent      : 374897, 714297, 2947832, 3127132  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=90 +k=0.9996 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : C:\Users\aking\Documents\RF.1.1.1\data\BTN\Landcover\Derived\landcover.tif 
names       : landcover 
min values  :        11 
max values  :       210 

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  : 1793, 3394, 6085442, 1  (nrow, ncol, ncell, nlayers)
resolution  : 100, 100  (x, y)
extent      : 374897, 714297, 2947832, 3127132  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=90 +k=0.9996 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : C:\Users\aking\Documents\RF.1.1.1\data\BTN\NPP\Derived\npp.tif 
names       :   npp 
min values  :     0 
max values  : 19151 

plot of chunk covariate_reports


Suomi NPP VIIRS-Derived 2012 Lights at Night, 15 arc-second

Folder: Lights
File Name: DEFAULT: VIIRS 2012
Source: http://ngdc.noaa.gov/eog/viirs/download_viirs_ntl.html
Description: These 'Lights at Night' data were derived from imagery collected by the Suomi National Polar-orbiting Partnership (NPP) Visible Infrared Imaging Radiometer Suite (VIIRS) sensor. Data were collected in 2012 on moonless nights and though background noise associated with fires, gas-flares, volcanoes or aurora have not been removed it represents the best-available data for night-time light production.
Class: raster
Derived Covariates:
,

class       : RasterBrick 
dimensions  : 1793, 3394, 6085442, 1  (nrow, ncol, ncell, nlayers)
resolution  : 100, 100  (x, y)
extent      : 374897, 714297, 2947832, 3127132  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=90 +k=0.9996 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : C:\Users\aking\Documents\RF.1.1.1\data\BTN\Lights\Derived\lights.tif 
names       : lights 
min values  :  -0.19 
max values  :     30 

plot of chunk covariate_reports


WorldClim/BioClim Mean Annual Temperature 1950-2000, 30 arc-second

Folder: Temp
File Name: DEFAULT: BIO1
Source: http://www.worldclim.org/current
Description: WorldClim/BioClim 1950-2000 mean annual precipitation (BIO12) and mean annual temperature (BIO1) estimates (Hijmans et al., 2005) were downloaded, mosaicked and subset to match the extent of our land cover data for the mapping of this region.
Class: raster
Derived Covariates:
,

class       : RasterBrick 
dimensions  : 1961, 3620, 7098820, 1  (nrow, ncol, ncell, nlayers)
resolution  : 100, 100  (x, y)
extent      : 362997, 724997, 2941032, 3137132  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=90 +k=0.9996 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : C:\Users\aking\Documents\RF.1.1.1\data\BTN\Temp\Derived\temp.tif 
names       : temp 
min values  : -117 
max values  :  245 

plot of chunk covariate_reports


WorldClim/BioClim Mean Annual Precipitation 1950-2000, 30 arc-second

Folder: Precip
File Name: DEFAULT: BIO12
Source: http://www.worldclim.org/current
Description: WorldClim/BioClim 1950-2000 mean annual precipitation (BIO12) and mean annual temperature (BIO1) estimates (Hijmans et al., 2005) were downloaded, mosaicked and subset to match the extent of our land cover data for the mapping of this region.
Class: raster
Derived Covariates:
,

class       : RasterBrick 
dimensions  : 1961, 3620, 7098820, 1  (nrow, ncol, ncell, nlayers)
resolution  : 100, 100  (x, y)
extent      : 362997, 724997, 2941032, 3137132  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=90 +k=0.9996 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : C:\Users\aking\Documents\RF.1.1.1\data\BTN\Precip\Derived\precip.tif 
names       : precip 
min values  :    244 
max values  :   4670 

plot of chunk covariate_reports


Roads (OSM), 2015

Folder: Roads
File Name: roads.shp
Source: Open Street Map, Downloaded 2013-09-16, http://extract.bbbike.org/
Description: These data were downloaded as part of a per-country package of data layers made availalble as shapefiles through the http://extract.bbbike.org website, extracted from the Open Street Map (OSM) database.
Class: linear
Derived Covariates:
cls, dst,

class       : SpatialLinesDataFrame 
nfeatures   : 3279 
extent      : 366403, 720305, 2943856, 3121376  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
nvariables  : 7

plot of chunk covariate_reports


Rivers (OSM), 2015

Folder: Rivers
File Name: waterways.shp
Source: Open Street Map, Downloaded 2013-09-16, http://extract.bbbike.org/
Description: These data were downloaded as part of a per-country package of data layers made availalble as shapefiles through the http://extract.bbbike.org website, extracted from the Open Street Map (OSM) database.
Class: linear
Derived Covariates:
cls, dst,

class       : SpatialLinesDataFrame 
nfeatures   : 463 
extent      : 368061, 719767, 2943293, 3126574  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
nvariables  : 4

plot of chunk covariate_reports


Populated Places

Folder: Populated
File Name: DEFAULT: Merged pop/builtupp, pop/builtupa, pop/mispopp
Source: National Geospatial-Intelligence Agency (NGA), http://geoengine.nga.mil/geospatial/SW_TOOLS/NIMAMUSE/webinter/rast_roam.html
Description: The VMAP0 data area downloaded as separate files, grouped roughly by continent, and merged into individual shapefiles for subsetting and further processing for population mapping efforts. These data were obtained directly from the original VMAP0 data sources provided by the NGA and pre-processed using Military Analyst in ArcGIS 10.0. Point data sources are buffered to 100 m and then all polygon data sources are merged to a single shapefile prior to processing.
Class: polygon
Derived Covariates:
cls, dst,

class       : SpatialPolygonsDataFrame 
nfeatures   : 66 
extent      : 386650, 712556, 2948612, 3114893  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
nvariables  : 13

plot of chunk covariate_reports


Inland Waterbodies

Folder: Waterbodies
File Name: DEFAULT: hydro/watrcrsl
Source: National Geospatial-Intelligence Agency (NGA), http://geoengine.nga.mil/geospatial/SW_TOOLS/NIMAMUSE/webinter/rast_roam.html
Description: The VMAP0 data area downloaded as separate files, grouped roughly by continent, and merged into individual shapefiles for subsetting and further processing for population mapping efforts. These data were obtained directly from the original VMAP0 data sources provided by the NGA and pre-processed using Military Analyst in ArcGIS 10.0.
Class: polygon
Derived Covariates:
cls, dst,

class       : SpatialPolygonsDataFrame 
nfeatures   : 23 
extent      : 383420, 643595, 2952597, 3126437  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
nvariables  : 10

plot of chunk covariate_reports


Protected Areas

Folder: Protected
File Name: DEFAULT: WDPAfgdb_Sept2012.gdb
Source: World Database on Protected Areas, Downloaded September, 2012, UNEP, http://www.wdpa.org, http://protectedplanet.net
Description: These data are compiled by UNEP and distributed via the Protected Planet website. All protected areas were downloaded regardless of International Union for Conservation of Nature (IUCN) or any other designation, so they include sanctuaries, national parks, game reserves, World Heritage Sites, etc.
Class: polygon
Derived Covariates:
cls, dst,

class       : SpatialPolygonsDataFrame 
nfeatures   : 16 
extent      : 371942, 719823, 2952093, 3124784  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
nvariables  : 26

plot of chunk covariate_reports


Urban Extents

Folder: Urban
File Name: DEFAULT: schneider-urban.shp
Source: Schneider, et al., United Nations
Description: These data were constructed from MODIS-derived imagery and provided to WorldPop researchers by Schneider, et al. as part of a global urban extents datasets.
Class: polygon
Derived Covariates:
cls, dst,

class       : SpatialPolygonsDataFrame 
nfeatures   : 35 
extent      : 381912, 711369, 2955538, 3115098  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
nvariables  : 2

plot of chunk covariate_reports


Elevation and Derived Slope, 3 second

Folder: Elevation
File Name: ele_btn_mos.tif File Name: ele_btn_mos_prj.tif
Source: HydroSHEDS Void-Filled DEM (Lehnert, et al., 2006), http://hydrosheds.cr.usgs.gov/dataavail.php
Description: The HydroSHEDS data are the result of an effort to provide a globally consistent dataset consisting of NASA's Shuttle Radar Topography Mission (SRTM) data and have been processed, void-filled and corrected for use at large scales.
Class: raster
Derived Covariates:
, slope,

class       : RasterBrick 
dimensions  : 1956, 3607, 7055292, 1  (nrow, ncol, ncell, nlayers)
resolution  : 100, 100  (x, y)
extent      : 363497, 724197, 2941232, 3136832  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=90 +k=0.9996 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0 
data source : C:\Users\aking\Documents\RF.1.1.1\data\BTN\Elevation\Derived\elevation.tif 
names       : elevation 
min values  :        43 
max values  :      7074 

plot of chunk covariate_reports


Digitized Building Locations (OSM), 2015

Folder: Buildings
File Name: buildings.shp
Source: Open Street Map, Downloaded 2013-09-16, http://extract.bbbike.org/
Description: These data were downloaded as part of a per-country package of data layers made availalble as shapefiles through the http://extract.bbbike.org website, extracted from the Open Street Map (OSM) database.
Class: polygon
Derived Covariates:
cls, dst,

class       : SpatialPolygonsDataFrame 
nfeatures   : 3263 
extent      : 390507, 709813, 2953918, 3095290  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
nvariables  : 3

plot of chunk covariate_reports


Points of Interest Locations (OSM), 2015

Folder: Points
File Name: points.shp
Source: Open Street Map, Downloaded 2013-09-16, http://extract.bbbike.org/
Description: These data were downloaded as part of a per-country package of data layers made availalble as shapefiles through the http://extract.bbbike.org website, extracted from the Open Street Map (OSM) database.
Class: point
Derived Covariates:
cls, dst,

class       : SpatialPointsDataFrame 
nfeatures   : 204 
extent      : 387010, 654609, 2957156, 3073146  (xmin, xmax, ymin, ymax)
coord. ref. : NA 
nvariables  : 4

plot of chunk covariate_reports