Ejemplo n.º 1
0
def main():

    file_name = '/media/DISK_02/databank/templates/esrf_ID19/tomo.edf'
    dark_file_name = '/media/DISK_02/databank/templates/esrf_ID19/dark.edf'
    white_file_name = '/media/DISK_02/databank/templates/esrf_ID19/flat.edf'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/ESRF.h5'

    sample_name = 'edf test'

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(file_name,
                                              white_file_name=white_file_name,
                                              dark_file_name=dark_file_name,
                                              data_type='edf',
                                              log='INFO')

    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         sample_name=sample_name,
                         data_exchange_type='tomography_raw_projections')
def main():

    file_name = '/media/DISK_02/databank/templates/xradia_dtu/sample_name.txrm'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/DTU.h5'
    sample_name = 'halvmaane_150kV-HE6-20X-60s'

    experimenter_name = "Martin Skovgaard Andersen"
    experimenter_affiliation = "Technical University of Denmark"
    experimenter_email = "*****@*****.**"
    instrument_comment = "Xradia Versa micro CT scanner"

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(file_name,
                                              data_type='xradia',
                                              log='INFO')
    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         experimenter_name=experimenter_name,
                         experimenter_affiliation=experimenter_affiliation,
                         experimenter_email=experimenter_email,
                         instrument_comment=instrument_comment,
                         sample_name=sample_name,
                         data_exchange_type='tomography_raw_projections')
def main():

    file_name = '/media/DISK_02/databank/templates/sls_tomcat/Hornby_SLS/Hornby_b.tif'
    log_file = '/media/DISK_02/databank/templates/sls_tomcat/Hornby_SLS/Hornby.log'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/SLS.h5'

    sample_name = 'Hornby_b'
    
    #Read SLS log file data
    file = open(log_file, 'r')
    for line in file:
        linelist=line.split()
        if len(linelist)>1:
            if (linelist[0]=="Number" and linelist[2]=="darks"):
                number_of_darks = int(linelist[4])
            elif (linelist[0]=="Number" and linelist[2]=="flats"):
                number_of_flats = int(linelist[4])
            elif (linelist[0]=="Number" and linelist[2]=="projections"):
                number_of_projections = int(linelist[4])
            elif (linelist[0]=="Rot" and linelist[2]=="min"):
                rotation_min = float(linelist[6])
            elif (linelist[0]=="Rot" and linelist[2]=="max"):
                rotation_max = float(linelist[6])
            elif (linelist[0]=="Angular" and linelist[1]=="step"):
                angular_step = float(linelist[4])
    file.close()

    dark_start = 1
    dark_end = number_of_darks + 1
    white_start = dark_end
    white_end = white_start + number_of_flats
    projections_start = white_end
    projections_end = projections_start + number_of_projections
    projections_angle_end = 180 + angular_step

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(file_name,
                                                       projections_start = projections_start,
                                                       projections_end = projections_end,
						       projections_angle_end = projections_angle_end,
                                                       projections_digits=4,
                                                       white_start = white_start,
                                                       white_end = white_end,
                                                       dark_start = dark_start,
                                                       dark_end = dark_end,
                                                       log='INFO'
                                                       )
    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data = data,
                          data_white = white,
                          data_dark = dark,
                          theta = theta,
                          hdf5_file_name = hdf5_file_name,
                          sample_name = sample_name,
                          data_exchange_type = 'tomography_raw_projections'
                          )
Ejemplo n.º 4
0
def main():

    file_name = '/media/DISK_02/databank/templates/elettra_SYRMEP/tomo_.tif'
    dark_file_name = '/media/DISK_02/databank/templates/elettra_SYRMEP/dark_.tif'
    white_file_name = '/media/DISK_02/databank/templates/elettra_SYRMEP/flat_.tif'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/Elettra.h5'

    projections_start = 1
    projections_end = 1441
    white_start = 1
    white_end = 11
    white_step = 1
    dark_start = 1
    dark_end = 11
    dark_step = 1

    sample_name = 'Volcanic_rock'

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(
        file_name,
        projections_start=projections_start,
        projections_end=projections_end,
        projections_digits=4,
        white_file_name=white_file_name,
        white_start=white_start,
        white_end=white_end,
        white_step=white_step,
        dark_file_name=dark_file_name,
        dark_start=dark_start,
        dark_end=dark_end,
        dark_step=dark_step,
        data_type='compressed_tiff',
        projections_zeros=True,
        white_zeros=False,
        dark_zeros=False,
        log='INFO')
    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         sample_name=sample_name,
                         data_exchange_type='tomography_raw_projections')
