如何解决三次高斯而不是椭圆高斯
我想绘制一个由高斯调制的椭圆,如以下示例所示:https://www.l3harrisgeospatial.com/docs/GAUSS2DFIT.html(请参见示例,然后创建一个椭圆部分)
from matplotlib.colors import LogNorm
import numpy as np
import matplotlib.pyplot as plt
B_max = 130
B_max2 = 120
theta = 80
x = np.linspace(-B_max*2,B_max*2,100)
XX,YY = np.meshgrid(x,x)
x_p = XX*np.cos(2*np.pi-np.deg2rad(theta))-YY*np.sin(2*np.pi-np.deg2rad(theta))
y_p = XX*np.sin(2*np.pi-np.deg2rad(theta))+YY*np.cos(2*np.pi-np.deg2rad(theta))
gaussian = np.exp(-((x_p/(B_max))**6+(y_p/(B_max2))**6)/2)
plt.figure()
plt.imshow(gaussian,extent =[-B_max*2,-B_max*2,B_max*2])
plt.figure()
plt.imshow(gaussian,B_max*2],norm=LogNorm(vmin=0.01,vmax=1))
但是这不是期望的形状:我更喜欢高斯立方体形状,而不是椭圆调制的高斯形状。看起来长轴和短轴领先于立方体的对角线,而不是我的椭球形状的限制。 meshgrid
我在做错什么吗?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。