<div dir="ltr">Hello Michele,<div><br></div><div>I haven't tried myself, but have you seen Markus's blog from some time ago?? </div><div><br></div><div><a href="http://courses.neteler.org/processing-landsat-8-data-in-grass-gis-7-rgb-composites-and-pan-sharpening/">http://courses.neteler.org/processing-landsat-8-data-in-grass-gis-7-rgb-composites-and-pan-sharpening/</a></div><div><br></div><div>Maybe useful :)</div><div><br></div><div>Best, </div><div>Vero</div></div><div class="gmail_extra"><br><div class="gmail_quote">2016-07-18 12:00 GMT-03:00 Michele Toma <span dir="ltr"><<a href="mailto:michele.toma@airborne.aero" target="_blank">michele.toma@airborne.aero</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div lang="EN-US" link="blue" vlink="purple">
<div>
<p class="MsoNormal"><span style="color:#1f497d">Hi everyone,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="color:#1f497d">I was wondering if anyone had any suggestions or workarounds to the issue that I am having with pansharpening Landsat imagery? My email below has the steps that I took to attempt to pansharpen, but once I create
a natural color image and run the pansharpen tool, I end up with distorted colors. Are there any other methods that I should try?<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="color:#1f497d">Thanks in advance,<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="color:#1f497d">Michele<u></u><u></u></span></p>
<p class="MsoNormal"><span style="color:#1f497d"><u></u> <u></u></span></p>
<div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">From:</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> grass-user [mailto:<a href="mailto:grass-user-bounces@lists.osgeo.org" target="_blank">grass-user-bounces@lists.osgeo.org</a>]
<b>On Behalf Of </b>Michele Toma<br>
<b>Sent:</b> Wednesday, July 13, 2016 5:39 PM<br>
<b>To:</b> <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<b>Subject:</b> [GRASS-user] Pansharpening Landsat 8 Imagery<u></u><u></u></span></p>
</div>
</div><div><div class="h5">
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Hi everyone,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I am having difficulty with pansharpening Landsat 8 imagery using this tutorial:
<a href="http://planet.qgis.org/planet/tag/landsat/" target="_blank">http://planet.qgis.org/planet/tag/landsat/</a>. I am able to create a natural color image, but the pansharpening results in images with distorted colors.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Here are the steps that I have followed:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p><u></u><span>1.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>Import all Landsat 8 bands to GRASS 7.0.4<u></u><u></u></p>
<p style="margin-left:1.0in">
<u></u><span>a.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>r.in.gdal input=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_B1.TIF output=LC81100362016082LGN00_B1 -k -e<u></u><u></u></p>
<p><u></u><span>2.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>Converted digital number data to reflectance using i.landsat.toar<u></u><u></u></p>
<p style="margin-left:1.0in">
<u></u><span>a.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>i.landsat.toar input=LC81100362016082LGN00_B output=LC81100362016082LGN00_refl metfile=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_MTL.txt sensor=oli8 date=2016-03-22<u></u><u></u></p>
<p><u></u><span>3.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>Color balanced using i.colors.enhance for bands 4, 3, 2 as RGB respectively<u></u><u></u></p>
<p style="margin-left:1.0in">
<u></u><span>a.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>i.colors.enhance red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka<u></u><u></u></p>
<p><u></u><span>4.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>Displayed RGB map using bands 4, 3, 2 (I get a really great natural color image)<u></u><u></u></p>
<p><u></u><span>5.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>Pansharpened bands 4, 3, 2 using band 8 as the panchromatic band with the IHS method<u></u><u></u></p>
<p style="margin-left:1.0in">
<u></u><span>a.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>i.pansharpen red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka pan=LC81100362016082LGN00_refl8@Osaka output=LC81100362016082LGN00_ihs method=ihs<u></u><u></u></p>
<p><u></u><span>6.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>Displayed RGB map with the “pansharpened” bands for 4, 3, 2<u></u><u></u></p>
<p><u></u><span>7.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>Result: An image that is entirely gray (I found this link that says to use i.fusion.hpf, but following those instructions, I get “'make' is not recognized as an internal or external command, operable program or batch file.”
<a href="http://gis.stackexchange.com/questions/111763/pan-sharpen-landsat-8-using-grass-7" target="_blank">
http://gis.stackexchange.com/questions/111763/pan-sharpen-landsat-8-using-grass-7</a> Am I supposed to use the GUI for this?)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">From searching online, I have found that the i.pansharpen module doesn’t take 16 bit imagery (which Landsat imagery is) but rather 8 bit. I found this link (<a href="http://gis.stackexchange.com/questions/39174/how-to-color-balance-false-composites-of-pan-sharpened-landsat-images-in-grass/67319#67319" target="_blank">http://gis.stackexchange.com/questions/39174/how-to-color-balance-false-composites-of-pan-sharpened-landsat-images-in-grass/67319#67319</a>)
that says to rescale the image. I attempted this twice:<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p><u></u><span>1.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>When I run it on the bands produced after running the i.landsat.toar module, I get numbers such as min=-0.00813278402972343 and max=1.58301384751301. After rescaling the images to 0-255 and proceeding with step 3 I end up with an entirely
white image. <u></u><u></u></p>
<p><u></u><span>2.<span style="font:7.0pt "Times New Roman"">
</span></span><u></u>However, if I rescale the original bands which all have ranges of 0-65535, to 0-255, then pansharpen, the result is a washed out version of the natural color image. I believe this has been pansharpened, but the colors are off. I also
attempted to color balance the pansharpened bands but was unsuccessful as well.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Am I rescaling the images incorrectly or is there another step I am getting wrong? Any guidance would be greatly appreciated.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thank you,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Michele Toma<u></u><u></u></p>
<p class="MsoNormal"><span style="color:#1f497d"><a href="mailto:michele.toma@airborne.aero" target="_blank">michele.toma@airborne.aero</a></span><u></u><u></u></p>
<p class="MsoNormal"><span style="color:#1f497d"><a href="http://www.airborne.aero" target="_blank">www.airborne.aero</a></span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:8.0pt;color:#1f497d"> </span><u></u><u></u></p>
<p class="MsoNormal"><span style="font-size:14.0pt;color:#1f497d"><img border="0" width="133" height="36" src="cid:image001.jpg@01D1E0C9.E5979560" alt="New Small2"></span><u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div></div></div>
</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" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br></blockquote></div><br></div>