Seperated inner phase example

Introduction

This example is taken from Listing S6 in the 2013 JBNMR nmrglue paper. In this example a pseudo-3D NMRPipe data set is seperated into 2D data sets using the script seperate.py where there is an innermost quadrature phase loop.

Instructions

Execute python seperate.py to seperate the pseudo-3D data set into 2D data sets. Ten directories with names nredor_XX.fid will be created.

The data used in this example is available for download.

Listing S6

[seperate.py]

import nmrglue as ng
import numpy as np

# read the NMR data, forcing the data to be two dimensional
dic, data = ng.varian.read('arrayed_data.fid', as_2d=True)

# set the new size of the separated data
array_size = len(dic['procpar']['nredor']['values'])
out_shape = int(data.shape[0] / array_size), data.shape[1]
dic['nblocks'] = out_shape[0]

# loop over the redor multiples, separating and saving each 2D
for i, nredor in enumerate(dic['procpar']['nredor']['values']):
    dir_name = 'nredor_' + nredor + '.fid'
    print "Creating directory:", dir_name
    sdata = np.empty(out_shape, dtype=data.dtype)
    sdata[::2] = data[2 * i::2 * array_size]
    sdata[1::2] = data[2 * i + 1::2 * array_size]
    ng.varian.write(dir_name, dic, sdata, overwrite=True)

Table Of Contents

Previous topic

Seperated interleaved example

Next topic

Processing S3E filtered data example

This Page