Ejemplo n.º 5
0
def main():

    file_name = '/media/DISK_02/databank/templates/australian_micro-tomography/SAMPLE_T_.tif'
    dark_file_name = '/media/DISK_02/databank/templates/australian_micro-tomography/DF__BEFORE_.tif'
    white_file_name = '/media/DISK_02/databank/templates/australian_micro-tomography/BG__BEFORE_.tif'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/Australian.h5'

    sample_name = 'Teeth'

    projections_start = 0
    projections_end = 18
    white_start = 0
    white_end = 10
    white_step = 1
    dark_start = 0
    dark_end = 10
    dark_step = 1

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(
        file_name,
        projections_start=projections_start,
        projections_end=projections_end,
        white_file_name=white_file_name,
        white_start=white_start,
        white_end=white_end,
        white_step=white_step,
        dark_file_name=dark_file_name,
        dark_start=dark_start,
        dark_end=dark_end,
        dark_step=dark_step,
        projections_digits=4,
        white_digits=2,
        dark_digits=2,
        projections_zeros=True,
        log='INFO')
    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         sample_name=sample_name,
                         data_exchange_type='tomography_raw_projections')
Ejemplo n.º 6
0
def main():

    # oster: pj: from 0 -> 1440; bf from 0 -> 19; df from 0 -> 19
    file_name = '/media/DISK_02/databank/templates/petraIII_P05/sample_name00_0000/scan_0002/ccd/pco01/ccd_.tif'
    dark_file_name = '/media/DISK_02/databank/templates/petraIII_P05/sample_name00_0000/scan_0000/ccd/pco01/ccd_.tif'
    white_file_name = '/media/DISK_02/databank/templates/petraIII_P05/sample_name00_0000/scan_0001/ccd/pco01/ccd_.tif'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/PetraIII.h5'

    projections_start = 0
    projections_end = 1441
    white_start = 0
    white_end = 20
    white_step = 1
    dark_start = 0
    dark_end = 20
    dark_step = 1

    sample_name = 'sample_name00_0000'

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(
        file_name,
        projections_start=projections_start,
        projections_end=projections_end,
        white_file_name=white_file_name,
        white_start=white_start,
        white_end=white_end,
        white_step=white_step,
        dark_file_name=dark_file_name,
        dark_start=dark_start,
        dark_end=dark_end,
        dark_step=dark_step,
        projections_digits=4,
        projections_zeros=True,
        log='INFO')
    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         sample_name=sample_name,
                         data_exchange_type='tomography_raw_projections')
Ejemplo n.º 7
0
def main():

    file_name = '/media/DISK_02/databank/templates/anka_topo-tomo/radios/image_.tif'
    dark_file_name = '/media/DISK_02/databank/templates/anka_topo-tomo/darks/image_.tif'
    white_file_name = '/media/DISK_02/databank/templates/anka_topo-tomo/flats/image_.tif'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/Anka.h5'

    projections_start = 0
    projections_end = 3167
    white_start = 0
    white_end = 100
    dark_start = 0
    dark_end = 100

    sample_name = 'insect'

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(
        file_name,
        projections_start=projections_start,
        projections_end=projections_end,
        white_file_name=white_file_name,
        white_start=white_start,
        white_end=white_end,
        dark_file_name=dark_file_name,
        dark_start=dark_start,
        dark_end=dark_end,
        projections_digits=5,
        log='INFO')

    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         sample_name=sample_name,
                         data_exchange_type='tomography_raw_projections')
Ejemplo n.º 8
0
def main():

    file_name = '/media/DISK_02/databank/templates/aps_13-BM/SPE/sample_name_.SPE'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/APS_13_BM_spe.h5'

    white_start = 1
    white_end = 8
    white_step = 2
    projections_start = 2
    projections_end = 7
    projections_step = 2

    sample_name = 'run2_soln1_2'

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(
        file_name,
        projections_start=projections_start,
        projections_end=projections_end,
        projections_step=projections_step,
        white_start=white_start,
        white_end=white_end,
        white_step=white_step,
        projections_zeros=False,
        white_zeros=False,
        dark_zeros=False,
        projections_digits=1,
        data_type='spe',
        log='INFO')
    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         sample_name=sample_name,
                         data_exchange_type='tomography_raw_projections')
