收容站

你一定见过Facebook那张著名的友邻图：

图中两个地点之间的连线不是直线，而是测地线（Geodesic），或者说大圆路径、航线等等。其实就是球体上两点之间最短路径投影后的表示。用测地线做图需要有比较多的数据，这样可以表现的比较“繁忙”。而且地理上必须是大范围的，不然看上去还是直线，就缺少了曲线的美感。

在ArcGIS中，也有计算测地线的工具，XY To Line：

它需要一张包含起始点地理坐标的表，就能输出两点之间的测地线了。使用起来很简单，当然真要自己算应该比较麻烦的。

下面举个不恰当的例子：美驻华使领馆与省会城市之间的距离关系。

美国在华有7个使领馆：北京、沈阳、上海、武汉、成都、广州和香港，不考虑香港。

大陆有31个省会城市，能享受“休假式治疗”的起码得是地方大员，处级干部就不考虑了，当然也不包括港澳台。

数据来自国家基础地理信息中心。

使邻馆是点数据，省会城市也是点数据，那么只需计算点和点之间的距离关系就可以了。

具体步骤见代码：

简要说明下，首先为点数据添加坐标属性，之后计算两两之间的距离，然后为距离表连接起始点的坐标属性，最后根据距离表计算测地线。

结果如下（点击可看大图）：

仅供参考。

上周汤在课上讲了三维阶梯图的制作方法，非常受启发。虽然我以前看到过类似的图，但是一直不清楚是怎么制作的。汤的思路是这样的：阶梯图一般用来表示地形，所以处理的是栅格数据。先将栅格数据转成矢量点，即每个栅格用一个点来表示。但是要注意栅格大小，栅格太多要先重采样，否则生成的矢量点太多，效果会不好。然后通过矢量点来创建Voronoi图，因为是矢量点是规则排列的，所以创建好的Voronoi多边形其实都是相同的正方形。之后，将多边形与点的值关联起来。最后，在三维场景中，使用设置图层的“Extrusion”属性，把多边形“挤”出来。

陕西某地的地形表示

这个方法在实验课的PPT上有说，虽然我没去上实验课，但是我看了PPT后，觉得还是有些地方可以改进的：

1. 多边形与点关联的那一步，PPT中使用的是“Attribute Join”。即先对点进行统计，再为多边形添加字段，然后通过“ID”进性关联，最后将点的统计值赋给多边形的新字段。在这里，属性关联不是最好的方法。因为每个点是包含于多边形内的，有很明显的位置关系，使用“Spatial Join”可以更方便的进性关联。
2. 处理创建后Voronoi图范围时，PPT中是将创建Voronoi图中的环境变量“Extent”值，设为输入栅格的范围。由于“Extent”是一个矩形范围，如果处处理的栅格中，有值的部分是矩形，那这样可以满足要求。而如果处理的是非矩形的，例如一个省的边界，那这就不能满足要求了。这个问题解决起来也不麻烦，只需对多边形使用“Clip”操作就行。其中，“Clip Features”参数由输入栅格通过“Raster To Polygon”获得，这样就能很好地除掉无效区域。

我觉得这种方法比较有用，可以用来处理许多数据，所以我按我说的方法做了一个模型。流程图如下，应该比较清晰明了：

这个模型还有很多能改进的地方。比如，并不是每个数据都需要重采样的，所以“Resample Cell Size”应该是一个可选参数。我去试着做过，要将这个参数改为可选，相关的改动要很多，要制作选择结构的脚本以及使用“Merge Branch”工具等。另外，还应添加一个重采样方法的参数，而这个参数也应是可选参数，这就比较麻烦了。当然，还有一些没发现的不足。

另外，ModelBuilder这东西做简单顺序结构的模型还挺方便的，但如果模型稍微复杂点，多一些选择、循环结构，用它实现起来就比较困难了。对于复杂的模型，我觉得还是直接用ArcPy写为好，方便自如多了。那么既然ArcPy可以更快速的解决所有问题，那还有必要掌握ModelBuilder吗？这个问题我也没答案。不过，对于ModelBuilder的高级建模技术，我现在掌握的也不多，如果都会了，或许这些问题就简单了，慢慢学吧。

最后，如果你有兴趣，可以下载这个工具试试，ArcGIS9.3及以前版本的请无视：(

地址：http://goo.gl/XIfzd