GIS系统中最短路径查找算法优化之一：双向广度优先搜索

     去年毕业的时候进入一家GIS软件开发公司，正式接手的第一个项目是GIS系统的拓扑分析，在公司原有的算法基础上进行维护和优化。在优化过程有了以下心得，本该去年就更新在blog上的，但是一直没有时间去弄，现更新之。

     在GIS系统中，最常用的一个功能就是在拓扑网络结构中进行最短路径分析。

    

     注：最短路径分析公司取名为FindAPath.

 

     公司前辈最初使用的算法是广度优先算法。因为广搜一旦搜索到解即为最小的步数，这是深搜所忘尘莫及的。

     代码流程大致如下：

     list<long> MiddleNodes; //用于保存待搜索子节点 list<long>::iterator listIter; map<long,long> PreMap; //用于保存搜索路径 map<long,long>::iterator mapIter; set<long> SearchNodes; //已搜索过的节点eid，防止环路，死循环 set<long>::iterator setIter; MiddleNodes.push_back(lBeginNode);//初始化分析起点 PreMap[lBeginNode] = -1;//用于在回溯路径的时候中止条件 while(!MiddleNodes.empty()) { listIter = MiddleNodes.begin(); MiddleNodes.erase(listIter); //取要搜索的，在队列中删除之 SearchNodes.insert(*listIter); //将搜索过的节点记录下来 //获取到待分析点的所有子节点,函数设计业务不列出 list<long> lChildNods = GetChildNodes(*ListIter); list<long>::iterator tempIter(lChildNods.begin()); list<long>::iterator EndIter(lChildNods.end()); while (tempIter != EndIter) { setIter = SearchNodes.find(*tempIter); //防止回环 if (setIter == SearchNodes.end()) { PreMap[*setIter] = *tempIter;//记录每条路径 MiddleNodes.push_back(*setIter);//将子节点插入到待分析队