基于osgeo创建渔网（不使用arcpy）

John H.

已于 2022-09-29 17:16:14 修改

阅读量476

点赞数 1

分类专栏： python OSGEO Gdal 文章标签： python

于 2022-05-03 11:58:22 首次发布

本文链接：https://blog.csdn.net/weixin_42778485/article/details/124553191

版权

本文介绍如何利用Python的osgeo模块，而非arcpy，来创建渔网图形。通过实例代码详细讲解了整个过程，适合地理信息系统（GIS）开发者参考。

摘要由CSDN通过智能技术生成

def shapefile_fishnet_creat(
        _input_boundary_shapefile, _output_grid_fishnet,
        flag=1, _rows=5, _cols=5,
        _gridHeight=100, _gridWidth=100):
    """
    根据矢量文件边界创建渔网
    _input_boundary_shapefile需要创建渔网的矢量边界，坐标系一定要是投影坐标系
    _output_grid_fishnet渔网单元格矢量文件目标绝对路径
    flag是制作渔网标志，如果flag=1是按照行列数做渔网，就是最后生成的渔网有几行几列
                       如果flag=0是按照渔网单元格高度、宽度生成渔网的
    rows生成的渔网的行数
    cols生成的渔网的列数
    _gridHeight渔网单元格高度，单位：米
    _gridWidth渔网单元格宽度，单位：米
    """
    _ds = ogr.Open(_input_boundary_shapefile)
    _lyr = _ds.GetLayer()
    x_min, x_max, y_min, y_max = _lyr.GetExtent()

    # get rows and columns
    if flag == 1:
        # 按照行列数生产渔网
        _gridHeight = (x_max - x_min)/_rows
        _gridWidth = (y_max - y_min)/_cols

        rows = math.ceil((y_max-y_min)/_gridHeight)
        cols = math.cei