matlab中如何提取等高线,在Python或MATLAB中从等高线图中提取数据

这是一个小型的Matlab脚本,可以完成这项工作(使用一些GUI,在图的斜角处读取guidlines):

%// Import the data:

imdata = importdata('your_picture_file');

Gray = rgb2gray(imdata.cdata);

colorLim = [-1 1]; %// this should be set manually

%// Get the area of the data:

f = figure('Position',get(0,'ScreenSize'));

imshow(imdata.cdata,'Parent',axes('Parent',f),'InitialMagnification','fit');

%// Get the area of the data:

title('Click with the cross on the most top left area of the *data*')

da_tp_lft = round(getPosition(impoint));

title('Click with the cross on the most bottom right area of the *data*')

da_btm_rgt = round(getPosition(impoint));

dat_area = double(Gray(da_tp_lft(2):da_btm_rgt(2),da_tp_lft(1):da_btm_rgt(1)));

%// Get the area of the colorbar:

title('Click with the cross within the upper most color of the *colorbar*')

ca_tp_lft = round(getPosition(impoint));

title('Click with the cross within the bottom most color of the *colorbar*')

ca_btm_rgt = round(getPosition(impoint));

cmap_area = double(Gray(ca_tp_lft(2):ca_btm_rgt(2),ca_tp_lft(1):ca_btm_rgt(1)));

close(f)

%// Convert the colormap to data:

data = dat_area./max(cmap_area(:)).*range(colorLim)-abs(min(colorLim));

现在, data 就是你要找的 .

以下是使用问题中的图形输出的说明:

4439e398-f122-467f-b042-61960eddeb49.png

插图的代码:

figure('Position',[100 200 1200 400]);

subplot 121

imshow(imdata.cdata)

hold on

plot(da_tp_lft(1),da_tp_lft(2),'m+','MarkerSize',7,'LineWidth',2)

plot(da_btm_rgt(1),da_btm_rgt(2),'m+','MarkerSize',7,'LineWidth',2)

plot(ca_tp_lft(1),ca_tp_lft(2),'r+','MarkerSize',7,'LineWidth',2)

plot(ca_btm_rgt(1),ca_btm_rgt(2),'r+','MarkerSize',7,'LineWidth',2)

hold off

title('The original image')

subplot 122

surf(data)

shading interp

view(50,40)

colorbar

caxis([-1 1])

title('Illusration of the data')

axis tight

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

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

相关文章

【Android Studio】分类整理res/Layouts中的布局文件(创建子目录)

好吧,这个功能需求也是因为作为一个重度强迫症患者,非常想要把layouts中的文件分类整理下。否则看到一堆Activity和Fragment等文件混在一起特别难受。更重要的原因是,如果你的工程交给别人维护,别人看到一堆布局文件压根不知道从哪…

从零开始学python项目_从零开始可以这样学python

从零开始学python第二章 面向对象第一篇 太原理工大学机器人团队20日打卡day9 ​ 当我们讨论python时不得不讨论到面向对象,而我身为机器人团队的一名成员,在平时设计的时候可以通过把机器人的各个部分设计成不同的对象以此来简化重复的函数引用。并且py…

GDI绘图基本步骤总结(经典)

一、获得绘图的窗口句柄 方法(详细参数及其调用可以看考MSDN): 1、 HWND FindWindow(LPCTSTR lpClassName, LPCTSTR lpWindowName) HWND FindWindowEx(HWND hwndParent, HWND hwndChildAfter,LPCTSTR lpClassName, LPCTSTR lpWindowName) 2、 HWN…

SQL排序查询

1.建一张客户表,包含客户名称,联系人,电话号码,公司地址. 2.建一张订单表,包含订单编号,客户ID,供应商ID,订单金额,客户支付状态,供应商支付状态,订单日期. 3.建一张供应商表&#x…

php smarty框架案例,php封装的smarty类案例

