Function extract_band

pub(crate) fn extract_band(source: &Path, band: usize) -> PathBuf

Extracts a single band from a raster source into a VRT.

Uses GDAL’s VRT driver in-process instead of spawning a gdal_translate subprocess, which reduces overhead from ~100ms to ~0.1ms per band.