예제 #1
0
filename_SOM1 = "stuff1.dat"

SOM1 = SOM()
SOM1.attr_list["filename"] = filename_SOM1
SOM1.attr_list["epoch"] = time()
SOM1.attr_list["timestamp"] = DST.make_ISO8601(SOM1.attr_list["epoch"])
SOM1.attr_list["username"] = "******"
SOM1.setAllAxisLabels(["Q", "E"])
SOM1.setAllAxisUnits(["A-1", "meV"])
SOM1.setYLabel("Intensity")
SOM1.setYUnits("Counts/(meV A-1))")

SO1 = SO(2)
SO1.id = 0
SO1.axis[0].val.extend(range(5))
SO1.axis[1].val.extend(range(10))

y_len = (len(SO1.axis[0].val) - 1) * (len(SO1.axis[1].val) - 1)
y = range(y_len)
SO1.y.extend(y)
SO1.var_y.extend(y)

SOM1.append(SO1)

file = open(filename_SOM1, "w")

d2d = DST.Dave2dDST(file)
d2d.writeSOM(SOM1)
d2d.release_resource()
예제 #2
0
# warranty, express or implied, or assumes any legal liability or
# responsibility for the accuracy, completeness, or usefulness of any
# information, apparatus, product, or process disclosed, or represents that
# its use would not infringe privately owned rights.
#

# $Id$

import sys
import DST

filename = sys.argv[1]

ifile = open(filename, "r")

d2d = DST.Dave2dDST(ifile)

som = d2d.getSOM()

d2d.release_resource()
x_index = 1
y_index = 31
N_y = len(som[0].axis[1].val)
print "N_y:", N_y
channel = y_index + (x_index * N_y)
print "channel:", channel
print som.attr_list
print "Y Label:", som.getYLabel()
print "Y Units:", som.getYUnits()
print "X Labels:", som.getAllAxisLabels()
print "X Units:", som.getAllAxisUnits()