.. _pylab_examples-image_demo2:

pylab_examples example code: image_demo2.py
===========================================



.. plot:: /home/mandrake/rpm/BUILD/matplotlib-1.1.1/doc/mpl_examples/pylab_examples/image_demo2.py

::

    #!/usr/bin/env python
    from pylab import *
    import matplotlib.cbook as cbook
    
    w, h = 512, 512
    
    datafile = cbook.get_sample_data('ct.raw', asfileobj=False)
    print 'loading', datafile
    s = file(datafile, 'rb').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 :ref:`how-to-search-examples`)