Python遥感开发之解决TIF数据之间行列不一致的问题
前言:主要解决在同一分辨率的情况下,遥感数据之间行和列数据不一致的问题。
1.问题如图所示
我们发现这两个TIF的分辨率是一样的,都是0.01x0.01,但是tasmax数据比lucc多出一个行和列,我们通过对比发现是多出了第一行和最后一列,只需要把tasmax的数据的第一行和最后一列删除即可。可以使用这两行代码解决问题。
data = np.delete(data, 0, axis=0)#行
data = np.delete(data, 554, axis=1)#列
2.完整代码如下所示
import os
import numpy as np
from osgeo import gdal, gdalnumeric
def read_tif02(filepath):
data = gdalnumeric.LoadFile(filepath)
data = data.astype(np.float32)
a = data[0][0]
data[data == a