Blame view

stim/matlab/rtsEnviSaveImage.m 1.02 KB
1ff9af85   David Mayerich   added MATLAB/OCTA...
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
  function rtsEnviSaveImage(S, fileName, headerName, W)
  
  %  S = 3D hyperspectral data
  %  W = list of wavelength values
  
  %save the image
  rtsSaveRAW(S, fileName, 'float32', 'l');
  
  %create the header file
  headerFile = [fileName '.hdr'];
  if nargin == 3
      headerFile = headerName;
  end
  
  fid = fopen(headerFile, 'w');
  
  %write the header and description
  fprintf(fid, 'ENVI\n');
  fprintf(fid, 'description = {\n');
  fprintf(fid, '  File Saved using Matlab.}\n');
  fprintf(fid, 'samples = %d\n', size(S, 1));
  fprintf(fid, 'lines   = %d\n', size(S, 2));
  fprintf(fid, 'bands   = %d\n', size(S, 3));
  fprintf(fid, 'header offset = 0\n');
  fprintf(fid, 'file type = ENVI Standard\n');
  fprintf(fid, 'data type = 4\n');
  fprintf(fid, 'interleave = bsq\n');
  fprintf(fid, 'sensor type = Unknown\n');
  fprintf(fid, 'byte order = 0\n');
  
  if nargin == 4
      %output wavelength units
      fprintf(fid, 'wavelength = {\n');
      for i = 1:length(W)-1
          fprintf(fid, ' %f,\n', W(i));
      end
      fprintf(fid, ' %f\n', W(length(W)));
      fprintf(fid, '}\n');
  end
  
  fclose(fid);