1. Introduction
China’s economic growth has drawn much attention over recent decades [
1,
2]. There are a number of studies focusing on China’s economic development. Economic statistics, especially long-term time series data, are the basis for various economic studies [
3,
4,
5]. However, incompleteness regarding GDP statistics in both spatial and temporal domains still exists [
6], especially on smaller scales (e.g., at the county level, including county cities and districts), and in the years before 2000. Economic statistics are not complete in both underdeveloped and developed provinces in China. Therefore, obtaining a spatially complete and temporally continuous indicator is crucial for analyzing Chinese economic phenomena.
Remotely sensed data provide a promising option to address the aforementioned shortcomings. With the development of technology, remotely sensed data has gradually become one of the most reliable data sources in the field of economics. The nighttime light data acquired from the Operational Linescan System (OLS) sensor of the Defense Meteorological Satellite Program (DMSP) is one of the most important representatives and has been widely used in socioeconomic studies [
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23]. Because of the unique optical amplification capability, these data can detect weak visible and near-infrared radiation at night, eliminate cloud reflection of the moonlight, and enhance the light generated by human activities, such as those involved in urban and industrial enterprises [
9]. In particular, the DMSP/OLS data have been accumulating since 1992, and so represent a solution to the lack of economic data in China over that time period.
Extracting statistical metrics from nighttime light data is widely used to assess economic development on local, regional, and global scales [
7,
8,
9,
10]. These statistical metrics, such as the sum or mean values of lighting data extracted by administrative divisions, are used to establish a regression relationship with economic data, such as GDP [
11]. Previous studies have attempted to quantify the correlation between nighttime light and China’s economic statistical data [
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22]. However, there are still several limitations. First and foremost, previous studies have not paid enough attention to the saturation problem of OLS lighting data in China. The saturation problem is that the DMSP/OLS stable light data range from 0 to 63, consequently, any areas with a brightness greater than 63 are only expressed as 63. Part of the reason for this is that the impact of saturation on provincial and higher-level data, which is the most commonly used by most researchers, is small [
4,
16,
23]. However, the saturation issue cannot be overlooked when applied on a smaller scale, such as at the county-level [
15]. With China’s economic development, the prosperity of certain cities is close to that of developed countries—these include Beijing, Shanghai, Guangzhou, and Shenzhen (
Appendix Figure A1). Saturation of nighttime light data can bias the results when assessing regional economic development, especially in small units, leading to unreliable economical conclusions for policy makers. Secondly, most previous studies focus on the provincial and prefectural levels [
12,
13,
14,
16,
17,
18,
19,
20], with only a few studies directed at finer spatial level data in China [
15]. Finally, most studies only use nighttime light data from a certain year to compare with GDP [
13,
14,
15] and pay little attention to long-term sequence data, which is crucial for studying China’s long-term sustained economic growth. Several studies have pointed out that the data from the Visible Infrared Imaging Radiometer Suite (VIIRS), carried by the Suomi National Polar-Orbiting Partnership (NPP) satellite, are more accurate than DMSP/OLS data in terms of modeling economic indicators [
13,
14,
15]. However, NPP-VIIRS began collecting data in early 2013, thus limiting the time span that can be studied, which is an important consideration in economic studies.
There are many existing methods to correct the NSL (nighttime stable lights) saturation problem, which generally fall into three categories. The first method involves utilizing other satellite data. For example, Zhang et al. [
24] proposed to use the Normalized Vegetation Index (NDVI) as an auxiliary dataset to remove the saturation phenomenon. Zhuo et al. [
25] suggested removing saturation using Enhanced Vegetation Index (EVI) data. Jing et al. [
26] used 1 km
2 GDP grid data to calibrate DMSP/OLS NSL data. The second method uses nighttime light calibration data. Letu et al. [
27] used 1996–1997 DMSP/OLS radiation calibration data to eliminate the saturation of NSL in 1999. Wu et al. [
28] used 2006 radiation calibration data and a power function regression method to perform saturation correction on global multi-period NSL data. Using the same calibration data and regression model, Cao et al. [
17] corrected 33 years of NSL data in China using Hegang in Heilongjiang Province as the calibration area. The final method uses statistical censoring approaches. Bluhm and Krause [
29] used a Pareto Distribution to address the saturation problem at city level. While these methods corrected the NSL data successfully for certain applications, there are some shortcomings. The use of other satellite data for saturation correction requires other remotely sensed data. In addition, long-term time series of image calibration need to match the reference data each year. The availability of reference data has become an important factor affecting the accuracy of saturation correction [
24]. In contrast, the method of saturation correction using DMSP/OLS radiation calibration data can remove the saturation and does not depend on other auxiliary data. However, these methods, which use only one year of non-saturated nighttime light data, assume that the intensity in saturated areas did not change over a long period of time, which is not realistic in developing countries, such as India and China [
24]. Statistical censoring approaches normally require several assumptions to be made with regards to cities and populations.
To address these limitations, we used GDP growth rate to correct the saturated pixels and so to solve the saturation problem of the nighttime light data. To verify the reliability of the new method, we compared statistics of corrected NSL data generated by our method with two other widely used methods (i.e., the second-order function [
30,
31] and power function [
17,
28]). We then further modeled the relationship between economic statistics and the calibrated nighttime light statistics at the county level over the whole of China, and estimated national county-level GDP using corrected NSL data from 1992 to 2013 and compared this with authoritative GDP statistics. In doing so, we attempted to provide a robust method to address the saturation problem in DMSP/OLS data, and subsequently improve the GDP estimation at county level, thus enhancing understanding in this field of study.