Hello Bikram,

Converting a .raw image to a .jpg before calculating the intensity profile is not advisable because JPEG is a lossy format that compresses the image, losing detail and accuracy. It's best to work directly with the raw data to maintain accuracy.

To read a .raw image you can simply use the fopen and fread and to calculate the intensity profile you can use the improfile function (https://www.mathworks.com/help/images/ref/improfile.html).

In the end, a simple code would look like this:

% Define parameters and read the .raw file

fid = fopen('image.raw', 'r');

width = 512; height = 512; channels = 1;

image = fread(fid, width * height * channels, 'uint8');

fclose(fid);

image = reshape(image, [width, height, channels])';

% Calculate and plot the intensity profile along a horizontal line

imshow(image, []);

x = [1 width]; y = [256 256];

intensityProfile = improfile(image, x, y);

figure; plot(intensityProfile);

title('Intensity Profile'); xlabel('Pixel Position'); ylabel('Intensity');

To calculate the intensity profile of a .raw image in MATLAB, first read the image data using fopen and fread. Reshape the data into a matrix and then sum the pixel values along the desired axis (e.g., rows for horizontal profile). Use plot to visualize the intensity profile.

% Define image parameters (you should know these in advance)

width = 512; % Image width (in pixels)

height = 512; % Image height (in pixels)

bitDepth = 8; % Image bit depth (e.g., 8, 16)

% Open the .raw file

fileID = fopen('image.raw', 'rb');

% Read the image data

rawData = fread(fileID, width * height, 'uint8'); % 'uint8' for 8-bit images

fclose(fileID);

% Reshape the raw data into a matrix form

imageData = reshape(rawData, [width, height]);

% Display the image

imshow(imageData, []);

% Display the image to allow user to select a line

imshow(imageData, []);

% Use improfile to extract intensity profile along a line

intensityProfile = improfile(imageData);

% Plot the intensity profile

figure;

plot(intensityProfile);

title('Intensity Profile');

xlabel('Position along the line');

ylabel('Intensity');

Marcos Moura i need your guidance regarding COMSOl boundary conditions