当前位置：首页 > news >正文

【图像融合】基于非负矩阵分解分解 CNMF的高光谱和多光谱数据融合附MATLAB代码

news 来源：原创 2025/5/22 20:53:04

基于CNMF的高光谱与多光谱数据融合技术详解

一、非负矩阵分解（NMF）与约束非负矩阵分解（CNMF）的核心原理

NMF的基本概念
非负矩阵分解（NMF）是一种通过将非负矩阵分解为两个非负矩阵乘积的降维方法。给定非负矩阵 $ V \in \mathbb{R}^{m \times n} $，NMF的目标是找到基矩阵 $ W \in \mathbb{R}^{m \times r} $ 和系数矩阵 $ H \in \mathbb{R}^{r \times n} $，使得 $ V \approx WH $（$ r $ 为潜在特征数）。其核心优势在于分解结果的非负性和可解释性，例如图像处理中“局部构成整体”的特性。
CNMF的改进与约束条件
CNMF在NMF基础上引入约束条件（如类别信息、判别性约束），增强分解的物理意义和分类精度。例如：
- 半监督学习：通过部分样本的类别信息约束分解过程。
- 判别性约束：结合Fisher判别分析，提升人脸识别等任务的分类性能。
- 线性混合模型（LMM） ：在遥感数据融合中，CNMF假设高光谱图像由端元（光谱特征）和丰度（空间分布）的线性组合构成，通过交替解混高光谱（HSI）和多光谱（MSI）数据，实现融合。

二、高光谱与多光谱数据的核心区别

特征	高光谱数据（HSI）	多光谱数据（MSI）
波段数量	数百至数千个连续窄波段（1-10 nm宽度）	3-10个离散宽波段（如可见光、近红外）
光谱分辨率	高（可捕捉细微光谱差异）	低（适用于一般分类任务）
空间分辨率	通常较低（受传感器限制）	较高（适用于细节分析）
应用场景	矿物识别、精准农业、化学分析	土地利用分类、植被监测

三、基于CNMF的高光谱与多光谱融合方法

融合框架与数学模型
- 观测模型：假设HSI和MSI分别由高分辨率HSI通过空间降采样和光谱降采样生成。
- 目标函数：通过最小化重建误差，交替更新端元矩阵 $ W $ 和丰度矩阵 $ H $：
  $min_{W,H} \|X - W H\|_F^2 + \|Y - W R H\|_F^2$

其中，$ X $ 为HSI，$ Y $ 为MSI，$ R $ 为光谱响应函数。

乘性迭代更新规则：通过乘法更新保证非负性，如 $ W \leftarrow W \odot \frac{X H^T}{W H H^T} $ 。

算法流程
- 步骤1：数据预处理
  对HSI进行空间下采样生成低分辨率HSI（LR-HSI），对原始HSI进行光谱下采样生成高分辨率MSI（HR-MSI）。
- 步骤2：交替解混
  交替分解HSI和MSI，优化端元和丰度矩阵，例如：

HSI解混：固定MSI的丰度矩阵，更新HSI的端元矩阵。
MSI解混：固定HSI的端元矩阵，更新MSI的丰度矩阵。
- 步骤3：融合重构
  利用优化后的端元和丰度矩阵，重构高空间-光谱分辨率的HSI。

性能优化与改进方法
- 图正则化与稀疏约束：引入图正则项（如GLNMF）增强局部几何结构保留能力。
- 自适应参数选择：通过正则化参数（如MINC-CNMF中的 $ \alpha $）平衡光谱与空间误差。
- 深度学习结合：如CUCaNet网络自适应学习光谱响应函数，提升泛化能力。
评价指标
- 光谱保真度：光谱角制图（SAM）、均方根误差（RMSE）。
- 空间质量：峰值信噪比（PSNR）、结构相似性（SSIM）。
- 分类精度：支持向量机（SVM）分类结果对比。

四、MATLAB代码实现框架

以下为基于CaImAn-MATLAB项目和CNMF融合原理的简化代码框架：

