2000~2023 年各省市区县夜间灯光最亮点的经纬度

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

2000～2023 年各省市区县乡镇类 NPP-VIIRS 夜间灯光面板数据：https://rstata.duanshu.com/#/brief/course/44c8eb23c9f74676ba7dcd23f886e5b1

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

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

• 2000~2023年各省夜间灯光最亮点的经纬度.xlsx
• 2000~2023年各区县夜间灯光最亮点的经纬度.xlsx
• 2000~2023年各城市夜间灯光最亮点的经纬度.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()

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

绘图展示

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

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

点击这里跳转到 RStata 短书平台获取附件：2000~2023 年各省市区县夜间灯光最亮点的经纬度