def convert_to_latlon(geoimg, lines):
srs = osr.SpatialReference(geoimg.srs()).ExportToProj4()
projin = Proj(srs)
projout = Proj(init='epsg:4326')
newlines = []
for line in lines:
l = []
for point in line:
pt = transform(projin, projout, point[0], point[1])
l.append(pt)
newlines.append(l)
return antimeridian_linesplit(newlines)
评论列表
文章目录