# 设置driver
driver = ogr.GetDriverByName('ESRI Shapefile')
# 打开输入的矢量
inDs = driver.Open(r'E:/wang2/python/polygon_jingweidu.shp', 0)
if inDs is None:
print("Could not open", 'sites.shp')
sys.exit(1)
# 打开输入矢量的图册
inLayer = inDs.GetLayer()
srs = inLayer.GetSpatialRef()
outshp = r'E:/wang2/python
在Python编程环境中，结合GDAL库，我们可以进行各种地理数据处理任务，包括根据已有矢量数据创建新的矢量数据。这个过程通常涉及到读取、处理和写入地理空间信息。以下将详细阐述如何使用GDAL库来实现这个功能。
我们需要导入必要的库，如`os`用于文件路径操作，`ogr`是GDAL库中的核心模块，用于处理矢量数据：
```python
import os
from osgeo import ogr, gdal
```
在给定的代码段中，`driver`被设置为`ESRI Shapefile`，这是GDAL支持的矢量数据格式之一。`driver.Open()`方法用于打开一个现有的Shapefile，这里的路径是`E:/wang2/python/polygon_jingweidu.shp`。如果文件无法打开，程序会输出错误信息并退出：
```python
driver = ogr.GetDriverByName('ESRI Shapefile')
inDs = driver.Open(r'E:/wang2/python/polygon_jingweidu.shp', 0)
if inDs is None:
print("Could not open", 'sites.shp')
sys.exit(1)
```
接着，`inLayer`获取了输入矢量数据集的图层，`srs`则获取了该图层的空间参考系统（Spatial Reference System）：
```python
inLayer = inDs.GetLayer()
srs = inLayer.GetSpatialRef()
```
然后，我们检查新矢量文件`polygon1.shp`是否已经存在，如果存在，使用`driver.DeleteDataSource()`删除它，以确保不会覆盖已有数据：
```python
outshp = r'E:/wang2/python/polygon1.shp'
if os.path.exists(outshp):
driver.DeleteDataSource(outshp)
```
接下来，我们创建一个新的矢量数据集`outDs`，并创建一个图层`outLayer`。这里，`geom_type=ogr.wkbPolygon`表示我们要创建的图层包含的是多边形几何对象：
```python
outDs = driver.CreateDataSource(outshp)
if outDs is None:
print("Could not create file")
sys.exit(1)
outLayer = outDs.CreateLayer('test', geom_type=ogr.wkbPolygon)
```
为了保持新矢量数据与原矢量数据的属性一致性，我们需要获取原图层的第一个特征（特征ID为0）的属性字段定义，并在新图层中创建相同的字段：
```python
fieldDefn = inLayer.GetFeature(0).GetFieldDefnRef('id')
outLayer.CreateField(fieldDefn)
```
遍历输入图层的所有特征，创建相应的输出特征，设置其几何和属性信息，然后添加到输出图层中：
```python
featureDefn = outLayer.GetLayerDefn()
cnt = 0
inFeature = inLayer.GetNextFeature()
while inFeature:
outFeature = ogr.Feature(featureDefn)
outFeature.SetGeometry(inFeature.GetGeometryRef())
# 如果需要复制属性，可以使用以下代码：
for i in range(inFeature.GetFieldCount()):
outFeature.SetField(i, inFeature.GetField(i))
outLayer.CreateFeature(outFeature)
inFeature.Destroy()
outFeature.Destroy()
cnt = cnt + 1
if cnt < 10:
inFeature = inLayer.GetNextFeature()
else:
break
inDs.Destroy()
outDs.Destroy()
```
这段代码实现了将现有Shapefile中的所有特征（多边形）复制到新的Shapefile中，同时保留了原始文件的属性字段。需要注意的是，如果需要复制所有属性，可以遍历每个特征的字段并将其设置到新特征中。在这个例子中，没有显示这部分代码，但可以通过上面注释的代码实现。
总结来说，通过GDAL库，我们可以方便地处理GIS矢量数据，包括读取、创建、修改和保存。这在进行地理信息系统的开发、数据分析和地图制图时非常有用。在实际应用中，可能还需要处理更多复杂的情况，比如处理不同数据格式、进行空间分析、投影转换等。
