欧美一级特黄大片做受成人-亚洲成人一区二区电影-激情熟女一区二区三区-日韩专区欧美专区国产专区

Python3+gdal如何讀取tiff格式數(shù)據(jù)-創(chuàng)新互聯(lián)

這篇文章將為大家詳細(xì)講解有關(guān)Python3+gdal如何讀取tiff格式數(shù)據(jù),小編覺得挺實(shí)用的,因此分享給大家做個(gè)參考,希望大家閱讀完這篇文章后可以有所收獲。

創(chuàng)新互聯(lián)建站是專業(yè)的晉源網(wǎng)站建設(shè)公司,晉源接單;提供成都做網(wǎng)站、網(wǎng)站設(shè)計(jì),網(wǎng)頁(yè)設(shè)計(jì),網(wǎng)站設(shè)計(jì),建網(wǎng)站,PHP網(wǎng)站建設(shè)等專業(yè)做網(wǎng)站服務(wù);采用PHP框架,可快速的進(jìn)行晉源網(wǎng)站開發(fā)網(wǎng)頁(yè)制作和功能擴(kuò)展;專業(yè)做搜索引擎喜愛的網(wǎng)站,專業(yè)的做網(wǎng)站團(tuán)隊(duì),希望更多企業(yè)前來(lái)合作!

1、遇到的問(wèn)題:numpy版本

im_data = dataset.ReadAsArray(0,0,im_width,im_height)#獲取數(shù)據(jù) 這句報(bào)錯(cuò)

升級(jí)numpy:pip install -U numpy 但是提示已經(jīng)是最新版本

解決:卸載numpy 重新安裝

2.直接從壓縮包中讀取tiff圖像

參考:http://gdal.org/gdal_virtual_file_systems.html#gdal_virtual_file_systems_vsizip

當(dāng)前情況是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+ "掩膜失敗,文件無(wú)法打開")
		return
	im_width = dataset.RasterXSize #柵格矩陣的列數(shù)
	print('im_width:', im_width) 
 
	im_height = dataset.RasterYSize #柵格矩陣的行數(shù)
	print('im_height:', im_height) 
	im_bands = dataset.RasterCount #波段數(shù)
	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) #獲取數(shù)據(jù)
		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)#獲取數(shù)據(jù)
	# 	# im_blueBand = im_data[0,0:im_width,0:im_height] #獲取藍(lán)波段
	# 	# 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 #獲取數(shù)據(jù)
		data2=band2.ReadAsArray(0,0,im_width,im_height).astype(np.uint16) #g #獲取數(shù)據(jù)
		data3=band3.ReadAsArray(0,0,im_width,im_height).astype(np.uint16) #b #獲取數(shù)據(jù)
		data4=band4.ReadAsArray(0,0,im_width,im_height).astype(np.uint16) #R #獲取數(shù)據(jù)
	# 	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

關(guān)于“Python3+gdal如何讀取tiff格式數(shù)據(jù)”這篇文章就分享到這里了,希望以上內(nèi)容可以對(duì)大家有一定的幫助,使各位可以學(xué)到更多知識(shí),如果覺得文章不錯(cuò),請(qǐng)把它分享出去讓更多的人看到。

另外有需要云服務(wù)器可以了解下創(chuàng)新互聯(lián)scvps.cn,海內(nèi)外云服務(wù)器15元起步,三天無(wú)理由+7*72小時(shí)售后在線,公司持有idc許可證,提供“云服務(wù)器、裸金屬服務(wù)器、高防服務(wù)器、香港服務(wù)器、美國(guó)服務(wù)器、虛擬主機(jī)、免備案服務(wù)器”等云主機(jī)租用服務(wù)以及企業(yè)上云的綜合解決方案,具有“安全穩(wěn)定、簡(jiǎn)單易用、服務(wù)可用性高、性價(jià)比高”等特點(diǎn)與優(yōu)勢(shì),專為企業(yè)上云打造定制,能夠滿足用戶豐富、多元化的應(yīng)用場(chǎng)景需求。

本文標(biāo)題:Python3+gdal如何讀取tiff格式數(shù)據(jù)-創(chuàng)新互聯(lián)
分享地址:http://aaarwkj.com/article40/ccoeeo.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供關(guān)鍵詞優(yōu)化、網(wǎng)站收錄搜索引擎優(yōu)化、網(wǎng)站制作、網(wǎng)站內(nèi)鏈、網(wǎng)站維護(hù)

廣告

聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請(qǐng)盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場(chǎng),如需處理請(qǐng)聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來(lái)源: 創(chuàng)新互聯(lián)

外貿(mào)網(wǎng)站制作
午夜激情在线观看国产| 国产一区精品在线免费看| 国产成人亚洲合色婷婷| 最新人妻少妇精品中文字幕视频| 亚洲国产精品自拍视频| 在线看黄视频网站永久免费| 亚洲精品成人午夜av| 中文字幕不卡一区在线| 啄木乌法国一区二区三区| 国产成十人十综合十亚洲| 久久精品一区二区三区乱码| 久久婷婷国产综合色啪| 日本免费播放一区二区视频| 亚州无吗一区二区三区| 国产视频在线一区二区| 人妻乱人伦中文字幕在线| 一区二区亚洲国产精品| 久久午夜视频在线观看| 中文字幕一区二区三天| 国产成人综合亚洲一区| av免费在线观看麻豆| 亚洲综合久久国产一区二区| 精品国产一区二区日韩91| 国产美女口爆吞精久久| sedoge在线播放免费有码| av毛片高清在线观看| 亚洲最色最黄大片在线视频| 亚洲午夜一区二区三区精品影院| 麻豆视传媒官网免费观看| 激情婷婷亚洲五月综合网| 91色综合久久久久婷婷| 国产精品色呦呦一区二区| 亚洲精品在线免费av| 亚洲av毛片一区二区三区网| 日韩在线国产精品视频| 色呦呦视频在线免费看| 国产三级视频在线观看视频 | 欧美在线免费一级黄片| 成人激情视频在线网页| 亚洲天堂成人av蜜臀| 国产亚洲中文久久网久久|