filters = tiledb.FilterList([tiledb.LZ4Filter()])
dom = tiledb.Domain(tiledb.Dim(name='scan', domain=(1, 506), tile=20, dtype=float),
                    tiledb.Dim(name='mz', domain=(0, 2000), tile=10, dtype=float))
schema = tiledb.ArraySchema(domain=dom, sparse=True, capacity=1024,
                            attrs=[tiledb.Attr(name='intensity', dtype=float, filters=filters)])
spec = tiledb.SparseArray.create('spec', schema)

with tiledb.SparseArray('spec', mode='w') as spec:
    #npoints = 27_464_448
    #scan_arr = np.zeros(npoints)
    #mz_arr = np.zeros(npoints)
    #intens_arr = np.zeros(npoints)

    tiledb.stats_enable()

    n_scans = 1_000_000
    scan_arr = np.repeat(np.arange(1,507), 1000)
    mz_arr = np.tile(np.arange(0,1000,step=1), 506)
    intens_arr = np.random.rand(1000*506)

    #i = 0
    #scan = 0
    #for line in open('nefazodone.raw'):
    #    line = line.strip()
    #    if line.startswith('RetTime='):
    #        scan += 1
    #    elif line.startswith('Mz='):
    #        mz, intens = line.split('=')[1].split(' ')
    #        scan_arr[i] = scan
    #        mz_arr[i] = float(mz)
    #        intens_arr[i] = float(intens)
    #        i += 1

    spec[scan_arr, mz_arr] = {'intensity': intens_arr}
    #assert i == npoints
    tiledb.stats_dump()
    tiledb.stats_disable()

tiledb.stats_enable()
data = spec[468.5:469.5, :]
tiledb.stats_dump()
tiledb.stats_disable()
print(data['intensity'])

</details>