%% 数据加载与预处理
% 加载高光谱数据（HSI）和多光谱数据（MSI）
load('paviaU.mat'); % 示例数据，包含HSI和MSI
HSI = double(paviaU_HSI); % 低分辨率HSI
MSI = double(paviaU_MSI); % 高分辨率MSI%% CNMF参数设置
params = struct();
params.n_components = 20;     % 端元数量（根据数据特性调整）
params.max_iter = 500;         % 最大迭代次数
params.tol = 1e-5;            % 收敛阈值
params.verbose = 1;            % 显示迭代信息%% CNMF融合核心步骤
% 初始化端元矩阵W和丰度矩阵H
[W_init, H_init] = initialize_nmf(HSI, params.n_components);% 交替优化W和H
for iter = 1:params.max_iter% 更新HSI端元矩阵WW = W .* (HSI * H_init') ./ (W * (H_init * H_init') + eps);% 更新MSI丰度矩阵HH = H_init .* (W' * MSI) ./ ((W' * W) * H_init + eps);% 检查收敛条件if norm(H - H_init, 'fro') < params.tolbreak;endH_init = H;
end%% 重构高分辨率HSI
HSI_fused = W * H;%% 结果可视化
figure;
subplot(1,3,1); imagesc(HSI(:,:,30)); title('原始HSI'); axis off;
subplot(1,3,2); imagesc(MSI(:,:,3)); title('输入MSI'); axis off;
subplot(1,3,3); imagesc(HSI_fused(:,:,30)); title('融合结果'); axis off;
colormap('jet'); colorbar;

五、关键函数与参数说明

initialize_nmf：初始化端元和丰度矩阵，常用方法包括随机初始化或基于SVD的初始化。
乘法更新规则：通过元素乘除保证非负性，避免负值出现。
数据分块与并行计算：使用run_CNMF_patches分块处理大规模数据，结合MATLAB并行工具箱加速。
可视化函数：imagesc显示图像矩阵，colormap调整颜色映射。

六、挑战与未来方向

计算效率：CNMF迭代过程耗时，需结合GPU加速或分布式计算。
端元数量选择：$ r $ 值需通过交叉验证或基于信息准则（如AIC）优化。
多模态融合扩展：结合深度学习（如卷积神经网络）提升复杂场景下的融合鲁棒性。

【图像融合】基于非负矩阵分解分解 CNMF的高光谱和多光谱数据融合附MATLAB代码

基于CNMF的高光谱与多光谱数据融合技术详解一、非负矩阵分解（NMF）与约束非负矩阵分解（CNMF）的核心原理 NMF的基本概念非负矩阵分解（NMF）是一种通过将非负矩阵分解为两个非负矩阵乘积的降维方法。给定非负…...

编程日记 2025/5/22 20:53:04

HarmonyOS NEXT 诗词元服务项目开发上架全流程实战（一、项目介绍及实现效果）

在当今数字化时代，如何让传统文化与现代科技相结合，成为了一个值得思考的问题。诗词作为中国传统文化的重要组成部分，承载着丰富的历史信息和文化内涵。为了让更多人了解和欣赏诗词的魅力，我们决定开发一款基于HarmonyOS NEXT的诗…...

编程日记 2025/5/15 17:57:49

线性代数与数据学习

The Functions of Deep Learning (essay from SIAM News, December 2018) Deep Learning and Neural Nets...

编程日记 2025/5/22 20:27:06

Linux中的计划任务

一次性任务功能介绍： 如果我们希望在将来的某个时间点去执行某件事件，这个事件执行完后任务就结束，那么我们就可以使用一性计划任务。而要实现这种功能，我们需要任务 atd 服务。我们先查询一下系统是否存在这个服务。查看是…...

编程日记 2025/5/22 20:04:28

【C++】线程池

C 线程池介绍什么是线程池？ 线程池（Thread Pool） 是一种并发编程模型，用于管理和复用多个线程，避免频繁创建/销毁线程的开销。它通过预创建一组线程，并将任务提交到队列中，由空闲线程自动执行…...

编程日记 2025/5/22 19:28:53

美团社招一面

美团社招一面做题 1、面试题 <style> .outer{width: 100px;background: red;height: 100px; }.inner {width: 50px;height: 50px;background: green; }</style> <div class"outer"><div class"inner"></div> </div>…...

编程日记 2025/5/22 19:27:08

C++:BST、AVL、红黑树

C:BST、AVL、红黑树二叉搜索树（BST）二叉平衡搜索树（AVL）红黑树（RBT）三者对比什么情况下使用？C 标准库中的使用总结二叉搜索树（BST） 二叉搜索树（Binary Sea…...

编程日记 2025/5/19 18:34:20

算法笔记.染色法判断二分图

题目：（来自AcWing） 给定一个 n 个点 m 条边的无向图，图中可能存在重边和自环。请你判断这个图是否是二分图。输入格式第一行包含两个整数 n 和 m。接下来 m 行，每行包含两个整数 u 和 v，表示点 u …...

编程日记 2025/4/29 19:26:38

在 IDEA 中写 Spark 程序：从入门到实践

