I'm sending the code I used to run the i.atcorr module (including the i.landsat.toar transformation) <br><br>(I also used the code sent by Yann Chemin (May 23) with the same results)<br><br>Of two files attached one has the code and other one the univariate stats output of the band 1 TOAR raster for one of the 1986 scenes <br>
<br>A copy of the original landsat used have been placed temporarily at the following location " <a href="http://www.briolat.org/assets">assets</a> "<br><br>hope somebody can help<br><br><br>thanks <br><br>juan<br>
<br><br>I'm using a desktop running Ubuntu 10.10 with a 64bit Intel(R) Core(TM) i7 CPU 860 @ 2.80GHz<br><br><br><br><br><div class="gmail_quote">2011/5/27 Gaspar Reyes Póndigo <span dir="ltr"><<a href="mailto:gos47@hotmail.com">gos47@hotmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div style="word-wrap: break-word;">I had same problem, i.atcorr not work with landsat 5TM band 1, 2 and sometimes with band 3<br>
I tried i.atcorr in different versions of grass, this is my results:<br><br>1. GRASS-6.4.1-2 under Mac Snow-Leopard 10.6.7: Not work band 1, 2, and sometimes with band 3 (for years 1992, 1993, 1995 and 1997).<br>2. GRASS-6.5. revision 46428 compiled in Mac Snow-Leopard 10.6.7: Not work band 1, 2, and sometimes with band 3 (for years 1992, 1993, 1995 and 1997)<br>
3. GRASS-6.5svn45719(2011) under Microsoft Windows Xp sp3: All the bands work,THIS IS VERY RARE <br><br>My images landsat 5 TM are from Tonameca, Oaxaca, México of years 1985 to 1999, and 2010 to 2011 (one image series for year), i am work with AOD and visibility value, and convert manually DN to radiance (a).<br>
<br>(a) Lλ = [(LMAXλ − LMINλ) / (Qcalmax − Qcalmin)] (Qcal − Qcalmin) + LMINλ <br><br>Gaspar Reyes Póndigo<br>Universidad del Mar, Oaxaca, México<br>Campus Puerto-Ángel<br>Carretera principal Puerto Ángel Zipolite km 1.5 <br>
C.P. 70902<br><br><blockquote type="cite">-------------------------------------------------------------------<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Message: 1<br></blockquote><blockquote type="cite">
Date: Fri, 27 May 2011 08:53:32 +0530<br></blockquote><blockquote type="cite">From: "Chemin, Yann (IWMI)" <<a href="mailto:Y.Chemin@cgiar.org" target="_blank">Y.Chemin@cgiar.org</a>><br></blockquote><blockquote type="cite">
Subject: RE: [GRASS-user] Re: i.atcorr returns all NULL values<br></blockquote><blockquote type="cite">To: "Juan Benavides Duque" <<a href="mailto:jbenavid@siu.edu" target="_blank">jbenavid@siu.edu</a>>,<br>
</blockquote><div class="im"><blockquote type="cite"><span style="white-space: pre-wrap;"> </span><<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>><br></blockquote></div><blockquote type="cite">
Cc: Markus Metz <<a href="mailto:markus.metz.giswork@googlemail.com" target="_blank">markus.metz.giswork@googlemail.com</a>>,<span style="white-space: pre-wrap;"> </span>Elena Mezzini<br></blockquote><blockquote type="cite">
<span style="white-space: pre-wrap;"> </span><<a href="mailto:elena.mezzini@gmail.com" target="_blank">elena.mezzini@gmail.com</a>><br></blockquote><blockquote type="cite">Message-ID:<br></blockquote><blockquote type="cite">
<span style="white-space: pre-wrap;"> </span><<a href="mailto:001F233294EA894387CFDAADAB87DE25085063E9@iwmix.IWMI.CGIARAD.ORG" target="_blank">001F233294EA894387CFDAADAB87DE25085063E9@iwmix.IWMI.CGIARAD.ORG</a>><br>
</blockquote><blockquote type="cite">Content-Type: text/plain; charset="us-ascii"<br></blockquote><div><div></div><div class="h5"><blockquote type="cite"><br></blockquote><blockquote type="cite">Hi Juan,<br></blockquote>
<blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Can you please give the download link for your images,<br></blockquote><blockquote type="cite">
<br></blockquote><blockquote type="cite">Could you also check the output of i.landsat.toar for them and tell if<br></blockquote><blockquote type="cite">reflectance ranges are valid/logic.<br></blockquote><blockquote type="cite">
<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Thank you,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Yann<br>
</blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">From: Juan Benavides Duque [mailto:<a href="mailto:jbenavid@siu.edu" target="_blank">jbenavid@siu.edu</a>] <br>
</blockquote><blockquote type="cite">Sent: Thursday, May 26, 2011 11:38 PM<br></blockquote><blockquote type="cite">To: Chemin, Yann (IWMI)<br></blockquote><blockquote type="cite">Cc: Markus Metz; Elena Mezzini<br></blockquote>
<blockquote type="cite">Subject: Re: [GRASS-user] Re: i.atcorr returns all NULL values<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote>
<blockquote type="cite"><br></blockquote><blockquote type="cite">I tried Yann's code for the i.atcorr on some Landsat TM5 images<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">It worked very nice for images taken after 1995 but before that still<br>
</blockquote><blockquote type="cite">returns null values when running i.atcorr for bands 1 and 2 <br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I used TM5 scenes from the same location for 1986, 1989, and 1996 and<br>
</blockquote><blockquote type="cite">the first two (1986 and 1989) had the same problem with the null<br></blockquote><blockquote type="cite">values... <br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">
sorry I can provide a way to fix it but I hope somebody will<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">juan c <br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">
<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Tue, May 24, 2011 at 3:49 AM, Chemin, Yann (IWMI)<br>
</blockquote><blockquote type="cite"><<a href="mailto:Y.Chemin@cgiar.org" target="_blank">Y.Chemin@cgiar.org</a>> wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><blockquote type="cite">
Each created visibility map holds only one value r.mapcalc<br></blockquote></blockquote><blockquote type="cite">expression="visibility=${vis_list[$i]}" --overwrite so why not just use<br></blockquote><blockquote type="cite">
this visibility value as input in the >parameter file instead of<br></blockquote><blockquote type="cite">creating a map where all cells have the same value?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">
Yes, this is a valid point, the script is designed to (hopefully soon)<br></blockquote><blockquote type="cite">encompass raster maps of visibility data (maybe interpolated...), so<br></blockquote><blockquote type="cite">this is a placeholder for the next steps.<br>
</blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">-----Original Message-----<br></blockquote><blockquote type="cite">From: Markus Metz [mailto:<a href="mailto:markus.metz.giswork@googlemail.com" target="_blank">markus.metz.giswork@googlemail.com</a>]<br>
</blockquote><blockquote type="cite">Sent: Tuesday, May 24, 2011 12:52 PM<br></blockquote><blockquote type="cite">To: Chemin, Yann (IWMI)<br></blockquote><blockquote type="cite">Cc: <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>; Elena Mezzini; Markus Neteler<br>
</blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Subject: Re: [GRASS-user] Re: i.atcorr returns all NULL values<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">
On Mon, May 23, 2011 at 9:36 PM, Markus Neteler <<a href="mailto:neteler@osgeo.org" target="_blank">neteler@osgeo.org</a>><br></blockquote><blockquote type="cite">wrote:<br></blockquote><blockquote type="cite"><blockquote type="cite">
Hi Yann,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">On Mon, May 23, 2011 at 5:07 AM, Chemin, Yann (IWMI)<br></blockquote>
</blockquote><blockquote type="cite"><<a href="mailto:Y.Chemin@cgiar.org" target="_blank">Y.Chemin@cgiar.org</a>> wrote:<br></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">p192 r030 image of July 2003 of Italy (L5TM) successfully corrected<br>
</blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">for band 1 and 2.<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">
The following script works well in GRASS GIS Trunk SVN.<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">
(better send as attachment since a series of line breaks got broken)<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">
...<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite"># For i.atcorr scripting<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">
#-----------------------------------------------------<br></blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">vis_list=(10 10 8 9.7 15 8 7 10 10 9.7 12 9.7 7 12 12 12 3 15 12 9.7<br>
</blockquote></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><blockquote type="cite">6 15) vis_len=${#vis_list[*]} echo $vis_len<br></blockquote></blockquote></blockquote><blockquote type="cite">
<blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">... this is not clear to me: why so many values and where do they<br></blockquote></blockquote><blockquote type="cite">
<blockquote type="cite">originate from? Since we have only 6 channels to work with...<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite">Moreover, how can it be possible that visibility is different for<br>
</blockquote><blockquote type="cite">different channels of the same scene? From the documentation it seems<br></blockquote><blockquote type="cite">that visibility is related to aerosol concentration at 550nm which can<br>
</blockquote><blockquote type="cite">not possibility be different for different bands of the same scene.<br></blockquote><blockquote type="cite">Each created visibility map holds only one value r.mapcalc<br></blockquote><blockquote type="cite">
expression="visibility=${vis_list[$i]}" --overwrite so why not just use<br></blockquote><blockquote type="cite">this visibility value as input in the parameter file instead of creating<br></blockquote><blockquote type="cite">
a map where all cells have the same value?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Markus M<br></blockquote><blockquote type="cite">_______________________________________________<br>
</blockquote><blockquote type="cite">grass-user mailing list<br></blockquote><blockquote type="cite"><a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br></blockquote><blockquote type="cite">
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote>
<blockquote type="cite"><br></blockquote></div></div><blockquote type="cite">-------------- next part --------------<br></blockquote><blockquote type="cite">An HTML attachment was scrubbed...<br></blockquote><blockquote type="cite">
URL: <a href="http://lists.osgeo.org/pipermail/grass-user/attachments/20110527/211bcf66/attachment-0001.html" target="_blank">http://lists.osgeo.org/pipermail/grass-user/attachments/20110527/211bcf66/attachment-0001.html</a></blockquote>
</div><br>_______________________________________________<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>
<br></blockquote></div><br>