<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">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">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">jbenavid@siu.edu</a>>,<br></blockquote><blockquote type="cite"><span class="Apple-tab-span" style="white-space: pre; "> </span><<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>><br></blockquote><blockquote type="cite">Cc: Markus Metz <<a href="mailto:markus.metz.giswork@googlemail.com">markus.metz.giswork@googlemail.com</a>>,<span class="Apple-tab-span" style="white-space: pre; "> </span>Elena Mezzini<br></blockquote><blockquote type="cite"><span class="Apple-tab-span" style="white-space: pre; "> </span><<a href="mailto:elena.mezzini@gmail.com">elena.mezzini@gmail.com</a>><br></blockquote><blockquote type="cite">Message-ID:<br></blockquote><blockquote type="cite"><span class="Apple-tab-span" style="white-space: pre; "> </span><<a href="mailto:001F233294EA894387CFDAADAB87DE25085063E9@iwmix.IWMI.CGIARAD.ORG">001F233294EA894387CFDAADAB87DE25085063E9@iwmix.IWMI.CGIARAD.ORG</a>><br></blockquote><blockquote type="cite">Content-Type: text/plain; charset="us-ascii"<br></blockquote><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:jbenavid@siu.edu] <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">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:markus.metz.giswork@googlemail.com]<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">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">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">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">grass-user@lists.osgeo.org</a><br></blockquote><blockquote type="cite"><a href="http://lists.osgeo.org/mailman/listinfo/grass-user">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><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">http://lists.osgeo.org/pipermail/grass-user/attachments/20110527/211bcf66/attachment-0001.html</a></blockquote></body></html>