#============================================================================== # %% 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')