Satellite image processing in Matlab

I am very new to Matlab and ı want to process my satellite images in it. I have 96 different monthly satellite images (2003-2010) and I want to make some calculations for them. I want to subtract respective months from each other to calculate changes. Let say m(j) is the jth month and m(j+1) is the j+1 month. Δm= m(j+1)-m(j). I want to apply this formula for all respective months and obtain 96 different results separately how can I do this ?

[Cyrus]

Do you have one image per month or separated bands? What type of satellite images, are you using? (Landsat-8? Landsat 7?. ,,,)

[Gokhan Kayan]

I have one GRACE satellite image for every month (96 image) and it only includes one band showing water thickness. To calculate water thickness change ı should subtract respective bands from each other.

Answer 1

[Image Analyst]

Just have a loop, like from the FAQ: http://matlab.wikia.com/wiki/FAQ#How_can_I_process_a_sequence_of_files.3F

% Specify the folder where the files live.
myFolder = 'C:\Users\yourUserName\Documents\My Pictures';
% Check to make sure that folder actually exists.  Warn user if it doesn't.
if ~isdir(myFolder)
  errorMessage = sprintf('Error: The following folder does not exist:\n%s', myFolder);
  uiwait(warndlg(errorMessage));
  return;
end
% Get a list of all files in the folder with the desired file name pattern.
filePattern = fullfile(myFolder, '*.jpg'); % Change to whatever pattern you need.
theFiles = dir(filePattern);
for k = 1 : length(theFiles)
  baseFileName = theFiles(k).name;
  fullFileName = fullfile(myFolder, baseFileName);
  fprintf(1, 'Now reading %s\n', fullFileName);
  % Now do whatever you want with this file name,
  % such as reading it in as an image array with imread()
  imageArray = imread(fullFileName);
  if k == 1
        priorImage = imageArray;
        continue; % Skip to bottom of loop.
  end
  imshow(imageArray);  % Display image.
  drawnow; % Force display to update immediately.
  differenceImage{k} = double(imageArray) - double(priorImage);
  % Now do whatever you want with differenceImage{k}.....
end

[Gokhan Kayan]

thank you frien ı try to run it.

Answer 2

[BANDARU UMAMADHURI]

Hii,Can someone give the insights of satellite image processing using matlab.Where we take a remote sensing data and need to find lime stone in a perticular area