使用python将Quickbird 影像作为参考影像, 用二次多项式实现高景(SV1)影像的相对配准

时间: 2024-03-04 22:48:52 浏览: 16
好的,下面是使用Python实现将Quickbird影像作为参考影像,用二次多项式实现高景(SV1)影像的相对配准的代码示例: ```python import cv2 import numpy as np # 读取Quickbird影像和高景(SV1)影像 quickbird = cv2.imread('quickbird.tif', cv2.IMREAD_GRAYSCALE) sv1 = cv2.imread('sv1.tif', cv2.IMREAD_GRAYSCALE) # 提取Quickbird影像和高景(SV1)影像的特征点 orb = cv2.ORB_create() kp1, des1 = orb.detectAndCompute(quickbird, None) kp2, des2 = orb.detectAndCompute(sv1, None) # 匹配特征点 bf = cv2.BFMatcher(cv2.NORM_HAMMING, crossCheck=True) matches = bf.match(des1, des2) # 排序匹配结果 matches = sorted(matches, key=lambda x: x.distance) # 取前50个匹配结果 good_matches = matches[:50] # 提取匹配点的坐标 src_pts = np.float32([kp1[m.queryIdx].pt for m in good_matches]).reshape(-1, 1, 2) dst_pts = np.float32([kp2[m.trainIdx].pt for m in good_matches]).reshape(-1, 1, 2) # 使用二次多项式变换模型进行配准 M, mask = cv2.findHomography(src_pts, dst_pts, cv2.RANSAC, 5.0) # 对高景(SV1)影像进行变换 result = cv2.warpPerspective(sv1, M, (quickbird.shape[1], quickbird.shape[0])) # 显示配准结果 cv2.imshow('Quickbird', quickbird) cv2.imshow('SV1', sv1) cv2.imshow('Result', result) cv2.waitKey(0) cv2.destroyAllWindows() ``` 在这个示例中,我们使用了OpenCV库中的ORB特征检测器和BFMatcher特征匹配器来提取Quickbird影像和高景(SV1)影像的特征点,并进行匹配。然后使用cv2.findHomography函数计算二次多项式变换参数,将高景(SV1)影像变换到与Quickbird影像相对应的位置,并显示配准结果。

相关推荐

优化这段代码% 读取第一组数据 imgfilename1 = 'C:\Users\86182\Desktop\tif\QB2013.dat'; data1 = read_data(imgfilename1); % 读取第二组数据 imgfilename2 = 'C:\Users\86182\Desktop\tif\SV2018.dat'; data2 = read_data(imgfilename2); % 显示Quickbird影像 im1 = data1(:, :, 2:4); im1 = uint8(im1); show_image(im1, 'Quickbird影像432波段显示', 1); % 显示高景影像 im2 = data2(:, :, 2:4); im2 = uint8(im2); show_image(im2, '高景影像432波段显示', 2); %NDVI计算结果图 X1=data1; im3=ndvi(X1); show_image(im3, '2013年(QB)NDVI专题图', 3); X2=data2; im4=ndvi(X2); show_image(im4, '2018年(SV)NDVI专题图', 4); function data = read_data(filename) %读取数据 data = multibandread(filename, [1989, 2126, 4], 'int16', 0, 'bsq', 'ieee-le', { 'Band', 'Direct', [1 4 3 2]}); %lines,samples,bands,hdr文件里查看 % 调整波段排列 %B=data(:,:,1); %G=data(:,:,2); %R=data(:,:,3); %C=data(:,:,4); %data=cat(4,B,C,R,G); % 将数据转换为0-255的整型用于显示 data_unit8 = uint8(data); for k = 1:4 data_k = double(data(:, :, k)); min_val = min(data_k, [], 'all'); max_val = max(data_k, [], 'all'); data_unit8(:, :, k) = uint8((data_k - min_val) / (max_val - min_val) * 255); end % 返回处理后的数据 data = data_unit8; end %ndvi计算 function data_ndvi(X) NIR = double(X(:,:,2));%近红外光谮带 red = double(X(:,:,3));%可见光红色光谱带 data_ndvi(:, :, k)=ndvi((NIR - red) ./ (NIR + red));%归一化 X=data_ndvi; end % 影像显示 function show_image(im, title_str, fig_num) figure(fig_num); imshow(im, [min(im(:)), max(im(:))]); title(title_str); end

