WebGIS最佳实践4 使用pgRouting实现最佳路径搜索

tongyi007

　　  下载代码
　　  我在WebGIS最佳实践-2 在WebGIS程序中实现路径分析中提到过pgRouting，现在来看看这头小象能给我们带来些什么。

　　先上效果图给大家鼓鼓劲。：》
　　为了实现以上效果我用jsp和openlayers编写了服务端和客户端，但是本文的重点是介绍pgRouting所以不在这里详细介绍这些这些代码了，大家可以从下载包中找到这两个文件。要使用pgRouting得先安装，我假设大家已经安装了Postgresql和PostGIS。
　　从这里下载合适的二进制包，我下的是Windows版本，解压缩后把Lib下面的动态库复制到Postgresql的bin目录下，把Share\Contrib下面的sql文件复制到Postgresql的share\contrib目录下，然后重新启动Postgresql（这一步我不确定是否必要）。然后就可以进行路径搜素了，下面我们从数据准备开始，一步步得创建出完整可用的路网数据表。Let·s rock！
　　我们先要创建一个新数据库，就叫routing吧，我们会先让它支持pgRouting的函数然后给他增加一张表，一张路网表。随后我们会需要执行一些sql脚本，你可以用postgresql提供的小程序psql.exe来执行也可以用图形化工具pgAdmin来完成。我选择用psql。postgresql提供的程序createdb.exe可以帮助我们创建数据库，我们假设用户是postgres：
createdb -U postgres routing
　　  然后让数据库支持PostGIS和pgRouting的函数和基础表

psql –U postgres –d routing –f postgis.sql
psql –U postgres –d routing –f spatial_ref_sys.sql
psql –U postgres –d routing –f routing_core.sql
psql –U postgres –d routing –f routing_core_wrappers.sql
psql –U postgres –d routing –f routing_topology.sql
　　 其中，&#8220;postgis.sql，spatial_ref_sys.sql&#8221;可以在&#8220;<Postgresql的安装路径>\<version>\share\contrib\postgis<-1.4>&#8221;下面找到，到此数据库创建完成。
　　  为了创建路网表，我们需要一些道路数据。可以从OpenStreetMap上免费获得XML格式的道路数据。由于我们平时更容易获得Shapefile格式的数据，所以我决定用shp来完成数据导入。从Geobase下载道路数据。这个网站提供加拿大全国的道路数据，选择British Columbia , Edition 8.0的Shapefile格式（nrn_rrn_bc_shp_en.zip）。解压缩，我们要使用的数据就在&#8220;NRN_BC_8_0_ROADSEG&#8221;中。这个表的字段很多，大家可以用熟悉的GIS数据处理工具删掉一些字段，当然不处理也没关系。
　　先用shp2psql.exe程序把shp文件装换成sql脚本

shp2pgsql -s 4326 NRN_BC_8_0_ROADSEG public.street_bc > street_bc.sql
　　如果一切顺利的话，我们会得到文件&#8220;street_bc.sql&#8221;，这是一个很大的文件。然后我们执行这个文件：

psql -U postgres -d routing -f street_bc.sql
　　下面创建路网的拓扑结构

psql &#8211;U postgres &#8211;d routing &#8211;c &#8220;alter table street_bc add &#8220;source&#8221; integer&#8221;
psql &#8211;U postgres &#8211;d routing &#8211;c &#8220;alter table street_bc add &#8220;target&#8221; integer&#8221;
psql &#8211;U postgres &#8211;d routing &#8211;c &#8220;select assign_vertex_id(&#8220;street_bc&#8221;,0.001,&#8221;the_geom&#8221;,&#8221;gid&#8221;)&#8221;
psql &#8211;U postgres &#8211;d routing &#8211;c &#8220;alter table street_bc add &#8220;length&#8221; double precision&#8221;
psql &#8211;U postgres &#8211;d routing &#8211;c &#8220;update routing set length=st_length(the_geom)&#8221;
　　到这里路网数据就做好了。来测试一下：

SELECT astext(st_transform(the_geom,900913)),l_stname_c as name FROM shortest_path('SELECT gid as id,source::integer,target::integer,length::double precision as cost FROM street_bc',158390, 63634, false, false),street_bc where edge_id = gid;