获取经纬度
library ( sf ) library ( terra ) library ( tidyverse ) library ( readxl ) library ( duckdb ) library ( httr ) # 必须加载,提供GET()函数 library ( jsonlite ) # 解析JSON返回结果 library ( arrow ) # 读取接种单位数据 jzdw <- read_xlsx ( './接种单位列表.xlsx' ) |> select ( 接种单位名称 , 接种单位所属区县名称 ) |> distinct () get_jzdw_location_amap <- function ( jzdw_df , amap_key , max_retry = 3 , base_timeout = 15 ) { # 步骤1:拼接地址 jzdw_df <- jzdw_df |> mutate ( 接种单位所属区县名称 = ifelse ( is.na ( 接种单位所属区县名称 ) | 接种单位所属区县名称 == "" , "" , 接种单位所属区县名称 ) , 完整地址 = paste0 ( "甘肃省" , 接种单位所属区县名称 , 接种单位名称 ) |> str_squish () , 地址编码 = URLencode ( 完整地址 , reserved = TRUE ) , 行号 = row_number () # 添加行号,方便跟踪进度 )...