利用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 相关文章推荐
在Linux下调试Python代码的各种方法
Apr 17 Python
Python中tell()方法的使用详解
May 24 Python
python实现字符串连接的三种方法及其效率、适用场景详解
Jan 13 Python
python正则表达式爬取猫眼电影top100
Feb 24 Python
Python的条件表达式和lambda表达式实例
Jan 31 Python
Python爬取智联招聘数据分析师岗位相关信息的方法
Aug 13 Python
Python控制台输出时刷新当前行内容而不是输出新行的实现
Feb 21 Python
利用 Python ElementTree 生成 xml的实例
Mar 06 Python
Python在终端通过pip安装好包以后在Pycharm中依然无法使用的问题(三种解决方案)
Mar 10 Python
python实现opencv+scoket网络实时图传
Mar 20 Python
python是怎么被发明的
Jun 15 Python
python 制作网站筛选工具(附源码)
Jan 21 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调用mysql存储过程
2007/02/14 PHP
php Ajax乱码
2008/04/09 PHP
PHP中用接口、抽象类、普通基类实现“面向接口编程”与“耦合方法”简述
2011/03/23 PHP
PHP操作MongoDB时的整数问题及对策说明
2011/05/02 PHP
php操作redis缓存方法分享
2015/06/03 PHP
PHP获取数据库表中的数据插入新的表再原删除数据方法
2018/10/12 PHP
JQuery扩展插件Validate—4设置错误提示的样式
2011/09/05 Javascript
js中取得变量绝对值的方法
2015/01/03 Javascript
jQuery实现仿腾讯微博滑出效果报告每日天气的方法
2015/05/11 Javascript
简单讲解AngularJS的Routing路由的定义与使用
2016/03/05 Javascript
Web性能优化系列 10个提升JavaScript性能的技巧
2016/09/27 Javascript
JavaScript设计模式之工厂模式和抽象工厂模式定义与用法分析
2018/07/26 Javascript
微信小程序框架wepy之动态控制类名
2018/09/14 Javascript
vue 使某个组件不被 keep-alive 缓存的方法
2018/09/21 Javascript
javascript实现遮罩层动态效果实例
2019/05/14 Javascript
Vue 设置axios请求格式为form-data的操作步骤
2019/10/29 Javascript
vue倒计时刷新页面不会从头开始的解决方法
2020/03/03 Javascript
vue style width a href动态拼接问题的解决
2020/08/07 Javascript
js实现筛选功能
2020/11/24 Javascript
Python可跨平台实现获取按键的方法
2015/03/05 Python
讲解Python中的递归函数
2015/04/27 Python
利用Pandas 创建空的DataFrame方法
2018/04/08 Python
Python设计模式之模板方法模式实例详解
2019/01/17 Python
Python多线程多进程实例对比解析
2020/03/12 Python
CSS3实现文本垂直排列的方法
2018/07/10 HTML / CSS
使用HTML5中的contentEditable来将多行文本自动增高
2016/03/01 HTML / CSS
欧洲当代手工玻璃和瓷器的领先品牌:LSA International
2018/06/03 全球购物
机关办公室岗位职责
2014/04/16 职场文书
分公司任命书
2014/06/06 职场文书
公司优秀员工获奖感言
2014/08/14 职场文书
公司给客户的感谢信
2015/01/23 职场文书
幼儿园小班个人工作总结
2015/02/12 职场文书
《抽屉原理》教学反思
2016/02/20 职场文书
opencv用VS2013调试时用Image Watch插件查看图片
2021/07/26 Python
《王国之心》迎来了发售的20周年, 野村哲发布贺图
2022/04/11 其他游戏
iOS 16进一步确认,一共支持16款iPhone
2022/04/28 数码科技