def main():

    file_name = '/media/DISK_02/databank/templates/diamond_JEEP/tiff/sample_name/im_.tif'
    white_file_name = '/media/DISK_02/databank/templates/diamond_JEEP/tiff/sample_name/flat_.tif'

    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/Diamond.h5'

    sample_name = 'waterflow'

    white_start = 0
    white_end = 1
    projections_start = 1000
    projections_end = 2440

    #print "Dark", dark_start, dark_end
    print "White", white_start, white_end
    print "Projection", projections_start, projections_end

    read = dataexchange.Import()
    # Read series of images
    data, white, dark, theta = read.xtomo_raw(
        file_name,
        projections_start=projections_start,
        projections_end=projections_end,
        projections_angle_end=360,
        projections_digits=6,
        white_file_name=white_file_name,
        white_start=white_start,
        white_end=white_end,
        log='INFO')
    write = dataexchange.Export()
    # Create minimal data exchange hdf5 file
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         sample_name=sample_name,
                         hdf5_file_name=hdf5_file_name,
                         data_exchange_type='tomography_raw_projections')
def main():

    file_name = '/local/dataraid/databank/templates/diamond_JEEP/sample_name_subx.nxs'
    hdf5_file_name = '/local/dataraid/databank/dataExchange/tmp/Diamond.h5'

    # Read raw data
    # set slice_start and slice_end to  full size. 1600-1610 is for testing only
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(file_name,
                                              data_type='nxs',
                                              slices_start=1600,
                                              slices_end=1610,
                                              slices_step=1,
                                              log='INFO')

    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         data_exchange_type='tomography_raw_projections')
Ejemplo n.º 11
0
def main():

    file_name = '/media/DISK_02/databank/templates/aps_13-BM/NC/sample_name_.nc'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/APS_13_BM_netCDF.h5'

    white_start = 1
    white_end = 4
    white_step = 2
    projections_start = 2
    projections_end = 3
    projections_step = 1

    sample_name = 'Dorthe_F'

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(
        file_name,
        projections_start=projections_start,
        projections_end=projections_end,
        projections_step=projections_step,
        white_start=white_start,
        white_end=white_end,
        white_step=white_step,
        projections_digits=3,
        data_type='nc',
        log='INFO')
    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         sample_name=sample_name,
                         data_exchange_type='tomography_raw_projections')
def main():

    file_name = '/media/DISK_02/databank/templates/aps_1-ID/data_.tif'
    log_file = '/media/DISK_02/databank/templates/aps_1-ID/TomoStillScan.dat'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/APS_1_ID.h5'

    #Read APS 1-ID log file data
    file = open(log_file, 'r')
    for line in file:
        linelist=line.split()
        if len(linelist)>1:
            if (linelist[0]=="First" and linelist[1]=="image"):
                projections_start = int(linelist[4])
            elif (linelist[0]=="Last" and linelist[1]=="image"):
                projections_end = int(linelist[4])
            elif (linelist[0]=="Dark" and linelist[1]=="field"):
                dark_start = int(linelist[6])
            elif (linelist[0]=="Number" and linelist[2]=="dark"):
                number_of_dark = int(linelist[5])
            elif (linelist[0]=="White" and linelist[1]=="field"):
                white_start = int(linelist[6])
            elif (linelist[0]=="Number" and linelist[2]=="white"):
                number_of_white = int(linelist[5])
    file.close()
    
    dark_end = dark_start + number_of_dark
    white_end = white_start + number_of_white

    # to fix a data collection looging bug ? 
    white_start = white_start + 1
    dark_start = dark_start +1
    projections_start = projections_start + 11
    projections_end = projections_end - 9
    
    
##    # these are correct per Peter discussion
##    projections_start = 943
##    projections_end = 1853
##    white_start = 1844
##    white_end = 1853
##    dark_start = 1854
##    dark_end = 1863
   
    sample_name = 'CAT4B_2'

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(file_name,
                                                       projections_start = projections_start,
                                                       projections_end = projections_end,
                                                       white_start = white_start,
                                                       white_end = white_end,
                                                       dark_start = dark_start,
                                                       dark_end = dark_end,
                                                       projections_digits = 6,
                                                       log='INFO'
                                                       )
    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data = data,
                          data_white = white,
                          data_dark = dark,
                          theta = theta,
                          hdf5_file_name = hdf5_file_name,
                          sample_name = sample_name,
                          data_exchange_type = 'tomography_raw_projections'
                          )
