标题如何给定任意坐标利用相关库进行写图?
-
使用的库有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:直线属性展示图