（Python）Gdal与Opencv遥感影像Siticher拼接

最新推荐文章于 2024-04-19 09:38:41 发布

卫少东

最新推荐文章于 2024-04-19 09:38:41 发布

阅读量1.5k

点赞数 1

分类专栏：摄影测量文章标签： python opencv 图像拼接遥感影像

本文链接：https://blog.csdn.net/weishaodong/article/details/119538309

版权

摄影测量专栏收录该内容

11 篇文章 10 订阅

订阅专栏

文章目录

前言
使用Opencv的Stitcher(傻瓜式)
TODO：自己实现
运行结果

前言

“ 很多时候，我们需要研究较大空间范围的研究区域时，单幅遥感图无法覆盖整个研究区域，这时就需要对遥感图像进行镶嵌（或拼接）。”
本文使用Python+Opencv+GDAL实现相邻两个遥感图像之间重叠区域的拼接。
基本步骤为：

对每幅图进行特征点提取
对特征点进行匹配
进行图像配准
把图像拷贝到另一幅图像的特定位置
对重叠边界进行特殊处理（颜色融合）
最近使用python进行一些的遥感影像处理，发现摄影测量方面相比于CV领域资料非常匮乏，特此记录一些小工具的实现，也为遥感领域的开源做出自己的一些贡献，水平有限，能用而已。

使用Opencv的Stitcher(傻瓜式)

import numpy as np
import cv2
from time import *
def normalization(data):
    _range = np.max(data) - np.min(data)
    return (data - np.min(data)) / _range
def Tiff16to8bit(img_16):
    if (np.max(img_16) - np.min(img_16) != 0):
        img_nrm = normalization(img_16)
        img_8 = np.uint8(255 * img_nrm)
    return img_8

#读入图像
img_b = cv2.imread('D:\\SongshanZY3\\ortho\\BWDSC.tif', cv2.IMREAD_LOAD_GDAL)
img_f = cv2.imread('D:\\SongshanZY3\\ortho\\FWDSC.tif', cv2.IMREAD_LOAD_GDAL)
#取一小块出来
x_offset,y_offset = 5000,6000
block_size_x = 3000
block_size_y = 3000
block_B = img_b[y_offset:y_offset + block_size_y,x_offset:x_offset + block_size_x]
block_F = img_f[y_offset+1000:y_offset+1000 + block_size_y,x_offset:x_offset + block_size_x]
block_B = Tiff16to8bit(block_B)
block_F = Tiff16to8bit(block_F)

#构建拼接器
stitcher = cv2.Stitcher.create(cv2.Stitcher_PANORAMA)#根据不同的OpenCV版本来调用
#需要RGB格式图
i1 = cv2.cvtColor(block_B,cv2.COLOR_GRAY2RGB)
i2 = cv2.cvtColor(block_F,cv2.COLOR_GRAY2RGB)

begin_time = time()
#key CODE
(_result, pano) = stitcher.stitch((i1, i2))
end_time = time()
run_time = end_time-begin_time
print ('运行时间：',run_time,'s') #该循环程序运行时间： 1.4201874732

cv2.namedWindow('pano',0)
cv2.imshow('pano', pano)
cv2.imwrite('pano.png', pano)
cv2.waitKey(0)