使用Python怎么修復(fù)遙感影像條帶-創(chuàng)新互聯(lián)

今天就跟大家聊聊有關(guān)使用Python怎么修復(fù)遙感影像條帶,可能很多人都不太了解,為了讓大家更加了解,小編給大家總結(jié)了以下內(nèi)容,希望大家根據(jù)這篇文章可以有所收獲。

創(chuàng)新互聯(lián)公司主要從事成都網(wǎng)站制作、成都做網(wǎng)站、網(wǎng)頁(yè)設(shè)計(jì)、企業(yè)做網(wǎng)站、公司建網(wǎng)站等業(yè)務(wù)。立足成都服務(wù)達(dá)拉特,十年網(wǎng)站建設(shè)經(jīng)驗(yàn),價(jià)格優(yōu)惠、服務(wù)專(zhuān)業(yè),歡迎來(lái)電咨詢(xún)建站服務(wù):13518219792

python的數(shù)據(jù)類(lèi)型有哪些?

python的數(shù)據(jù)類(lèi)型:1. 數(shù)字類(lèi)型,包括int(整型)、long(長(zhǎng)整型)和float(浮點(diǎn)型)。2.字符串,分別是str類(lèi)型和unicode類(lèi)型。3.布爾型,Python布爾類(lèi)型也是用于邏輯運(yùn)算,有兩個(gè)值:True(真)和False(假)。4.列表,列表是Python中使用最頻繁的數(shù)據(jù)類(lèi)型,集合中可以放任何數(shù)據(jù)類(lèi)型。5. 元組,元組用”()”標(biāo)識(shí),內(nèi)部元素用逗號(hào)隔開(kāi)。6. 字典,字典是一種鍵值對(duì)的集合。7. 集合,集合是一個(gè)無(wú)序的、不重復(fù)的數(shù)據(jù)組合。

使用GDAL修復(fù)影像條帶的代碼如下:

def gdal_repair(tif_name, out_name, bands):
  """
    tif_name(string): 源影像名
    out_name(string): 輸出影像名
    bands(integer): 影像波段數(shù)
  """
  # 打開(kāi)影像文件
  tif = gdal.Open(tif_name)
  
  # 根據(jù)文件類(lèi)型獲取對(duì)應(yīng)的驅(qū)動(dòng)程序
  driver = gdal.GetDriverByName('GTiff')
  
  # 根據(jù)指定文件的驅(qū)動(dòng)程序,使用現(xiàn)有數(shù)據(jù)集創(chuàng)建新的可寫(xiě)數(shù)據(jù)集
  # 所有支持創(chuàng)建新文件的驅(qū)動(dòng)程序都支持該`CreateCopy()`方法,   # 但僅`Create()`部分支持該方法
  # CreateCopy的第一個(gè)參數(shù)為目標(biāo)文件名,第二個(gè)參數(shù)為源數(shù)據(jù)集
  # 第三個(gè)參數(shù)的值是`0`或`1`,值是`0`。即使無(wú)法將原始數(shù)據(jù)準(zhǔn)確地轉(zhuǎn)換為目標(biāo)數(shù)據(jù),程序仍將執(zhí)行
  new_img = driver.CreateCopy(out_name, tif, 0)
 
  for i in tqdm(range(1, bands)):
    # 分別對(duì)每個(gè)波段處理
    band = new_img.GetRasterBand(i)
    
    # 使用FillNodata對(duì)條帶部分進(jìn)行插值
    gdal.FillNodata(targetBand = band, maskBand = band, maxSearchDist = 15, smoothingIterations=0)
    
    # 將修復(fù)好的波段寫(xiě)入新數(shù)據(jù)集中
    new_img.GetRasterBand(i).WriteArray(band.ReadAsArray())

修復(fù)之后的效果圖如下所示:

使用Python怎么修復(fù)遙感影像條帶

Opencv修復(fù)Landsat ETM+影像條帶

使用opencv修復(fù)影像的代碼如下:

def cv2_repair(tif_name):
  # 讀取tif影像
  tif_data = gdal_array.LoadFile(tif_name).astype('float32')

  # 獲取掩膜
  mask = tif_data.sum(axis=0)
  mask = (mask == 0).astype(np.uint8)
  
  bands = tif_data.shape[0]

  res = []
  for i in tqdm(range(bands)):
    # cv.Inpaint(src, inpaintMask, dst, inpaintRadius, flags)
    # src:源圖像,可以是8位、16位無(wú)符號(hào)整型和32位浮點(diǎn)型1通道或者8位無(wú)符號(hào)3通道
    # inpaintMask:掩膜,8位無(wú)符號(hào)整型
    # dst:和源圖像具有一樣大小的輸出
    # inpaintRadius:算法考慮的每個(gè)已修復(fù)點(diǎn)的圓形鄰域的半徑     # flags:修復(fù)算法類(lèi)型,可選cv2.INPAINT_NS和cv2.INPAINT_TELEA
    
    repaired = cv2.inpaint(tif_data[i], mask, 3, flags=cv2.INPAINT_TELEA)
    res.append(repaired)

  return np.array(res)

修復(fù)之后的結(jié)果圖:

使用Python怎么修復(fù)遙感影像條帶

看完上述內(nèi)容,你們對(duì)使用Python怎么修復(fù)遙感影像條帶有進(jìn)一步的了解嗎?如果還想了解更多知識(shí)或者相關(guān)內(nèi)容,請(qǐng)關(guān)注創(chuàng)新互聯(lián)成都網(wǎng)站設(shè)計(jì)公司行業(yè)資訊頻道,感謝大家的支持。

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

當(dāng)前名稱(chēng):使用Python怎么修復(fù)遙感影像條帶-創(chuàng)新互聯(lián)
文章轉(zhuǎn)載:http://muchs.cn/article6/pdiog.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供定制網(wǎng)站、用戶(hù)體驗(yàn)、App開(kāi)發(fā)、網(wǎng)站制作、商城網(wǎng)站品牌網(wǎng)站制作

廣告

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

手機(jī)網(wǎng)站建設(shè)