Next Article in Journal
Accumulative Errors Optimization for Visual Odometry of ORB-SLAM2 Based on RGB-D Cameras
Previous Article in Journal
Fire Risk Assessment in Dense Urban Areas Using Information Fusion Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Mapping of GDP at 1 km2 Using VIIRS Nighttime Satellite Imagery

1
Department of Geography and the Environment, University of Denver, Denver, CO 80208, USA
2
Morgridge College of Education, University of Denver, Denver, CO 80210, USA
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2019, 8(12), 580; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8120580
Submission received: 5 November 2019 / Revised: 6 December 2019 / Accepted: 7 December 2019 / Published: 11 December 2019

Abstract

:
Frequent and rapid spatially explicit assessment of socioeconomic development is critical for achieving the Sustainable Development Goals (SDGs) at both national and global levels. Over the past decades, scientists have proposed many methods for estimating human activity on the Earth’s surface at various spatiotemporal scales using Defense Meteorological Satellite Program Operational Line System (DMSP-OLS) nighttime light (NTL) data. However, the DMSP-OLS NTL data and the associated processing methods have limited their reliability and applicability for systematic measuring and mapping of socioeconomic development. This study utilized Visible Infrared Imaging Radiometer Suite (VIIRS) NTL and the Isolation Forest machine learning algorithm for more intelligent data processing to capture human activities. We used machine learning and NTL data to map gross domestic product (GDP) at 1 km2. We then used these data products to derive inequality indexes (e.g., Gini coefficients) at nationally aggregate levels. This flexible approach processes the data in an unsupervised manner at various spatial scales. Our assessments show that this method produces accurate subnational GDP data products for mapping and monitoring human development uniformly across the globe.

1. Introduction

The United Nations has established a set of sustainable development goals to achieve a better future for people and the planet. Building on the success of the Millennium Development Goals (MDGs), the 2030 Agenda for Sustainable Development aims to promote and stimulate a series of actions to transform our world. The 17 Sustainable Development Goals (SDGs) with 169 associated targets will unite and mobilize efforts from countries across the world to tackle and address urgent development issues like poverty, inequality, and climate change [1,2,3,4]. Although significant progress has been made towards the achievement of these goals, some of the actions and policies have not been implemented effectively because of the complexity of the Earth system and human–environment interactions. In other words, global climate change is progressing at a quick pace and many people are still living in poverty. Therefore, it is important to understand the global distribution of wealth, characterize socioeconomic well-being, and predict environmental change at appropriate spatiotemporal resolutions to facilitate the implementation of policies and the achievement of SDGs [5].
Measuring socioeconomic data in a timely and accurate manner is important for evaluating current socioeconomic status and assessing policy effectiveness. Doing this well helps countries achieve many of the SDGs including sustainable development, eradication of poverty, and reduction of inequality and exclusion. It also helps practitioners, scientists, and policymakers compare levels of development across the globe to inform efforts toward achieving the SDGs. However, collecting these data can be costly and challenging for many less-developed countries. In recent years, the availability of remotely sensed images has greatly helped scientists monitor human activity on the Earth [6]. For instance, nighttime light (NTL) data are widely used for estimating and evaluating socioeconomic activities since they can capture the artificial light on the Earth’s surface [7,8,9]. Remote sensing technology and satellite imagery have provided us with global and regional economic data to understand and evaluate the relationship between human development and nature [10]. There are many difficulties associated with collecting traditional census data for measuring human well-being. For example, accurate information about the size and distribution of the human population is not available for many regions of the world and sometimes these data are of poor quality [11]. Hence, remote sensing data can be an alternative way for scientists to study and monitor human activities in a timely, consistent, and affordable way. NTL data are different from other remote sensing data as they capture the artificial light on the Earth’s surface and offer a unique view of human activity [9,12,13,14,15,16]. For example, NTL imagery has been used to generate and demonstrate the quantitative relationships between the NTL and population and energy consumption in the USA [17,18].
The Visible Infrared Imaging Radiometer Suite (VIIRS) platform is a new and improved vehicle for developing global NTL data products. Prior to VIIRS, the Defense Meteorological Satellite Program Operational Line System (DMSP-OLS) was primarily designed and developed for cloud cover image detection. Researchers discovered that DMSP-OLS nighttime images of the visible and near-infrared (VNIR) band could help scientists observe and detect the VNIR emission sources (e.g., city lights, auroras, gas flares, and fires). Thus, the DMSP-OLS NTL data have been used in many fields including (1) the measuring of human settlements, (2) urban population and socioeconomic activity, (3) energy and electricity consumption, (4) the monitoring of gas flaring, (5) forest fires, and (6) the impacts of military actions and natural disasters [19]. In recent years, the Visible Infrared Imaging Radiometer Suite (VIIRS) sensor, which is equipped with the Day/Night Band (DNB), has outperformed its predecessor DMSP-OLS in many ways. In general, VIIRS exceeds DMSP-OLS including greater dynamic range, finer spatial resolution, and lower detection limits [20,21].
Over the past decades, scientists have proposed many methods to estimate GDP at national and subnational levels using NTL [19,22,23,24]. For example, Sutton et al. [25] estimated the marketed and non-marketed economic value on a global scale and discovered that the GDP was correlated with the amount of light energy emitted by that country based on DMSP-OLS NTL. Shi et al. [26] used VIIRS NTL to estimate GDP and electricity power consumption and concluded that it can be a strong tool to evaluate socioeconomic indicators. Nevertheless, some researchers found that using NTL alone is insufficient to capture the spatial heterogeneity of GDP at subnational levels. First of all, NTL is not a direct measurement of socioeconomic activities. In addition to that, in many developing regions like Sub-Saharan Africa, a large portion of the population are engaged in agricultural activities. Thus, using NTL alone cannot estimate GDP accurately in these regions [27]. Therefore, some researchers have started to estimate various socioeconomic indicators using NTL based on urban and rural regions separately in order to capture the different levels of productivity [28].
This paper presents an approach that utilizes the VIIRS NTL data for estimating the socioeconomic development metrics for the world. We mainly studied socioeconomic development with a focus on the measurement of inequality as these indicators can reflect the socioeconomic status and the current level of development. Rising levels of economic inequality can lead to a series of consequences including lower rates of economic growth, happiness, and higher rates of crime, health, and poverty problems [29]. Moreover, to improve the current data processing method to achieve better results, we separated GDP estimation based on rural and urban regions using land cover classification data. In addition to that, we also used a machine learning-based data processing method for filtering the NTL outliers that were not related to socioeconomic activities to capture spatial heterogeneity of GDP distribution at the pixel level. We adopted the unsupervised Isolation Forest (iForest) machine learning model to help us automatically detect and remove irrelevant NTL data so as to improve model accuracy [30]. We applied this method to develop two NTL-based indexes including (1) NTL-Gini, which is adapted from the Gini coefficient and developed based on the cumulative share of population and gross domestic product (GDP) and (2) NTL-2020, which is adapted from the 20:20 ratio and used to show how much richer the top 20% of populations are to the bottom 20% based on the cumulative distribution of GDP and population. We produced the two NTL-based indexes to investigate and estimate the current development progress for countries around the world [31].
This paper is organized as follows. Section 2 describes the data and methods we used for developing NTL-based indexes and the model that was developed for evaluating social and economic status. In Section 3 and Section 4, we present and evaluate the results and compare them with the actual data. Finally, we summarize the results and draw conclusions in Section 5.