在大数据处理领域，Apache Spark 凭借其出色的性能和丰富的功能受到广泛欢迎。而 IntelliJ IDEA 作为一款功能强大的 Java 集成开发环境，为编写 Spark 程序提供了极大的便利。本文将详细介绍如何在 IDEA 中搭建 Spark 开发环境并编写运行 Spark 程序&…...

编程日记 2025/4/29 19:20:05

[Spring] Sentinel详解

🌸个人主页:https://blog.csdn.net/2301_80050796?spm1000.2115.3001.5343 🏵️热门专栏: 🧊 Java基本语法(97平均质量分)https://blog.csdn.net/2301_80050796/category_12615970.html?spm1001.2014.3001.5482 🍕 Collection与…...

编程日记 2025/4/29 19:14:28

让数据优雅落地：用 serde::Deserialize 玩转结构体实体

前言想象一下，服务器突然飞来一堆 JSON 数据，就像一群无头苍蝇冲进办公室，嗡嗡作响，横冲直撞。此刻，你的任务，就是把这群“迷路数据”安置进正确的格子里，分门别类，秩序井然，不混不乱，不漏一只。好在 Rust 早就为我们备好瑞士军刀：serde::Deserialize。它不仅刀…...

编程日记 2025/5/15 13:43:23

【wpf】 WPF中实现动态加载图片浏览器（边滚动边加载）

WPF中实现动态加载图片浏览器（边滚动边加载） 在做图片浏览器程序时，遇到图片数量巨大的情况（如几百张、上千张），一次性加载所有图片会导致界面卡顿甚至程序崩溃。本文介绍一种 WPF Prism 实现动态分页加…...

编程日记 2025/5/2 7:43:02

Flow原理

