利用python GDAL库读写geotiff格式的遥感影像方法


Posted in Python onNovember 29, 2018

如下所示:

from osgeo import gdal
import numpy as np
def read_tiff(inpath):
  ds=gdal.Open(inpath)
  row=ds.RasterXSize
  col=ds.RasterYSize
  band=ds.RasterCount
  geoTransform=ds.GetTransform()
  proj=ds.GetTransform()
  data=np.zeros([row,col,band])
  for i in range(band):
   dt=ds.GetRasterBand(1)
   data[:,:,i]=dt.ReadAsArray(0,0,col,row)
  return data
 
def array2raster(outpath,array,geoTransform,proj):
 cols=array.shape[1]
 rows=array.shape[0]
 driver=gdal.GetDriverByName('Gtiff')
 outRaster=driver.Create(newRasterfn,cols,rows,1,gdal.GDT_Byte)
 outRaster.SetGeoTransform(geoTransform)#参数2,6为水平垂直分辨率,参数3,5表示图片是指北的
 outband=outRaster.GetRasterBand(1)
 outband.WriteArray(array)
 outRaster.SetProjection(proj)#将几何对象的数据导出为wkt格式
 outRaster.FlushCache()
 
if _name=="_main_":
 
 data,geoTransform,proj=read_tiff('d:/a.tif')
 
 array2raster("d:/b.tif",np.zeros[2400,2400],geoTransform,proj)

利用python GDAL库读写geotiff格式的遥感影像,并生成与原影像具有相同地理坐标和投影坐标的geotiff格式图片。

以上这篇利用python GDAL库读写geotiff格式的遥感影像方法就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
详解Django缓存处理中Vary头部的使用
Jul 24 Python
使用Python的urllib2模块处理url和图片的技巧两则
Feb 18 Python
基于Django的python验证码(实例讲解)
Oct 23 Python
同时安装Python2 & Python3 cmd下版本自由选择的方法
Dec 09 Python
Python爬取当当、京东、亚马逊图书信息代码实例
Dec 09 Python
Python数据结构与算法之常见的分配排序法示例【桶排序与基数排序】
Dec 15 Python
Python模块搜索路径代码详解
Jan 29 Python
python的pandas工具包,保存.csv文件时不要表头的实例
Jun 14 Python
python遍历文件目录、批量处理同类文件
Aug 31 Python
python用Tkinter做自己的中文代码编辑器
Sep 07 Python
python编写函数注意事项总结
Mar 29 Python
pytest实现多进程与多线程运行超好用的插件
Jul 15 Python
在python中利用GDAL对tif文件进行读写的方法
Nov 29 #Python
使用python判断你是青少年还是老年人
Nov 29 #Python
pycham查看程序执行的时间方法
Nov 29 #Python
pycharm在调试python时执行其他语句的方法
Nov 29 #Python
Selenium chrome配置代理Python版的方法
Nov 29 #Python
在PyCharm中实现关闭一个死循环程序的方法
Nov 29 #Python
selenium使用chrome浏览器测试(附chromedriver与chrome的对应关系表)
Nov 29 #Python
You might like
php实现上传图片保存到数据库的方法
2015/02/11 PHP
php遍历目录下文件并按修改时间排序操作示例
2019/07/12 PHP
Javascript select下拉框操作常用方法
2009/11/09 Javascript
用document.documentElement取代document.body的原因分析
2009/11/12 Javascript
为JavaScript添加重载函数的辅助方法
2010/07/04 Javascript
纯js分页代码(简洁实用)
2013/11/05 Javascript
javascript实例--教你实现扑克牌洗牌功能
2014/05/15 Javascript
node.js中的dns.getServers方法使用说明
2014/12/08 Javascript
JavaScript中的fontsize()方法使用详解
2015/06/08 Javascript
简述JavaScript的正则表达式中test()方法的使用
2015/06/16 Javascript
JavaScript判断是否是微信浏览器
2016/06/13 Javascript
深入理解Angular中的依赖注入
2017/06/26 Javascript
详解jQuery同步Ajax带来的UI线程阻塞问题及解决办法
2017/08/09 jQuery
基于bootstrap写的一点localStorage本地储存
2017/11/21 Javascript
vuex提交state&&实时监听state数据的改变方法
2018/09/16 Javascript
Vue 配合eiement动态路由,权限验证的方法
2018/09/26 Javascript
layui table单元格事件修改值的方法
2019/09/24 Javascript
nodejs实现UDP组播示例方法
2019/11/04 NodeJs
[47:53]DOTA2上海特级锦标赛主赛事日 - 1 败者组第一轮#2COL VS Spirit
2016/03/02 DOTA
python re正则表达式模块(Regular Expression)
2014/07/16 Python
在Python的列表中利用remove()方法删除元素的教程
2015/05/21 Python
Python批量重命名同一文件夹下文件的方法
2015/05/25 Python
Python中super函数的用法
2017/11/17 Python
python中sort和sorted排序的实例方法
2019/08/26 Python
django执行数据库查询之后实现返回的结果集转json
2020/03/31 Python
python实现感知机模型的示例
2020/09/30 Python
美国正宗设计师眼镜在线零售商:EYEZZ
2019/03/23 全球购物
介绍一下linux的文件系统
2012/03/20 面试题
《周恩来的四个昼夜》观后思想汇报范文两篇
2014/09/10 职场文书
市场督导岗位职责
2015/04/10 职场文书
社区节水倡议书
2015/04/29 职场文书
家长意见和建议怎么写
2015/06/04 职场文书
大学生先进个人主要事迹材料
2015/11/04 职场文书
创业计划书之健康营养产业
2019/10/15 职场文书
python 用递归实现通用爬虫解析器
2021/04/16 Python
剖析后OpLog订阅MongoDB的数据变更就没那么难了
2022/02/24 MongoDB