掩膜行列和空间分辨率不一致的数据

本文主要是介绍掩膜行列和空间分辨率不一致的数据，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

欢迎关注我的个人公众号：小Rser

掩膜行列和空间分辨率不一致的数据https://mp.weixin.qq.com/s?__biz=MzkyNjMzNTQ2Mw==&mid=2247483929&idx=1&sn=c16db2a848b6013e3526fbfc79e517b8&chksm=c239aa91f54e2387c15f5dc48f7ef368e6a745112a98f22a0302d9093c534512603efa69f0b3&token=2086761678&lang=zh_CN#rd

import gdal#基准数据
ref_file = r"D:\微信公众号\掩膜\基准数据.tif"
#待匹配数据
pend_file = r"D:\微信公众号\掩膜\待匹配数据.tif"
# 结果文件
outfile = r"D:\微信公众号\掩膜\匹配后数据.tif"
ref_ds = gdal.Open(ref_file)
base_xsize = ref_ds.RasterXSize
base_ysize = ref_ds.RasterYSize
base_proj = ref_ds.GetProjection()
base_gt = ref_ds.GetGeoTransform()
xmin = base_gt[0]
ymax = base_gt[3]
xmax = xmin + base_xsize * base_gt[1]
ymin = ymax + base_ysize * base_gt[5]
#待处理数据
pend_ds = gdal.Open(pend_file)
#数据匹配
dataset = gdal.Warp(outfile,pend_ds,width=base_xsize,height=base_ysize,srcNodata=-999 ,dstNodata=-999,dstSRS=str(base_proj),outputBounds=(xmin, ymin, xmax, ymax),format="GTiff",resampleAlg=gdal.GRA_Bilinear)

import gdal#匹配后文件
infile = r"D:\微信公众号\掩膜\匹配后数据.tif"
#掩膜文件
maskfile = r"D:\微信公众号\掩膜\基准数据.tif"
#结果文件
outfile = r"D:\微信公众号\掩膜\掩膜后数据.tif"#读取匹配后文件
inds = gdal.Open(infile)
indata = inds.ReadAsArray()
cols = inds.RasterXSize
rows = inds.RasterYSize
geo = inds.GetGeoTransform()
proj = inds.GetProjection()
band = inds.GetRasterBand(1)
#读取掩膜文件
maskds = gdal.Open(maskfile)
maskdata = maskds.ReadAsArray()
band_mask = maskds.GetRasterBand(1)
#掩膜值
mask_nodata = band_mask.GetNoDataValue()
#掩膜
indata[maskdata == mask_nodata] = -999
#输出结果，无效值设置为-999
driver = gdal.GetDriverByName("Gtiff")
out_ds = driver.Create(outfile, cols, rows, 1, band.DataType)
out_ds.SetGeoTransform(geo)
out_ds.SetProjection(proj)
out_band = out_ds.GetRasterBand(1)
out_band.WriteArray(indata)
out_band.SetNoDataValue(-999)