From c90c04249fb195b33ee9bf3ce2a815419a90df86 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Mon, 26 May 2025 16:56:02 +0200 Subject: [PATCH] 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. --- Data-Analyzer/plotImages.m | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/Data-Analyzer/plotImages.m b/Data-Analyzer/plotImages.m index 9e97009..c5b5551 100644 --- a/Data-Analyzer/plotImages.m +++ b/Data-Analyzer/plotImages.m @@ -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; angle = 0; -center = [1630, 1700]; -span = [500, 500]; +center = [1375, 2020]; +span = [150, 150]; fraction = [0.1, 0.1]; pixel_size = 5.86e-6; @@ -53,9 +53,9 @@ for k = 1 : length(files) fprintf(1, 'Now reading %s\n', fullFileName); - atm_img = im2double(imrotate(h5read(fullFileName, append(groupList(cam), "/atoms")), angle)); - bkg_img = im2double(imrotate(h5read(fullFileName, append(groupList(cam), "/background")), angle)); - dark_img = im2double(imrotate(h5read(fullFileName, append(groupList(cam), "/dark")), angle)); + atm_img = double(imrotate(h5read(fullFileName, append(groupList(cam), "/atoms")), angle)); % im2double rescales values to between [0, 1], use double instead + bkg_img = double(imrotate(h5read(fullFileName, append(groupList(cam), "/background")), angle)); + dark_img = double(imrotate(h5read(fullFileName, append(groupList(cam), "/dark")), angle)); refimages(:,:,k) = subtractBackgroundOffset(cropODImage(bkg_img, center, span), fraction)'; absimages(:,:,k) = subtractBackgroundOffset(cropODImage(calculateODImage(atm_img, bkg_img, dark_img), center, span), fraction)';