如何将多个 gpx 文件加载到 PostGIS 中？

Question

我有一堆来自 GPSLogger for Android 应用程序的 gpx 文件。

文件看起来像：

<?xml version="1.0" encoding="UTF-8"?>
<gpx  version="1.0" creator="GPSLogger - http://gpslogger.mendhak.com/"  
  xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"  
  xmlns="http://www.topografix.com/GPX/1/0" 
  xsi:schemaLocation="http://www.topografix.com/GPX/1/0 
    http://www.topografix.com/GPX/1/0/gpx.xsd" >
  <time>2011-08-26T06:25:20Z</time>
  <bounds></bounds>
  <trk>
    <trkseg>
      <trkpt  lat="46.94681501102746"  lon="7.398453755309032" >
        <ele>634.0</ele>
        <speed>0.0</speed>
        <src>gps</src>
        <sat>6</sat>
        <time>2011-08-26T06:25:20Z</time>
      </trkpt>
      <trkpt  lat="46.94758878281887"  lon="7.398622951942811" >
        <ele>748.0</ele>
        <speed>0.0</speed>
        <src>gps</src>
        <sat>5</sat>
        <time>2011-08-26T06:30:56Z</time>
      </trkpt>

...   ...   ...

    </trkseg>
  </trk>
</gpx>

是否可以使用 SQL 或 Python 遍历包含这些文件的目录并将它们加载到一个 PostGIS 表中？

我偶然发现这篇博客文章提到：

我不知道有什么可以直接从 GPX 转换为 PostGIS

这篇文章 给出了使用 SQL 来执行此操作的示例，但我无法理解代码：/

Answer 1

ogr2ogr （GDAL 的一部分）是一个简单直接的 Unix shell 工具，用于将 GPX 文件加载到 PostGIS 中。

ogr2ogr -append -f PostgreSQL PG:dbname=walks walk.gpx

ogr2ogr 在 PostGIS 中使用自己的模式创建自己的数据库表。 tracks 表中每个 GPS 轨迹占一行； tracks.wkb_geometry 包含 GPS 轨迹本身作为 MultiLineString。 track_points 表包含各个位置修正（带有时间戳）。

这是导入之前数据库 walks 的样子：

walks=# \d
               List of relations
 Schema |       Name        | Type  |  Owner   
--------+-------------------+-------+----------
 public | geography_columns | view  | postgres
 public | geometry_columns  | view  | postgres
 public | raster_columns    | view  | postgres
 public | raster_overviews  | view  | postgres
 public | spatial_ref_sys   | table | postgres
(5 rows)

...和导入之后：

walks=# \d
                    List of relations
 Schema |           Name           |   Type   |  Owner   
--------+--------------------------+----------+----------
 public | geography_columns        | view     | postgres
 public | geometry_columns         | view     | postgres
 public | raster_columns           | view     | postgres
 public | raster_overviews         | view     | postgres
 public | route_points             | table    | postgres
 public | route_points_ogc_fid_seq | sequence | postgres
 public | routes                   | table    | postgres
 public | routes_ogc_fid_seq       | sequence | postgres
 public | spatial_ref_sys          | table    | postgres
 public | track_points             | table    | postgres
 public | track_points_ogc_fid_seq | sequence | postgres
 public | tracks                   | table    | postgres
 public | tracks_ogc_fid_seq       | sequence | postgres
 public | waypoints                | table    | postgres
 public | waypoints_ogc_fid_seq    | sequence | postgres
(15 rows)

Answer 2

如果您使用的是 Linux，您可以尝试以下操作：

1. 使用程序将 GPX 转换为 SHP：gpx2shp

   sudo apt-get install gpx2shp
   ...
   gpx2shp -o 输出文件.shp infile.gpx
   

2. 然后使用 shp2pgsql

   将该文件加载到支持 postgis 的数据库中

   sudo apt-get install postgis
   ...
   shp2pgsql 输出文件.shp gis_table
   

您当然可以使用管道并在一个命令行中进行所有操作

有关详细信息，请参阅联机帮助页。

编辑
如果您仍然需要 python 脚本，可以在这里找到帮助
http://pypi.python.org/pypi/gpxtools