Python根据数据模板创建shapefile的实现-创新互联
废话不多说,我就直接上代码让大家看看吧!
公司专注于为企业提供成都网站制作、成都做网站、微信公众号开发、商城开发,小程序制作,软件定制设计等一站式互联网企业服务。凭借多年丰富的经验,我们会仔细了解各客户的需求而做出多方面的分析、设计、整合,为客户设计出具风格及创意性的商业解决方案,创新互联更提供一系列网站制作和网站推广的服务。#!/usr/bin/env python # -*- coding: utf-8 -*- # @File : copyShapefile.py # @Author: huifer # @Date : 2018-4-28 from os.path import exists import gdal from osgeo import ogr from os import remove gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "YES") # 路径中文 gdal.SetConfigOption("SHAPE_ENCODING", "GBK") # 属性中文 in_shapefile = "dataSample/wang_point.shp"# 数据模板 out_shapefile = "shapefileAa.shp" # 输出数据集 in_ds = ogr.Open(in_shapefile) # 读取模板数据 in_lyr = in_ds.GetLayerByIndex(0) if exists(out_shapefile): remove(out_shapefile) drv = ogr.GetDriverByName("ESRI Shapefile") # 指定数据驱动 out_ds = drv.CreateDataSource(out_shapefile) # 创建数据源 proj = in_lyr.GetSpatialRef() # 获取模板坐标系 out_lyr = out_ds.CreateLayer(out_shapefile.split(".")[0], proj, ogr.wkbPoint) # copy the schema of the original shapefile to the destination shapefile lyr_def = in_lyr.GetLayerDefn() for i in range(lyr_def.GetFieldCount()): # 获取字段长度 out_lyr.CreateField(lyr_def.GetFieldDefn(i)) # 创建字段 feature = ogr.Feature(lyr_def) wkt = "POINT(88615.730000 75345.486000)" point = ogr.CreateGeometryFromWkt(wkt) feature.SetGeometry(point) # 添加点 out_lyr.CreateFeature(feature) # 关闭 特征 feature = None # 关闭数据 data_source = None
另外有需要云服务器可以了解下创新互联scvps.cn,海内外云服务器15元起步,三天无理由+7*72小时售后在线,公司持有idc许可证,提供“云服务器、裸金属服务器、高防服务器、香港服务器、美国服务器、虚拟主机、免备案服务器”等云主机租用服务以及企业上云的综合解决方案,具有“安全稳定、简单易用、服务可用性高、性价比高”等特点与优势,专为企业上云打造定制,能够满足用户丰富、多元化的应用场景需求。
标题名称:Python根据数据模板创建shapefile的实现-创新互联
文章来源:http://myzitong.com/article/ijpjd.html