开发者社区

文档建议反馈控制台

文章/答案/技术大牛

发布

社区首页 >问答首页 >利用scipy和matplotlib实现散乱三维数据等值线的插值和绘制

问利用scipy和matplotlib实现散乱三维数据等值线的插值和绘制
EN

Stack Overflow用户

提问于 2017-11-29 17:26:17

回答 1查看 1.4K关注 0票数 0

我有一个散乱网格上的三维有限元计算结果。我想用python在一个水平面上插值这些数据，并使用matplotlib在那个平面上绘制结果的等高线。输出文件results.dat包含6列，即x、y、z坐标和三个位移分量。我对第三部分感兴趣。我想在一个100 x 100的矩形网格上进行插值，网格跨度从(-50,10)到(50,110)，位置为z=-5。我想过以下几种可能性：

import numpy as np
from scipy.interpolate import grid data

res = np.loadtxt('results.dat')
points = res[:,0:3]
values = res[:,-1]
xi,yi,zi = np.mgrid[-50.:50.:101j,10.:100.:101j,-5:-5:1j]

gridres = griddata(points, values, (xi, yi, zi))

为了检查结果是否正确，我应该能够绘制等高线，如下所示

plt.contour(xi,yi,griddata)

但是我不能以plt.contour要求的格式处理xi、yi和griddata

multidimensional-array

EN

回答 1

Stack Overflow用户

发布于 2017-11-29 17:48:08

经过反复试验，我终于发现，为了得到我想要的结果，我必须这样做：

import matplotlib.pyplot as plt

x_contour = np.mgrid[-50.:50.:101j]
y_contour = np.mgrid[10.:100.:101j]

plt.contourf(x_contour, y_contour, gridres[:,:,0].T)
plt.show()

如果有人能帮助我理解为什么会这样，我将不胜感激！有没有其他建议的插值方法？如果我可以使用任何高阶插值方法而不是线性插值方法，那就太好了，这似乎是这种3D情况下唯一的选择。

票数 0

EN

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/47558392

复制

相关文章

Postgresql源码（86）varchar的创建与插入分析

一般PG数据类型不会带括号，varchar和numeric是比较特殊的基本类型，后面可以带括号指定长度或精度。

mingjie

2022/10/28

5470

PostgreSQL 列转行问题

java https 网络安全

SELECT relative_label_content FROM frk_s.label_cor_gene

全栈程序员站长

2022/07/04

1.1K0

postgresql 触发器简介（转）

– 把before for each row的触发器删掉, 再测试插入 : postgres=# drop trigger tg02 on t_ret; DROP TRIGGER postgres=# drop trigger tg2 on t_ret; DROP TRIGGER postgres=# insert into t_ret values(1,’digoal’,now()); NOTICE: 00000: tg01 LOCATION: exec_stmt_raise, pl_exec.c:2840 NOTICE: 00000: tg1 LOCATION: exec_stmt_raise, pl_exec.c:2840 NOTICE: 00000: tg03, after for each row 的触发器函数返回空, 不影响后续的触发器是否被调用. 因为只要表上面发生了真正的行操作, after for each row就会被触发, 除非when条件不满足. (这个后面会讲到) LOCATION: exec_stmt_raise, pl_exec.c:2840 NOTICE: 00000: tg3 LOCATION: exec_stmt_raise, pl_exec.c:2840 NOTICE: 00000: tg04 LOCATION: exec_stmt_raise, pl_exec.c:2840 NOTICE: 00000: tg4 LOCATION: exec_stmt_raise, pl_exec.c:2840 INSERT 0 1 – 有数据插入. 这也说明了before for each statement的返回值为空并不会影响数据库对行的操作. 只有before for each row的返回值会影响数据库对行的操作. postgres=# select * from t_ret ; id | info | crt_time —-+——–+—————————- 1 | digoal | 2013-03-10 16:50:39.551481 (1 row)

qubianzhong

2019/07/01

4K0

Postgresql 强制修改列的类型

stackoverflow table using

当列的类型为字符类型，当我们想修改为数值类型时，是无法成功的，这个时候我们可以通过以下方法进行修改。

飞奔去旅行

2019/06/13

2.9K0

Excel快速插入空白列

如何快速插入空白列呢？在开始菜单栏选择插入，插入工作表列（即可在左边创建一个新的空白列）

magize

2023/07/28

2690

hashmap 的重新散列和装载因子

hashmap 编程算法

HashMap 的装载因子是 0.75，用人话说就是当 HashMap 的容量达到定义容量的 75% 的时候，HashMap 会进行扩容，当 HashMap 进行扩容的时候就会重新散列（rehashing）。

HoneyMoose

2022/11/30

5940

hashmap 的重新散列和装载因子

【说站】python哈希散列的映射

编程算法 python

以上就是python哈希散列的映射，希望对大家有所帮助。更多Python学习指路：python基础教程

很酷的站长

2022/11/23

7470

【说站】python哈希散列的映射

hashmap 的重新散列和装载因子

hashmap 编程算法

