使用pyshp创建直线段的shapefile?

时间:2017-08-22 14:33:50

标签: python shapefile pyshp

使用Python和pyshp库,我试图从下面的数据创建一个shapefile(存储在列表中):

edge_list = [
    [-40.5, -20.666],
    [-39.849998, -18.700001],
    [-39.816002, -19.6],
    [-40.071999, -19.391001],
    [-40.150002, -19.933001],
    [-39.733002, -18.533001],
    [-39.833, -18.733],
    [-39.708, -18.419001],
    [-39.370998, -17.891001],
    [-39.200001, -17.417],
    [-39.216999, -17.299999],
    [-39.167, -17.083],
    [-39.049999, -16.433001],
    [-38.932999, -13.967],
    [-39.083, -16.583],
    [-39.0, -13.916],
    [-38.900002, -13.6],
]

以下是我的代码的一部分(其中edge_list是上面的列表):

w = shapefile.Writer()
w.line(parts=[edge_list])
w.field("COMMON_ID", 'C')
w.save("test")

我明白了:

enter image description here

但我希望得到这个:

enter image description here

任何提示?

编辑:这是完整的测试代码,但没有太多内容。文件" temp.csv"只包含上面显示的两列点,用逗号分隔,并为标题(x,y)添加一行。

import csv
import shapefile

data = csv.reader(open("test.csv", "rb"), delimiter = ',')
data.next() # skip header line
edge_list = []
for row in data:
    edge_list.append([float(row[0]), float(row[1])])

for e in range(len(edge_list)):
    print "x=", edge_list[e][0], "y=", edge_list[e][1]

w = shapefile.Writer()
w.line(parts=[edge_list])
w.field("COMMON_ID", 'C')
w.save("test")

2 个答案:

答案 0 :(得分:1)

免责声明:我没有使用shapefile或pyshp。但我知道我的画线方式。

我所看到的是它按照您输入点的顺序绘制线条。它连接点,这是你的点的顺序。您需要做的是重新排序edge_list中的点数。

在您的情况下,如果订购了y变量,您的圆点看起来会很好。

所以,试着替换这一行:

w.line(parts=[edge_list])

用这个:

w.line(parts=sorted(edge_list, key=lambda point: point[1]))

这将按y变量对点进行排序,并应按照您想要的方式绘制线条。

答案 1 :(得分:1)

我通过简单地将列表排序>将其传递给shapefile.line()来获得解决方案。

以下是完整,完整的代码:

import csv
import shapefile

data = csv.reader(open("test.csv", "rb"), delimiter = ',')
data.next() # skip header line
edge_list = []
for row in data:
    edge_list.append([float(row[0]), float(row[1])])

## sort list before passing it to shapefile.line function
edge_list = list(sorted(edge_list, key=lambda point: point[1]))

for e in range(len(edge_list)):
    print "x=", edge_list[e][0], "y=", edge_list[e][1]

w = shapefile.Writer()
w.line(parts=[edge_list]) ## this now works OK!
w.field("COMMON_ID", 'C')
w.save("test")