最新推荐

recommend-type

ENVI4.2 DEM Extraction 模块简介

DEM模型是表示地形表面高程的规则栅格...ENVI4.2具有从立体像对进行DEM提取的最新工具,能够从ASTER,IKONOS,OrbView-3,QuickBird,以及SPOT卫星影像提取DEM。 DEM Extraction 向导指导您从含有RPC的立体像对中提取DEM。
recommend-type

高级色系PPT11.pptx

高级色系PPT11.pptx
recommend-type

node-v7.9.0-linux-x86.tar.xz

Node.js,简称Node,是一个开源且跨平台的JavaScript运行时环境,它允许在浏览器外运行JavaScript代码。Node.js于2009年由Ryan Dahl创立,旨在创建高性能的Web服务器和网络应用程序。它基于Google Chrome的V8 JavaScript引擎,可以在Windows、Linux、Unix、Mac OS X等操作系统上运行。 Node.js的特点之一是事件驱动和非阻塞I/O模型,这使得它非常适合处理大量并发连接,从而在构建实时应用程序如在线游戏、聊天应用以及实时通讯服务时表现卓越。此外,Node.js使用了模块化的架构,通过npm(Node package manager,Node包管理器),社区成员可以共享和复用代码,极大地促进了Node.js生态系统的发展和扩张。 Node.js不仅用于服务器端开发。随着技术的发展,它也被用于构建工具链、开发桌面应用程序、物联网设备等。Node.js能够处理文件系统、操作数据库、处理网络请求等,因此,开发者可以用JavaScript编写全栈应用程序,这一点大大提高了开发效率和便捷性。 在实践中,许多大型企业和组织已经采用Node.js作为其Web应用程序的开发平台,如Netflix、PayPal和Walmart等。它们利用Node.js提高了应用性能,简化了开发流程,并且能更快地响应市场需求。
recommend-type

基于tensorflow的的cnn卷积神经网络的图像识别分类

【作品名称】:基于tensorflow的的cnn卷积神经网络的图像识别分类 【适用人群】:适用于希望学习不同技术领域的小白或进阶学习者。可作为毕设项目、课程设计、大作业、工程实训或初期项目立项。
recommend-type

### 数据分析概念、使用技巧、优缺点的文章

数据分析是指通过收集、清洗、处理和解释数据,以发现其中的模式、趋势和关联,从而提供决策支持或洞察见解的过程。它在各行各业中都扮演着至关重要的角色,从市场营销到科学研究,从金融领域到医疗保健,都有广泛的应用。
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

实现实时数据湖架构:Kafka与Hive集成

![实现实时数据湖架构:Kafka与Hive集成](https://img-blog.csdnimg.cn/img_convert/10eb2e6972b3b6086286fc64c0b3ee41.jpeg) # 1. 实时数据湖架构概述** 实时数据湖是一种现代数据管理架构,它允许企业以低延迟的方式收集、存储和处理大量数据。与传统数据仓库不同,实时数据湖不依赖于预先定义的模式,而是采用灵活的架构,可以处理各种数据类型和格式。这种架构为企业提供了以下优势: - **实时洞察:**实时数据湖允许企业访问最新的数据,从而做出更明智的决策。 - **数据民主化:**实时数据湖使各种利益相关者都可
recommend-type

如何用python编写api接口

在Python中编写API接口可以使用多种框架,其中比较流行的有Flask和Django。这里以Flask框架为例,简单介绍如何编写API接口。 1. 安装Flask框架 使用pip命令安装Flask框架: ``` pip install flask ``` 2. 编写API接口 创建一个Python文件,例如app.py,编写以下代码: ```python from flask import Flask, jsonify app = Flask(__name__) @app.route('/api/hello', methods=['GET']) def hello():
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。