2. Materials and Methods

2.1. Data Collections

The datasets used in this study are described in Table 1. This study used multisource geospatial data in tandem with aggregate national socioeconomic data to develop NTL-based inequality indexes. We used the stable, cloud-free VIIRS NTL (vcm-orm-ntl) product which is produced by the National Oceanic Atmospheric Administration (NOAA) and the National Aeronautics and Space Administration (NASA). We selected the “vcm-orm-ntl” data to calculate the sum of total NTL intensity in the urban region for each of the administrative units to estimate the development in the urban regions. The NTL can be used to estimate economic activities contributed by commercial and industrial activities. This product improves our estimates of socioeconomic data, since it contains cloud-free average radiance values and the outliers caused by fires and other ephemeral light have been removed. The administrative boundary file was collected from the Database of Global Administrative Areas (GADM) [32]. We mainly used the country (level 0) and subdivision level (levels 1 and 2) data for national and subnational data processing. Additionally, population distribution and settlement data at the 1 km2 level were obtained from the Global Human Settlement (GHS) datasets, which contain location, population, and urban extent information for the human presence on the planet from 1975 to 2015. Human settlement information was obtained from the GHS Settlement Model grid (SMOD), which contains urban center (densely populated areas), urban clusters (towns and suburbs), and rural grid cells (rural areas) that can help us separate urban regions from rural regions. The population data of GHS contain the distribution and density of population and they are expressed as the number of people in each cell. For example, Figure 1 shows the NTL, SMOD, and population data distribution in mainland China and Afghanistan. Socioeconomic statistics were obtained from the World Bank and United Nations Development Programme (UNDP) databases. We used income classification data to group countries and agriculture, forestry, and fishing with value added (% of GDP) for calculating the proportion of GDP in rural regions to estimate the agriculture activity.

2.2. Data Pre-Processing

Results from previous studies show that using NTL alone is insufficient to accurately measure the GDP at subnational levels [27]. In order to assess the different levels of socioeconomic development in the country, we separated the data into urban and rural data to capture the different levels of industrial, commercial, and agricultural activities on a subnational scale using level 0, 1, and 2 administrative districts from GADM. The data preparation consisted of three steps (Figure 2). In step 1, we converted GHS population raster data to points in ArcGIS Pro 2.4.1 and used the converted population data points to sample and extract the values of human settlement information and NTL intensity values based on SMOD data and VIIRS NTL. Therefore, all data points contained attribute information of population, nighttime light value, SMOD, administrative districts information at different levels, income classification (based on country), and unique point identifier. In step 2, we separated population data into urban and rural data based on the SMOD value. In step 3, we applied the iForest method to identify NTL outliers in urban regions that were not related to socioeconomic activities and reclassified them as 0.
Although we selected the “vcm-orm-ntl”, irrelevant light sources that are not filtered by the product processing algorithm still remain due to its sensitivity [28]. Moreover, NTL is not directly measuring socioeconomic activities and contains irrelevant data that can greatly affect the GDP estimates especially at the 1 km2 level [38,39]. Therefore, we adopted a series of measures to remove these irrelevant data in step 3. Since we were only using NTL to measure GDP in the urban regions, we first filtered the data by selecting the data points in the urban regions based on their attributes (SMOD value >20). Then, we used the unsupervised machine learning model of iForest to detect the anomalies for the urban pixels based on the population and NTL attributes of urban data points. Over the past decades, many anomaly detection models have been developed based on classification, clustering, and statistical methods. Some researchers have used the DMSP data as a mask to extract the VIIRS data or build a normal profile for the NTL data in order to remove outliers based on the range of NTL intensity value distribution. However, it is very difficult to apply this method globally due to the variations among countries. Furthermore, some of the methods can only be applied to data of low dimensionality and smaller size. By contrast, the iForest machine learning model can automatically detect irrelevant NTL outliers [30] because it is (1) not reliant on the distribution profile of NTL data, (2) more suitable for processing large datasets, and (3) specifically designed to detect anomalies. Studies have demonstrated that the iForest method can outperform many other existing model-based, distance-based, and density-based methods [30]. It is also more suitable for processing large datasets compared to traditional methods like the density-based spatial clustering of applications with noise (DBSCAN) [40]. We used data points for all countries as input so that iForest could detect outliers based on the population and NTL values’ patterns. For instance, points with extremely high NTL value and low population were identified as anomalies (Figure 3). The identified NTL outlier values were changed to 0. Since we had already used a filtered VIIRS NTL product, only a small proportion (about 0.1%) of NTL outliers were detected.

