<div dir="ltr"><div class="gmail_default" style="font-family:tahoma,sans-serif"><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Wed, Sep 4, 2013 at 10:17 AM, Nikos Alexandris <span dir="ltr"><<a href="mailto:nik@nikosalexandris.net" target="_blank">nik@nikosalexandris.net</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">On Friday 30 of August 2013 02:56:14 Beatrice wrote:<br>
<br>
> Dear All,<br>
<br>
Good Day Beatrice,<br>
<div><div class="h5"><br>
<br>
> I have been running i.landsat.toar but have just noted this problem. For me<br>
> it shows zeros both in the command output (verbose) and output metadata<br>
> (<a href="http://r.info" target="_blank">r.info</a>). An example for one of the bands:<br>
<br>
<br>
> From command output:<br>
><br>
>  LANDSAT: 5 SENSOR: TM<br>
>  ACQUISITION DATE 2011-07-05 [production date 2011-11-29]<br>
>    Earth-sun distance    = 0.00000000<br>
>    Solar elevation angle = 0.00000000<br>
>    Atmospheric correction = dos4<br>
>    Percent of solar irradiance in path radiance = 0.0000<br>
> -------------------<br>
>  BAND 1  (code 1)<br>
>    calibrated digital number (DN): 0.0 to 0.0<br>
>    calibration constants (L): 0.000 to 0.000<br>
>    at-surface radiance = 0.00000000 * DN + -0.000<br>
>    mean solar exoatmospheric irradiance (ESUN): -0.000<br>
>    at-surface reflectance = radiance / 0.00000<br>
>    the darkness DN with a least 1000 pixels is 1<br>
>    the DN mode is 53<br>
> -------------------<br>
<br>
<br>
> From <a href="http://r.info" target="_blank">r.info</a> for the same band:<br>
>  |     Earth-sun distance (d) ................ -0.00000000<br>
>  |<br>
>  |     Digital number (DN) range ............. 0 to 0<br>
>  |<br>
>  |     Calibration constants (Lmin to Lmax) .. +0.000 to +0.000<br>
>  |<br>
>  |     DN to Radiance (gain and bias) ........ +0.00000 and +0.00000<br>
>  |<br>
>  |     Mean solar irradiance (ESUN) .......... 0.000<br>
>  |<br>
>  |     Reflectance = Radiance divided by ..... 0.00000<br>
>  |<br>
>  |<br>
>  |<br>
>  |     Dark object (1000 pixels) DN = ........ 1<br>
>  |<br>
>  |     Mode in reflectance histogram ......... 0.00000<br>
<br>
<br>
</div></div>This doesn't look healthy.  However, I think it is a past problem which is<br>
fixed.<br>
<div><div class="h5"><br>
<br>
> Running i.landsat.toar -p results to<br>
><br>
> number=5<br>
> creation=2011-11-29<br>
> date=2011-07-05<br>
> sun_elev=45.585310<br>
> sensor=TM<br>
> bands=7<br>
> sunaz=44.259879<br>
><br>
> Which I think confirms that i.landsat.toar can read input values from my<br>
> metafile. The output bands have values in the range of surface<br>
> reflectance/temperature (for band 6). For example<br>
><br>
> Band 1 output - Range of data:    min = 0.01  max = 0.461925184292249, Data<br>
> Units:   unitless<br>
> Band 2 output - Range of data:    min = 0.01  max = 0.925906812458597, Data<br>
> Units:   unitless<br>
> Band 3 output - Range of data:    min = 0.01  max = 0.785829845850446, Data<br>
> Units:   unitless<br>
> .<br>
> .<br>
> Band 6 output - Range of data:    min = 203.371307953013  max =<br>
> 327.669907753679, Data Units:   Kelvin<br>
><br>
> I have same results when using MTL or MTLOld as source of metadata, on TM or<br>
> ETM+ sensors and different atmospheric correction methods (uncorrected,<br>
> DOS). I am running winGRASS7 revision 57466.<br>
><br>
> Is there a problem with the printed output from <a href="http://r.info/command" target="_blank">r.info/command</a> output or are<br>
> these zeros the actual values that i.landsat.toar uses as input to<br>
> computations?<br>
<br>
</div></div>I doubt that <a href="http://r.info" target="_blank">r.info</a> has a problem reading the metadata that pertain to a<br>
queried raster map.  If I am not wrong, it merely reads simple text files that<br>
reside inside the "hist" directory inside the corresponding Mapset.<br>
<br>
My guess is that you use an older version of i.landsat.toar whicj, indeed,<br>
produces empty raster maps.  Unfortunately, I do not use (nor have any access<br>
to a Windows machine) so as to be in position to help directly.  We need<br>
feedback from someone with a Windows OS.<br>
<br>
(<br>
On/Off-Topic:  is there a way to run and test GRASS in an emulated Windows<br>
environment for those who don't have a legally licensed Windows OS?<br>
)<br>
<br></blockquote><div><div class="gmail_default" style="font-family:tahoma,sans-serif">​There are two options that I am aware of:</div><div class="gmail_default" style="font-family:tahoma,sans-serif"><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif">
1) Create a virtual machine and install windows on it (meh..)<br><br></div><div class="gmail_default" style="font-family:tahoma,sans-serif">2) Use WINE, which as they state in their site  they are "<span style="color:rgb(0,0,0);font-family:'bitstream vera sans',verdana,arial,helvetica,sans-serif;font-size:13px">a compatibility layer capable of running Windows applications on several POSIX-compliant operating systems". <br>
<br></span></div><div class="gmail_default" style="font-family:tahoma,sans-serif"><span style="color:rgb(0,0,0);font-family:'bitstream vera sans',verdana,arial,helvetica,sans-serif;font-size:13px"><a href="http://www.winehq.org">http://www.winehq.org</a></span></div>
<br><div class="gmail_default" style="font-family:tahoma,sans-serif">​Nikos.​</div><br><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">

Best, Nikos<br>
<div class=""><div class="h5">_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</div></div></blockquote></div><br></div></div>