1、GDAL进行图像一块区域的裁剪和大图像裁剪为多个小图像

最新推荐文章于 2023-08-15 17:12:40 发布

原创

最新推荐文章于 2023-08-15 17:12:40 发布 · 2.6k 阅读

26 ·

CC 4.0 BY-SA版权

文章标签：

#计算机视觉 #图像处理 #opencv #python

本文介绍了一种将大型遥感影像切割成多个较小图像的方法，包括保持地理坐标信息的切割方式及忽略坐标信息的裁剪技术。通过Python编程实现，使用GDAL库处理图像，并提供完整的代码示例。

将大的遥感影像拆建为多个

1 代码实现（一小块儿区域裁剪）
2 将一大块区域裁剪为小块图像（损失了坐标）
参考资料

1 代码实现（一小块儿区域裁剪）

# 将一幅遥感影像拆建为多张
import gdal

# 读取要切的原图
in_ds = gdal.Open("zx1.jpg")
print("open tif file succeed")
width = in_ds.RasterXSize  # 获取数据宽度
height = in_ds.RasterYSize  # 获取数据高度
outbandsize = in_ds.RasterCount  # 获取数据波段数
im_geotrans = in_ds.GetGeoTransform()  # 获取仿射矩阵信息
im_proj = in_ds.GetProjection()  # 获取投影信息
datatype = in_ds.GetRasterBand(1).DataType
im_data = in_ds.ReadAsArray()  # 获取数据

# 读取原图中的每个波段
in_band1 = in_ds.GetRasterBand(1)
in_band2 = in_ds.GetRasterBand(2)
in_band3 = in_ds.GetRasterBand(3)

# 定义切图的起始点坐标
offset_x = 0
offset_y = 0

# offset_x = width/2  # 这里是随便选取的，可根据自己的实际需要设置
# offset_y = height/2

# 定义切图的大小（矩形框）
block_xsize = 60  # 行
block_ysize = 60  # 列

k = 0
for i in range(width // block_xsize):
    for j in range(height // block_xsize):
        out_band1 = in_band1.ReadAsArray(i * block_xsize, j * block_xsize, block_xsize, block_ysize)
        out_band2 = in_band2.ReadAsArray(i * block_xsize, j * block_xsize, block_xsize, block_ysize)
        out_band3 = in_band3.ReadAsArray(i * block_xsize, j * block_xsize, block_xsize, block_ysize)
        print(out_band3)
        k += 1

        ## 从每个波段中切需要的矩形框内的数据(注意读取的矩形框不能超过原图大小)
        # out_band1 = in_band1.ReadAsArray(offset_x, offset_y, block_xsize, block_ysize)
        # out_band2 = in_band2.ReadAsArray(offset_x, offset_y, block_xsize, block_ysize)
        # out_band3 = in_band3.ReadAsArray(offset_x, offset_y, block_xsize, block_ysize)

        # 获取Tif的驱动，为创建切出来的图文件做准备
        gtif_driver = gdal.GetDriverByName("GTiff")

        # 创建切出来的要存的文件（3代表3个不都按，最后一个参数为数据类型，跟原文件一致）
        out_ds = gtif_driver.Create(str(k) + 'clip4.tif', block_xsize, block_ysize, outbandsize, datatype)
        # print("create new tif file succeed")

        # 获取原图的原点坐标信息，# 获取仿射矩阵信息
        ori_transform = in_ds.GetGeoTransform()
        if ori_transform:
            print(ori_transform)
            print("Origin = ({}, {})".format(ori_transform[0], ori_transform[3]))
            print("Pixel Size = ({}, {})".format(ori_transform