fun main() {runBlocking {launch {flow4.collect{println("---collect-4")}println("---flow4")}}val flow4 flow<Boolean>{delay(5000)emit(false) } 我们分析下整个流程 1.flow为什么之后在collect之后才会发送数据 2.collect的调用流程我…...

编程日记 2025/5/4 13:39:30

业绩回暖、股价承压，三只松鼠赴港上市能否重构价值锚点？

在营收重返百亿俱乐部后，三只松鼠再度向资本市场发起冲击。 4月25日，这家坚果零食巨头正式向港交所递交上市申请书，若成功登陆港股，将成为国内首个实现“AH”双上市的零食品牌。其赴港背后的支撑力，显然来自近期披露…...

编程日记 2025/5/22 19:46:33

基于大模型的胆总管结石全流程预测与临床应用研究报告

目录一、引言 1.1 研究背景 1.2 研究目的与意义 1.3 研究方法和创新点二、大模型在胆总管结石预测中的应用原理 2.1 大模型概述 2.2 模型构建的数据来源与处理 2.3 模型训练与优化三、术前预测与准备 3.1 术前胆总管结石存在的预测 3.2 基于预测结果的术前检查方…...

编程日记 2025/5/22 20:07:40

QT—布局管理器之BoxLayout篇

1.布局管理器的概述在Qt中，使用布局管理器的主要原因是它能够自动管理组件的大小和位置，从而实现灵活且动态的界面布局。布局管理器可以自动调整组件以适应窗口大小的变化，确保界面在不同分辨率和设备上都能保持良好的显示效果。这不仅减少了…...

编程日记 2025/5/15 4:30:58

如何在 IntelliJ IDEA 中编写 Speak 程序

在当今数字化时代，语音交互技术越来越受到开发者的关注。如果你想在 IntelliJ IDEA（一个强大的集成开发环境）中编写一个语音交互（Speak）程序，那么本文将为你提供详细的步骤和指南。一、环境准备在开始编…...

编程日记 2025/5/2 16:47:25

湖北理元理律师事务所：债务优化的法律机制与民生实践

在债务纠纷日益增多的社会背景下，合法、规范的债务管理服务成为民生需求的重要环节。湖北理元理律师事务所作为经国家司法局注册登记的债事服务机构，以法律为工具，探索出一套覆盖债务咨询、规划与风险防控的服务体系。 1.法律服务的专业化框…...

编程日记 2025/5/21 3:51:10

练习普通话，说话更有节奏

玲珑塔，塔玲珑，玲珑宝塔第一层，一张高桌四条腿， 一个和尚一本经。一个铙钹一口磬，一个木鱼一盏灯。一个金玲，整四两，风儿一刮响哗愣。玲珑塔，隔过两层数三层，三张高桌十…...

编程日记 2025/5/22 20:16:46

链表相关——Python实现

一、链表的创建及数据插入示例代码： #1.定义一个结点类 class ListNode():def __init__(self,x,nextNone):self.valxself.nextnext #2.定义单链表 class LinkList():#2.1 创建一个头指针为空的链表def __init__(self,headNone):self.headNone#2.2 将数据插入链表…...

编程日记 2025/5/22 19:51:35

[OS_9] C 标准库和实现 | musl libc | offset

在你感觉有困难的时候，计算机一定有解决办法操作系统为我们提供了对象和操作它们的 API：我们学习了进程管理的 fork, execve, exit, waitpid；内存管理的 mmap；文件 (对象) 管理的 open, read, write, dup, close, pipe, …… 大…...

编程日记 2025/5/13 10:25:33

【氮化镓】质子辐照对 GaN-on-GaN PiN 二极管电导调制的影响

2025 年，中国科学技术大学的 Xuan Xie 等人采用实验研究的方法，深入探究了 10-MeV 和 50-MeV 质子辐照（ fluence 最高达 11014 cm−2）对 kV 级垂直 GaN-on-GaN PiN 二极管的导电调制影响。研究背景在于空间应用中的功率电子电路易受质子、α 粒子和重离子等影响而降级甚至失…...

编程日记 2025/5/22 19:28:26

深入浅出限流算法（一）：简单但有“坑”的固定窗口计数器

在现代分布式系统和 API 设计中，限流 (Rate Limiting) 是一个不可或缺的环节。它像一个尽职的门卫，保护着我们的服务资源，防止因突发流量或恶意攻击导致系统过载，同时也能确保资源的公平分配。实现限流的算法多种多样&#xff0…...

编程日记 2025/5/20 11:57:23

项目上线流程梳理（Linux宝塔面板）

项目部署（Linux宝塔面板） 一、准备工作 1、后端项目梳理配置添加application-prod.yml使用maven进行打包生成jar包 2、前端项目修改request.ts请求的后端端口（服务器地址）打包二、服务器 1、服务器环境安装 2、初始化数…...

编程日记 2025/5/18 8:12:59

MCP Servers快速了解一文搞懂 MCP Servers mcp server网站 https://mcpservers.orghttps://mcp.sohttps://glama.ai/mcp/servershttps://www.pulsemcp.com 一、安装&配置mcp clients mcp clients可以使用客户端（常见claude desktop，但claude需…...

编程日记 2025/5/22 1:02:49

永磁同步电机控制算法-转速环电流环SMC控制器

一、原理介绍为改进传统PI转速环电流环控制器易超调、抗干扰性能差的问题，转速环采用一阶滑模控制器，电流环采用二阶滑模控制器。二、仿真验证在MATLAB/simulink里面验证所提算法，采用和实验中一致的控制周期1e-4，电机部分计…...

编程日记 2025/5/21 11:17:31

Nginx支持HTTP2/HTTP3的并用CURL测试

对比HTTP2和HTTP3 HTTP/3 相比 HTTP/2 的主要优势： 项目HTTP/2HTTP/3传输层基于 TCP（ TLS）基于 QUIC（内置 TLS）建立连接速度慢，至少 2~3 次握手（TCP TLS）快，只需 1 次握…...

编程日记 2025/5/22 16:23:06

阅读MySQL实战45讲第11天

目录引言： 基本原理排序方式索引排序文件排序（File Sort） 优化建议一、全字段排序二、rowid 排序如果 MySQL 认为排序的单行长度太大会怎么做呢？ 1. 使用磁盘临时表 2. 分阶段排序 3. 产生警告或错误三、全字段排序 VS ro…...

编程日记 2025/5/21 10:59:44

linux 使用nginx部署vue、react项目

前言本文基于：操作系统 CentOS Stream 8 使用工具：Xshell8、Xftp8 1.安装依赖安装gcc、gcc-c yum install gcc gcc-c -y安装pcre、pcre-devel yum install pcre pcre-devel -y安装zlib、zlib-devel yum install zlib zlib-devel -y安装openssl、…...

编程日记 2025/4/29 17:15:03

逆向设计——CWDM_splitter

一、创建结构 switchtolayout; selectall; delete;## SIM PARAMS num_wg = 4; wg_width = 0.5e-6; out_wg_dist = 1e-6; mode_width = 3*wg_width; total_wg_h = num_wg*wg_width + (num_wg-1)*out_wg_dist;opt_size_x=6e-6; opt_size_y=6e-6;size_x=opt_size_x+1e-6; size_y=o…...

编程日记 2025/4/29 17:13:57