Esempio n. 1
0
#==============================================================================
# %% LOAD GRAIN AND EXPERIMENT DATA
#==============================================================================

experiment, nf_to_ff_id_map  = nfutil.gen_trial_exp_data(grain_out_file,det_file,mat_file, x_ray_energy, mat_name, max_tth, comp_thresh, chi2_thresh, misorientation_bnd, \
                       misorientation_spacing,ome_range_deg, num_imgs, beam_stop_width)

#==============================================================================
# %% NEAR FIELD - MAKE MEDIAN DARK
#==============================================================================
print '>>>>>>>>>>>>>>>>>>loading images>>>>>>>>>>>>>>>>>>'
dark = nfutil.gen_nf_dark(data_folder,
                          img_nums,
                          num_for_dark,
                          experiment.nrows,
                          experiment.ncols,
                          dark_type='median',
                          num_digits=6)

#==============================================================================
# %% NEAR FIELD - LOAD IMAGE DATA AND PROCESS
#==============================================================================

image_stack = gen_nf_cleaned_image_stack(
    data_folder,
    img_nums,
    dark,
    ome_dilation_iter,
    threshold,
    experiment.nrows,
#==============================================================================
# %%  TOMO PROCESSING - VIEW TOMO_MASK FOR SAMPLE BOUNDS
#==============================================================================
plt.close('all')
plt.imshow(tomo_mask,interpolation='none')

#==============================================================================
# %%  TOMO PROCESSING - CONSTRUCT DATA GRID
#==============================================================================

test_crds, n_crds, Xs, Ys, Zs = nfutil.gen_nf_test_grid_tomo(tomo_mask.shape[1], tomo_mask.shape[0], v_bnds, voxel_spacing)

#==============================================================================
# %% NEAR FIELD - MAKE MEDIAN DARK
#==============================================================================
dark=nfutil.gen_nf_dark(data_folder,img_nums,num_for_dark,experiment.nrows,experiment.ncols)


#==============================================================================
# %% NEAR FIELD - LOAD IMAGE DATA AND PROCESS
#==============================================================================

image_stack=nfutil.gen_nf_image_stack(data_folder,img_nums,dark,num_erosions,num_dilations,ome_dilation_iter,threshold,experiment.nrows,experiment.ncols)


#==============================================================================
# %% VIEW IMAGES FOR DEBUGGING TO LOOK AT IMAGE PROCESSING PARAMETERS
#==============================================================================
plt.close('all')
img_to_view=0
plt.imshow(image_stack[img_to_view,:,:],interpolation='none')