zoukankan      html  css  js  c++  java
  • 给SHP文件定义投影

    #!/usr/bin/env python
    # -*- coding: utf-8 -*-
    
    import urllib.request
    import os
    
    def get_epsg_code(epsg):
       """
       Get the ESRI formatted .prj definition
       usage get_epsg_code(4326)
    
       We use the http://spatialreference.org/ref/epsg/4326/esriwkt/
       """
      
       f=urllib.request.urlopen("http://spatialreference.org/ref/epsg/{0}/esriwkt/".format(epsg))
       return (f.read())
    
    # Shapefile filename must equal the new .prj filename
    shp_filename = "UTM_Zone_Boundaries"
    
    # path to where Shapefile with no .prj is located
    os.chdir("../geodata/")
    
    # Here we write out a new .prj file with the same name
    # as our Shapefile named "schools" in this example
    
    with open("%s.prj" % shp_filename, "wb") as prj:
        epsg_code = get_epsg_code(4326)
        print(epsg_code)
        prj.write(epsg_code)
        print("done writing projection definition " + prj.name + "to EPSG:" + str(epsg_code))
    
    # now change abck to our code directory
    os.chdir("../code/")
  • 相关阅读:
    2021.4.11(每周总结)
    2021.4.9
    2021.4.8
    2021.4.7
    k8s集群安装(一)安装方案介绍
    k8s映射外部服务Endpoints
    Linux服务器配置上网代理squid
    SonarQube学习笔记
    yum安装指定版本mysql
    linux目录没有颜色的处理
  • 原文地址:https://www.cnblogs.com/gispathfinder/p/5743682.html
Copyright © 2011-2022 走看看