HashMap 的装载因子是 0.75，用人话说就是当 HashMap 的容量达到定义容量的 75% 的时候，HashMap 会进行扩容，当 HashMap 进行扩容的时候就会重新散列（rehashing）。

HoneyMoose

2022/12/03

4650

hashmap 的重新散列和装载因子

Postgresql插入或更新操作upsert

云推荐引擎 postgresql

幂等性的一个要求是多次操作的结果一致。对于update操作，多次直接的结果都是最后update的值，是满足需求的。但对于insert，如果已经插入，第二次会报错，duplicate error，主键重复或者unique key duplicate。所以需要做一下处理。

Ryan-Miao

2019/01/24

2K0

重新认识 Java 中的内存映射（mmap）

缓存编程算法自动驾驶 java

mmap 是一种内存映射文件的方法，即将一个文件映射到进程的地址空间，实现文件磁盘地址和一段进程虚拟地址的映射。实现这样的映射关系后，进程就可以采用指针的方式读写操作这一段内存，而系统会自动回写脏页到对应的文件磁盘上，即完成了对文件的操作而不必再调用 read,write 等系统调用函数。相反，内核空间对这段区域的修改也直接反映用户空间，从而可以实现不同进程间的文件共享。

kirito-moe

2021/11/17

4.7K0

重新认识 Java 中的内存映射（mmap）

Postgresql在哪里使用列统计信息？

postgresql 遍历索引统计优化

对pg_statistic表的查询都是走syscache的，要找到所有使用列统计信息地方，遍历系统表索引即可

mingjie

2023/10/13

1580

Postgresql在哪里使用列统计信息？

我的 Serverless 实战 — 云函数与触发器的创建与使用 ( 开通腾讯云 “ 云开发 “ 服务 | 创建云函数 | 创建触发器 | 测试触发器 )

云函数云开发腾讯云服务函数

阿里云 , 腾讯云 , 都提供了相关 Serverless 服务 , 这里以腾讯云为例进行展示 ;

韩曙亮

2023/03/29

1.7K0

我的 Serverless 实战 — 云函数与触发器的创建与使用 ( 开通腾讯云 “ 云开发 “ 服务 | 创建云函数 | 创建触发器 | 测试触发器 )

PostgreSQL表用户列最大个数

postgresql 磁盘函数系统优化

有些业务可能有这么个需求：需要增加用户列，即通过ALTER TABLE ... ADD...来添加用户列。那么PG/GP中是否会有列个数的限制呢？

yzsDBA

2023/09/08

3250

PostgreSQL表用户列最大个数

触发器与视图的创建与使用

今日小知识点：DML触发器按出发时刻分两类：after触发器（在表中数据修改之后出发，是默认类型）、instead of 触发器（在表中数据修改以前触发）

ellipse

2019/08/16

1.4K0

触发器创建删除等操作

编程算法 sql server 数据库 linux sql

触发器是一种特殊的存储过程，类似于事件函数，SQL Server™ 允许为 INSERT、UPDATE、DELETE 创建触发器，即当在表中插入、更新、删除记录时，触发一个或一系列 T-SQL语句。

全栈程序员站长

2022/07/15

1.7K0

PostgreSQL创建表分析

脚本准备创建表的脚本 CREATE DATABASE sampledb OWNER perrynzhou; GRANT ALL PRIVILEGES ON DATABASE sampledb TO perrynzhou; 数据登录脚本 psql -h 127.0.0.1 -d sampledb sampledb=# CREATE TABLE stu_xx_01(NAME TEXT NOT NULL,AGE INT NOT NULL); 表创建分析表创建过程概述服务进程接受SQ

用户4700054

2022/08/17

1.7K0

PostgreSQL 创建表格的方法

云推荐引擎 postgresql 数据库 sql

CREATE TABLE 在当前数据库创建一个新的空白表，该表将由发出此命令的用户所拥有。

用户4988085

2021/09/17

2K0

PostgreSQL 教程

postgresql 教程

本 PostgreSQL 教程可帮助您快速了解 PostgreSQL。您将通过许多实际示例快速掌握 PostgreSQL，并将这些知识应用于使用 PostgreSQL 开发应用程序。

postgres

2023/10/22

6390

数据库 PostgreSQL 常用命令

https 网络安全数据库 postgresql nosql

在我们开始学习 PostgreSQL 数据库前，让我们先了解下 ORDBMS 的一些术语：

为为为什么

2023/02/21

2.3K0

数据库 PostgreSQL 常用命令

数据库 PostgreSQL 常用命令

https 网络安全数据库 postgresql nosql

在我们开始学习 PostgreSQL 数据库前，让我们先了解下 ORDBMS 的一些术语：

为为为什么

2023/01/16

2.2K0

相似问题

如何在PostgreSQL中编写重新插入触发器？

22

为postgresql中的每个新插入创建触发器

13

Postgresql插入触发器连接

25

Postgresql插入规则/触发器

22

用于插入的postgresql触发器

13

添加站长进交流群

领取专属 10元无门槛券

AI混元助手 在线答疑

关注 腾讯云开发者公众号

洞察 腾讯核心技术

剖析业界实践案例