pylab_examples example code: image_demo2.py

(Source code, png)

#!/usr/bin/env python

from __future__ import print_function
from pylab import *
import matplotlib.cbook as cbook

w, h = 512, 512

datafile = cbook.get_sample_data('ct.raw.gz', asfileobj=True)
s = datafile.read()
A = fromstring(s, uint16).astype(float)
A *= 1.0/max(A)
A.shape = w, h

extent = (0, 25, 0, 25)
im = imshow(A, cmap=cm.hot, origin='upper', extent=extent)

markers = [(15.9, 14.5), (16.8, 15)]
x, y = zip(*markers)
plot(x, y, 'o')
#axis([0,25,0,25])


#axis('off')
title('CT density')

if 0:
    x = asum(A, 0)
    subplot(212)
    bar(arange(w), x)
    xlim(0, h - 1)
    ylabel('density')
    setp(gca(), 'xticklabels', [])

show()

Keywords: python, matplotlib, pylab, example, codex (see Search examples)