python - 无法使用 Rectbivariatespline 对网格网格上的二维数组进行插值
问题描述
我正在尝试S
在精细的网格上绘图。S
是一个二维数组,其值位于SR
和处Temp
。我想制作 和 的精细网格SR
并Temp
绘制 的等高线图S
。我正在尝试使用RectBivariateSpline
,但在查找时出现错误Z2
。
我得到的错误是:
in line Z2 = zspline(TT,SR)
raise ValueError("Error code returned by bispev: %s" % ier)
ValueError: Error code returned by bispev: 10
我的代码:
import numpy as np
import matplotlib.pyplot as plt
from scipy.interpolate import RectBivariateSpline
T = np.array([850,800,750])
strain_rate = np.array([0.0003,0.001,0.01,0.1])
flow_stress1 = np.array([[28.16, 49.46, 95.96],
[46.45, 80.51, 126.62],
[151.46, 107.94, 235.14],
[165.63, 228.77, 319.15]])
Temp = 1000/(273+T)
TT = np. linspace(max(Temp), min(Temp),50)
S = np.log10(flow_stress1)
SR= np.log10(strain_rate)
srate= np.linspace(min(SR), max(SR),100)
zspline = RectBivariateSpline(SR,Temp,S , kx=2, ky=2)
Z2 = zspline(TT,SR)
fig = plt.figure()
ax = fig.add_subplot(111)
x,y = np.meshgrid(Temp,SR)
ax.contourf(x,y,Z2)
plt.show()
我试过寻找但没有找到解决这个问题的令人满意的方法。
解决方案
推荐阅读
- c - C奇怪行为中的简单字符串
- .net-core - 无法将 switch 表达式与 dotnet-script 一起使用
- java - 如何用我的流中的空列表替换任何 null 实例?
- vb.net - 从 SQL 读取到 Button
- ssl - IBM-MQ 遥测通道中的刷新 SSL 安全性未按预期工作
- pandas - How to mark or select rows in one dataframe, where value lies between any of the ranges in another dataframe featuring additional identifier
- testing - 转换 Quarkus 测试的响应正文
- git - git rebase 在已经推送到远程的功能分支上
- json - 如何将分号分隔的参数转换为json?
- xsd - 使用 XSD (1.1) 将属性设置为某个值时使元素成为必需