home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 6 / 06.iso / b / b002 / 3.ddi / LANDSAT.PAS < prev    next >
Encoding:
Pascal/Delphi Source File  |  1988-03-17  |  3.1 KB  |  144 lines
  1. program landsat(input,output,inf,outf);
  2.  
  3.     (* (c) 1988 J.Ronald Eastman *)
  4.  
  5. const
  6.  
  7.     rec_length    = 4319;  (* record length - 1 [4319 for TM/3595 for MSS] *)
  8.     header_recs   = 1;     (* the number of header records *)
  9.     record_header = 32;    (* documentation on the left side of each record *)
  10.  
  11. type
  12.  
  13.     byte     = char;      (* defines byte data type -- one pixel per byte *)
  14.     scanline = record     (* this record holds a full scaline *)
  15.                   pixel : array[0..rec_length] of byte;
  16.                end;
  17.  
  18. var
  19.  
  20.     inf            : file of scanline;           (* the LANDSAT input file *)
  21.     outf           : file of byte;               (* the output file *)
  22.     infile,outfile : packed array[1..80] of char;
  23.     col1,col2      : integer;
  24.     row1,row2      : integer;
  25.     cols,rows      : integer;
  26.     min,max        : integer;
  27.  
  28. procedure get_parameters;
  29.  
  30. begin
  31.  
  32.     write(chr(27),'[2J');         (* these two lines clear the screen *)
  33.     writeln(chr(27),'[1;1H');
  34.  
  35.     writeln('              IDRISI Landsat CCT-X Format Extraction Program');
  36.     writeln('              ----------------------------------------------');
  37.     writeln;
  38.     writeln;
  39.  
  40.     write('Enter the name of the raw tape file      : ');
  41.     readln(infile);
  42.     writeln;
  43.     writeln;
  44.  
  45.     write('Enter a new name for the extracted image : ');
  46.     readln(outfile);
  47.     writeln;
  48.     writeln;
  49.  
  50.     write('Enter the start row                      : ');
  51.     readln(row1);
  52.     writeln;
  53.     writeln;
  54.  
  55.     write('Enter the finish row                     : ');
  56.     readln(row2);
  57.     writeln;
  58.     writeln;
  59.  
  60.     write('Enter the start column                   : ');
  61.     readln(col1);
  62.     col1:=col1+record_header;
  63.     writeln;
  64.     writeln;
  65.  
  66.     write('Enter the finish column                  : ');
  67.     readln(col2);
  68.     col2:=col2+record_header;
  69.     writeln;
  70.     writeln;
  71.  
  72.     rows:=(row2-row1)+1;
  73.     cols:=(col2-col1)+1;
  74.  
  75. end;
  76.  
  77. procedure extract;
  78.  
  79. var i,j    : integer;
  80.     raw    : scanline;
  81.  
  82. begin
  83.  
  84.     for i:=1 to header_recs do read(inf,raw);    (* by-pass the header *)
  85.  
  86.     max:=(0);
  87.     min:=(255);
  88.  
  89.     for i:=0 to row2 do begin
  90.         read(inf,raw);
  91.     if i>=row1 then begin
  92.     for j:=(col1) to (col2) do begin
  93.         if ord(raw.pixel[j])>max then max:=ord(raw.pixel[j]);
  94.         if ord(raw.pixel[j])<min then min:=ord(raw.pixel[j]);
  95.         write(outf,raw.pixel[j]);
  96.      end; (* j *)
  97.      end; (* if *)
  98.      end; (* i *)
  99.  
  100. end;
  101.  
  102. procedure do_it;
  103.  
  104. begin
  105.  
  106.     writeln;
  107.     writeln('Extraction in Progress ...');
  108.     writeln;
  109.  
  110.     open(inf,infile,old);
  111.     open(outf,outfile,new);
  112.     reset(inf);
  113.     rewrite(outf);
  114.  
  115.     extract;
  116.  
  117.     close(inf);
  118.     close(outf);
  119.  
  120. end;
  121.  
  122. procedure print_results;
  123.  
  124. begin
  125.  
  126.     writeln;
  127.     writeln;
  128.  
  129.     writeln('The number of rows is : ',rows);
  130.     writeln('The number of cols is : ',cols);
  131.     writeln('The minimum value  is : ',(min));
  132.     writeln('The maximum value  is : ',(max));
  133.  
  134. end;
  135.  
  136. begin
  137.  
  138.     get_parameters;
  139.     do_it;
  140.     print_results;
  141.  
  142. end.
  143.  
  144.