python 之 shapely写图

标题如何给定任意坐标利用相关库进行写图？

• 使用的库有shapely，fiona

• 代码如下：

import fiona
from shapely.geometry import mapping,Point,Polygon,LineString
out_dir="F:/03-code/Test/L1.shp"#写图从名字直接在路径中体现，路径要加双引号‘’
coor1=(502411.0023114022, 212402.7633508926)
coor2=(499831.0023114023, 212402.7633508926)
pt1=LineString([coor1,coor2])#这里要用括号括来，不管是(),[];即[coor1,coor2]，(coor1,coor2)
#print(p1.distance(p2))#如果要用到计算面积，距离什么，要用Point(),LineString(),Polygon()转换下
schema={'geometry':'LineString',
	'properties': {'id':'int','name':'str'}}
#schema是一个字典结构，指定了geometry及其它属性结构      
#第一行格式要注意/第一行后的‘，’不要忘记写/标识几何的有'Polygon''LineString','Point'。
#第二行'id','name':代表所写图里面包含这两个图层；后面是设置了图层的具体类型。
with fiona.open (out_dir,'w','ESRI Shapefile',schema=schema) as od:
#fiona.open():打开读写文件;'w'：可写状态； shema：主题格式，要提前设置 ；'ESRI Shapefile':这个注意大小写，固定使用，不能写错。			
	line1=mapping(pt1)#将pt1对象转为GeoJSON格式
 	feature={'geometry':line1,
		'properties':{'id':'1','name':'lin1'}}
	#对直线的设置：'id'图层为1，名字为'lin1'
	od.write(feature)#将上面设置的feature写到图od中

---------------我是分割线------------------------下面看运行结果

图一：该条直线写图处理后在gis软件上的展示图

图2：直线属性展示图