
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-7">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body dir="ltr">

<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Arial,Helvetica,sans-serif;" dir="ltr">

<p><span style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;">I am trying to apply atmospheric correction to Landsat 8 scenes using i.atcorr. I have read lots of discussion about it, but I cannot seem to figure out the right way to use

 it. After tons of tests I am still struggling with the correct input and output range.</span><br>

</p>

<div style="color: rgb(0, 0, 0);">

<div>

<div id="divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">

<p>The results deriving from radiance as input and from reflectance as input are different! The correct flags are used and the radiance and reflectance values are from i.landsat.toar and their validity is checked through the mapcalculator by computing them

 manually from the .mtl file data.</p>

<p>So, I have radiance and reflectance and try to get the i.atcorr results. I try several combination of input and output range and the only ones that seem to produce a nice RGB are:</p>

<p>Input: the full range of radiance (Lmin, Lmax)</p>

<p>Output: 0,1 or the full range of reflectance (computed by dividing Lmin and Lmax by the factor that converts radiance to reflectance and is computed by i.landsat.toar)</p>

<p>Or in case of reflectance</p>

<p>Input: Lmin, Lmax divided by the convertion factor mentioned above</p>

<p>Output: the same as radiance</p>

<p><br>

</p>

<p>In both cases the resulting values are different which cannot be correct (Radiance or reflectance should produce the same outcome) and when I examine the resulting RGB, in QGIS for example, the result is almost exactly the same as the digital number image

 before the atmospheric correction.</p>

<p>The only thing that i.atcorr is doing is shifting the histogram (that is why the output values are different) and of course a histogram shift produces the same visual result.</p>

<p><br>

</p>

<p>Does anyone really uses i.atcorr to perform atmospheric correction to Landsat? And if so, which is the right way?</p>

<p><br>

</p>

<p>Thanks in advance,</p>

<p>Mat</p>

</div>

</div>

</div>

</div>

</body>

</html>