2.3. GDP and Inequality

The Gini coefficient is a statistical measure of economic inequality based on income distribution [41,42]. A higher Gini coefficient value indicates a higher degree of income inequality, whereas a lower value indicates a lower degree of income inequality. Elvidge et al. [12] developed the night light development index (NLDI) based on the Lorenz curve analysis to analyze the co-distribution of NTL and population by sorting the NTL values in an ascending order. Nevertheless, this index may not be sufficient to capture the distribution of economic activity as it cannot accurately represent the spatiotemporal variation of income distribution since NTL is not a direct measure of income or wealth. Therefore, we combined the NTL values (nanoWatts/cm2/sr) with the actual agricultural production ratios and population distribution in order to improve this characterization of economic activity. We used (1) the urban data points (SMOD value >20) with NTL values to measure the distribution of economic activity in the urban regions and (2) the rural data points (SMOD value <20) with population density to measure the distribution of economic activities in the rural regions. We obtained the aggregate GDP in each district based on the sum of rural and urban GDP (in constant 2011 U.S. dollars) for that district. We defined the urban ( UV j ) and rural pixel value ( RV j ) of GDP as follows:
UV i = SnNTL ( 1 AgRatio ) GDP PopV i TotNTL TotUrPop
RV i = PopV i AgRatio GDP TotRuPop
where i is the unique identification of population pixel (derived from the GHS population layer), SnNTL is the total NTL in the district, TotNTL is the total NTL for the country, PopVi is the population count of the corresponding pixel, TotRuPop and TotUrPop are the total rural and urban population for the country, AgRatio is the proportion of agriculture production of the total GDP, and GDP is the national GDP at purchasing power parity (in 2011 constant U.S. dollars) data obtained from the World Bank database. Based on the procedures described above, we produced a gridded GDP product at the 1 km2 level for countries around the world (Figure 4). The subnational GDP calculation was based on the aggregate NTL and population using level 2 GADM districts. In addition, for countries without AgRatio data, we used the total rural population divided by the total national population to calculate the estimated AgRatio.
NTL-Gini and NTL-2020 ratios were calculated based on the accumulative distribution of aggregate GDP at level 1 and 2 districts for countries around the world. We sorted the aggregate GDP data at district level in an ascending order to construct a GDP distribution profile for each country based on the fraction of population and the cumulative share of GDP at the subnational levels and to plot the Lorenz curve (Figure 5). The NTL-Gini coefficient is equal to the area marked A divided by the sum of the areas marked A and B in Figure 5a (Gini index = Area A/(Area A + Area B)).
Moreover, we also calculated the 20:20 ratio based on the distribution of NTL GDP at district level to measure inequality (Appendix A). Higher 20:20 ratios indicated higher income inequality [43,44,45]. The 20:20 ratio can be more revealing than the Gini coefficient since it compares how much wealthier the top 20% of the population is to the bottom 20% of the population. Many studies have shown that this can be a more useful measure to evaluate other development issues like health and social problems. In order to calculate the 20:20 ratio, we used the same distribution profile for each country and calculated the ratio between the total GDP for the top 20% of the population and total GDP for the bottom 20% of the population.

3. Results

3.1. Subnational GDP Validation

The performance of NTL-based development indexes was evaluated using the actual GDP data from the Organisation for Economic Co-operation and Development (OECD) Regional Statistics and Indicators, the Gini coefficient data from the World Bank databank, and the 20:20 Ratios from UNDP. We first compared and validated the subnational GDP products by using the 249 regional GDP (Large regions TL2) data from OECD administrative units that matched the level 1 districts from GADM. Regional total GDP results of the NTL-based GDP were aggregated using the zonal statistics tool in ArcGIS Pro based on the 1 km2 gridded NTL GDP product. We produced the cross-sectional fit comparing the NTL-based GDP against the actual GDP from OECD regions. In Figure 6, the overall result shows that NTL-based subnational GDP has a high coefficient of determination (R2 = 0.761). In addition, since many researchers have studied the relationship between total NTL values within regions with GDP [26,46] based on simple linear regression, we also produced the cross-sectional fit comparing the sum of NTL within districts against the actual GDP from OECD regions (R2 = 0.684). Results in Figure 6 show that NTL GDP can better reflect the actual GDP values. Nevertheless, due to the small size of validation data (n = 246), it is difficult for us to evaluate the results’ accuracy globally at various spatial scales.
We also compared our data based on the electricity accessibility [47] and the Gridded GDP datasets [48]. Figure 7 shows that because the Gridded GDP datasets only contain GDP per capita information for Uganda at the national level, the GDP data (Figure 7c) is mainly dependent on population density (Figure 7a). Therefore, it fails to show the subnational variation of economic activities in urban and rural regions. Figure 7b shows the electricity access estimation (distribution of people without access to electricity) near Kampala, Uganda [47]. Both the NTL GDP data and the electricity access rate show that the Gridded GDP dataset [48] overestimates GDP in many regions outside Kampala despite the fact that these regions have a low electricity access rate and are less developed. This is possibly because the Gridded GDP product is incapable of differentiating levels of productivity within districts and fails to capture the spatial heterogeneity of GDP at various spatial scales.

3.2. Inequality Validation

