阿杞

摘要： 方法一： from shapely.geometry import shapefrom shapely import wkt shape() geojson -> geometryshape.wkt geometry-> wktshape.__geo_interface__ geometry-> g 阅读全文

摘要： 前一阵子开始接触PostGis，在做空间数据导入的时候感觉不是很方便，一种方法是postgis自带的一个数据导入的工具shp2pgsql， 但是只能导入shp文件，其他格式的数据基本上都是转成shp格式，然后再导入；另一种是使用ogr的导入命令ogr2ogr，支持多种数据的导入。我们拿到的数据主要是mapinfo的tab数据，可以用ogr2ogr导入，为了其他人的使用以及过程中的控制，自己做了一个tab导入postgis的工具。tab文件的读取使用的是开源的读写工具mitab（mitab使用的ogr library）， 基本流程就是把读到的空间数据转换成wkt的字符串，生成sql语句，然后执行 阅读全文

摘要： 原来使用超图软件进行道路长度的统计，后来决定在项目中使用开源的工具，两点之间距离量算的方法，网上很容易找到，基本上是把地球作为一个标准球体，使用固定半径来进行计算，这样的每一个计算结果与实际长度差距不是很大，一般是可以接受的，但是计算对象较多时统计的结果与超图的计算结果对比就会发现有较大的差距，这是因为超图的长度计算中地球是作为椭球体处理的。 下面提供一种椭球体表面两点间的距离量算，“Vincenty's formulae”，代码非原创，有趣的朋友可以去原作者页面研究一下，源码下载页面http://www.codeproject.com/Articles/19939/GPS-Recei 阅读全文