zoukankan      html  css  js  c++  java
  • 基于pyproj地理坐标和投影坐标转换

    from shapely.geometry import Point, LineString
    from pyproj import Transformer
    from pyproj import CRS
    

    class wgs842proj(): def __init__(self, fcrs=None, tcrs=None): if fcrs is None: fcrs = CRS("EPSG:4326") if tcrs is None: tcrs = CRS("EPSG:4527") self.transformer = Transformer.from_crs(crs_4326, crs_proj) def pointProj(self, p): # p Point(113, 34) lng, lat = p.x, p.y y, x = self.transformer.transform(lat, lng) return Point([x, y]) def lineProj(self, line): # line LineString() lngs, lats = line.xy pts = [self.transformer.transform(lat, lng)for lat, lng in zip(lats, lngs)] pts = [_[::-1]for _ in pts] return LineString(pts) if __name__ == '__main__': p0 = [113.0, 34.0] p1 = [113.5, 34.0] ls = LineString([p0, p1]) p0 = Point(p0) p1 = Point(p1) op = wgs842proj() p0_ = op.pointProj(p0) p1_ = op.pointProj(p1) ls_ = op.lineProj(ls)
  • 相关阅读:
    【CodeForces 788B】奇妙的一笔画问题
    数论day2——离散对数、元根
    学习阶段总结(1)
    Flask特殊装饰器
    Flask蓝图
    Flask对象配置
    Flask实例化配置
    Flask路由
    Flask Session
    Flask jinja2
  • 原文地址:https://www.cnblogs.com/ddzhen/p/15341955.html
Copyright © 2011-2022 走看看