libgeotiff库的示例程序

libgeotiff是一个开源的地理空间图像文件读写库，主要用于处理TIFF（Tagged Image File Format）格式的数据，特别是那些包含地理元数据（如坐标系统、投影等）的图像。它支持多种GIS相关的操作，如读取、修改和写入地理TIFF文件。

下面是一个简单的Python示例，使用GDAL（Python版本的封装库，底层依赖libgeotiff）来打开并显示一个TIFF文件：

from osgeo import gdal

# 打开TIFF文件
ds = gdal.Open("path_to_your_geotiff_file.tif")

# 检查文件是否成功打开
if ds is None:
    print("无法打开文件")
else:
    # 获取GeoTransform信息（用于描述栅格的几何变换）
    geotransform = ds.GetGeoTransform()
    
    # 获取第一幅图层（bands数目等于1的单波段文件）
    band = ds.GetRasterBand(1)
    
    # 计算并打印图像宽度和高度
    width = ds.RasterXSize
    height = ds.RasterYSize
    print(f"Image size: {width}x{height}")
    
    # 读取并显示部分像素值
    x, y = 0, 0
    pixel_value = band.ReadAsArray(x, y, 5, 5)  # 读取5x5像素块
    print(f"Pixel value at ({x}, {y}):", pixel_value)

    # 关闭文件
    ds = None