1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
|
import pycbf
object = pycbf.cbf_handle_struct() # FIXME
object.read_file("../img2cif_packed.cif",pycbf.MSG_DIGEST)
object.rewind_datablock()
print "Found",object.count_datablocks(),"blocks"
object.select_datablock(0)
print "Zeroth is named",object.datablock_name()
object.rewind_category()
categories = object.count_categories()
for i in range(categories):
print "Category:",i,
object.select_category(i)
category_name = object.category_name()
print "Name:",category_name,
rows=object.count_rows()
print "Rows:",rows,
cols = object.count_columns()
print "Cols:",cols
loop=1
object.rewind_column()
while loop is not 0:
column_name = object.column_name()
print "column name \"",column_name,"\"",
try:
object.next_column()
except:
break
print
for j in range(rows):
object.select_row(j)
object.rewind_column()
print "row:",j
for k in range(cols):
name=object.column_name()
print "col:",name,
object.select_column(k)
typeofvalue=object.get_typeofvalue()
print "type:",typeofvalue
if typeofvalue.find("bnry") > -1:
print "Found the binary!!",
s=object.get_integerarray_as_string()
print type(s)
print dir(s)
print len(s)
try:
import numpy
d = numpy.frombuffer(s,numpy.uint32)
# Hard wired Unsigned Int32
print d.shape
print d[0:10],d[d.shape[0]/2],d[-1]
print d[d.shape[0]/3:d.shape[0]/3+20]
d=numpy.reshape(d,(2300,2300))
# from matplotlib import pylab
# pylab.imshow(d,vmin=0,vmax=1000)
# pylab.show()
except ImportError:
print "You need to get numpy and matplotlib to see the data"
else:
value=object.get_value()
print "Val:",value,i
print
del(object)
#
print dir()
#object.free_handle(handle)
|