To evaluate whether NTL-based GDP distribution can predict inequality accurately, we compared the NTL-based inequality indexes against the actual Gini index and 20:20 ratios data using root mean square error (RMSE) and mean absolute error (MAE) (Figure 8). We normalized all inequality data into the range of 0 to 1. We also compared the data by categorizing the countries based on the income level classification. The overall RMSE and MAE for all countries without using income classification were also compared. The inequality validation results show that there is an overall smaller deviation between the NTL-2020 ratios and the actual data from UNDP, indicating that using NTL and population data can better capture the differences of wealth distribution for the top and bottom 20% of the population in urban and rural regions. Both of the NTL-based Gini coefficient and NTL-2020 have similar RMSE and MAE for high-income and low-income countries, whereas the NTL-2020 ratios have smaller RMSE and MAE for upper-middle and lower-middle countries. This shows that the overall GDP distribution profile may be more accurate for developed and less-developed countries as they tend to rely more on tertiary industry (that can be captured by NTL) and primary industry (captured by population density). For many developing countries (with upper-middle and lower middle incomes), where there tends to be greater socioeconomic inequality, the NTL-2020 ratios can better capture this unequal distribution of income and opportunity. For instance, studies [27] have shown that the correlation between light intensity values and economic activity is much weaker for countries that are dependent on agriculture. This is probable since most of these countries are developed countries or industrialized countries that have advanced their technology infrastructures and developed their economies. Therefore, it is harder to measure the socioeconomic development in these countries.

4. Discussion

The NTL global GDP data at 1 km2 can be aggregated into different subnational levels to support analysis at multiple spatial scales. In general, the NTL data collected by the VIIRS can help us not only monitor the light sources but also study various human activities. By using the multisource data, we developed an NTL-based index to estimate different levels of socioeconomic development in urban and rural regions in an efficient and accurate manner. The VIIRS data were capable of capturing commercial and industrial activities more accurately than DMSP-OLS. Although DMSP-OLS NTL data have been widely used due to their detection of anthropogenic lighting sources to study human activities, they still present many significant problems. For instance, the data have deficiencies such as coarse spatial resolution, saturation, lack of in-flight calibration, and lack of low-light imaging spectral bands suitable for discriminating lighting types [8]. Elvidge et al. [20] compared capabilities of DMSP-OLS and VIIRS and concluded that VIIRS is superior to DMSP-OLS in many ways. Therefore, as more VIIRS products are released (monthly and annual), there is a great potential to capture human development at various spatiotemporal scales. Furthermore, the NTL imagery can potentially become an alternative method for scientists to measure and assess socioeconomic development to achieve SDGs. First, there are many limitations for collecting and calculating traditional inequality data. NTL can help us generate reliable estimates to evaluate if cities have achieved the sustainable goals on a global scale. Second, NTL estimates can be combined with multisource data to help people understand the current water, energy, and food security nexus to evaluate and manage the capacity of our growth. Third, it is important to develop different measures based on various sources of data and evaluate the nation’s sustainable development based on multiple indexes. The current method is also limited by the availability of accurate data for model optimization and validation. For example, in our model, we assumed that the rural region was mainly dependent on agricultural activities. However, there are also other labor-intensive activities that can contribute to the rural GDP like mining, oil extraction, and refinery. Therefore, it is also important to collect more accurate socioeconomic data on various spatial scales to improve the accuracy of GDP estimates.

5. Conclusions

Our approach is suitable for measuring the distribution of GDP on both subnational and national levels. The NTL-based GDP estimation using urban and rural separation helped us capture the spatial heterogeneity of GDP distribution compared to the simple linear regression method based on NTL values only. By utilizing the iForest machine learning solution, it was easier for us to detect outliers from the urban NTL data to better estimate GDP distribution at the pixel level (1 km2). Furthermore, the NTL-based indexes were useful for estimating a variety of inequality indicators. Nevertheless, due to the different levels of development, the performance of NTL-based indexes was also affected. In the future, several options can take this research to another level: (a) incorporating an advanced machine learning model or hybrid model to improve the model performance; (b) collecting more historical socioeconomic data to analyze development changes based on the trend and make forecasts for the future; (c) estimating the inequality at subnational levels and validating the results using ground-truth data; and (d) incorporating more variables to train the model so that it can be customized and adjusted for different inequality evaluation purposes. Moreover, monthly VIIRS NTL products are now also available. There is a great potential for us to understand the dynamics of human population changes within cities, assess our ecological footprints, estimate the demand of resources, and evaluate the limit of our growth [49,50,51].

Author Contributions

Conceptualization, Paul C. Sutton and Xuantong Wang; methodology, all of the authors; validation, all of the authors.

Funding

This research received no external funding.

Acknowledgments

