Python与开源GIS:使用OGR中拷贝方法创建新的Shapefile

Python与开源GIS:使用OGR中拷贝方法创建新的Shapefile


发布日期: 1970-01-01 更新日期: 2019-10-28 编辑:xuzhiping 浏览次数: 5537

标签:

摘要: Python与开源GIS教程的内容,开发了单独的内容,请打开 https://www.osgeo.cn/pygis/ 查看。 本页面的内容不再更新。 矢量数据集创建也并没有非常难,就是把读取的时候的顺序 重复以便就好了。总体来说这个过程就是构...

Python与开源GIS教程的内容,开发了单独的内容,请打开 https://www.osgeo.cn/pygis/ 查看。
本页面的内容不再更新。



矢量数据集创建也并没有非常难,就是把读取的时候的顺序 重复以便就好了。总体来说这个过程就是构建数据源->构建 层->构建要素->构建形状->关闭数据源。 我们首先来过看一 下如何使用已有的数据来构建我们的数据,也就是复制/粘贴 方法。在OGR的矢量数据模型中,矢量数据分为datasource,layer,与feature三个层次,在这三个层次上,OGR有不同的 复制数据的方法。

在datasource层次创建数据

我们使用世界各国的边界数据作为示例,利用这个数据,来创建 与之内容完全相同的另外一套数据。

    from osgeo import ogr
    import os,math
    inshp = '/gdata/world_borders.shp'
    ds = ogr.Open(inshp)
    driver = ogr.GetDriverByName("ESRI Shapefile")
    outputfile = '/gdata/world_borders_copy.shp'
    if os.access( outputfile, os.F_OK ):
        driver.DeleteDataSource( outputfile )
    pt_cp = driver.CopyDataSource(ds, outputfile)
    pt_cp.Release()

在这里,我们使用了CopyDataSource()的方法来对数据进行复制的 操作。这个函数的第1个参数是要进行复制的数据源,第2个参数是 要生成的数据的路径,并且返回一个指针。为了将数据写入到硬盘,必须释放这个指标,使用Release()方法。

在layer层次拷贝数据

下面我们在layer层次来完成同样的工作。

    from osgeo import ogr
    import os,math
    inshp = '/gdata/world_borders.shp'
    ds = ogr.Open(inshp)
    driver = ogr.GetDriverByName("ESRI Shapefile")
    outputfile = '/gdata/world_borders_copy2.shp'
    if os.access( outputfile, os.F_OK ):
        driver.DeleteDataSource( outputfile )
    layer = ds.GetLayer()
    newds = driver.CreateDataSource(outputfile)
    # layernew = newds.CreateLayer('linescopy',None,ogr.wkbLineString)
    pt_layer  = newds.CopyLayer(layer, 'abcd')
    newds.Destroy()

为了对layer进行拷贝操作,首先得有数据源。在这里首先建立了 newds这一数据源,然后使用newds的CopyLayer()方法来对图层进 行拷贝操作。值得注意的是,当图层拷贝完成之后,需要对newds 使用Destroy()操作,才能将数据写入到磁盘。 CopyLayer()方法的 第1个参数是OGR的Layer对象,第2个参数是要生成的图层的名称。对于Shapefile来说,这个名称是没有用的,但是无论如何,都要 给出一个字符串型的变量。

在feature层次拷贝数据

下面我们来看一下最复杂的在feature层次进行数据拷贝的使用 方法。抽象来看的话,这个方法只是再多进行了一步操作。

    from osgeo import ogr
    import os,math
    inshp = '/gdata/world_borders.shp'
    ds = ogr.Open(inshp)
    driver = ogr.GetDriverByName("ESRI Shapefile")
    outputfile = '/gdata/world_borders_copy3.shp'
    if os.access( outputfile, os.F_OK ):
        driver.DeleteDataSource( outputfile )
    newds = driver.CreateDataSource(outputfile)
    layernew = newds.CreateLayer('worldcopy',None,ogr.wkbLineString)
    layer = ds.GetLayer()
    extent = layer.GetExtent()
    print (extent)
    feature = layer.GetNextFeature()
    while feature is not None:
        layernew.CreateFeature(feature)
        feature = layer.GetNextFeature()
    newds.Destroy()

特别注意:Destroy()不能省略,Destroy()除了销毁数据还担 负了数据flush到磁盘的任务。如果没有这一句,那么刚才创 建的一系列Feature都不会写入磁盘的文件中。

返回目录:Python与开源GIS

关注公众号
获取免费资源

随机推荐


Copyright © Since 2014. 开源地理空间基金会中文分会 吉ICP备05002032号

Powered by TorCMS

OSGeo 中国中心 邮件列表

问题讨论 : 要订阅或者退订列表,请点击 订阅

发言 : 请写信给: osgeo-china@lists.osgeo.org