karthik/Calculations
1
0
Fork 0
Code Issues Pull Requests repo.project_board Releases Wiki Activity

Changed im2double to double to get original intensity counts without it being rescaled to within [0,1] - necessary for accurate OD determination at high intensity.

Browse Source
This commit is contained in:
Karthik 2025-05-26 16:56:02 +02:00
parent dc5ccb67ba
commit c90c04249f
1 changed files with 8 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
Data-Analyzer/plotImages.m
View File

@ -4,17 +4,17 @@ groupList = ["/images/MOT_3D_Camera/in_situ_absorption", "/images/ODT_1_Axis_Ca
folderPath = "C:/Users/Karthik/Documents/GitRepositories/Calculations/Data-Analyzer/"; folderPath = "D:/Data - Experiment/2025/05/22/";
run = '0060'; run = '0078';
folderPath = strcat(folderPath, run); folderPath = strcat(folderPath, run);
cam = 5; cam = 5;
angle = 0; angle = 0;
center = [1630, 1700]; center = [1375, 2020];
span = [500, 500]; span = [150, 150];
fraction = [0.1, 0.1]; fraction = [0.1, 0.1];
pixel_size = 5.86e-6; pixel_size = 5.86e-6;
@ -53,9 +53,9 @@ for k = 1 : length(files)
fprintf(1, 'Now reading %s\n', fullFileName); fprintf(1, 'Now reading %s\n', fullFileName);
atm_img = im2double(imrotate(h5read(fullFileName, append(groupList(cam), "/atoms")), angle)); atm_img = double(imrotate(h5read(fullFileName, append(groupList(cam), "/atoms")), angle)); % im2double rescales values to between [0, 1], use double instead
bkg_img = im2double(imrotate(h5read(fullFileName, append(groupList(cam), "/background")), angle)); bkg_img = double(imrotate(h5read(fullFileName, append(groupList(cam), "/background")), angle));
dark_img = im2double(imrotate(h5read(fullFileName, append(groupList(cam), "/dark")), angle)); dark_img = double(imrotate(h5read(fullFileName, append(groupList(cam), "/dark")), angle));
refimages(:,:,k) = subtractBackgroundOffset(cropODImage(bkg_img, center, span), fraction)'; refimages(:,:,k) = subtractBackgroundOffset(cropODImage(bkg_img, center, span), fraction)';
absimages(:,:,k) = subtractBackgroundOffset(cropODImage(calculateODImage(atm_img, bkg_img, dark_img), center, span), fraction)'; absimages(:,:,k) = subtractBackgroundOffset(cropODImage(calculateODImage(atm_img, bkg_img, dark_img), center, span), fraction)';