Integrating Nighttime Lights with Multisource Geospatial Indicators for County-Level GDP Spatialization: A Geographically Weighted Regression Approach in Mountainous Sichuan, China

Precise, spatially explicit sub-provincial GDP estimates are essential for regional planning, especially in mountainous areas where official economic data remain spatially coarse and unevenly distributed. This study develops a multisource county-level GDP spatialization framework for Sichuan Province, China, integrating corrected NPP/VIIRS nighttime-light (NTL) data with Points of Interest (POIs), land-use structure indicators (proportion of farmland (PFL); proportion of construction land (PCL)), elevation, precipitation, accessibility and population density within a unified indicator system. Two regression approaches—Ordinary Least Squares (OLS) as a global benchmark and Geographically Weighted Regression (GWR) as the spatially adaptive primary model—are calibrated on county-level cross-sectional data for 2020 (n = 183) and evaluated using R2, adjusted R2, AICc and residual spatial diagnostics. The multisource GWR model achieves R2 = 0.882 (adjusted R2 = 0.872, AICc = 5712.26), substantially outperforming both the global OLS benchmark (R2 = 0.801) and NTL-only GWR baseline (R2 = 0.662), confirming that spatial nonstationarity is an intrinsic feature of the GDP–proxy relationship and that integrating complementary geospatial proxies is the primary pathway to improved estimation accuracy in topographically heterogeneous regions. The GWR-based GDP surface exhibits a pronounced basin–plateau contrast: high-value clusters concentrate along the Chengdu Plain and adjacent city corridors, while extensive low-value zones prevail across the western highlands (global Moran’s I = 0.33, Z = 14.26, p < 0.001). Spatially varying GWR coefficients reveal that elevation and precipitation constrain GDP most strongly in high-altitude counties, construction land exerts a consistently positive but spatially graded effect, and the influences of accessibility and population density are context-dependent and locally differentiated. These findings support differentiated territorial development policies: plateau counties require accessibility-first strategies; hill counties benefit from targeted small-city industrialization; and basin cores need managed growth to balance agglomeration advantages against congestion pressures. The framework relies exclusively on globally or nationally available data and is portable to other mountainous regions, though cross-regional validation and extension to multi-year panels using geographically weighted panel regression remain important directions for future work.