[原]导入OpenStreetMap海图数据，并在GeoServer上发布

转载请注明作者think8848和出处(http://think8848.cnblogs.com)

上回我们说到如何《在GeoServer中为OpenStreetMap数据设置OSM样式》，文中本来是要发布海图的，因为时间太晚了，所以今天补上吧。

1. 下载OpenStreetMap海图数据

请打开这个网站下载http://openstreetmapdata.com/data/water-polygons OpenStreetMap海图数据(是否称为水图更好点:))

这里面有三个文件，本文将用到第二个文件

因为我们下载的OpenStreetMap的中国数据就是Mercator投影坐标系,SRID为3857，而上文参考的这个例子使用的是WGS84坐标系，SRID为4326，下载第一个文件你还需要自已转一下，所以为了避免麻烦，直接下第二个就好了。

2. 将shp文件导入到PostGis中

先将下载好的数据弄到GeoServer所以服务器，然后解压，然后使用shp2pgsql工具导入到PostGis中，命令如下：

su postgres

shp2pgsql -s 3857 -I -D /tmp/water-polygons-split-3857/water_polygons.shp ocean_all | psql -d chinaosmgisdb -U think8848

这里注意一点，在GeoServer中用了哪个Postgresql用户，这个 -U 参数后就用哪个用户，免得导了数据后没有访问权限

3. 将我们发布的地图范围所在区域(bounds)海图从完整的海图中切出来

psql -U think8848 -d chinaosmgisdb -W

CREATE TABLE ocean AS
WITH bounds AS (
  SELECT ST_SetSRID(ST_Extent(way)::geometry,4326) AS geom
  FROM planet_osm_line
  )
SELECT 1 AS id, ST_Intersection(b.geom, o.geom) AS geom
FROM bounds b, ocean_all o
WHERE ST_Intersects(b.geom, o.geom);