Python 데이터 템 플 릿 에 따라 shapefile 을 만 드 는 실현

1677 단어 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
이상 의 Python 은 데이터 템 플 릿 에 따라 shapefile 을 만 드 는 것 이 바로 작은 편집 이 여러분 에 게 공유 하 는 모든 내용 입 니 다.참고 가 되 고 저 희 를 많이 사랑 해 주 셨 으 면 좋 겠 습 니 다.

좋은 웹페이지 즐겨찾기