This work was supported in part by the Microsoft AI for Earth grant.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. NTL-Gini and NTL-2020 results.
Table A1. NTL-Gini and NTL-2020 results.
CountryNTL-GiniNTL-2020
American Samoa0.0101.057
Solomon Islands0.0181.100
San Marino0.0351.161
Cyprus0.0391.202
New Caledonia0.0531.291
Belize0.0541.353
Guam0.0781.498
Bermuda0.0791.409
Tonga0.0821.789
Qatar0.1031.580
Spain0.1091.747
Cayman Islands0.1202.619
Virgin Islands (U.S.)0.1211.788
Libya0.1231.854
Trinidad and Tobago0.1231.901
Italy0.1312.010
Greece0.1371.844
Israel0.1411.994
Liechtenstein0.1492.013
Belgium0.1522.268
Saudi Arabia0.1562.141
Singapore0.1582.099
Bosnia and Herzegovina0.1672.466
Finland0.1672.420
Iceland0.1684.659
Malta0.1682.175
Chile0.1712.431
Bahamas, The0.1732.877
Kuwait0.1812.563
Albania0.1892.698
Bahrain0.1902.540
Barbados0.1904.441
Uruguay0.1972.965
Argentina0.1982.878
St. Vincent and the Grenadines0.2023.381
Kyrgyz Republic0.2023.158
Jamaica0.2042.934
France0.2063.050
Hong Kong SAR, China0.2122.718
Sierra Leone0.2142.793
Nepal0.2152.863
Jordan0.2162.728
Japan0.2203.224
Korea, Rep.0.2223.103
Dominican Republic0.2233.647
United Kingdom0.2243.161
Latvia0.2263.747
Armenia0.2272.915
Puerto Rico0.2283.541
New Zealand0.2334.454
Morocco0.2353.296
Serbia0.2394.004
Malaysia0.2453.873
Turkmenistan0.2543.671
Togo0.2543.840
Mongolia0.2564.181
Montenegro0.2595.175
Iran, Islamic Rep.0.2593.299
Czech Republic0.2623.669
Egypt, Arab Rep.0.2623.541
Canada0.2634.361
Bangladesh0.2654.117
Switzerland0.2664.482
Belarus0.2684.785
Ireland0.2694.082
Brazil0.2765.263
Germany0.2774.293
Australia0.2794.875
Peru0.2796.108
Pakistan0.2804.334
Lebanon0.2803.699
Tajikistan0.2824.011
Ecuador0.2885.204
Portugal0.2899.392
Hungary0.2915.646
Tunisia0.2944.615
Turkey0.2964.884
United States0.2976.125
Andorra0.299113.727
Costa Rica0.2995.939
Algeria0.3014.461
Antigua and Barbuda0.30235.927
Luxembourg0.30418.170
Bolivia0.3076.176
North Macedonia0.3084.890
Comoros0.3117.412
South Africa0.3116.693
Colombia0.3117.124
Côte d’Ivoire0.3124.120
China0.3144.654
Uzbekistan0.3145.460
Oman0.3164.982
Venezuela, RB0.3164.780
West Bank and Gaza0.3176.818
Mexico0.3176.849
Mauritius0.3185.428
United Arab Emirates0.3195.260
Cuba0.3205.631
Sweden0.3207.993
Mali0.3215.199
Guyana0.3225.355
Paraguay0.3245.810
Indonesia0.3265.033
Guinea-Bissau0.3265.975
Bulgaria0.33214.342
Georgia0.3357.281
Lesotho0.3367.861
Liberia0.3364.940
Austria0.3377.776
Poland0.3436.219
Isle of Man0.34594.870
Panama0.35317.050
Dominica0.3565.079
Lithuania0.36211.281
Honduras0.3667.189
Iraq0.3679.034
Suriname0.3688.129
Denmark0.3697.777
El Salvador0.3717.590
Ghana0.3716.551
Slovak Republic0.3746.704
Myanmar0.3756.744
India0.3817.201
Syrian Arab Republic0.3839.901
Gambia, The0.3887.404
Ukraine0.3996.813
Russian Federation0.39910.015
Azerbaijan0.4028.961
Croatia0.40927.197
Vanuatu0.4095.782
St. Lucia0.41115.134
Grenada0.4167.682
Nicaragua0.41714.450
Brunei Darussalam0.42319.773
Thailand0.4289.299
Burkina Faso0.4288.286
Benin0.4298.287
Haiti0.43910.196
Fiji0.44615.726
Ethiopia0.4568.087
Congo, Rep.0.45633.783
Cameroon0.45712.810
Senegal0.46011.864
Equatorial Guinea0.46333.041
Angola0.46621.985
Moldova0.46810.032
Djibouti0.46881.186
Niger0.4729.449
Botswana0.47312.043
Rwanda0.4779.348
Norway0.47769.563
Vietnam0.47710.579
Slovenia0.47955.003
Central African Republic0.4809.520
Philippines0.48111.231
Kazakhstan0.48714.438
Madagascar0.4909.834
Sudan0.49917.698
Mauritania0.51114.726
Tanzania0.51413.423
Samoa0.51416.265
Kenya0.51610.632
São Tomé and Príncipe0.52213.227
Romania0.52326.556
Zambia0.52622.972
Guinea0.53214.223
Cambodia0.53412.629
Estonia0.53738.572
Mozambique0.53718.714
Guatemala0.54117.937
Sri Lanka0.54719.815
Gabon0.54817.021
Chad0.55023.489
Netherlands0.55012.488
Zimbabwe0.55334.817
Malawi0.56010.926
Somalia0.57736.880
Nigeria0.59523.318
Afghanistan0.61840.357
Uganda0.62421.898
Seychelles0.62492.991
Lao PDR0.63317.154
Namibia0.63633.687
Eswatini0.63927.596
Burundi0.64241.806
Korea, Dem. Rep.0.64221.546
Bhutan0.64716.960
Eritrea0.68072.123
Congo, Dem. Rep.0.70864.175
Timor-Leste0.73635.984
Cabo Verde0.746166.356
St. Kitts and Nevis0.755351.285
Papua New Guinea0.790651.035
Yemen, Rep.0.804562.975
South Sudan0.831336.357

