2

我查看了 Nominatim 的源代码，看起来它的大部分功能只是准备 SQL 查询并格式化它的结果。

并且可能的其他地理（解码）编码器以相同的方式工作。我可以使用直接 SQL 查询从数据库中获取数据吗？

我在哪里可以找到它们的例子？我想从坐标中获取地名。

4

1 回答 1

3

对的，这是可能的。

几年前我尝试过，想出了一个简单粗暴的地理编码器，你可以在 Github 上找到：https ://github.com/plechi/mini-geocoder

您需要一个安装了 PostGIS 和 HSTORE 扩展的 PostgreSQL-Server。

创建一个新的 Postgres 数据库。

启用 PostGIS 和 HSTORE 扩展：

CREATE EXTENSION postgis;  
CREATE EXTENSION hstore;

之后，您需要Osmosis和 Openstreetmap 的数据集，请查看 OSM-Wiki 了解详细信息：http ://wiki.openstreetmap.org/wiki/Downloading_data

从 Osmosis 创建数据库模式（在终端中运行命令）：

psql -d <database> -f <osmosis-folder>/script/pgsnapshot_schema_0.6.sql

导入数据（也在终端中执行）：

osmosis --read-xml file="<downloaded_datafile>.osm" --write-apidb host="<dbhost>" database="<dbname>" user="<dbuser>" password="<dbpassword>"

从理论上讲，您可以查询生成的数据库（检查架构以获取详细信息）。

对于我的地理编码器，我创建了一个“优化表”以便于查询：

CREATE TABLE geocode_optimized
AS SELECT
     w.tags -> 'addr:street'                            AS street,
     w.tags -> 'addr:housenumber'                       AS housenumber,
     w.tags -> 'addr:postcode'                          AS postcode,
     w.tags -> 'addr:city'                              AS city,
     w.tags -> 'addr:country'                           AS country,
     AVG(ST_X(n.geom))                                  AS longitude,
     AVG(ST_Y(n.geom))                                  AS latitude,
     to_tsvector(concat_ws(' ', w.tags -> 'addr:street',
                           w.tags -> 'addr:housenumber',
                           w.tags -> 'addr:postcode',
                           w.tags -> 'addr:city',
                           w.tags -> 'addr:country'
                 ))                                     AS full_text,
     st_makepoint(AVG(ST_X(n.geom)), AVG(ST_Y(n.geom))) AS point
   FROM ways w
     INNER JOIN way_nodes wn ON w.id = wn.way_id
     INNER JOIN nodes n ON n.id = wn.node_id
   WHERE exist(w.tags, 'addr:housenumber') AND exist(w.tags, 'addr:street')
   GROUP BY housenumber, street, postcode, city, country;

CREATE INDEX idx_geocode_full_text ON geocode_optimized USING GIN (full_text);

地理编码（地址到坐标）：

SELECT 
   street, 
   housenumber, 
   postcode, 
   city, 
   country, 
   longitude, 
   latitude 
FROM geocode_optimized 
WHERE full_text @@ plainto_tsquery('YOUR ADDRESS')

反向地理编码（地址坐标）

SELECT 
    street, 
    housenumber, 
    postcode, 
    city, 
    country, 
    longitude, 
    latitude, 
    CAST (st_distance_sphere(st_makepoint(longitude,latitude), st_makepoint('<longitude>','<latitude>')) AS FLOAT) as distance 
FROM geocode_optimized;

如前所述，这是相当粗糙的，它可能不是最高效的解决方案。

于 2017-03-23T19:04:50.980 回答