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

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

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

python的數(shù)據(jù)類型：1. 數(shù)字類型，包括int（整型）、long（長整型）和float（浮點型）。2.字符串，分別是str類型和unicode類型。3.布爾型，Python布爾類型也是用于邏輯運算，有兩個值：True（真）和False（假）。4.列表，列表是Python中使用最頻繁的數(shù)據(jù)類型，集合中可以放任何數(shù)據(jù)類型。5. 元組，元組用”()”標識，內(nèi)部元素用逗號隔開。6. 字典，字典是一種鍵值對的集合。7. 集合，集合是一個無序的、不重復的數(shù)據(jù)組合。

使用GDAL修復影像條帶的代碼如下：

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

修復之后的效果圖如下所示：

Opencv修復Landsat ETM+影像條帶

使用opencv修復影像的代碼如下：

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位無符號整型和32位浮點型1通道或者8位無符號3通道
    # inpaintMask：掩膜，8位無符號整型
    # dst：和源圖像具有一樣大小的輸出
    # inpaintRadius：算法考慮的每個已修復點的圓形鄰域的半徑     # flags：修復算法類型，可選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)

修復之后的結果圖：

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

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

當前名稱：使用Python怎么修復遙感影像條帶-創(chuàng)新互聯(lián)
文章轉載：http://chinadenli.net/article6/pdiog.html

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

廣告

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