mrcfile_create_populate_read_basic_example.py

1import mrcfile
2import numpy as np
3
4# Create a 3D numpy array of data
5data = np.arange(12, dtype=np.int16).reshape(3, 2, 2)
6
7# Create a new MRC file and save the data
8with mrcfile.new('test.mrc', overwrite=True) as mrc:
9    mrc.set_data(data)
10
11# Open an existing MRC file for reading
12with mrcfile.open('test.mrc') as mrc:
13    print(mrc.data)