1992～2024 年各省市区县夜间灯光最亮点的经纬度（余柏蒗版本2）

前不久给大家分享了一份 1992～2024 年各省市区县乡镇类 NPP-VIIRS 夜间灯光面板数据：

1992～2024 年各省市区县乡镇类 NPP-VIIRS 夜间灯光面板数据 & 栅格数据（余柏蒗版本2）：https://rstata.duanshu.com/#/course/5ee752e4c1d84ab48b2235366dfb83ec

然后有小伙伴想要使用各城市夜间灯光最亮处的坐标作为该城市 CBD 的坐标。

于是我就顺便提取了各年、各省市区县夜光最亮处的坐标数据：

• 1992~2024年各城市夜间灯光最亮点的经纬度.xlsx
• 1992~2024年各区县夜间灯光最亮点的经纬度.xlsx
• 1992~2024年各省份夜间灯光最亮点的经纬度.xlsx

数据预览如下：

需要注意的是，有些年份有些区域的最亮处可能不止一处，特别是对于区县数据，有很多区县最亮处有很多个点，因此我只保留了最亮点少于 5 个的。

处理方法

该数据的处理方法不是很复杂，这里简要介绍下：

首先从夜间灯光栅格数据数据中裁剪出某个省/市/区/县的：

regiondf %>%
  slice(1) -> tempsf

rast("rasterdata/cn2006.tif") -> rst

rst %>%
  terra::crop(vect(tempsf)) %>%
  terra::mask(vect(tempsf)) -> temprst

然后小于最大值的都替换成 NA：

temprst[temprst < max(temprst[], na.rm = T)] <- NA
temprst

再把栅格数据转换成矢量数据：

temprst %>%
  raster::raster() %>%
  raster::rasterToPoints(spatial = T) %>%
  st_as_sf()

这样就获取了最亮点的坐标。需要注意有可能一个区域有多个最亮点。

绘图展示

为了便于大家直观感受这个数据，下图展示了 2024 年各省份夜间灯光最亮点的位置：

下图展示了 2024 年各城市夜间灯光最亮点的位置：

数据引用格式

由于该数据包含较多 RStata 处理的内容，在研究中使用该数据请使用清晰的方式注明数据来源于 RStata 或者 RStata 数据中心，并使用如下格式引用：

RStata 数据中心: 1992～2024 年各省市区县夜间灯光最亮点的经纬度（余柏蒗版本2）. 2025. https://tidyfriday.cn/rsdb2/

英文文献可以使用下面的格式引用：

RStata Data Center: Geographic Coordinates of Brightest Nighttime Light Pixels by Province/City/District/County, 1992–2024 (Version 2 by BaiLang YU). 2025. https://tidyfriday.cn/rsdb2/

点击这里跳转到 RStata 短书平台获取附件：1992～2024 年各省市区县夜间灯光最亮点的经纬度（余柏蒗版本2）