
<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Because I am doing different pan sharpening algorithms, I think I got the equation backward previously. Now I am using<div><br></div><div><div class="BodyFragment"><font size="2"><span style="font-size: 10pt;"><div class="PlainText">r.mapcalc "histmatched= ( pan - pan_mean ) / pan_sd * colorband_sd + colorband_mean</div><div class="PlainText"><br></div><div class="PlainText">where pan=high res panchromatic band and colorband=lower resolution restricted frequency bands (e.g., red, green, blue)</div><div class="PlainText"><br></div><div class="PlainText">This gives me a sharpened image. But I don't seem much difference from matching by the more standard cumulative distribution function except perhaps the mean/stdev method is a bit faster (and maybe especially so for bands with many grey levels).</div><div class="PlainText"><br></div><div class="PlainText">Here is the new code. You can just drop this into the scripts directory for $GRASS_ADDON_BASE to test.</div><div class="PlainText"><br></div><div class="PlainText"></div></span></font></div></div></body></html>