昨日の続き

axis() コマンドを使うだけだった。実に単純。

## ipython -pylab で、すでにもろもろ import ずみ
delta = 0.5
extent = (-3.25,4.25,-4.25,3.25)

# データの生成
x = arange(-3.0, 4.001, delta)
y = arange(-4.0, 3.001, delta)
X, Y = meshgrid(x, y)
Z1 = bivariate_normal(X, Y, 1.0, 1.0, 0.0, 0.0)
Z2 = bivariate_normal(X, Y, 1.5, 0.5, 1, 1)
Z = (Z1 - Z2) * 10

levels = arange(-2.0, 1.601, 0.4) 
fig = figure()
im = imshow(Z, interpolation='nearest', extent=extent)
show()
contour(Z, levels, hold='on', colors = 'k', origin='image', extent=extent)
axis((-3., 3., -4., 2.))          # (x 軸最小, x 軸最大, y 軸最小, y 軸最大)
ylim = get(gca(), 'ylim')
setp(gca(), ylim=ylim[::-1])
title("Image, clipped, reversed y-axis")
colorbar(im)