Python叠加两幅栅格图像的实现方法


Posted in Python onJuly 05, 2019

目的

现有两幅栅格图像,一个是某地区道路栅格图,一个是某地区土地利用类型图,需要将道路叠加到土地利用类型图中,即叠加后,重合的像元值以道路图为准,其余的像元值仍是土地利用类型图原有的像元值。

图1 道路信息图

Python叠加两幅栅格图像的实现方法

图2 土地利用类型图

Python叠加两幅栅格图像的实现方法

图3 结果图

Python叠加两幅栅格图像的实现方法

具体实现

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叠加两幅栅格图像的实现方法就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
详解python分布式进程
Oct 08 Python
详解python-图像处理(映射变换)
Mar 22 Python
python将类似json的数据存储到MySQL中的实例
Jul 12 Python
Django ORM 自定义 char 类型字段解析
Aug 09 Python
python飞机大战pygame游戏背景设计详解
Dec 17 Python
python各层级目录下import方法代码实例
Jan 20 Python
Python环境管理virtualenv&amp;virtualenvwrapper的配置详解
Jul 01 Python
python对输出的奇数偶数排序实例代码
Dec 04 Python
Django前后端分离csrf token获取方式
Dec 25 Python
python中if嵌套命令实例讲解
Feb 25 Python
详解Django中的FBV和CBV对比分析
Mar 01 Python
浅谈Python协程asyncio
Jun 20 Python
用vue.js组件模拟v-model指令实例方法
Jul 05 #Python
python买卖股票的最佳时机(基于贪心/蛮力算法)
Jul 05 #Python
Python 计算任意两向量之间的夹角方法
Jul 05 #Python
python实现两个经纬度点之间的距离和方位角的方法
Jul 05 #Python
Python3+Appium实现多台移动设备操作的方法
Jul 05 #Python
Python PIL读取的图像发生自动旋转的实现方法
Jul 05 #Python
python读出当前时间精度到秒的代码
Jul 05 #Python
You might like
PHP array_flip() 删除重复数组元素专用函数
2010/05/16 PHP
有关phpmailer的详细介绍及使用方法
2013/01/28 PHP
PHP下通过QRCode类库创建中间带网站LOGO的二维码
2014/07/12 PHP
全面解读PHP的Yii框架中的日志功能
2016/03/17 PHP
PHP设计模式之装饰器模式定义与用法详解
2018/04/02 PHP
动态加载js文件 document.createElement
2006/10/14 Javascript
javascript 数据类型转换(parseInt,parseFloat)
2010/07/20 Javascript
原生JS实现表单checkbook获取已选择的值
2013/07/21 Javascript
javascript实现自动输出文本(打字特效)
2015/08/27 Javascript
JavaScript语言精粹经典实例(整理篇)
2016/06/07 Javascript
详解js中的apply与call的用法
2016/07/30 Javascript
Bootstrap模态框禁用空白处点击关闭
2016/10/20 Javascript
jQuery基于ajax操作json数据简单示例
2017/01/05 Javascript
浅谈angular.copy() 深拷贝
2017/09/14 Javascript
vue移动端轻量级的轮播组件实现代码
2018/07/12 Javascript
解决vue中虚拟dom,无法实时更新的问题
2018/09/15 Javascript
vue构建动态表单的方法示例
2018/09/22 Javascript
JS实现li标签的删除
2019/04/12 Javascript
pm2发布node配置文件ecosystem.json详解
2019/05/15 Javascript
在Vue项目中用fullcalendar制作日程表的示例代码
2019/08/04 Javascript
浅谈layui使用模板引擎动态渲染元素要注意的问题
2019/09/14 Javascript
简单实现python爬虫功能
2015/12/31 Python
wx.CheckBox创建复选框控件并响应鼠标点击事件
2018/04/25 Python
python3.X 抓取火车票信息【修正版】
2018/06/19 Python
python爬虫之自动登录与验证码识别
2020/06/15 Python
python读出当前时间精度到秒的代码
2019/07/05 Python
python使用 __init__初始化操作简单示例
2019/09/26 Python
python实现最短路径的实例方法
2020/07/19 Python
基于python判断字符串括号是否闭合{}[]()
2020/09/21 Python
pytorch学习教程之自定义数据集
2020/11/10 Python
css3动画效果小结(推荐)
2016/07/25 HTML / CSS
iHerb台湾:维生素、保健品和健康产品
2018/01/31 全球购物
学校社会实践活动总结
2014/07/03 职场文书
党员评议自我评价
2015/03/03 职场文书
大学体育课感想
2015/08/10 职场文书
男方家长婚礼答谢词
2015/09/29 职场文书