<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">Hey!</span><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">I am re-sending this email, I think it got blocked because of attached images jpg examples, if you want to see it, I will send it personally. </div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">This is my first post, unfortunate it's not to help somebody but to ask for some help.. I do have to say, GRASS seems to be great so far and I feel my problems are probably user related. </div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">I am pre-processing Landsat data (at the moment only TM and ETM) to use for a classification later, doesn't matter. </div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Since all the data I'm acquiring comes as L1T (204/032), then the pre-processing workflow was done as follows:</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">DN -> TOA reflectance(&temperature) -> Cloud detection and masking (if necessary) -> Topographic correction using ASTERM 30m DEM. </div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">All data projection is as delivered by USGS (I don't remember if I reprojected the DEM but most probably I did).</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Imagery dates - 204/32:</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div>
<div style="font-family:arial,sans-serif;font-size:13px">
TM4 - 1989001 (January 1, 1989)</div><div style="font-family:arial,sans-serif;font-size:13px">TM5 - 2007027 (January 27,2007)</div><div style="font-family:arial,sans-serif;font-size:13px">ETM+ - 2003023 (January 23,2003)</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">All until here seems fine I guess, but I am having the following issues:</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px"><ul><li style="margin-left:15px">Band seems saturated or affected by some "multiplier". Such as TM4 composites look "normal" or as expected but TM5 and ETM+ seem pinkish or saturated as if one of the bands is saturated. </li>
<li style="margin-left:15px">Since TM4 looks the most "normal", all the others seem affected by some kind of noise (white noise?). Should I apply a Low pass filter of some kind? </li><li style="margin-left:15px">
Both i.landsat.rgb (runs endlessly) and i.landsat.dehaze (says it can't find the GUI, I've re-installed it without success)</li></ul></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">
If you need added information, please ask and I'll send you the commands. Here are the ones which i have been using.</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">
Added note: I'm using Windows 8! (don't know how important that can be)</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">-------------------</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Calling the file</div><div style="font-family:arial,sans-serif;font-size:13px">r.in.gdal input=C:\Invader_B\GRASS_Workfolder_01\01_LT5\LT52040322007027MPS00\LT52040322007027MPS00_B1.TIF output=RAW_2007027_LT5_B.1</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">(...)</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">
(i used DOS2 in this latest attempt since both uncorrected and DOS1 yield similar problems)</div><div style="font-family:arial,sans-serif;font-size:13px">i.landsat.toar sensor=tm5 method=dos2 --overwrite --verbose input_prefix=RAW_2007027_LT5_B. output_prefix=2007027_LT5_B. metfile=C:\Invader_B\GRASS_Workfolder_01\01_LT5\LT52040322007027MPS00\LT52040322007027MPS00_MTL.txt</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">i.landsat.acca --overwrite --verbose -5 input_prefix=2007027_LT5_B. output=2007027_CloudMask<br>
</div><div style="font-family:arial,sans-serif;font-size:13px">(The mask is not used by me since i visually see there are no clouds [there is snow in Serra da Estrela, famous place here])<br></div><div style="font-family:arial,sans-serif;font-size:13px">
<br></div><div style="font-family:arial,sans-serif;font-size:13px">Getting the illumination raster</div><div style="font-family:arial,sans-serif;font-size:13px">i.topo.corr -i --overwrite base=ASTER_30m zenith=<a href="tel:63.09559142" value="+16309559142" target="_blank">63.09559142</a> azimuth=153.85758094 out=ASTER.illu.prj<br>
</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px"><div>i.topo.corr base=ASTER.illu.prj zenith=<a href="tel:63.09559142" value="+16309559142" target="_blank">63.09559142</a> method=c-factor --verbose --overwrite input=2007027_LT5_B.1,2007027_LT5_B.2,2007027_LT5_B.3,2007027_LT5_B.4,2007027_LT5_B.5,2007027_LT5_B.6,2007027_LT5_B.7 out=npca_topo_CM</div>
<div><br></div><div>r.colors -e map=npca_topo_CM.2007027_LT5_B.1@01_LT5 color=grey<br></div><div>(...)</div><div><br></div><div>d.rgb to display the image in attachment.</div></div><div><br></div><div>--------------------</div>
<div><br></div><div>Pop-up window box text for i.landsat.dehaze:</div><div><br></div><div>"unable to fetch interface description for command i.landsat.dehaze.bat". </div><div>Details: (...) ./scripts/i.landsat.dehaze: no such file or directory</div>
<div><br></div>-- <br><div dir="ltr"><br>Nuno César de Sá<br>+351 91 961 90 37<br><br></div>
</div>