对Python3+gdal 读取tiff格式数据的实例讲解
作者:竹子熊猫 发布时间:2023-03-10 10:21:53
标签:Python3,gdal,tiff
1、遇到的问题:numpy版本
im_data = dataset.ReadAsArray(0,0,im_width,im_height)#获取数据 这句报错
升级numpy:pip install -U numpy 但是提示已经是最新版本
解决:卸载numpy 重新安装
2.直接从压缩包中读取tiff图像
参考:http://gdal.org/gdal_virtual_file_systems.html#gdal_virtual_file_systems_vsizip
当前情况是2层压缩: /'/vsitar/C:/Users/summer/Desktop/a_PAN1.tiff'
3.读tiff
def readTif(fileName):
merge_img = 0
driver = gdal.GetDriverByName('GTiff')
driver.Register()
dataset = gdal.Open(fileName)
if dataset == None:
print(fileName+ "掩膜失败,文件无法打开")
return
im_width = dataset.RasterXSize #栅格矩阵的列数
print('im_width:', im_width)
im_height = dataset.RasterYSize #栅格矩阵的行数
print('im_height:', im_height)
im_bands = dataset.RasterCount #波段数
im_geotrans = dataset.GetGeoTransform()#获取仿射矩阵信息
im_proj = dataset.GetProjection()#获取投影信息
if im_bands == 1:
band = dataset.GetRasterBand(1)
im_data = dataset.ReadAsArray(0,0,im_width,im_height) #获取数据
cdata = im_data.astype(np.uint8)
merge_img = cv2.merge([cdata,cdata,cdata])
cv2.imwrite('C:/Users/summer/Desktop/a.jpg', merge_img)
#
elif im_bands == 4:
# # im_data = dataset.ReadAsArray(0,0,im_width,im_height)#获取数据
# # im_blueBand = im_data[0,0:im_width,0:im_height] #获取蓝波段
# # im_greenBand = im_data[1,0:im_width,0:im_height] #获取绿波段
# # im_redBand = im_data[2,0:im_width,0:im_height] #获取红波段
# # # im_nirBand = im_data[3,0:im_width,0:im_height] #获取近红外波段
# # merge_img=cv2.merge([im_redBand,im_greenBand,im_blueBand])
# # zeros = np.zeros([im_height,im_width],dtype = "uint8")
# # data1 = im_redBand.ReadAsArray
# band1=dataset.GetRasterBand(1)
# band2=dataset.GetRasterBand(2)
# band3=dataset.GetRasterBand(3)
# band4=dataset.GetRasterBand(4)
data1=band1.ReadAsArray(0,0,im_width,im_height).astype(np.uint16) #r #获取数据
data2=band2.ReadAsArray(0,0,im_width,im_height).astype(np.uint16) #g #获取数据
data3=band3.ReadAsArray(0,0,im_width,im_height).astype(np.uint16) #b #获取数据
data4=band4.ReadAsArray(0,0,im_width,im_height).astype(np.uint16) #R #获取数据
# print(data1[1][45])
# output1= cv2.convertScaleAbs(data1, alpha=(255.0/65535.0))
# print(output1[1][45])
# output2= cv2.convertScaleAbs(data2, alpha=(255.0/65535.0))
# output3= cv2.convertScaleAbs(data3, alpha=(255.0/65535.0))
merge_img1 = cv2.merge([output3,output2,output1]) #B G R
cv2.imwrite('C:/Users/summer/Desktop/merge_img1.jpg', merge_img1)
4.图像裁剪:
import cv2
import numpy as np
import os
tiff_file = './try_img/2.tiff'
save_folder = './try_img_re/'
if not os.path.exists(save_folder):
os.makedirs(save_folder)
tif_img = cv2.imread(tiff_file)
width, height, channel = tif_img.shape
# print height, width, channel : 6908 7300 3
threshold = 1000
overlap = 100
step = threshold - overlap
x_num = width/step + 1
y_num = height/step + 1
print x_num, y_num
N = 0
yj = 0
for xi in range(x_num):
for yj in range(y_num):
# print xi
if yj <= y_num:
print yj
x = step*xi
y = step*yj
wi = min(width,x+threshold)
hi = min(height,y+threshold)
# print wi , hi
if wi-x < 1000 and hi-y < 1000:
im_block = tif_img[wi-1000:wi, hi-1000:hi]
elif wi-x > 1000 and hi-y < 1000:
im_block = tif_img[x:wi, hi-1000:hi]
elif wi-x < 1000 and hi-y > 1000:
im_block = tif_img[wi-1000:wi, y:hi]
else:
im_block = tif_img[x:wi,y:hi]
cv2.imwrite(save_folder + 'try' + str(N) + '.jpg', im_block)
N += 1
来源:https://blog.csdn.net/summermaoz/article/details/78346929
0
投稿
猜你喜欢
- 是否曾经有过这样的经历:把一个元素置于另一个元素之上,而希望下面的那个元素成为可点击的?现在,利用css的pointer-events属性即
- 第一章:日志管理 1.forcing log switchessql> alter system switch logfile;2.f
- 前言接着上一篇:AI识别照片是谁,人脸识别face_recognition开源项目安装使用根据项目提供的demo代码,调整了一下功能,自己写
- 本文实例讲述了Python删除windows垃圾文件的方法。分享给大家供大家参考。具体如下:#coding:utf-8import os#f
- 在ASP中,直接使用“Insert into” 语句与使用ADO中AddNew方法有什么区别?哪一种更好呢?AddNew方法的实质就是封装了
- 客户端从服务端下载文件的流程分析: 浏览器发送一个请求,请求访问服务器中的某个网页(如:down.php),该网页的代码如下。 服务器接受到
- 以XML格式查看查询结果通过使用传统—xml 选项调用MySQL命令行客户程序,你可以以XML格式(而不是传统的列表形式
- 在项目中,尤其是pc端的时候,我们在用户登录后会给前端返回一个标识,来判断用户是否登录,这个标识大多数都是用户的id  
- 这次我们来说一些很小的东西,相当小。不过先说一个故事:这天晚上,你打算出去下一次馆子,就行动了,找到了一家新开的小店进去坐了下来,服务员很殷
- 本文实例讲述了PHP函数按引用传递参数及函数可选参数用法。分享给大家供大家参考,具体如下:一、函数按引用传递参数1. 代码<!DOCT
- 前言通常执行 python 程序要有相应的 Python 环境,但某些特定场景下,我们可能并不愿意这么麻烦的去配置这些环境(比如将写好的脚本
- 内容摘要:一堆数据摆面前,数据背后有什么样的事情在发生,这些数据里面暗藏着什么样的用户需求,什么样的商业机会?看懂这些,将为未来产品设计的方
- 效果图实例代码如下:<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transit
- JS代码:function showFlash(src,w,h){ html&nbs
- Adobe AIR 在我们国庆节的时候发布了beta2版。 同时发布的还有: FLEX 3.0 beta2 官方说这一版本已经基本上接近正式
- 为了防止机器人频繁登陆网站或者破坏分子恶意登陆,很多用户登录和注册系统都提供了图形验证码功能。验证码(CAPTCHA)是“Completel
- 这一段CSS代码相当简单,目的就是想用CSS来控制某段文字的显示与隐藏。起初我采用了下面的代码,令人不可思议的是,它们在我的IE6.0上居然
- 本文实例讲述了Python 操作mysql数据库查询之fetchone(), fetchmany(), fetchall()用法。分享给大家
- 后台数据库用是Access,客户用了一年后说打开界面非常慢,查看了数据库后发现数据表中的记录已有五万多条,自己试过将记录复制到10 万条,打
- 本文实例讲述了Yii2框架实现登陆添加验证码功能。分享给大家供大家参考,具体如下:models中LoginForm.phppublic $v