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()
# 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()