Ejemplo n.º 13
0
def main():

    file_name = '/media/DISK_02/databank/templates/als_beamline_8.3.2/sample_name_0000_.tif'
    dark_file_name = '/media/DISK_02/databank/templates/als_beamline_8.3.2/sample_namedrk_.tif'
    white_file_name = '/media/DISK_02/databank/templates/als_beamline_8.3.2/sample_namebak_.tif'
    log_file = '/media/DISK_02/databank/templates/als_beamline_8.3.2/sample_name.sct'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/ALS.h5'

    verbose = True

    # Read ALS log file data
    file = open(log_file, 'r')
    for line in file:
        if '-scanner' in line:
            Source = re.sub(r'-scanner ', "", line)
            if verbose: print 'Facility', Source

        if '-object' in line:
            Sample = re.sub(r'-object ', "", line)
            if verbose: print 'Sample', Sample

        if '-senergy' in line:
            Energy = re.findall(r'\d+.\d+', line)
            if verbose: print 'Energy', Energy[0]

        if '-scurrent' in line:
            Current = re.findall(r'\d+.\d+', line)
            if verbose: print 'Current', Current[0]

        if '-nangles' in line:
            Angles = re.findall(r'\d+', line)
            if verbose: print 'Angles', Angles[0]

        if '-num_bright_field' in line:
            WhiteEnd = re.findall(r'\d+', line)

        if '-num_dark_fields' in line:
            DarkEnd = re.findall(r'\d+', line)

    file.close()

    dark_start = 0
    dark_end = int(DarkEnd[0])
    dark_step = 1
    white_start = 0
    white_end = int(WhiteEnd[0])
    white_step = 1
    projections_start = 0
    projections_end = int(Angles[0])

    print dark_end, white_end, projections_end

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(
        file_name=file_name,
        projections_start=projections_start,
        projections_end=projections_end,
        white_file_name=white_file_name,
        white_start=white_start,
        white_end=white_end,
        dark_file_name=dark_file_name,
        dark_start=dark_start,
        dark_end=dark_end,
        projections_digits=4,
        log='INFO')

    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         sample_name=Sample,
                         data_exchange_type='tomography_raw_projections')
def main():

    log_file = '/media/DISK_02/databank/templates/aps_2-BM/Sam01_exp.hdf'
    hdf5_file_name = '/media/DISK_02/databank/templates/dataExchange/tmp/APS_2_BM_hdf4.h5'

    #Read APS 2-BM log file data
    f = SD.SD(log_file)

    sds = f.select('base_name')
    data = sds.get()
    base_name = ''.join(data)
    print base_name

    file_name = os.path.split(
        log_file)[0] + "/" + "raw" + "/" + base_name + "_.hdf"

    sds = f.select('start_angle')
    start_angle = sds.get()[0]

    sds = f.select('end_angle')
    end_angle = sds.get()[0]

    sds = f.select('angle_interval')
    angle_interval = sds.get()[0]

    sds = f.select('num_dark_fields')
    num_dark_fields = sds.get()[0]

    f.end()

    white_start = 1
    white_end = 2
    projections_start = 2
    projections_end = projections_start + (int)(
        (end_angle - start_angle) / angle_interval) + 1
    dark_start = projections_end + 1
    dark_end = dark_start + num_dark_fields

    sample_name = base_name

    # Read raw data
    read = dataexchange.Import()
    data, white, dark, theta = read.xtomo_raw(
        file_name,
        projections_start=projections_start,
        projections_end=projections_end,
        white_start=white_start,
        white_end=white_end,
        dark_start=dark_start,
        dark_end=dark_end,
        projections_digits=5,
        data_type='hdf4',
        log='INFO')
    # Save data as dataExchange
    write = dataexchange.Export()
    write.xtomo_exchange(data=data,
                         data_white=white,
                         data_dark=dark,
                         theta=theta,
                         hdf5_file_name=hdf5_file_name,
                         sample_name=sample_name,
                         data_exchange_type='tomography_raw_projections')