python遥感影像按影像行列号批量刺点

这里提供一种使用Python和GDAL库对遥感影像进行按影像行列号批量刺点的方法：

1. 安装GDAL库

使用pip或conda安装GDAL库，安装命令如下：

pip install gdal

或

conda install gdal

2. 编写Python脚本

以下是一个示例脚本，其中包含了对影像进行读取、获取行列号、刺点等操作：

import gdal
import ogr

# 打开影像文件
dataset = gdal.Open('image.tif')

# 获取影像信息
x_size = dataset.RasterXSize
y_size = dataset.RasterYSize
geo_transform = dataset.GetGeoTransform()
srs = dataset.GetProjection()

# 创建矢量数据源
driver = ogr.GetDriverByName('ESRI Shapefile')
out_file = driver.CreateDataSource('points.shp')
out_layer = out_file.CreateLayer('points', srs, ogr.wkbPoint)

# 创建属性表字段
field_defn = ogr.FieldDefn('id', ogr.OFTInteger)
out_layer.CreateField(field_defn)

# 循环遍历影像行列号，刺点
for y in range(y_size):
    for x in range(x_size):
        # 获取像素中心点坐标
        x_coord = geo_transform[0] + x * geo_transform[1] + y * geo_transform[2]
        y_coord = geo_transform[3] + x * geo_transform[4] + y * geo_transform[5]
        # 创建点要素
        point = ogr.Geometry(ogr.wkbPoint)
        point.AddPoint(x_coord, y_coord)
        # 创建属性表记录
        feature_defn = out_layer.GetLayerDefn()
        feature = ogr.Feature(feature_defn)
        feature.SetField('id', y * x_size + x + 1)
        feature.SetGeometry(point)
        # 写入要素
        out_layer.CreateFeature(feature)
        # 释放要素
        feature = None

# 关闭数据源
out_file = None

在脚本中，首先使用GDAL库打开影像文件，获取影像的行列数、地理参考等信息。然后使用OGR库创建一个矢量数据源和一个点图层，并设置属性表字段。接着循环遍历影像的行列号，计算每个像素中心点的坐标，并创建一个点要素和对应的属性表记录，写入到矢量数据源中。最后关闭数据源。

3. 运行脚本

将脚本保存为.py文件，使用Python运行即可。在运行过程中，需要将'image.tif'替换成实际的影像文件路径，'points.shp'替换成实际的输出矢量文件路径。运行完成后，可以在指定的输出矢量文件中查看刺点结果。

需要注意的是，该方法只适用于影像行列号与实际地理坐标存在线性关系的情况，对于存在旋转、错位等变形的影像可能需要进行坐标转换或其他处理。