服务器之家

服务器之家 > 正文

对Python3+gdal 读取tiff格式数据的实例讲解

时间:2021-04-25 00:23     来源/作者:竹子熊猫

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

?
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
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.图像裁剪:

?
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
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

以上这篇对Python3+gdal 读取tiff格式数据的实例讲解就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持服务器之家。

原文链接:https://blog.csdn.net/summermaoz/article/details/78346929

标签:

相关文章

热门资讯

2020微信伤感网名听哭了 让对方看到心疼的伤感网名大全
2020微信伤感网名听哭了 让对方看到心疼的伤感网名大全 2019-12-26
yue是什么意思 网络流行语yue了是什么梗
yue是什么意思 网络流行语yue了是什么梗 2020-10-11
背刺什么意思 网络词语背刺是什么梗
背刺什么意思 网络词语背刺是什么梗 2020-05-22
苹果12mini价格表官网报价 iPhone12mini全版本价格汇总
苹果12mini价格表官网报价 iPhone12mini全版本价格汇总 2020-11-13
Intellij idea2020永久破解,亲测可用!!!
Intellij idea2020永久破解,亲测可用!!! 2020-07-29
返回顶部