Skip to content

Commit 7f071da

Browse files
zvezdochiotzvezdochiot
committed
1.0: init version
0 parents  commit 7f071da

8 files changed

+518
-0
lines changed

‎.gitmodules

+3
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,3 @@
1+
[submodule "src/stb"]
2+
path = src/stb
3+
url = https://github.com/nothings/stb

‎LICENSE

+24
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,24 @@
1+
This is free and unencumbered software released into the public domain.
2+
3+
Anyone is free to copy, modify, publish, use, compile, sell, or
4+
distribute this software, either in source code form or as a compiled
5+
binary, for any purpose, commercial or non-commercial, and by any
6+
means.
7+
8+
In jurisdictions that recognize copyright laws, the author or authors
9+
of this software dedicate any and all copyright interest in the
10+
software to the public domain. We make this dedication for the benefit
11+
of the public at large and to the detriment of our heirs and
12+
successors. We intend this dedication to be an overt act of
13+
relinquishment in perpetuity of all present and future rights to this
14+
software under copyright law.
15+
16+
THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
17+
EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
18+
MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
19+
IN NO EVENT SHALL THE AUTHORS BE LIABLE FOR ANY CLAIM, DAMAGES OR
20+
OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
21+
ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
22+
OTHER DEALINGS IN THE SOFTWARE.
23+
24+
For more information, please refer to <https://unlicense.org>

‎Makefile

+12
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,12 @@
1+
PNAME = stbicomix
2+
CFLAGS = -std=c99 -O2 -Wall -Wextra -Wno-unused-but-set-variable -Wno-unused-parameter -Werror
3+
LDLIBS = -lm -s
4+
SRCS = src/comix.c
5+
6+
all: $(PNAME)
7+
8+
$(PNAME): $(SRCS)
9+
$(CC) $(CFLAGS) $^ $(LDLIBS) -o $@
10+
11+
clean:
12+
rm -f $(PNAME)

‎README.md

+46
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,46 @@
1+
# stbicomix
2+
3+
Comix filter based of Gauss blur.
4+
Used:
5+
* Gauss blur: [iir_gauss_blur](https://github.com/arkanis/iir_gauss_blur).
6+
* STB: [stb](https://github.com/nothings/stb).
7+
8+
The Comix filter works on the difference between the unsigned difference
9+
of the original and blurry image and the blurry version of the this unsigned difference.
10+
11+
This filter was first applied in [STEX](https://github.com/ImageProcessing-ElectronicPublications/scantailor-experimental) (2025)
12+
in a single-component version (Y, the values of the color components were aligned
13+
in accordance with the brightness values before and after the filter).
14+
15+
Here this filter is implemented in a full-color version.
16+
17+
## Usage
18+
19+
`./stbicomix [-h] [-s sigma] [-m mixed] input-file output-file`
20+
21+
`-s sigma` The sigma of the gauss normal distribution (number >= 0.5).
22+
Larger values result in a stronger blur.
23+
24+
`-m mixed` The mixed coefficient.
25+
26+
`-h` display this help and exit.
27+
28+
You can use either sigma to specify the strengh of the blur.
29+
30+
The performance is independent of the blur strengh (sigma). This tool is an
31+
implementation of the paper "Recursive implementaion of the Gaussian filter"
32+
by Ian T. Young and Lucas J. van Vliet.
33+
34+
stb_image and stb_image_write by Sean Barrett and others is used to read and
35+
write images.
36+
37+
## Installation
38+
39+
- Clone the repo or download the source `git clone --recurse-submodules https://github.com/ImageProcessing-ElectronicPublications/stbicomix`
40+
- Execute `make`
41+
- Done. Either use the `stbicomix` executable directly or copy it somewhere in your PATH.
42+
43+
## Links
44+
45+
* STB: [stb](https://github.com/nothings/stb).
46+
* Gauss blur: [iir_gauss_blur](https://github.com/arkanis/iir_gauss_blur).

‎src/comix.c

+128
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,128 @@
1+
#define _GNU_SOURCE
2+
#include <stdio.h>
3+
#include <unistd.h>
4+
5+
#define STB_IMAGE_IMPLEMENTATION
6+
#define STBI_FAILURE_USERMSG
7+
#include "stb/stb_image.h"
8+
#define STB_IMAGE_WRITE_IMPLEMENTATION
9+
#include "stb/stb_image_write.h"
10+
#define IIR_GAUSS_BLUR_IMPLEMENTATION
11+
#include "iir_gauss_blur.h"
12+
#define FILTER_COMIX_IMPLEMENTATION
13+
#include "filter_comix.h"
14+
15+
void usage(char* progname)
16+
{
17+
fprintf(stderr,
18+
"%s %s %s\n",
19+
"Usage:", progname, "[-h] [-s sigma] [-m mixed] input-file output.png\n"
20+
"Comix filter an image and save it as PNG.\n"
21+
);
22+
}
23+
24+
void help(float sigma, float mix)
25+
{
26+
fprintf(stderr,
27+
"%s %f %s %f %s\n",
28+
" -s sigma The sigma of the gauss normal distribution (number >= 0.5, default =", sigma, ").\n"
29+
" Larger values result in a stronger blur.\n"
30+
" -m mixed The mixed coefficient (number, default =", mix, ").\n"
31+
" -h display this help and exit.\n"
32+
"\n"
33+
"You can use either sigma to specify the strengh of the blur.\n"
34+
"\n"
35+
"The performance is independent of the blur strengh (sigma). This tool is an\n"
36+
"implementation of the paper \"Recursive implementaion of the Gaussian filter\"\n"
37+
"by Ian T. Young and Lucas J. van Vliet.\n"
38+
"\n"
39+
"stb_image and stb_image_write by Sean Barrett and others is used to read and\n"
40+
"write images.\n"
41+
);
42+
}
43+
44+
unsigned char* image_copy(unsigned int width, unsigned int height, unsigned char components, unsigned char* image)
45+
{
46+
size_t image_size = height * width * components;
47+
unsigned char* dest = (unsigned char*)malloc(image_size * sizeof(unsigned char));
48+
if (dest != NULL)
49+
{
50+
for (size_t i = 0; i < image_size; i++)
51+
{
52+
dest[i] = image[i];
53+
}
54+
}
55+
return dest;
56+
}
57+
58+
int main(int argc, char** argv)
59+
{
60+
float sigma = 10.0f;
61+
float mix = 1.0f;
62+
63+
int opt;
64+
while ( (opt = getopt(argc, argv, "s:m:h")) != -1 )
65+
{
66+
switch(opt)
67+
{
68+
case 's':
69+
sigma = strtof(optarg, NULL);
70+
break;
71+
case 'm':
72+
mix = strtof(optarg, NULL);
73+
break;
74+
case 'h':
75+
usage(argv[0]);
76+
help(sigma, mix);
77+
return 0;
78+
default:
79+
usage(argv[0]);
80+
return 1;
81+
}
82+
}
83+
84+
// Need at least two filenames after the last option
85+
if (argc < optind + 2)
86+
{
87+
usage(argv[0]);
88+
return 1;
89+
}
90+
91+
int width = 0, height = 0, components = 1;
92+
unsigned char* image = stbi_load(argv[optind], &width, &height, &components, 0);
93+
if (image == NULL)
94+
{
95+
fprintf(stderr, "Failed to load %s: %s.\n", argv[optind], stbi_failure_reason());
96+
return 2;
97+
}
98+
99+
unsigned char* blur = image_copy(width, height, components, image);
100+
if (blur == NULL)
101+
{
102+
fprintf(stderr, "ERROR: not use memmory\n");
103+
return 3;
104+
}
105+
106+
unsigned char* delta = image_copy(width, height, components, image);
107+
if (delta == NULL)
108+
{
109+
fprintf(stderr, "ERROR: not use memmory\n");
110+
return 3;
111+
}
112+
113+
iir_gauss_blur(width, height, components, blur, sigma);
114+
115+
image_deltaabs(width, height, components, image, blur, delta);
116+
117+
iir_gauss_blur(width, height, components, delta, sigma);
118+
119+
image_filter_comix(width, height, components, mix, image, blur, delta);
120+
121+
if ( stbi_write_png(argv[optind+1], width, height, components, image, 0) == 0 )
122+
{
123+
fprintf(stderr, "Failed to save %s.\n", argv[optind+1]);
124+
return 4;
125+
}
126+
127+
return 0;
128+
}

‎src/filter_comix.h

+90
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,90 @@
1+
/**
2+
3+
Comix filter an image v1.0
4+
By zvezdochiot <[email protected]>
5+
This is free and unencumbered software released into the public domain.
6+
7+
QUICK START
8+
9+
#include ...
10+
#include ...
11+
#define FILTER_COMIX_IMPLEMENTATION
12+
#include "filter_comix.h"
13+
...
14+
unsigned int width = 0, height = 0;
15+
unsigned int widthb = 0, heightb = 0;
16+
unsigned int widthd = 0, heightd = 0;
17+
unsigned char components = 1;
18+
unsigned char componentsb = 1, componentsd = 1;
19+
float mix = 1.0f;
20+
21+
unsigned char* image = stbi_load("foo.png", &width, &height, &components, 0);
22+
unsigned char* blur = stbi_load("foo_blur.png", &widthb, &heightb, &componentsb, 0);
23+
unsigned char* delta = stbi_load("foo_delta.png", &widthd, &heightd, &componentsd, 0);
24+
25+
if ((width == widthb) && (height == heightb) && (components == componentsb) && (width == widthd) && (height == heightd) && (components == componentsd))
26+
{
27+
image_filter_comix(width, height, components, mix, image, blur, delta);
28+
stbi_write_png("foo.comix.png", width, height, components, image, 0);
29+
}
30+
31+
VERSION HISTORY
32+
33+
1.0 2025-01-30 "init" Initial release.
34+
35+
**/
36+
37+
#ifndef FILTER_COMIX_H
38+
#define FILTER_COMIX_H
39+
#ifdef __cplusplus
40+
extern "C" {
41+
#endif
42+
43+
void image_deltaabs(unsigned int width, unsigned int height, unsigned char components, unsigned char* src, unsigned char* blur, unsigned char* delta);
44+
void image_filter_comix(unsigned int width, unsigned int height, unsigned char components, float mix, unsigned char* src, unsigned char* blur, unsigned char* delta);
45+
46+
#ifdef __cplusplus
47+
}
48+
#endif
49+
50+
#ifdef FILTER_COMIX_IMPLEMENTATION
51+
#include <stdlib.h>
52+
#include <math.h>
53+
54+
void image_deltaabs(unsigned int width, unsigned int height, unsigned char components, unsigned char* src, unsigned char* blur, unsigned char* delta)
55+
{
56+
size_t image_size = height * width * components;
57+
if ((src != NULL) && (blur != NULL) && (delta != NULL))
58+
{
59+
for (size_t i = 0; i < image_size; i++)
60+
{
61+
unsigned char origin = src[i];
62+
unsigned char mean = blur[i];
63+
delta[i] = (origin < mean) ? (mean - origin) : (origin - mean);
64+
}
65+
}
66+
}
67+
68+
void image_filter_comix(unsigned int width, unsigned int height, unsigned char components, float mix, unsigned char* src, unsigned char* blur, unsigned char* delta)
69+
{
70+
size_t image_size = height * width * components;
71+
if ((src != NULL) && (blur != NULL) && (delta != NULL))
72+
{
73+
for (size_t i = 0; i < image_size; i++)
74+
{
75+
float origin = src[i];
76+
float mean = blur[i];
77+
float rebound = delta[i];
78+
float dlt = mean - origin;
79+
float cre = (dlt < 0.0f) ? (dlt + rebound) : (dlt - rebound);
80+
float fre = origin + cre;
81+
82+
float retval = mix * fre + (1.0f - mix) * origin + 0.5f;
83+
retval = (retval < 0.0f) ? 0.0f : (retval < 255.0f) ? retval : 255.0f;
84+
src[i] = (unsigned char) retval;
85+
}
86+
}
87+
}
88+
89+
#endif /* FILTER_COMIX_IMPLEMENTATION */
90+
#endif /* FILTER_COMIX_H */

0 commit comments

Comments
 (0)