使用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")
我明白了:
但我希望得到这个:
任何提示?
编辑:这是完整的测试代码,但没有太多内容。文件" 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")
答案 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")