def test_merge_Zmap():
reader = SeqReader()
reader.read_seq(resource('seq/au30_a1_.SEQ'))
Z1, _, N1 = reader.get_Zmap('orsnr___')
reader.read_seq(resource('seq/au30_m1_.SEQ'))
Z2, _, N2 = reader.get_Zmap('orsnr___')
Z, N = SeqReader.merge_Zmap(Z1, Z2, N1, N2)
def __init__(self, seqfiles):
Z_merged, N_merged = None, None
readers = map(SeqReader, seqfiles)
for reader in readers:
Z, _, N = reader.get_Zmap('orsnr___', thres=10)
if Z_merged is None:
Z_merged, N_merged = Z, N
else:
Z_merged, N_merged = SeqReader.merge_Zmap(Z, Z_merged, N, N_merged)
self.__Z = Z_merged
self.__NX = len(self.__Z[0])
self.__NY = len(self.__Z)
def test_read_seq():
reader = SeqReader(resource('seq/Quartz_500Mpa_.SEQ'))
reader.get_Om()
Z, _, N = reader.get_Zmap('orsnr___')