References

  1. Assembly, G. Sustainable Development Goals. (SDGs), Transforming Our World: The 2030; United Nations: New York, NY, USA, 2015. [Google Scholar]
  2. Griggs, D.; Stafford-Smith, M.; Gaffney, O.; Rockström, J.; Öhman, M.C.; Shyamsundar, P.; Steffen, W.; Glaser, G.; Kanie, N.; Noble, I. Policy: Sustainable development goals for people and planet. Nature 2013, 495, 305. [Google Scholar] [CrossRef] [PubMed]
  3. Robert, K.W.; Parris, T.M.; Leiserowitz, A.A. What is sustainable development? Goals, indicators, values, and practice. Environ. Sci. Policy Sustain. Dev. 2005, 47, 8–21. [Google Scholar] [CrossRef]
  4. Sachs, J.D. From millennium development goals to sustainable development goals. Lancet 2012, 379, 2206–2211. [Google Scholar] [CrossRef]
  5. Andries, A.; Morse, S.; Murphy, R.J.; Lynch, J.; Woolliams, E.R. Seeing Sustainability from Space: Using Earth Observation Data to Populate the UN Sustainable Development Goal Indicators. Sustainability 2019, 11, 5062. [Google Scholar] [CrossRef] [Green Version]
  6. Nordhaus, W.; Azam, Q.; Corderi, D.; Hood, K.; Victor, N.M.; Mohammed, M.; Miltner, A.; Weiss, J. The G-Econ Database on Gridded Output: Methods and Data; Yale University: New Haven, CT, USA, 2006; Volume 6. [Google Scholar]
  7. Zhang, Q.; Seto, K.C. Mapping urbanization dynamics at regional and global scales using multi-temporal DMSP/OLS nighttime light data. Remote Sens. Environ. 2011, 115, 2320–2329. [Google Scholar] [CrossRef]
  8. Baugh, K.; Elvidge, C.D.; Ghosh, T.; Ziskin, D. Development of a 2009 stable lights product using DMSP-OLS data. Proc. Asia Pac. Adv. Netw. 2010, 30, 114–130. [Google Scholar] [CrossRef]
  9. Elvidge, C.D.; Baugh, K.E.; Dietz, J.B.; Bland, T.; Sutton, P.C.; Kroehl, H.W. Radiance calibration of DMSP-OLS low-light imaging data of human settlements. Remote Sens. Environ. 1999, 68, 77–88. [Google Scholar] [CrossRef]
  10. Doll, C.H.; Muller, J.-P.; Elvidge, C.D. Night-time imagery as a tool for global mapping of socioeconomic parameters and greenhouse gas emissions. AMBIO J. Hum. Environ. 2000, 29, 157–163. [Google Scholar] [CrossRef]
  11. Skinner, C. Issues and challenges in census taking. Annu. Rev. Stat. Appl. 2018, 5, 49–63. [Google Scholar] [CrossRef]
  12. Elvidge, C.D.; Baugh, K.E.; Anderson, S.J.; Sutton, P.C.; Ghosh, T. The Night Light Development Index (NLDI): A spatially explicit measure of human development from satellite data. Soc. Geogr. 2012, 7, 23–35. [Google Scholar] [CrossRef]
  13. Elvidge, C.D.; Sutton, P.C.; Ghosh, T.; Tuttle, B.T.; Baugh, K.E.; Bhaduri, B.; Bright, E. A global poverty map derived from satellite data. Comput. Geosci. 2009, 35, 1652–1660. [Google Scholar] [CrossRef]
  14. Ghosh, T.; Anderson, S.J.; Elvidge, C.D.; Sutton, P.C. Using Nighttime Satellite Imagery as a Proxy Measure of Human Well-Being. Sustainability 2013, 5, 4988–5019. [Google Scholar] [CrossRef] [Green Version]
  15. Li, X.; Xu, H.; Chen, X.; Li, C. Potential of NPP-VIIRS nighttime light imagery for modeling the regional economy of China. Remote Sens. 2013, 5, 3057–3081. [Google Scholar] [CrossRef] [Green Version]
  16. Liu, Z.; He, C.; Zhang, Q.; Huang, Q.; Yang, Y. Extracting the dynamics of urban expansion in China using DMSP-OLS nighttime light data from 1992 to 2008. Landsc. Urban Plan. 2012, 106, 62–72. [Google Scholar] [CrossRef]
  17. Welch, R. Monitoring urban population and energy utilization patterns from satellite Data. Remote Sens. Environ. 1980, 9, 1–9. [Google Scholar] [CrossRef]
  18. Welch, R.; Zupko, S. Urbanized area energy-utilization patterns from Dmsp data. Photogramm. Eng. Remote Sens. 1980, 46, 201–207. [Google Scholar]
  19. Huang, Q.; Yang, X.; Gao, B.; Yang, Y.; Zhao, Y. Application of DMSP/OLS nighttime light images: A meta-analysis and a systematic literature review. Remote Sens. 2014, 6, 6844–6866. [Google Scholar] [CrossRef] [Green Version]
  20. Elvidge, C.D.; Baugh, K.E.; Zhizhin, M.; Hsu, F.-C. Why VIIRS data are superior to DMSP for mapping nighttime lights. Proc. Asia Pac. Adv. Netw. 2013, 35, 62–69. [Google Scholar] [CrossRef] [Green Version]
  21. Hillger, D.; Kopp, T.; Lee, T.; Lindsey, D.; Seaman, C.; Miller, S.; Solbrig, J.; Kidder, S.; Bachmeier, S.; Jasmin, T. First-light imagery from Suomi NPP VIIRS. Bull. Am. Meteorol. Soc. 2013, 94, 1019–1029. [Google Scholar] [CrossRef]
  22. Xu, T.; Ma, T.; Zhou, C.; Zhou, Y. Characterizing spatio-temporal dynamics of urbanization in China using time series of DMSP/OLS night light data. Remote Sens. 2014, 6, 7708–7731. [Google Scholar] [CrossRef] [Green Version]
  23. Yu, B.; Shi, K.; Hu, Y.; Huang, C.; Chen, Z.; Wu, J. Poverty Evaluation Using NPP-VIIRS Nighttime Light Composite Data at the County Level in China. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 1217–1229. [Google Scholar] [CrossRef]
  24. Hu, T.; Huang, X. A novel locally adaptive method for modeling the spatiotemporal dynamics of global electric power consumption based on DMSP-OLS nighttime stable light data. Appl. Energy 2019, 240, 778–792. [Google Scholar] [CrossRef]
  25. Sutton, P.C.; Elvidge, C.D.; Ghosh, T. Estimation of gross domestic product at sub-national scales using nighttime satellite imagery. Int. J. Ecol. Econ. Stat. 2007, 8, 5–21. [Google Scholar]
  26. Shi, K.; Yu, B.; Huang, Y.; Hu, Y.; Yin, B.; Chen, Z.; Chen, L.; Wu, J. Evaluating the ability of NPP-VIIRS nighttime light data to estimate the gross domestic product and the electric power consumption of China at multiple scales: A comparison with DMSP-OLS data. Remote Sens. 2014, 6, 1705–1724. [Google Scholar] [CrossRef] [Green Version]
  27. Bundervoet, T.; Maiyo, L.; Sanghi, A. Bright Lights, Big Cities: Measuring National and Subnational Economic Growth in Africa from Outer Space, with an Application to Kenya and Rwanda; The World Bank: Washington, DC, USA, 2015. [Google Scholar]
  28. Townsend, A.C.; Bruce, D.A. The Use of Night-time Lights Satellite Imagery as a Measure of Australia’s Regional Electricity Consumption and Population Distribution. Int. J. Remote Sens. 2010, 31, 4459–4480. [Google Scholar] [CrossRef]
  29. Birdsall, N. Rising inequality in the new global economy. Int. J. Dev. Issues 2006, 5, 1–9. [Google Scholar] [CrossRef] [Green Version]
  30. Liu, F.T.; Ting, K.M.; Zhou, Z.-H. Isolation forest. In Proceedings of the 2008 Eighth IEEE International Conference on Data Mining, Washington, DC, USA, 15–19 December 2008; pp. 413–422. [Google Scholar]
  31. United Nations. Inequality Measurement; UN: New York, NY, USA, 2015. [Google Scholar]
  32. Hijmans, R.; Garcia, N.; Wieczorek, J. Global Administrative Areas Database (GADM) Version 3.6; UN: New York, NY, USA, 2018. [Google Scholar]
  33. Schiavina, M.; Freire, S.; MacManus, K. GHS population grid multitemporal (1975, 1990, 2000, 2015) R2019A. Eur. Comm. JRC 2019. [Google Scholar] [CrossRef]
  34. Pesaresi, M.; Florczyk, A.; Schiavina, M.; Melchiorri, M.; Maffenini, L. GHS settlement grid, updated and refined REGIO model 2014 in application to GHS-BUILT R2018A and GHS-POP R2019A, multitemporal (1975-1990-2000-2015), R2019A. Eur. Comm. JRC 2019. Available online: https://data.europa.eu/euodp/en/data/dataset/42e8be89-54ff-464e-be7b-bf9e64da5218/resource/e8dfcd47-6ff6-4bc1-a1cd-9ffde5405736 (accessed on 1 January 2019).
  35. NOAA Version 1 VIIRS Day/Night Band Nighttime Lights. Available online: https://ngdc.noaa.gov/eog/viirs/download_dnb_composites.html (accessed on 01 August 2019).
  36. The World Bank. World Development Indicators; The World Bank: Washington, DC, USA, 2016. [Google Scholar]
  37. UNDP Income Quintile Ratio. Available online: http://hdr.undp.org/en/content/income-quintile-ratio (accessed on 10 May 2019).
  38. Elvidge, C.D.; Baugh, K.; Zhizhin, M.; Hsu, F.C.; Ghosh, T. VIIRS night-time lights. Int. J. Remote Sens. 2017, 38, 5860–5879. [Google Scholar] [CrossRef]
  39. Wang, R.; Wan, B.; Guo, Q.; Hu, M.; Zhou, S. Mapping regional urban extent using NPP-VIIRS DNB and MODIS NDVI data. Remote Sens. 2017, 9, 862. [Google Scholar] [CrossRef] [Green Version]
  40. Ester, M.; Kriegel, H.-P.; Sander, J.; Xu, X. A density-based algorithm for discovering clusters in large spatial databases with noise. KDD 1996, 96, 226–231. [Google Scholar]
  41. Gini, C. Measurement of inequality of incomes. Econ. J. 1921, 31, 124–126. [Google Scholar] [CrossRef]
  42. Gastwirth, J.L. The estimation of the Lorenz curve and Gini index. Rev. Econ. Stat. 1972, 54, 306–316. [Google Scholar] [CrossRef]
  43. Kawachi, I.; Kennedy, B.P. The relationship of income inequality to mortality: Does the choice of indicator matter? Soc. Sci. Med. 1997, 45, 1121–1127. [Google Scholar] [CrossRef]
  44. Mackenbach, J.P.; Kunst, A.E. Measuring the magnitude of socio-economic inequalities in health: An overview of available measures illustrated with two examples from Europe. Soc. Sci. Med. 1997, 44, 757–771. [Google Scholar] [CrossRef]
  45. De Maio, F.G. Income inequality measures. J. Epidemiol. Community Health 2007, 61, 849–852. [Google Scholar] [CrossRef]
  46. Sutton, P.C.; Costanza, R. Global estimates of market and non-market values derived from nighttime satellite imagery, land cover, and ecosystem service valuation. Ecol. Econ. 2002, 41, 509–527. [Google Scholar] [CrossRef]
  47. Falchetta, G.; Pachauri, S.; Parkinson, S.; Byers, E. A high-resolution gridded dataset to assess electrification in sub-Saharan Africa. Sci. Data 2019, 6, 110. [Google Scholar] [CrossRef]
  48. Kummu, M.; Taka, M.; Guillaume, J.H. Gridded global datasets for gross domestic product and Human Development Index over 1990–2015. Sci. Data 2018, 5, 180004. [Google Scholar] [CrossRef] [Green Version]
  49. Wackernagel, M.; Onisto, L.; Bello, P.; Linares, A.C.; Falfán, I.S.L.; García, J.M.; Guerrero, A.I.S.; Guerrero, M.G.S. National natural capital accounting with the ecological footprint concept. Ecol. Econ. 1999, 29, 375–390. [Google Scholar] [CrossRef]
  50. Costanza, R.; d’Arge, R.; De Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’neill, R.V.; Paruelo, J. The value of the world’s ecosystem services and natural capital. Nature 1997, 387, 253. [Google Scholar] [CrossRef]
  51. Imhoff, M.L.; Bounoua, L.; Zhang, P.; Wolfe, R. Assessing the Urban Heat Island Effect Across Biomes in the Continental USA Using Landsat and MODIS; NASA: Washington, DC, USA, 2011.