这篇文章主要介绍了php封装的smarty类,针对Smarty的基本操作技巧进行了封装整理,具有一定参考借鉴价值,需要的朋友可以参考下具体如下:/*** Project: Smarty: the PHP compiling template engine* File: Smarty.class.php* SVN: $Id: Smarty.class.php 4848 2014-06…

使用pip安装BeautifulSoup4模块

1.测试是否安装了BeautifulSoup4模块 import bs4 print bs4 执行报错说明没有安装该模块 Traceback (most recent call last):File "D:\work_java\python\src\baike_spider\test.py", line 7, in <module>import bs4 ImportError: No module named bs4 2.使用p…

mysql8创建用户并授权_新特性解读 | 从 wireshark 看 MySQL 8.0 加密连接

作者&#xff1a;秦福朗爱可生 DBA 团队成员&#xff0c;负责项目日常问题处理及公司平台问题排查。热爱 IT&#xff0c;喜欢在互联网里畅游&#xff0c;擅长摄影、厨艺&#xff0c;不会厨艺的 DBA 不是好司机&#xff0c;didi~本文来源&#xff1a;原创投稿*爱可生开源社区出品…

ACM进阶

ACM队不是为了一场比赛而存在的&#xff0c;为的是队员的整体提高。大学期间&#xff0c;ACM队队员必须要学好的课程有&#xff1a;l C/C两种语言l 高等数学l 线性代数l 数据结构l 离散数学l 数据库原理l 操作系统原理l 计算机组成原理l 人工智能l 编译原理l 算法设计与分析除此…

[转 TDD] 如何坚持TDD:使用者出现的问题以及解决方案

原文地址&#xff1a;http://www.infoq.com/cn/articles/levison-TDD-adoption-strategy 转载于:https://www.cnblogs.com/cly84920/archive/2010/03/24/4426820.html

DP备份任务失败原因解析

以JDC为例&#xff0c;DP的备份任务失败troubleshooting流程为&#xff1a;1&#xff09;/etc/opt/omni/server/datalist下都是备份的job code。如果要查看一个job code&#xff0c;可以more 此文件夹下的某个datalist,输入如下&#xff1a;最上面的两个脚本pre和post&#xff…

mysql求表中年龄同张三,mysql子查询与连接查询

表结构以及数据&#xff1a;CREATE TABLEstudent (idint(11) NOT NULLAUTO_INCREMENT,namevarchar(255) CHARACTER SET utf8 COLLATE utf8_general_ci NULL DEFAULT NULL,ageint(11) NULL DEFAULT NULL,PRIMARY KEY(id) USING BTREE) ENGINE InnoDB AUTO_INCREMENT 13 CHARACT…

jenkins war包_Jenkins的安装和部署(jenkins教程)

一、jenkins的下载、安装以及环境的搭建部署。1.什么是jenkins以及它的作用&#xff1a;Jenkins 是一个可扩展的持续集成&#xff08;CI&#xff09;平台。它只是一个平台&#xff0c;真正运作的都是插件。Jenkins的主要作用&#xff1a;持续的监控软件版本发布/测试项目。监控…

.net Web网站转换成Web应用程序 备忘

1.新建Web应用程序&#xff0c;把Web网站代码全部粘贴进来2.添加所有原来引用3.【*】Web应用程序不再自动预编译App_Code文件夹&#xff0c;但是会编译任何文件夹下的类文件(除App_Code),所有给App_Code文件夹换个其它名字。4.【*】把第3部改名后文件夹下的所有文件&#xff0c…

div没有设置高度,背景颜色却无法显示

在设计网页时&#xff0c;设置了div为自动高度或不给height值&#xff0c;当div设置背景时在IE6可以显示背景颜色&#xff0c;但是IE7,8却不显示。 解决办法如下&#xff1a; 1、把height:auto;改成height:100%; overflow:hidden; 原理&#xff1a;IE6可以再不给高度的情况下自…

内存恶鬼drawRect

转载于:http://bihongbo.com/内存恶鬼drawRect 标题有点吓人&#xff0c;但是对于drawRect的评价倒是一点都不过分。在平日的开发中&#xff0c;随意覆盖drawRect方法&#xff0c;稍有不慎就会让你的程序内存暴增。下面我们来看一个例子。 去年的某天午后&#xff0c;北京的雾霾…

c#和python同一主机直接udp_为什么Python 如此之慢

原文:https://hackernoon.com/why-is-python-so-slow-e5074b6fe55b​hackernoon.comPython 正在爆炸般流行起来&#xff0c;它被用于DevOps, 数据处理&#xff0c;web开发和安全领域。但是在速度方面却没有取得过什么胜利。 Java在速度方面和C/C/C#/Python比起来如何&#xff…

matlab编程ode45,matlab 中ode45的源代码如何看懂

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼求常微分方程的数值解 ode45方法的源代码怎么看懂呢&#xff1f;四百多行 如何理解这些代码的核心思想 以方便未来自己使用呢&#xff1f;求大神指点迷津 感激不尽function varargout ode45(ode,tspan,y0,options,varargin)%ODE45…

N2 API

非常详细的N2API文档&#xff0c;以及N2的Wiki和讨论区&#xff0c;特此贴出&#xff0c;以备参考。 参考:http://consultingblogs.emc.com/jamesbroome/archive/2009/04/24/why-i-like-n2-cms.aspx 转载于:https://www.cnblogs.com/chenjunsheep/archive/2010/03/25/1695300.…

C#第零篇

看看荒废的园子和无所作为的自己。真是坑了身边一大帮人。 C、C、C#都是知道语法&#xff0c;自己从来不敢认真地去编一下程啊。真是惭愧。不过现在觉得C#比较简单&#xff0c;就从这里从头开始。伤心啊。 现在没什么时间&#xff0c;大概一天可以有2个小时时间来学习。其它时间…

lldb 调试php,linux系统下如何在vscode中调试C++代码

本篇博客以一个简单的hello world程序&#xff0c;介绍在vscode中调试C代码的配置过程。1. 安装编译器vscode是一个轻量的代码编辑器&#xff0c;并不具备代码编译功能&#xff0c;代码编译需要交给编译器完成。linux下最常用的编译器是gcc&#xff0c;通过如下命令安装&#x…