<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
{mso-style-priority:99;
mso-style-link:"Balloon Text Char";
margin:0in;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Tahoma","sans-serif";}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
span.BalloonTextChar
{mso-style-name:"Balloon Text Char";
mso-style-priority:99;
mso-style-link:"Balloon Text";
font-family:"Tahoma","sans-serif";}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Hi Vero,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Thanks for the link. I believe the tutorial that I was originally following was the same as what you listed – just on a different site.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Paul Shapley was kind enough to share his tutorial (<a href="http://www.paulshapley.com/2014/04/processing-landsat-8-bands-in-grass-gis.html">http://www.paulshapley.com/2014/04/processing-landsat-8-bands-in-grass-gis.html</a>)
which has worked for me. (Thanks Paul!)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">For anyone that is interested, I followed Paul’s tutorial and was able to get a pansharpened image. I did end up color balancing (i.colors.enhance) as a final
step because urban areas were extremely washed out. This corrected the problem and I now have an image that has been pansharpened with natural colors.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">(Note: I used version 7.0.4 and installed the i.fusion.brovey extension as I believe this was replaced by i.pansharpen. To do this – Settings -> Addons extensions
-> Install extension from addons. Then search for i.fusion.brovey in the “imagery” drop down.)<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">Thanks again for all the help!<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D">-Michele<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D"><o:p> </o:p></span></p>
<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""> Veronica Andreo [mailto:veroandreo@gmail.com]
<br>
<b>Sent:</b> Monday, July 18, 2016 12:41 PM<br>
<b>To:</b> Michele Toma<br>
<b>Cc:</b> grass-user@lists.osgeo.org<br>
<b>Subject:</b> Re: [GRASS-user] Pansharpening Landsat 8 Imagery<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">Hello Michele,<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I haven't tried myself, but have you seen Markus's blog from some time ago?? <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><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><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Maybe useful :)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Best, <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Vero<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<p class="MsoNormal">2016-07-18 12:00 GMT-03:00 Michele Toma <<a href="mailto:michele.toma@airborne.aero" target="_blank">michele.toma@airborne.aero</a>>:<o:p></o:p></p>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#1F497D">Hi everyone,</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#1F497D"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><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?</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#1F497D"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#1F497D">Thanks in advance,</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#1F497D"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#1F497D">Michele</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#1F497D"> </span><o:p></o:p></p>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><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</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Hi everyone,<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">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.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Here are the steps that I have followed:<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p>1.<span style="font-size:7.0pt"> </span>Import all Landsat 8 bands to GRASS 7.0.4<o:p></o:p></p>
<p style="margin-left:1.0in">a.<span style="font-size:7.0pt"> </span>r.in.gdal input=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_B1.TIF output=LC81100362016082LGN00_B1 -k -e<o:p></o:p></p>
<p>2.<span style="font-size:7.0pt"> </span>Converted digital number data to reflectance using i.landsat.toar<o:p></o:p></p>
<p style="margin-left:1.0in">a.<span style="font-size:7.0pt"> </span>i.landsat.toar input=LC81100362016082LGN00_B output=LC81100362016082LGN00_refl metfile=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_MTL.txt sensor=oli8 date=2016-03-22<o:p></o:p></p>
<p>3.<span style="font-size:7.0pt"> </span>Color balanced using i.colors.enhance for bands 4, 3, 2 as RGB respectively<o:p></o:p></p>
<p style="margin-left:1.0in">a.<span style="font-size:7.0pt"> </span>i.colors.enhance red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka<o:p></o:p></p>
<p>4.<span style="font-size:7.0pt"> </span>Displayed RGB map using bands 4, 3, 2 (I get a really great natural color image)<o:p></o:p></p>
<p>5.<span style="font-size:7.0pt"> </span>Pansharpened bands 4, 3, 2 using band 8 as the panchromatic band with the IHS method<o:p></o:p></p>
<p style="margin-left:1.0in">a.<span style="font-size:7.0pt"> </span>i.pansharpen red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka pan=LC81100362016082LGN00_refl8@Osaka output=LC81100362016082LGN00_ihs
method=ihs<o:p></o:p></p>
<p>6.<span style="font-size:7.0pt"> </span>Displayed RGB map with the “pansharpened” bands for 4, 3, 2<o:p></o:p></p>
<p>7.<span style="font-size:7.0pt"> </span>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?)<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">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:<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p>1.<span style="font-size:7.0pt"> </span>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. <o:p></o:p></p>
<p>2.<span style="font-size:7.0pt"> </span>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.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Am I rescaling the images incorrectly or is there another step I am getting wrong? Any guidance would be greatly appreciated.<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Thank you,<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto">Michele Toma<o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#1F497D"><a href="mailto:michele.toma@airborne.aero" target="_blank">michele.toma@airborne.aero</a></span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="color:#1F497D"><a href="http://www.airborne.aero" target="_blank">www.airborne.aero</a></span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:8.0pt;color:#1F497D"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span style="font-size:14.0pt;color:#1F497D"><img border="0" width="133" height="36" id="_x0000_i1025" src="cid:image001.jpg@01D1E10B.B5E57A20" alt="New Small2"></span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"> <o:p></o:p></p>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><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><o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</body>
</html>