<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 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Aptos;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:11.0pt;
font-family:"Aptos",sans-serif;
color:#933131;
mso-ligatures:standardcontextual;
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0cm;
font-size:10.0pt;
font-family:"Courier New";
color:#933131;}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:"Consolas",serif;
color:#933131;
mso-ligatures:standardcontextual;
mso-fareast-language:EN-US;}
span.EmailStyle22
{mso-style-type:personal-reply;
font-family:"Aptos",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
mso-ligatures:none;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
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 bgcolor=white lang=EN-AU link=blue vlink=purple style='word-wrap:break-word'><div class=WordSection1><p class=MsoNormal><span style='color:windowtext'>Thanks Micha,<o:p></o:p></span></p><p class=MsoNormal><span style='color:windowtext'>That works but I don’t know why it is called “BundleToPerfectSensor”? I have searched for pansharpening algorithms previously and would never have guessed that one!<o:p></o:p></span></p><p class=MsoNormal><span style='color:windowtext'>Cheers Grant<o:p></o:p></span></p><p class=MsoNormal><span style='color:windowtext'>Perth WA<o:p></o:p></span></p><p class=MsoNormal><span style='color:windowtext'><o:p> </o:p></span></p><div><div style='border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal><b><span lang=EN-US style='font-family:"Calibri",sans-serif;color:windowtext;mso-ligatures:none;mso-fareast-language:EN-AU'>From:</span></b><span lang=EN-US style='font-family:"Calibri",sans-serif;color:windowtext;mso-ligatures:none;mso-fareast-language:EN-AU'> Micha Silver <tsvibar@gmail.com> <br><b>Sent:</b> Thursday, May 23, 2024 10:38 PM<br><b>To:</b> Grant Boxer <boxerg@iinet.net.au>; qgis-user@lists.osgeo.org<br><b>Subject:</b> Re: [Qgis-user] Pansharpening Images<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><span style='font-size:13.5pt'>In the past I've used the pan sharpening function in Orfeo Toolbox on WorldView 3 imagery. I chose the "bayes" method, not Gram-Schmidt that you want. Orfeo implements three sharpening algorithms:</span><span style='font-size:12.0pt;mso-ligatures:none;mso-fareast-language:EN-AU'><o:p></o:p></span></p><p class=MsoNormal><span style='font-size:13.5pt'><a href="https://www.orfeo-toolbox.org/CookBook/Applications/app_BundleToPerfectSensor.html">https://www.orfeo-toolbox.org/CookBook/Applications/app_BundleToPerfectSensor.html</a></span><o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><span style='font-size:13.5pt'>Here's the code I used (after preparing a *.vrt from both the PAN and the Multiband files)</span><o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><span style='font-size:13.5pt;font-family:"Courier New"'>WV_dir="....."<br>region="....."<br><br># Pan sharpening<br>otbcli_BundleToPerfectSensor -inp ${WV_dir}/${region}_PAN.vrt -inxs ${WV_dir}/${region}_Multi.vrt -out ${WV_dir}/${region}_pansharp_bayes.tif uint16 -method bayes -ram 8000 -progress true<br><br># Compress<br>gdal_translate -co COMPRESS=ZSTD -co PREDICTOR=2 ${WV_dir}/${region}_pansharp_bayes.tif ${WV_dir}/${region}_pansharp_bayes2.tif<br>mv ${WV_dir}/${region}_pansharp_bayes2.tif ${WV_dir}/${region}_pansharp_bayes.tif<br>gdaladdo -r average --config COMPRESS_OVERVIEW ZSTD --config PREDICTOR_OVERVIEW 2 -ro ${WV_dir}/${region}_pansharp_bayes.tif<br><br></span><span style='font-size:13.5pt'>HTH</span><o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><div><p class=MsoNormal>On 23/05/2024 6:43, Grant Boxer via QGIS-User wrote:<o:p></o:p></p></div><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><p class=MsoNormal>Does anyone know of a plug-in or algorithm for QGIS where I can pansharpen satellite imagery using the Gram-Schmidt method?<o:p></o:p></p><p class=MsoNormal>Cheers Grant<o:p></o:p></p><p class=MsoNormal>Perth WA<o:p></o:p></p><p class=MsoNormal> <o:p></o:p></p><p class=MsoNormal><span style='font-size:12.0pt;mso-ligatures:none;mso-fareast-language:EN-AU'><br><br><o:p></o:p></span></p><pre>_______________________________________________<o:p></o:p></pre><pre>QGIS-User mailing list<o:p></o:p></pre><pre><a href="mailto:QGIS-User@lists.osgeo.org">QGIS-User@lists.osgeo.org</a><o:p></o:p></pre><pre>List info: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><o:p></o:p></pre><pre>Unsubscribe: <a href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a><o:p></o:p></pre></blockquote><pre>-- <o:p></o:p></pre><pre>Micha Silver<o:p></o:p></pre><pre>Ben Gurion Univ.<o:p></o:p></pre><pre>Sde Boker, Remote Sensing Lab<o:p></o:p></pre><pre>cell: +972-523-665918<o:p></o:p></pre></div></body></html>