wenmin-wu/cv-brovey-pansharpening
skills/cv/brovey-pansharpening/SKILL.md
Fuse low-resolution multispectral bands with a high-resolution panchromatic band using the Brovey transform to produce sharp multi-band imagery
$ install --global
skillsauthnpx skillsauth add wenmin-wu/ds-skills cv-brovey-pansharpeningInstall this skill globally with one command. Works with Claude Code, Cursor, and Windsurf.
Security Scan Results
3 of 9 scanners reported clean
Some scanners were skipped, did not run, or reported a non-clean status. Review each row below.
3
Scanners Passed
9
Scanners in report
Clean
TrivyContainer and dependency vulnerability scanner
95%
Clean
SemgrepStatic code analysis for vulnerabilities
95%
Clean
mcp-scan (Snyk)Model Context Protocol security validation
95%
Skipped
Snyk (dep)Open source security scanning
50%
Skipped
Socket.devSupply chain security analysis
50%
Skipped
VirusTotalMulti-engine malware detection
50%
Skipped
CrowdStrikeAdvanced threat intelligence
50%
Skipped
OSV-ScannerOpen Source Vulnerability database check
50%
Skipped
OWASP Dep-Check
50%
Last scanned: Apr 23, 2026, 6:49 AM13.8s1 file scanned
SKILL.md
- name:
- cv-brovey-pansharpening
- description:
- Fuse low-resolution multispectral bands with a high-resolution panchromatic band using the Brovey transform to produce sharp multi-band imagery
Brovey Pansharpening
Overview
Satellite sensors capture high-resolution panchromatic (single-band) and lower-resolution multispectral images. Pansharpening fuses both to get high-resolution color imagery. The Brovey transform normalizes each band by a weighted sum, then multiplies by the panchromatic intensity. This preserves spectral ratios while injecting spatial detail.
Quick Start
import numpy as np
from skimage.transform import rescale
def brovey_pansharpen(ms_bands, pan, nir_weight=0.1):
R, G, B, NIR = ms_bands
dnf = (pan - nir_weight * NIR) / (
nir_weight * R + nir_weight * G + nir_weight * B + 1e-10)
sharp = np.stack([R * dnf, G * dnf, B * dnf, NIR * dnf], axis=-1)
return np.clip(sharp, 0, None)
ms_upscaled = np.stack([
rescale(ms[i], 4, order=3) for i in range(4)], axis=-1)
ms_upscaled = ms_upscaled[:pan.shape[0], :pan.shape[1]]
result = brovey_pansharpen(
[ms_upscaled[..., i] for i in range(4)], pan)
Workflow
- Upscale multispectral bands to match panchromatic resolution (typically 4x bicubic)
- Crop to align dimensions if sizes differ by a pixel
- Compute per-pixel normalization factor from pan and weighted band sum
- Multiply each band by the normalization factor
- Clip negatives to zero
Key Decisions
- NIR weight (W): 0.1 is typical; higher values reduce NIR influence on the ratio
- vs HSV pansharpening: Brovey preserves spectral ratios better; HSV preserves hue but can shift brightness
- Upscale method: bicubic (order=3) balances sharpness and artifacts
- Band alignment: always crop after rescale — off-by-one pixels cause misregistration
References
- Panchromatic sharpening
Related Skills
wenmin-wu/timeseries-scaled-pinball-loss
data-ai
VerifiedTrustedCommunity
Scaled Pinball Loss (SPL) metric for evaluating quantile forecasts, normalized by mean absolute successive differences of training data
31SKILL.mdUpdated Apr 23, 2026
wenmin-wu/timeseries-retroactive-outlier-rescaling
data-ai
VerifiedTrustedCommunity
Walk backward through a time series and multiplicatively rescale segments when jumps exceed a fraction of the running mean to correct data collection anomalies
31SKILL.mdUpdated Apr 23, 2026
wenmin-wu/timeseries-ratio-target-for-smape
testing
VerifiedTrustedCommunity
Transform forecasting target to next/current ratio minus one so that optimizing MAE or squared error implicitly minimizes SMAPE
31SKILL.mdUpdated Apr 23, 2026
wenmin-wu/timeseries-quantile-ratio-scaling
tools
VerifiedTrustedCommunity
Convert point forecasts to prediction intervals by scaling with logit-transformed quantile ratios passed through a Normal CDF
31SKILL.mdUpdated Apr 23, 2026