matlab二维曲面散点图插值方法

在 MATLAB 中,你可以使用以下函数进行二维曲面散点插值:

  1. griddata: 该函数可以在散点数据上进行二维插值,生成平滑的曲面。它支持多种插值方法,包括三次样条插值、最近邻插值、线性插值和自然邻近法插值。

  2. scatteredInterpolant: 这是一个基于散点数据的插值类,可以用于生成二维曲面。你可以通过指定插值方法(如三次样条插值、线性插值等)来创建一个 scatteredInterpolant 对象,并使用该对象进行插值计算。

  3. TriScatteredInterp: 这是一个基于三角剖分的插值类,用于在二维散点数据上进行插值。它将散点数据进行三角剖分,并使用插值方法对每个三角形进行插值计算,从而生成二维曲面。

griddata使用方法如下:

data = xlsread('程序参考2.xlsx','Sheet1');

x = data(:,1);
y = data(:,2);
z = data(:,3);

% 定义插值点的网格
n_points = 100; % 插值点个数

xi = linspace(min(x), max(x), n_points); % x 坐标范围
yi = linspace(min(y), max(y), n_points); % y 坐标范围
[XI, YI] = meshgrid(xi, yi); % 插值点的网格

% 二维插值
ZI = griddata(x, y, z, XI, YI);

% 绘制原始数据和插值结果
scatter3(x, y,  z, 'r');
hold on;
mesh(XI, YI, ZI);

使用scatteredInterpolant

% xyz 得到方法同前
% 创建插值函数对象
F = scatteredInterpolant(x, y, z);

% 对插值点进行插值计算
ZI = F(XI, YI);

% 绘制插值结果
figure(2)
mesh(XI, YI, ZI);
xlabel('X');
ylabel('Y');
zlabel('Z');
title('二维曲面散点插值');

 使用三角剖分

% xyz 得到方法同前
% 创建插值函数对象
F = TriScatteredInterp(x, y, z);

% 对插值点进行插值计算
ZI = F(XI, YI);

% 绘制插值结果
figure(3)
surf(XI, YI, ZI);
xlabel('X');
ylabel('Y');
zlabel('Z');
title('二维曲面散点插值');

结论:

griddata 是一个函数,而 scatteredInterpolant 是一个类。

  • griddata 函数需要提供散点数据的坐标和值,以及插值点的坐标,然后在插值点上生成插值结果。

  • scatteredInterpolant 类创建一个插值函数对象,需要提供散点数据的坐标和值,然后可以使用该对象对新的点进行插值计算。

  • griddata核心会调用scatteredInterpolant,本质应该没有区别

线性三角剖分插值法

https://blog.csdn.net/u012839396/article/details/48174901

这个博主解释的不错,主要流程是:

步骤:

1、根据离散坐标,参考凸包算法,获得离散点的凸包;

2、构件Delaunary三角网;

3、根据格网化方法,获取格网点的坐标:

 可以根据离散点控制的范围,然后纵横坐标的格网间隔进行格网化,也可以根据离散点的控制范围,然后纵横坐标的节点数进行格网化,根据需要来确定。

4、因为要求待插值点在Delaunary三角形内,所以,遍历插值区域内所有的格网点,根据插值点所在的三角形内的三个点,使用三次方程内插出待定点坐标上的值。
 

用matlab做的,有一个项目可以用来参考

https://www.mathworks.com/matlabcentral/fileexchange/38925-linearly-interpolate-triangulation

 这个地方可以注释掉,感觉不一定非要限 是三角形

简单写个调用例子

%x y z数据同前

% 构建三角剖分
DT = delaunayTriangulation(x, y, z);

% Get the connectivity table
tri = DT.ConnectivityList;
tri = tri(:, [1, 2, 3]);

ZI=interptri(tri,x,y,z,XI,YI);

 % 绘制插值结果
figure(1)
surf(XI, YI, ZI);
xlabel('X');
ylabel('Y');
zlabel('Z');
title('二维曲面散点插值');

效果还可以

https://blog.csdn.net/u012839396/article/details/48174901

