Skip to content

Commit

Permalink
Merge pull request #4826 from ericmjl/mep12_image_demo2.py
Browse files Browse the repository at this point in the history
mep12 on image_demo2.py
  • Loading branch information
WeatherGod committed Jul 30, 2015
2 parents afc0280 + 1df0074 commit 12efb7f
Showing 1 changed file with 9 additions and 20 deletions.
29 changes: 9 additions & 20 deletions examples/pylab_examples/image_demo2.py
@@ -1,35 +1,24 @@
#!/usr/bin/env python

from __future__ import print_function
from pylab import *
# from pylab import *
import matplotlib.pyplot as plt
import numpy as np
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 = np.fromstring(s, np.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)
im = plt.imshow(A, cmap=plt.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')
plt.plot(x, y, 'o')

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

show()
plt.show()

0 comments on commit 12efb7f

Please sign in to comment.