Figure 1. (a) NTL for mainland China, (b) GHS SMOD data for mainland China, (c) GHS population data for mainland China, (d) NTL for Afghanistan, (e) GHS SMOD data for Afghanistan, (f) GHS population data for Afghanistan.
Figure 1. (a) NTL for mainland China, (b) GHS SMOD data for mainland China, (c) GHS population data for mainland China, (d) NTL for Afghanistan, (e) GHS SMOD data for Afghanistan, (f) GHS population data for Afghanistan.
Ijgi 08 00580 g001
Figure 2. Procedures for processing population, SMOD, and NTL data.
Figure 2. Procedures for processing population, SMOD, and NTL data.
Ijgi 08 00580 g002
Figure 3. Sample result from iForest outlier detection. Outlier identified in Iran that there is (a) very high NTL value and (b) low population density in the coastal region.
Figure 3. Sample result from iForest outlier detection. Outlier identified in Iran that there is (a) very high NTL value and (b) low population density in the coastal region.
Ijgi 08 00580 g003
Figure 4. Gridded GDP product at 1 km2 level.
Figure 4. Gridded GDP product at 1 km2 level.
Ijgi 08 00580 g004
Figure 5. Lorenz curve for Gini and 20:20 ratios estimation. Calculation of (a) NTL-Gini and (b) NTL-2020 ratios based on the Lorenz curve, and (c) sample Lorenz curve for China based on the cumulative distribution of population and GDP.
Figure 5. Lorenz curve for Gini and 20:20 ratios estimation. Calculation of (a) NTL-Gini and (b) NTL-2020 ratios based on the Lorenz curve, and (c) sample Lorenz curve for China based on the cumulative distribution of population and GDP.
Ijgi 08 00580 g005
Figure 6. (a) Scatterplot of Organisation for Economic Co-operation and Development (OECD) regional GDP and sum of NTL values within districts, and (b) scatterplot of OECD regional GDP and district-level NTL GDP (n = 246 subnational districts).
Figure 6. (a) Scatterplot of Organisation for Economic Co-operation and Development (OECD) regional GDP and sum of NTL values within districts, and (b) scatterplot of OECD regional GDP and district-level NTL GDP (n = 246 subnational districts).
Ijgi 08 00580 g006
Figure 7. Comparisons of population distribution, electricity accessibility [47], Gridded GDP at the 1 km2 level [48], and NTL GDP data around Kampala, Uganda. (a) Population data, (b) electricity access data, (c) GDP data from Gridded global datasets, (d) NTL-based GDP.
Figure 7. Comparisons of population distribution, electricity accessibility [47], Gridded GDP at the 1 km2 level [48], and NTL GDP data around Kampala, Uganda. (a) Population data, (b) electricity access data, (c) GDP data from Gridded global datasets, (d) NTL-based GDP.
Ijgi 08 00580 g007
Figure 8. Validation of NTL-based Gini coefficient and 20:20 ratios based on root mean square error (RMSE) and mean absolute error (MAE) for all countries, low-income countries, lower-middle income countries, upper-middle income countries, and high-income countries.
Figure 8. Validation of NTL-based Gini coefficient and 20:20 ratios based on root mean square error (RMSE) and mean absolute error (MAE) for all countries, low-income countries, lower-middle income countries, upper-middle income countries, and high-income countries.
Ijgi 08 00580 g008
Table 1. Table of data inputs for estimating GDP and inequality indexes.
Table 1. Table of data inputs for estimating GDP and inequality indexes.
DatasetDescriptionSources
PopulationGlobal spatial information for the human presence on the planet in 2015 with 1 km2 spatial resolution.GHS [33]
Human SettlementGlobal spatial information for the human settlement (urban and rural) in 2015 with 1 km2 spatial resolution.GHS [34]
VIIRSVIIRS Cloud Mask-Outlier Removed-Night-Time Lights (vcm-orm-ntl) annual data from 2015 with a spatial resolution of 15 arc-second.NOAA/NASA [35]
Global Administrative AreasGlobal administrative areas of countries including the sub-divisions (v3.6).GADM [32]
Productivity RatiosAgriculture, forestry, and fishing with value added (% of GDP) from 2015 at national level.The World Bank [36]
National GDP 2015 National GDP at purchasing power parity in constant 2011 U.S. dollars. The World Bank [36]
Gini IndexGini index estimates based on household survey data from 2015 at national level.The World Bank [36]
Income Quintile RatioRatio of the average income between the richest 20% and the poorest 20% of the population. Only 2013 income quintile ratios are available.UNDP [37]

Share and Cite

MDPI and ACS Style

Wang, X.; Sutton, P.C.; Qi, B. Global Mapping of GDP at 1 km2 Using VIIRS Nighttime Satellite Imagery. ISPRS Int. J. Geo-Inf. 2019, 8, 580. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8120580

AMA Style

Wang X, Sutton PC, Qi B. Global Mapping of GDP at 1 km2 Using VIIRS Nighttime Satellite Imagery. ISPRS International Journal of Geo-Information. 2019; 8(12):580. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8120580

Chicago/Turabian Style

Wang, Xuantong, Paul C. Sutton, and Bingxin Qi. 2019. "Global Mapping of GDP at 1 km2 Using VIIRS Nighttime Satellite Imagery" ISPRS International Journal of Geo-Information 8, no. 12: 580. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi8120580

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop