retinex 0.1.0

Single-scale and multi-scale Retinex image enhancement with color restoration (MSRCR)
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148

use clap::{Parser, ValueEnum};
use retinex::{
    extract_illumination, multi_scale_retinex, multi_scale_retinex_color_restored,
    multi_scale_retinex_full, normalize_reflectance, single_scale_retinex,
    single_scale_retinex_color_restored, single_scale_retinex_full,
};
use std::path::PathBuf;

#[derive(Parser, Debug)]
#[command(name = "retinex")]
#[command(about = "Apply Retinex enhancement to images", long_about = None)]
struct Cli {
    #[arg(help = "Input image file path")]
    input: PathBuf,

    #[arg(help = "Output image file path (reflectance result)")]
    output: PathBuf,

    #[arg(short, long, value_enum, default_value = "single")]
    mode: Mode,

    #[arg(long, value_delimiter = ',', default_value = "15.0")]
    sigmas: Vec<f32>,

    /// Save the estimated illumination to a separate file
    #[arg(long)]
    illumination: Option<PathBuf>,

    /// Save the raw reflectance (before color restoration) to a separate file
    #[arg(long)]
    reflectance: Option<PathBuf>,

    /// Enable color restoration (MSRCR) to prevent grayscale output
    #[arg(long, default_value_t = false)]
    color_restore: bool,
}

#[derive(Copy, Clone, Debug, PartialEq, Eq, ValueEnum)]
enum Mode {
    Single,
    Multi,
}

fn main() {
    let cli = Cli::parse();

    if cli.sigmas.is_empty() {
        eprintln!("Sigmas cannot be empty");
        std::process::exit(1);
    }

    if let Some(&invalid) = cli.sigmas.iter().find(|&&sigma| sigma <= 0.0) {
        eprintln!("Sigmas must be positive, got {invalid}");
        std::process::exit(1);
    }

    let image = match image::open(&cli.input) {
        Ok(img) => img,
        Err(err) => {
            eprintln!("Failed to open {}: {err}", cli.input.display());
            std::process::exit(1);
        }
    };

    // Save illumination if requested
    if let Some(illum_path) = &cli.illumination {
        let illum_result = match cli.mode {
            Mode::Single => extract_illumination(&image, cli.sigmas[0]),
            Mode::Multi => {
                // For multi-scale, use the first sigma for illumination visualization
                // or we could average them - using first for simplicity
                extract_illumination(&image, cli.sigmas[0])
            }
        };

        match illum_result {
            Ok(illum) => {
                if let Err(err) = illum.save(illum_path) {
                    eprintln!(
                        "Failed to save illumination {}: {err}",
                        illum_path.display()
                    );
                    std::process::exit(1);
                }
                println!("Saved illumination to {}", illum_path.display());
            }
            Err(err) => {
                eprintln!("Failed to extract illumination: {err}");
                std::process::exit(1);
            }
        }
    }

    // Save raw reflectance if requested
    if let Some(refl_path) = &cli.reflectance {
        let refl_result = match cli.mode {
            Mode::Single => {
                let output = single_scale_retinex_full(&image, cli.sigmas[0]);
                output.map(|o| normalize_reflectance(&o.reflectance))
            }
            Mode::Multi => {
                let output = multi_scale_retinex_full(&image, &cli.sigmas);
                output.map(|o| normalize_reflectance(&o.reflectance))
            }
        };

        match refl_result {
            Ok(refl) => {
                if let Err(err) = refl.save(refl_path) {
                    eprintln!("Failed to save reflectance {}: {err}", refl_path.display());
                    std::process::exit(1);
                }
                println!("Saved raw reflectance to {}", refl_path.display());
            }
            Err(err) => {
                eprintln!("Failed to extract reflectance: {err}");
                std::process::exit(1);
            }
        }
    }

    // Process the main output
    let result = if cli.color_restore {
        match cli.mode {
            Mode::Single => single_scale_retinex_color_restored(&image, cli.sigmas[0]),
            Mode::Multi => multi_scale_retinex_color_restored(&image, &cli.sigmas),
        }
    } else {
        match cli.mode {
            Mode::Single => single_scale_retinex(&image, cli.sigmas[0]),
            Mode::Multi => multi_scale_retinex(&image, &cli.sigmas),
        }
    };

    match result {
        Ok(processed) => {
            if let Err(err) = processed.save(&cli.output) {
                eprintln!("Failed to save {}: {err}", cli.output.display());
                std::process::exit(1);
            }
            println!("Saved result to {}", cli.output.display());
        }
        Err(err) => {
            eprintln!("Processing failed: {err}");
            std::process::exit(1);
        }
    }
}