目的
现有两幅栅格图像,一个是某地区道路栅格图,一个是某地区土地利用类型图,需要将道路叠加到土地利用类型图中,即叠加后,重合的像元值以道路图为准,其余的像元值仍是土地利用类型图原有的像元值。
图1 道路信息图
图2 土地利用类型图
图3 结果图
具体实现
from gdalconst import * from osgeo import gdal import osr import sys import copy #叠加两个栅格图像(一个道路栅格图,一个土地利用类型图),两幅图像重叠的像元值都是第一个图像的值, #未重叠的像元值还是土地利用类型图上的值,最终结果便是土地利用类型图上面多了道路信息。 roadFile = 'E:\\Exercise\\test\\grasstest\\road_rastercalc.tif' landuseFile = 'E:\\Exercise\\test\\grasstest\\landuse.tif' roadDs = gdal.Open(roadFile, GA_ReadOnly) landuseDs = gdal.Open(landuseFile, GA_ReadOnly) if roadDs is None: print 'Can not open ', roadFile sys.exit(1) geotransform = roadDs.GetGeoTransform() projection=roadDs.GetProjection() cols = roadDs.RasterXSize rows = roadDs.RasterYSize roadBand = roadDs.GetRasterBand(1) roadData = roadBand.ReadAsArray(0,0,cols,rows) roadNoData = roadBand.GetNoDataValue() landuseBand = landuseDs.GetRasterBand(1) landuseData = landuseBand.ReadAsArray(0,0,cols,rows) landuseNoData = landuseBand.GetNoDataValue() result = landuseData for i in range(0,rows): for j in range(0,cols): if(abs(roadData[i,j] - 20) < 0.0001): result[i,j] = 20 if((abs(landuseData[i,j] - landuseNoData)>0.0001) and (abs(roadData[i,j] - roadNoData) < 0.0001)): result[i,j] = landuseData[i,j] if((abs(landuseData[i,j] - landuseNoData)<0.0001) and (abs(roadData[i,j] - roadNoData) < 0.0001)): result[i,j] = landuseNoData #write result to disk resultPath = 'E:\\Exercise\\test\\grasstest\\result_landuse.tif' format = "GTiff" driver = gdal.GetDriverByName(format) ds = driver.Create(resultPath, cols, rows, 1, GDT_Float32) ds.SetGeoTransform(geotransform) ds.SetProjection(projection) ds.GetRasterBand(1).SetNoDataValue(landuseNoData) ds.GetRasterBand(1).WriteArray(result) ds = None print 'ok---------'
以上这篇Python叠加两幅栅格图像的实现方法就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持。
华山资源网 Design By www.eoogi.com
广告合作:本站广告合作请联系QQ:858582 申请时备注:广告合作(否则不回)
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
免责声明:本站资源来自互联网收集,仅供用于学习和交流,请遵循相关法律法规,本站一切资源不代表本站立场,如有侵权、后门、不妥请联系本站删除!
华山资源网 Design By www.eoogi.com
暂无评论...
更新日志
2024年09月23日
2024年09月23日
- 华少翌.2005-华少翌【中唱艺能】【WAV+CUE】
- 群星《十二女声Ⅱ》[WAV/CUE/分轨][527MB]
- 大自然音乐系列《漫步在山裡》[WAV/CUE/分轨][502.8MB]
- 大自然音乐系列《星光夜语》[WAV/CUE/分轨][453.3MB]
- 史依弘《史依弘交响乐京剧演唱专辑HQ》头版限量编号[WAV+CUE]
- 曼丽2020《女人三十风继续吹》HQCD[WAV+CUE]
- 曼丽2021《女人三十2雾之恋》HQCD[WAV+CUE]
- 隔壁老樊.2019-我曾【好听音乐】【FLAC分轨】
- 高林生.1993-牵挂你的人是我【白天鹅】【WAV+CUE】
- 洪卓立.2012-找个懂我的女孩(国专)【英皇娱乐】【WAV+CUE】
- 曼丽2022《女人三十3晚秋》HQCD[WAV+CUE]
- 曼丽2018《女人三十4冬天》HQCD[WAV+CUE]
- 柏菲-《好歌30年壹》限量开盘母带ORMCD[低速原抓WAV+CUE]
- 大自然音乐系列《水舞》[WAV/CUE/分轨][424.3MB]
- 大自然音乐系列《星光夜语》[WAV/CUE/分轨][453.3MB]