C++要使用估计得先建立库,

https://jingyan.baidu.com/article/59a015e30c5ab3f794886520.html

https://blog.csdn.net/sikong00/article/details/108754058?spm=1001.2101.3001.6650.1&utm_medium=distribute.pc_relevant.none-task-blog-2%7Edefault%7ECTRLIST%7ERate-1-108754058-blog-102868479.235%5Ev38%5Epc_relevant_anti_vip&depth_1-utm_source=distribute.pc_relevant.none-task-blog-2%7Edefault%7ECTRLIST%7ERate-1-108754058-blog-102868479.235%5Ev38%5Epc_relevant_anti_vip&utm_relevant_index=2

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/145388.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

update与merge基本语法规则、异同点、使用场景

– UPDATE更新数据SQL UPDATE /append parallel(A,8) nologging/ CRSMAK.M_CM_GS_BILL_DISH_INFO A --票据贴现信息 SET A.YEAR_PAR_AMT (SELECT /parallel(B,8)/ YEAR_PAR_AMT FROM (SELECT BIZ_ID, SUM(NVL(PAR_AMT, 0)) YEAR_PAR_AMT --年累计票面金额 FROM CRSMAK.M_CM_G…

存钱虚拟计划,嘚

存钱计划—虚拟 2024年 (第一年) 1月 2月 3月 4月 5月 6 月 7月 8月 9月 10月 11月 12月 1w 1w 1w 1w 1w 1w 1w 1w 1w 1w 1w 1w 2025年(第二年) 1w1w*102.5% A 懒得写A代替 A A A A A A A A A A 2026年(第三年) (1w1w*10…

Centos7.9用rancher来快速部署K8S

什么是 Rancher? Rancher 是一个 Kubernetes 管理工具,让你能在任何地方和任何提供商上部署和运行集群。 Rancher 可以创建来自 Kubernetes 托管服务提供商的集群,创建节点并安装 Kubernetes,或者导入在任何地方运行的现有 Kube…

React路由跳转但是页面未发生改变,刷新后正常显示

react项目,配置好路由后,触发跳转发现url已经改变但页面内容未改变,刷新后才能正常显示。 发现是被react的严格模式影响,将BrwserRouter放到其外层即可 const root ReactDOM.createRoot(document.getElementById("root&qu…

OpenCV入门5——OpenCV的算术与位运算

文章目录 图像的加法运算图像的减法运算图像的乘除运算图像的融合OpenCV位运算-非操作OpenCV位操作-与运算OpenCV位操作-或与异或为图像添加水印 图像的加法运算 # -*- coding: utf-8 -*- import cv2 import numpy as npimg cv2.imread(E://pic//4.jpg)# 图的加法运算就是矩阵…

EasyCVR视频监控+AI智能分析网关如何助力木材厂安全生产?

旭帆科技有很多工厂的视频监管方案,小编也经常分享出来供大家参考。近期,又有伙伴后台私信我们想要关于木材厂的方案。针对木材厂的生产过程与特性以及安全风险等,我们来分享一下相关的监管方案: 1)温湿度监测&#xf…

技巧篇:Mac 环境PyCharm 配置 python Anaconda

Mac 中 PyCharm 配置 python Anaconda环境 在 python 开发中我们最常用的IDE就是PyCharm,有关PyCharm的优点这里就不在赘述。在项目开发中我们经常用到许多第三方库,用的最多的命令就是pip install 第三方库名 进行安装。现在你可以使用一个工具来帮你解…

three.js相机调用reset方法复原但无法完全复原

目录 一、问题 二、原因 三、总结 一、问题 需求:经过多次旋转、平移等变换后,希望恢复到初次渲染的角度、位置。结果发现: three.js oribtControls调用 一次reset方法 无法完全复原;需要多次调用后才能复原。 this.orbitControls.reset…

前端开发学习 (一) 搭建Vue基础环境

一、环境搭建 1、安装nodejs #下载地址 https://nodejs.org/dist/v20.9.0/node-v20.9.0-x64.msi 2、配置环境变量 上面下载完安装包后自行安装,安装完成后安装下图操作添加环境变量 #查看版本 node --version v20.9.0# npm --version 10.1.03、配置npm加速源 np…

2.c++基础语法

文章目录 1.c 程序结构关键字标识符、操作符、标点预处理指令注释main 主函数命名空间 2.c 变量和常量变量 3.c 数组和容器4.c 程序流程5.c字符和字符串 1.c 程序结构 关键字 关键字事程序保留的,程序员不能使用,c的常见关键字如下图: 标识…

机器人导航+OPENCV透视变换示例代码

透视变换又称四点变换,所以不能用于5边形这样的图形变换,不是真正的透视变换,但是这个方法可以把机器人看到的图像转换为俯视图,这样就可以建立地图,要不然怎么建立地图呢。 void CrelaxMyFriendDlg::OnBnClickedOk()…

Python如何调用ixchariot进行吞吐量测试

Python如何调用ixchariot进行吞吐量测试 要使用Python调用IxChariot进行吞吐量测试,您可以使用 subprocess 模块来执行IxChariot的TCL命令行。下面是一个简单的示例代码: import subprocess# 定义IxChariot的安装路径和测试脚本路径 ixchariot_path &q…

Ubuntu 18.04无网络连接的n种可能办法

文章目录 网络图标消失,Ubuntu无网络连接VMware上Ubuntu18.04,桥接了多个网卡,其中一个用来上网,均设置为静态ip网络桥接链路没有接对路由不对 网络图标消失,Ubuntu无网络连接 sudo service network-manager stop sud…

增删改查mysql

查询 -- 查询表结果-- 查看 当前数据库下的表show tables;-- 查看指定的表desc tb_emp; -- td_emp 是表名-- 查看 数据库的见表语句show create table tb_emp; 修改 -- 修改表结构 -- 修改 为表 tb_emp 添加字段 qq varchar(11) alter table tb_emp add qq varchar(11) …

【LeetCode刷题-滑动窗口】--340.至多包含K个不同字符的最长子串

340.至多包含K个不同字符的最长子串 class Solution {public int lengthOfLongestSubstringKDistinct(String s, int k) {int len s.length();if(len < k){return len;}//滑动窗口的左右指针int left 0,right 0;//定义一个哈希映射HashMap<Character,Integer> hash…

2023-11-18 mysql-sysbench压测TPS/QPS-记录

摘要: 2023-11-18 mysql-sysbench压测TPS/QPS sysbench压测TPS/QPS 使用 sysbench 压测 TPS/QPS 的基础测试数据是 6 张表,每张表写入 10 万行数据,然后在此基础测试数据上分别进行 300 秒的混合写、更新、删除压测。 sysbench --mysql-host=192.168.41.xxx --mysql-port=3…

java中的深度复制和浅复制的BUG

刷题刷到LeetCode回溯DFS的算法题39题的时候,碰见一个Arraylist里面的bug,其中dfs函数里面的第一个if判断里面的语句 paths.add(path); path.clear();其中path是添加了path,但是添加之后path.clear(),导致原来添加到paths的path置为空数组,因为ArrayList的add只是把一个引用指…

使用Postman进行压力测试

1.打开Postman新建测试接口 2.点击右边保存&#xff0c;选择一个文件集合&#xff0c;如果没有就创建&#xff0c;然后保存 就是这个东西&#xff0c;这里不便展示出来&#xff0c;压力测试需要在文件夹里面进行 3.选择要测试的接口&#xff0c;iterations 表示请求发起次数&a…

基于Vue+SpringBoot的大学计算机课程管理平台 开源项目

目录 一、摘要1.1 项目介绍1.2 项目录屏 二、功能模块2.1 实验课程档案模块2.2 实验资源模块2.3 学生实验模块 三、系统设计3.1 用例设计3.2 数据库设计3.2.1 实验课程档案表3.2.2 实验资源表3.2.3 学生实验表 四、系统展示五、核心代码5.1 一键生成实验5.2 提交实验5.3 批阅实…