在Python项目中，如何利用GDAL库进行栅格和矢量数据的读取、写入操作，并展示数据结构？

为了更好地理解和掌握GDAL在Python项目中的应用，特别是在处理栅格和矢量数据方面，本回答将提供一个详细的指南，内容涵盖读取、写入操作以及数据结构的展示。首先，确保你已经安装了GDAL库及其Python绑定。

参考资源链接：[GDAL-Python库详解：操作地理空间数据的利器](https://wenku.csdn.net/doc/mbr98i2c3e?spm=1055.2569.3001.10343)

开始之前，你需要导入必要的模块。对于栅格数据，你可以使用`gdal`模块，而对于矢量数据，则使用`ogr`模块。

    from osgeo import gdal
    from osgeo import ogr

接下来，我们来读取栅格数据。GDAL提供了`gdal.Open()`函数来打开栅格文件。例如，读取一个TIFF文件可以这样做：

    raster_dataset = gdal.Open('path_to_your_raster_file.tif')
    if raster_dataset is None:
        raise RuntimeError('无法打开栅格数据')

`raster_dataset`变量现在是栅格数据集的句柄，你可以从中获取栅格数据的结构信息，例如行列数、波段数和地理变换信息：

    print('行列数:', raster_dataset.RasterXSize, 'x', raster_dataset.RasterYSize)
    print('波段数:', raster_dataset.RasterCount)

对于矢量数据，使用`ogr`模块打开矢量文件：

    vector_dataset = ogr.Open('path_to_your_vector_file.shp')
    if vector_dataset is None:
        raise RuntimeError('无法打开矢量数据')

`vector_dataset`提供了对矢量数据集中图层的访问。你可以遍历图层来了解其内容：

    for layer in vector_dataset:
        print('图层名:', layer.GetName())
        print('要素数:', layer.GetFeatureCount())

读取完数据后，我们来看看如何写入数据。对于栅格数据，可以创建一个新的栅格数据集，写入数据，然后保存：

    driver = gdal.GetDriverByName('GTiff')  # 使用GTiff驱动创建TIFF文件
    out_dataset = driver.Create('output_raster.tif', cols, rows, bands, gdal.GDT_Byte)
    out_dataset.SetGeoTransform(geotransform)  # 设置地理变换信息
    out_dataset.SetProjection(projection)     # 设置坐标系统

对于矢量数据，你可以创建一个新的矢量文件，然后添加图层和要素：

    driver = ogr.GetDriverByName('ESRI Shapefile')
    out_dataset = driver.CreateDataSource('output_vector.shp')
    out_layer = out_dataset.CreateLayer('output_layer', srs=None, geom_type=ogr.wkbPoint)

最后，不要忘记在操作完成后释放资源：

    raster_dataset = None
    vector_dataset = None
    out_dataset = None

为了深入学习GDAL的高级操作和细节，建议阅读《GDAL-Python库详解：操作地理空间数据的利器》这本书。它不仅详细介绍了如何读取和写入数据，还包括了如何处理复杂的数据转换和分析任务。通过全面的学习，你将能够更有效地利用GDAL库来处理地理空间数据，为GIS项目提供强大的支持。

参考资源链接：[GDAL-Python库详解：操作地理空间数据的利器](https://wenku.csdn.net/doc/mbr98i2c3e?spm=1055.2569.3001.10343)