Sophie

Sophie

distrib > Mandriva > 8.2 > i586 > media > contrib > by-pkgid > 8b2b1fb157760a0d31e072e140388824 > files > 215

gri-2.8.0-1mdk.i586.rpm

<html>
<head>
<title>example6.gri</title>
</head>
<body bgcolor="#FFFFFF">
<pre>
<font color=#B22222># Example 6 -- Plot IR image of Gulf of Maine</font>

<font color=#B22222># Define characteristics of norda images</font>
<font color=#B22222># Note that the pixel to temperature conversion formula is</font>
<font color=#B22222>#</font>
<font color=#B22222>#   Temperature = 5C + pixel_value / 10</font>
<font color=#B22222>#</font>
<font color=#B22222># where pixel_value ranges from 0 to 255.  Thus, a pixel</font>
<font color=#B22222># value of 0 corresponds to a temperatuer of 5C, and</font>
<font color=#B22222># 255 corresponds to 30.5C; this is why the limits</font>
<font color=#B22222># \0val and \255val, for use by the `set image</font>
<font color=#B22222># range' command, take on these values.</font>
<font color=#0000EE>\0val =</font> <font color=#2F4F4F>"5"</font>			<font color=#B22222># 0 in image</font>
<font color=#0000EE>\255val =</font> <font color=#2F4F4F>"30.5"</font>		<font color=#B22222># 255 in image</font>
<font color=#0000EE>.rows. =</font> 128
<font color=#0000EE>.cols. =</font> 128
<font color=#0000EE>.pixel_width. =</font> 2
<font color=#0000EE>.km. =</font> {<font color=#9400D3>rpn</font><font color=#CDAD00> .cols.</font><font color=#CDAD00> .pixel_width.</font> *}

<font color=#B22222># get filenames</font>
query<font color=#0000EE> \filename</font> <font color=#2F4F4F>"Name image file" ("example6image.dat"</font>)
query<font color=#0000EE> \maskname</font> <font color=#2F4F4F>"Name mask  file" ("example6mask.dat"</font>)

<font color=#B22222># get data</font>
open<font color=#0000EE> \filename</font> binary uchar
set image range<font color=#0000EE> \0val</font><font color=#0000EE> \255val
</font>read image<font color=#CDAD00> .rows.</font><font color=#CDAD00> .cols.</font> box 0 0<font color=#CDAD00> .km.</font><font color=#CDAD00> .km.</font>
close
open<font color=#0000EE> \maskname</font> binary uchar
read image mask<font color=#CDAD00> .rows.</font><font color=#CDAD00> .cols.</font>
close

<font color=#B22222># find out what grayscale method to use</font>
query<font color=#0000EE> \histo</font> <font color=#2F4F4F>"Do histogram enhancement? (yes|no)" ("no"</font>)
query<font color=#0000EE> \minT</font>  <font color=#2F4F4F>"T/deg for white on page?          " ("10"</font>)
query<font color=#0000EE> \maxT</font>  <font color=#2F4F4F>"T/deg for black on page?          " ("15"</font>)
<font color=#0000EE>\incT =</font> <font color=#2F4F4F>"1"</font>

<font color=#B22222># set up scales. </font>
set x size 12.8
set y size 12.8
set x name <font color=#2F4F4F>"km"</font>
set y name <font color=#2F4F4F>"km"</font>
set x axis 0<font color=#CDAD00> .km.</font> 32
set y axis 0<font color=#CDAD00> .km.</font> 32

<font color=#B22222># plot image, grayscale, and histogram</font>
<font color=#9400D3>if</font> {<font color=#2F4F4F>"<font color=#0000EE>\histo" =</font>= "yes"</font>}
    set image grayscale using histogram black<font color=#0000EE> \maxT</font> white<font color=#0000EE> \minT
</font><font color=#9400D3>else</font>
    set image grayscale black<font color=#0000EE> \maxT</font> white<font color=#0000EE> \minT
</font><font color=#9400D3>end if</font>
draw image
draw image palette left<font color=#0000EE> \minT</font> right<font color=#0000EE> \maxT</font> increment<font color=#0000EE> \incT
</font>draw image histogram
<font color=#9400D3>if</font> {<font color=#2F4F4F>"<font color=#0000EE>\histo" =</font>= "yes"</font>}
    draw title <font color=#2F4F4F>"Example 6: grayscale histogram enhanced"</font>
<font color=#9400D3>else</font>
    draw title <font color=#2F4F4F>"Example 6: grayscale linear<font color=#0000EE> \minT</font> to<font color=#0000EE> \maxT"</font>
</font><font color=#9400D3>end if</font>

</pre>
</body>
</html>