<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<style id="bidiui-paragraph-margins" type="text/css">body p { margin-bottom: 0cm; margin-top: 0pt; } </style>
</head>
<body bidimailui-detected-decoding-type="UTF-8"
bidimailui-charset-is-forced="true" style="direction: ltr;"
text="#933131" bgcolor="#FFFFFF">
<p><font size="4">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:</font></p>
<p><font size="4"><a class="moz-txt-link-freetext" href="https://www.orfeo-toolbox.org/CookBook/Applications/app_BundleToPerfectSensor.html">https://www.orfeo-toolbox.org/CookBook/Applications/app_BundleToPerfectSensor.html</a></font></p>
<p><font size="4"><br>
</font></p>
<p><font size="4">Here's the code I used (after preparing a *.vrt
from both the PAN and the Multiband files)</font></p>
<p><font size="4"><br>
</font></p>
<p><font size="4" face="monospace">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>
</font><font size="4">HTH</font><font size="4" face="monospace"><br>
</font></p>
<p><font size="4"><br>
</font></p>
<div class="moz-cite-prefix">On 23/05/2024 6:43, Grant Boxer via
QGIS-User wrote:<br>
</div>
<blockquote type="cite"
cite="mid:00b601daacc3$6934f750$3b9ee5f0$@iinet.net.au">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<meta name="Generator"
content="Microsoft Word 15 (filtered medium)">
<style>@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}@font-face
{font-family:Aptos;}p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Aptos",sans-serif;
mso-ligatures:standardcontextual;
mso-fareast-language:EN-US;}span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Aptos",sans-serif;
color:windowtext;}.MsoChpDefault
{mso-style-type:export-only;
font-size:11.0pt;
mso-fareast-language:EN-US;}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]-->
<div class="WordSection1">
<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>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
QGIS-User mailing list
<a class="moz-txt-link-abbreviated" href="mailto:QGIS-User@lists.osgeo.org">QGIS-User@lists.osgeo.org</a>
List info: <a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a>
Unsubscribe: <a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/qgis-user">https://lists.osgeo.org/mailman/listinfo/qgis-user</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918</pre>
</body>
</html>