简单使用deno实现反向地理编码

反向地理编码（Reverse Geocoding）是将经纬度等坐标信息转换为人类可读的地址（如城市名、街道名等）的过程。实现反向地理编码有多种办法，可以调用百度、高德、谷歌等服务商提供的第三方API，也可自行去解析地图离线数据实现。如果需求不是要得到精确地址（例如只要得到省、市名），那么采用后者无疑是最优的。

下载省、市边界数据

点击打开：https://datav.aliyun.com/portal/school/atlas/area_selector，在页面右侧，就可以下载省或市的边界数据了：

或者终端里边运行：

# 市级粒度数据
wget -O china_full_city.geo.json https://geo.datav.aliyun.com/areas_v3/bound/100000_full_city.json

最后下载来的 china_full_city.geo.json 数据结构大概长这样：

interface FeatureCollection {
    features: {
        // 各个市级行政区域的轮廓数据
        geometry: {
            coordinates: [float, float][][];
            type: 'MultiPolygon';
        };
        properties: {
            name: string;
            parent: {
                areacode: int;
            };
            // …此处省略其他的一些属性
        };
        type: 'Feature';
    }[];
    type: 'FeatureCollection';
}

它是遵循 GeoJSON 的。所谓“GeoJSON”，就是一种用于编码地理数据的开源格式，专门用于描述几何对象（点、线、多边形等）这类数据的一种格式规范，被广泛应用于地理信息系统（GIS）、地图服务和Web应用。

有了全国所有市的轮廓数据，我们就知道全国每个市在地图上的板块范围。这样，通过遍历每个市的多边形区域，使用turf.js中几何的算法，判断待查询的GPS坐标是否在多边形的范围内，即可得出所在城市。

使用 Turf.js 来检索GeoJSON文件

在 china_full_city.geo.json 同级目录下新建 getCityName.ts 文件，编写TypeScript代码：

import { FeatureCollection, Feature, MultiPolygon } from 'https://esm.sh/v135/@types/geojson@7946.0.14/index.d.ts';

// # 首先是根据这个json文件的结构，可以写出它的类型`XFeatureCollection`：

// properties的属性
type XGeoJsonProperties = {
    acroutes: number[],
    adcode: number,
    center: number[],
    centroid: number[],
    childrenNum: 0,
    level: string,
    name: string, // 市级行政区名称
    parent: {adcode: number}, // 市级行政区所在的省级行政区编码
    subFeatureIndex: number,
}

type XFeatureCollection = FeatureCollection<MultiPolygon, XGeoJsonProperties>;

// # 然后遍历features列表的各个城市，使用turf.js封装的`booleanPointInPolygon`方法，
//   判断待查询经纬度的点是不是在市级行政区图块多边形之中，是：则表示坐标在该市范围之内

import chinaProvincesGeoJSON from "./china_full_city.geo.json" with {type: 'json'};
import { booleanPointInPolygon } from "https://esm.sh/@turf/boolean-point-in-polygon@7.1.0";
import * as turf from "https://esm.sh/@turf/turf@7.1.0";

// 定义一个待查询的国内经纬度坐标，形式是“度分秒”
const [[d0, m0, s0], [d1, m1, s1]] = [[112, 34, 56], [34, 56, 12]];
// 转化为“度”的坐标形式
const [lng, lat] = [d0+m0/60+s0/3600, d1+m1/60+s1/3600];
// 定义为Point坐标类型
const pt = turf.point([lng, lat]);

// 定义类型
const collections = chinaProvincesGeoJSON as XFeatureCollection;

let targetFeature!: Feature<MultiPolygon, XGeoJsonProperties>;

for (const f of collections.features) {
    // 使用turf的这个函数判断这个点是不是在（市级行政区）多边形区域内
    const ret = booleanPointInPolygon(pt, turf.multiPolygon(f.geometry.coordinates));
    if (ret) {
        targetFeature = f;
        break;
    }
}

console.log(targetFeature ? targetFeature.properties.name : '（超出查询范围）');

然后终端运行：

# 返回：焦作市
deno run -A getCityName.ts

如果想要知道城市所在的省级行政区名，可以通过查询父级行政区域码 targetFeature.properties.parent.adcode 得知：

const provinceAdcodeNameMapping: {[k: number|string]: string} = {
    "110000":"北京市","120000":"天津市","130000":"河北省","140000":"山西省","150000":"内蒙古自治区",
    "210000":"辽宁省","220000":"吉林省","230000":"黑龙江省",
    "310000":"上海市","320000":"江苏省","330000":"浙江省","340000":"安徽省","350000":"福建省","360000":"江西省","370000":"山东省",
    "410000":"河南省","420000":"湖北省","430000":"湖南省","440000":"广东省","450000":"广西壮族自治区","460000":"海南省",
    "500000":"重庆市","510000":"四川省","520000":"贵州省","530000":"云南省","540000":"西藏自治区",
    "610000":"陕西省","620000":"甘肃省","630000":"青海省","640000":"宁夏回族自治区","650000":"新疆维吾尔自治区",
    "710000":"台湾省",
    "810000":"香港特别行政区","820000":"澳门特别行政区",
    "100000_JD":""
};

// 返回：河南省
console.log(targetFeature ? provinceAdcodeNameMapping[targetFeature.properties.parent.adcode] : '');

最后

这种办法适合对解析精度要求不太高的场景，如果想要解析到区和街道或者更高级别，考虑到未来数据更新和维护的复杂度，建议还是使用各大第三方平台的API比较好。