2012-07-14 128 views
3

我试图用numpy.meshgrid和pylab.imshow()创建一个轮廓图;这个工作非常出色,在首先使用一些技巧我计算器上找到(感谢家伙!d)使用imshow修正坐标轴

f = scipy.linspace(1e5,1e6,100) 
A = scipy.linspace(1e3,1e5,100) 
ff,AA = numpy.meshgrid(f,A) 
SLP = calc_SLP2D(ff,AA) 
maxAmps = maxA(f) 
print maxAmps 
brezovich = calc_SLP2D(f,maxAmps) 
print brezovich 
pylab.imshow(SLP,origin='lower') 
pylab.plot(f,maxA(f)) 
pylab.colorbar() 
pylab.xlabel('Frequency [kHz]',{'fontsize':20}) 
pylab.ylabel('Field Amplitude [A/m]',{'fontsize':20}) 
pylab.title('Brezovich Criterion',{'fontsize':20}) 
pylab.grid() 
pylab.show() 

contour image with incorrect axes http://web.mit.edu/scottnla/Public/SLP_contour.pdf

但是,你会发现,轴由两个输入的大小编号矩阵而不是实际值。横坐标应该在10万到1,000,000纵坐标去从1000到5000我的计算器阅读该解决方案是使用“范围”选项,所以:

pylab.imshow(SLP,origin='lower',extent=(ff.min(),ff.max(),AA.min(),AA.max())) 

这并修复轴,但是尺度图像是一个非常奇怪的方式:

contour plot with correct axes but incorrect graphic http://web.mit.edu/scottnla/Public/SLP_contour2.pdf

而且我不知道是什么原因造成这一点。

有关如何重新调整坐标轴而不使图像看起来很奇怪的想法吗?

谢谢!!

弥敦道lachenmyer

回答

11

尝试将参数aspect='auto'添加到imshow。像这样:

pylab.imshow(SLP,aspect='auto',origin='lower',extent=(ff.min(),ff.max(),AA.min(),AA.max())) 
+0

非常感谢!这固定了一切! – 2012-07-14 08:20:24

+0

很高兴我可以帮助:) – 2012-07-14 16:48:43