For now, 1 km squares are sampled in all survey areas. These 1 kilometer squares are abstractions for a village or neighborhood and the aim is to interview 25 people in each. However, there is a need to get a more diverse picture of each 1 km square and hence, they will be sliced up into smaller 100 meter squares that in turn get drawn by PPS. The 100 meter squares are eyeballed by using Google Earth Pro, to check if those are inhabited.
The sampling of the 100 meter square units is done by the function add_100m_units. Here, missing or 0 values in WorldPop are set to 0.0001 to make sure that all 100 meter squares are taking into account that are lying in the sampling areas.
The function preapre_100m_eyeballing reduces the SpatialPolygonsDataFrame to a CSV file. The source files can be loaded from GitHub and the executable exe from here.
lusaka_1k_sample_complete_100m <- add_100m_units(lusaka_1k_sample_complete,lusaka_population)
nairobi_1k_sample_complete_100m <- add_100m_units(nairobi_1k_sample_complete,nairobi_population)
lilongwe_1k_sample_complete_100m <- add_100m_units(lilongwe_1k_sample_complete,lilongwe_population)
malawi_1k_sample_complete_100m <- add_100m_units(malawi_sample_200_complete,malawi_population)
zambia_1k_sample_complete_100m <- add_100m_units(zambia_sample_200_complete,zambia_population)## Warning in fld_names == attr(res, "ofld_nms"): longer object length is not
## a multiple of shorter object length
## Warning in fld_names == attr(res, "ofld_nms"): longer object length is not
## a multiple of shorter object length
## Warning in fld_names == attr(res, "ofld_nms"): longer object length is not
## a multiple of shorter object length
## Warning in fld_names == attr(res, "ofld_nms"): longer object length is not
## a multiple of shorter object length
## Warning in fld_names == attr(res, "ofld_nms"): longer object length is not
## a multiple of shorter object length
## Warning in fld_names == attr(res, "ofld_nms"): longer object length is not
## a multiple of shorter object length
## Warning in fld_names == attr(res, "ofld_nms"): longer object length is not
## a multiple of shorter object length
## Warning in fld_names == attr(res, "ofld_nms"): longer object length is not
## a multiple of shorter object length
## Warning in fld_names == attr(res, "ofld_nms"): longer object length is not
## a multiple of shorter object length
## Warning in fld_names == attr(res, "ofld_nms"): longer object length is not
## a multiple of shorter object length