遥感影像分割成多幅小图像方法

目的

本人需要将一幅Landsat8影像，切割成多幅256 * 256的小幅影像，用于深度学习相关处理。
经过在envi中的可视化的裁剪，得到了一幅5130 * 5230的影像，我经过计算，发现可以横向切出
256 * 20，纵向同样切出256*20 =5120， 共计400幅影像。

方法

利用python

• gdal模块读取影像并存入numpy数组
• 利用循环和opencv进行写入（实际用的时候要用gdal，可关注我后面有关landsat8处理小工具的博文）
  https://blog.youkuaiyun.com/RSstudent/article/details/108724556?spm=1001.2014.3001.5501
  其中代码的329-378行就是所需内容，还附带了地理编码和投影，地理信息不会丢失，可在实际中使用。
• 利用os模块得到路径

代码

import numpy as np
from osgeo import gdal
from osgeo import gdal_array
import cv2
import os

filename = 'tiffbigcut.tif'
dataset = gdal.Open(filename)

datatype = dataset.G