├── .gitignore ├── CMakeLists.txt ├── LICENSE ├── README.md ├── dirent ├── LICENSE └── dirent.h ├── examples ├── batch.bat ├── lenna.bat ├── lenna.png ├── lenna_1_0_0.png ├── lenna_osx.sh ├── lenna_test_16_1_1.png ├── lenna_test_1_0_0.png ├── lenna_test_1_1_1.png ├── lenna_test_2_1_1.png └── lenna_test_4_1_1.png ├── noi_image.h ├── stb ├── stb_image.h └── stb_image_write.h └── test ├── impl.c └── main.c /.gitignore: -------------------------------------------------------------------------------- 1 | # Prerequisites 2 | *.d 3 | 4 | # Object files 5 | *.o 6 | *.ko 7 | *.obj 8 | *.elf 9 | 10 | # Linker output 11 | *.ilk 12 | *.map 13 | *.exp 14 | 15 | # Precompiled Headers 16 | *.gch 17 | *.pch 18 | 19 | # Libraries 20 | *.lib 21 | *.a 22 | *.la 23 | *.lo 24 | 25 | # Shared objects (inc. Windows DLLs) 26 | *.dll 27 | *.so 28 | *.so.* 29 | *.dylib 30 | 31 | # Executables 32 | *.exe 33 | *.out 34 | *.app 35 | *.i*86 36 | *.x86_64 37 | *.hex 38 | 39 | # Debug files 40 | *.dSYM/ 41 | *.su 42 | *.idb 43 | *.pdb 44 | 45 | # Kernel Module Compile Results 46 | *.mod* 47 | *.cmd 48 | .tmp_versions/ 49 | modules.order 50 | Module.symvers 51 | Mkfile.old 52 | dkms.conf 53 | build/ 54 | -------------------------------------------------------------------------------- /CMakeLists.txt: -------------------------------------------------------------------------------- 1 | set(CMAKE_OSX_ARCHITECTURES "arm64") 2 | 3 | cmake_minimum_required(VERSION 3.16) 4 | 5 | # set the project name 6 | project(NotOkImageFormat) 7 | 8 | set(EXECUTABLE_OUTPUT_PATH ${PROJECT_SOURCE_DIR}/bin/) 9 | 10 | # add the executable 11 | add_executable(noi 12 | test/main.c 13 | test/impl.c 14 | noi_image.h 15 | stb/stb_image.h 16 | stb/stb_image_write.h 17 | ) 18 | target_compile_definitions(noi PUBLIC NOI_THREADS) 19 | target_include_directories(noi PUBLIC ${CMAKE_SOURCE_DIR}) 20 | 21 | if(UNIX AND NOT APPLE) 22 | set(THREADS_PREFER_PTHREAD_FLAG ON) 23 | find_package(Threads REQUIRED) 24 | target_link_libraries(noi PRIVATE Threads::Threads m) 25 | endif() 26 | -------------------------------------------------------------------------------- /LICENSE: -------------------------------------------------------------------------------- 1 | MIT License 2 | 3 | Copyright (c) 2021 Boris Batkin 4 | 5 | Permission is hereby granted, free of charge, to any person obtaining a copy 6 | of this software and associated documentation files (the "Software"), to deal 7 | in the Software without restriction, including without limitation the rights 8 | to use, copy, modify, merge, publish, distribute, sublicense, and/or sell 9 | copies of the Software, and to permit persons to whom the Software is 10 | furnished to do so, subject to the following conditions: 11 | 12 | The above copyright notice and this permission notice shall be included in all 13 | copies or substantial portions of the Software. 14 | 15 | THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 16 | IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 17 | FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE 18 | AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 19 | LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, 20 | OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE 21 | SOFTWARE. 22 | -------------------------------------------------------------------------------- /README.md: -------------------------------------------------------------------------------- 1 | # NotOkImageFormat 2 | Lossy fixed-rate GPU-friendly image compression\decompression. 3 | Supported profiles 4 | 5 | 16:1:1 2.8125 bpp yuv 6 | 4:1:1 3.75 bpp yuv 7 | 2:1:1 5.0 bpp yuv 8 | 1:1:1 7.5 bpp rgb 9 | 1:1:0 5 bpp rg (normal-maps) 10 | 1:0:0 2.5 bpp greyscale 11 | 12 | Tested on Windows (Windows 10, MSVC 2019 and Clang 12), Mac OSX (12.0 Monterey, Apple Clang 13), 13 | Linux (Ubuntu 20.04LTS, GCC-9). 14 | 15 | Currently 2-4.5 times faster decompression than STBI JPEG implementation, with a lot of potential to optimize. 16 | 17 | Some work in progress numbers from my M1 Max 2021 apple laptop: 18 | 19 | noi_compress 512 x 512 profile YUV_16_1_1 20 | 0 mb in 0.26 sec, 2.854mb/sec 21 | PSNR = -32.0 PSNR(YUV) = -37.6 22 | running noi_decompressing 100 times, 110600 bytes 23 | decompression speed 1609.4mb/sec 24 | 25 | noi_compress 512 x 512 profile YUV_4_1_1 26 | 0 mb in 0.35 sec, 2.148mb/sec 27 | PSNR = -33.3 PSNR(YUV) = -39.1 28 | running noi_decompressing 100 times, 141320 bytes 29 | decompression speed 1146.8mb/sec 30 | 31 | noi_compress 512 x 512 profile YUV_2_1_1 32 | 0 mb in 0.45 sec, 1.672mb/sec 33 | PSNR = -34.0 PSNR(YUV) = -39.8 34 | running noi_decompressing 100 times, 182280 bytes 35 | decompression speed 1117.7mb/sec 36 | 37 | noi_compress 512 x 512 profile RGB_1_1_1 38 | 0 mb in 0.65 sec, 1.148mb/sec 39 | PSNR = -36.6 PSNR(YUV) = -41.8 40 | running noi_decompressing 100 times, 264200 bytes 41 | decompression speed 1241.7mb/sec 42 | 43 | noi_compress 512 x 512 profile Y_1_0_0 44 | 0 mb in 0.35 sec, 2.162mb/sec 45 | PSNR = -37.7 PSNR(YUV) = -37.7 46 | running noi_decompressing 100 times, 100360 bytes 47 | decompression speed 2895.8mb/sec 48 | 49 | bash-3.2$ ../bin/noi -stbjpg lenna.png lenna.jpg 50 | running stbi_load_from_memory 100 times, 68593 bytes 51 | decompression speed 354.3mb/sec 52 | 53 | I finally got to implement this really old idea of mine, of combining a quantizer with Hadamard transform. 54 | 55 | This is how compression works: 56 | 57 | 1. RGB->YUV color conversion for the YUV profiles 58 | 2. 4x4 HDT 59 | 3. combined weight (0,0) is stored as is. there are 4 bits there which can be used for something 60 | 4. 'corners' of size 3, 5, and 7 of the 4x4 matrix are quantized with k-mean quantizer - down to 256 means 61 | 5. as the result we have 5 byte blocks - 2 bytes for weight, and 3 quantization indices 62 | 6. index pallet is stored as 256 entries of 15 2-byte coefficients 63 | 64 | This is what happens during decompression: 65 | 66 | 1. original blocks are restored from the 3 palette 67 | 2. 4x4 iHDT 68 | 3. YUV->RGB for the YUV profiles 69 | 70 | NOI is really fast to decompress, even on the CPU. GPU is probably fast enough to decompress as it textures. 71 | 72 | Compression can be speed-up significantly with better k-means implementation. However I would not want to waste any time on it. This really ought to be a shader. GPU implementation of k-means would be crazy fast and completely parallel. 73 | 74 | Future work (in no particular order) 75 | 76 | * GPU implementation 77 | * better PSNR by interpolating U, V - what's currently there is a nearest filter, which is horrible 78 | * expose number of passes for minor improvement in quality. At around 8 passes PSNR goes down 0.1db 79 | 80 | [Kodak dataset numbers](https://docs.google.com/spreadsheets/d/e/2PACX-1vROIuXdb9BQB0Gem7Pn0q9Y4heimPg6y8xvhhnJ1Cgaqr1qaJ4LmQsBXUk4pBaG7HcME4SPS2JNNUb2/pubhtml?gid=1381620930&single=true) 81 | 82 | YUV_16_1_1 PSNR=-32.3db PSNR(YUV)=-38.2db 83 | YUV_4_1_1 PSNR=-33.2db PSNR(YUV)=-39.3db 84 | YUV_2_1_1 PSNR=-33.6db PSNR(YUV)=-39.7db 85 | RGB_1_1_1 PSNR=-34.8db PSNR(YUV)=-40.4db 86 | 87 | top left - original, top right - 1:1:1, mid left 2:1:1, mid right 4:1:1, bottom left - 16:1:1, bottom right 1:0:0 88 | ![combine_images](https://user-images.githubusercontent.com/272689/144731527-4fc04c56-bed0-4f91-828e-4960c294b430.png) 89 | 90 | -------------------------------------------------------------------------------- /dirent/LICENSE: -------------------------------------------------------------------------------- 1 | The MIT License (MIT) 2 | 3 | Copyright (c) 1998-2019 Toni Ronkko 4 | 5 | Permission is hereby granted, free of charge, to any person obtaining a copy 6 | of this software and associated documentation files (the "Software"), to deal 7 | in the Software without restriction, including without limitation the rights 8 | to use, copy, modify, merge, publish, distribute, sublicense, and/or sell 9 | copies of the Software, and to permit persons to whom the Software is 10 | furnished to do so, subject to the following conditions: 11 | 12 | The above copyright notice and this permission notice shall be included in all 13 | copies or substantial portions of the Software. 14 | 15 | THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 16 | IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 17 | FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE 18 | AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 19 | LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, 20 | OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE 21 | SOFTWARE. 22 | -------------------------------------------------------------------------------- /dirent/dirent.h: -------------------------------------------------------------------------------- 1 | /* 2 | * Dirent interface for Microsoft Visual Studio 3 | * 4 | * Copyright (C) 1998-2019 Toni Ronkko 5 | * This file is part of dirent. Dirent may be freely distributed 6 | * under the MIT license. For all details and documentation, see 7 | * https://github.com/tronkko/dirent 8 | */ 9 | #ifndef DIRENT_H 10 | #define DIRENT_H 11 | 12 | /* Hide warnings about unreferenced local functions */ 13 | #if defined(__clang__) 14 | # pragma clang diagnostic ignored "-Wunused-function" 15 | #elif defined(_MSC_VER) 16 | # pragma warning(disable:4505) 17 | #elif defined(__GNUC__) 18 | # pragma GCC diagnostic ignored "-Wunused-function" 19 | #endif 20 | 21 | /* 22 | * Include windows.h without Windows Sockets 1.1 to prevent conflicts with 23 | * Windows Sockets 2.0. 24 | */ 25 | #ifndef WIN32_LEAN_AND_MEAN 26 | # define WIN32_LEAN_AND_MEAN 27 | #endif 28 | #include 29 | 30 | #include 31 | #include 32 | #include 33 | #include 34 | #include 35 | #include 36 | #include 37 | #include 38 | #include 39 | #include 40 | 41 | /* Indicates that d_type field is available in dirent structure */ 42 | #define _DIRENT_HAVE_D_TYPE 43 | 44 | /* Indicates that d_namlen field is available in dirent structure */ 45 | #define _DIRENT_HAVE_D_NAMLEN 46 | 47 | /* Entries missing from MSVC 6.0 */ 48 | #if !defined(FILE_ATTRIBUTE_DEVICE) 49 | # define FILE_ATTRIBUTE_DEVICE 0x40 50 | #endif 51 | 52 | /* File type and permission flags for stat(), general mask */ 53 | #if !defined(S_IFMT) 54 | # define S_IFMT _S_IFMT 55 | #endif 56 | 57 | /* Directory bit */ 58 | #if !defined(S_IFDIR) 59 | # define S_IFDIR _S_IFDIR 60 | #endif 61 | 62 | /* Character device bit */ 63 | #if !defined(S_IFCHR) 64 | # define S_IFCHR _S_IFCHR 65 | #endif 66 | 67 | /* Pipe bit */ 68 | #if !defined(S_IFFIFO) 69 | # define S_IFFIFO _S_IFFIFO 70 | #endif 71 | 72 | /* Regular file bit */ 73 | #if !defined(S_IFREG) 74 | # define S_IFREG _S_IFREG 75 | #endif 76 | 77 | /* Read permission */ 78 | #if !defined(S_IREAD) 79 | # define S_IREAD _S_IREAD 80 | #endif 81 | 82 | /* Write permission */ 83 | #if !defined(S_IWRITE) 84 | # define S_IWRITE _S_IWRITE 85 | #endif 86 | 87 | /* Execute permission */ 88 | #if !defined(S_IEXEC) 89 | # define S_IEXEC _S_IEXEC 90 | #endif 91 | 92 | /* Pipe */ 93 | #if !defined(S_IFIFO) 94 | # define S_IFIFO _S_IFIFO 95 | #endif 96 | 97 | /* Block device */ 98 | #if !defined(S_IFBLK) 99 | # define S_IFBLK 0 100 | #endif 101 | 102 | /* Link */ 103 | #if !defined(S_IFLNK) 104 | # define S_IFLNK 0 105 | #endif 106 | 107 | /* Socket */ 108 | #if !defined(S_IFSOCK) 109 | # define S_IFSOCK 0 110 | #endif 111 | 112 | /* Read user permission */ 113 | #if !defined(S_IRUSR) 114 | # define S_IRUSR S_IREAD 115 | #endif 116 | 117 | /* Write user permission */ 118 | #if !defined(S_IWUSR) 119 | # define S_IWUSR S_IWRITE 120 | #endif 121 | 122 | /* Execute user permission */ 123 | #if !defined(S_IXUSR) 124 | # define S_IXUSR 0 125 | #endif 126 | 127 | /* Read group permission */ 128 | #if !defined(S_IRGRP) 129 | # define S_IRGRP 0 130 | #endif 131 | 132 | /* Write group permission */ 133 | #if !defined(S_IWGRP) 134 | # define S_IWGRP 0 135 | #endif 136 | 137 | /* Execute group permission */ 138 | #if !defined(S_IXGRP) 139 | # define S_IXGRP 0 140 | #endif 141 | 142 | /* Read others permission */ 143 | #if !defined(S_IROTH) 144 | # define S_IROTH 0 145 | #endif 146 | 147 | /* Write others permission */ 148 | #if !defined(S_IWOTH) 149 | # define S_IWOTH 0 150 | #endif 151 | 152 | /* Execute others permission */ 153 | #if !defined(S_IXOTH) 154 | # define S_IXOTH 0 155 | #endif 156 | 157 | /* Maximum length of file name */ 158 | #if !defined(PATH_MAX) 159 | # define PATH_MAX MAX_PATH 160 | #endif 161 | #if !defined(FILENAME_MAX) 162 | # define FILENAME_MAX MAX_PATH 163 | #endif 164 | #if !defined(NAME_MAX) 165 | # define NAME_MAX FILENAME_MAX 166 | #endif 167 | 168 | /* File type flags for d_type */ 169 | #define DT_UNKNOWN 0 170 | #define DT_REG S_IFREG 171 | #define DT_DIR S_IFDIR 172 | #define DT_FIFO S_IFIFO 173 | #define DT_SOCK S_IFSOCK 174 | #define DT_CHR S_IFCHR 175 | #define DT_BLK S_IFBLK 176 | #define DT_LNK S_IFLNK 177 | 178 | /* Macros for converting between st_mode and d_type */ 179 | #define IFTODT(mode) ((mode) & S_IFMT) 180 | #define DTTOIF(type) (type) 181 | 182 | /* 183 | * File type macros. Note that block devices, sockets and links cannot be 184 | * distinguished on Windows and the macros S_ISBLK, S_ISSOCK and S_ISLNK are 185 | * only defined for compatibility. These macros should always return false 186 | * on Windows. 187 | */ 188 | #if !defined(S_ISFIFO) 189 | # define S_ISFIFO(mode) (((mode) & S_IFMT) == S_IFIFO) 190 | #endif 191 | #if !defined(S_ISDIR) 192 | # define S_ISDIR(mode) (((mode) & S_IFMT) == S_IFDIR) 193 | #endif 194 | #if !defined(S_ISREG) 195 | # define S_ISREG(mode) (((mode) & S_IFMT) == S_IFREG) 196 | #endif 197 | #if !defined(S_ISLNK) 198 | # define S_ISLNK(mode) (((mode) & S_IFMT) == S_IFLNK) 199 | #endif 200 | #if !defined(S_ISSOCK) 201 | # define S_ISSOCK(mode) (((mode) & S_IFMT) == S_IFSOCK) 202 | #endif 203 | #if !defined(S_ISCHR) 204 | # define S_ISCHR(mode) (((mode) & S_IFMT) == S_IFCHR) 205 | #endif 206 | #if !defined(S_ISBLK) 207 | # define S_ISBLK(mode) (((mode) & S_IFMT) == S_IFBLK) 208 | #endif 209 | 210 | /* Return the exact length of the file name without zero terminator */ 211 | #define _D_EXACT_NAMLEN(p) ((p)->d_namlen) 212 | 213 | /* Return the maximum size of a file name */ 214 | #define _D_ALLOC_NAMLEN(p) ((PATH_MAX)+1) 215 | 216 | 217 | #ifdef __cplusplus 218 | extern "C" { 219 | #endif 220 | 221 | 222 | /* Wide-character version */ 223 | struct _wdirent { 224 | /* Always zero */ 225 | long d_ino; 226 | 227 | /* File position within stream */ 228 | long d_off; 229 | 230 | /* Structure size */ 231 | unsigned short d_reclen; 232 | 233 | /* Length of name without \0 */ 234 | size_t d_namlen; 235 | 236 | /* File type */ 237 | int d_type; 238 | 239 | /* File name */ 240 | wchar_t d_name[PATH_MAX+1]; 241 | }; 242 | typedef struct _wdirent _wdirent; 243 | 244 | struct _WDIR { 245 | /* Current directory entry */ 246 | struct _wdirent ent; 247 | 248 | /* Private file data */ 249 | WIN32_FIND_DATAW data; 250 | 251 | /* True if data is valid */ 252 | int cached; 253 | 254 | /* Win32 search handle */ 255 | HANDLE handle; 256 | 257 | /* Initial directory name */ 258 | wchar_t *patt; 259 | }; 260 | typedef struct _WDIR _WDIR; 261 | 262 | /* Multi-byte character version */ 263 | struct dirent { 264 | /* Always zero */ 265 | long d_ino; 266 | 267 | /* File position within stream */ 268 | long d_off; 269 | 270 | /* Structure size */ 271 | unsigned short d_reclen; 272 | 273 | /* Length of name without \0 */ 274 | size_t d_namlen; 275 | 276 | /* File type */ 277 | int d_type; 278 | 279 | /* File name */ 280 | char d_name[PATH_MAX+1]; 281 | }; 282 | typedef struct dirent dirent; 283 | 284 | struct DIR { 285 | struct dirent ent; 286 | struct _WDIR *wdirp; 287 | }; 288 | typedef struct DIR DIR; 289 | 290 | 291 | /* Dirent functions */ 292 | static DIR *opendir(const char *dirname); 293 | static _WDIR *_wopendir(const wchar_t *dirname); 294 | 295 | static struct dirent *readdir(DIR *dirp); 296 | static struct _wdirent *_wreaddir(_WDIR *dirp); 297 | 298 | static int readdir_r( 299 | DIR *dirp, struct dirent *entry, struct dirent **result); 300 | static int _wreaddir_r( 301 | _WDIR *dirp, struct _wdirent *entry, struct _wdirent **result); 302 | 303 | static int closedir(DIR *dirp); 304 | static int _wclosedir(_WDIR *dirp); 305 | 306 | static void rewinddir(DIR* dirp); 307 | static void _wrewinddir(_WDIR* dirp); 308 | 309 | static int scandir(const char *dirname, struct dirent ***namelist, 310 | int (*filter)(const struct dirent*), 311 | int (*compare)(const struct dirent**, const struct dirent**)); 312 | 313 | static int alphasort(const struct dirent **a, const struct dirent **b); 314 | 315 | static int versionsort(const struct dirent **a, const struct dirent **b); 316 | 317 | static int strverscmp(const char *a, const char *b); 318 | 319 | /* For compatibility with Symbian */ 320 | #define wdirent _wdirent 321 | #define WDIR _WDIR 322 | #define wopendir _wopendir 323 | #define wreaddir _wreaddir 324 | #define wclosedir _wclosedir 325 | #define wrewinddir _wrewinddir 326 | 327 | /* Compatibility with older Microsoft compilers and non-Microsoft compilers */ 328 | #if !defined(_MSC_VER) || _MSC_VER < 1400 329 | # define wcstombs_s dirent_wcstombs_s 330 | # define mbstowcs_s dirent_mbstowcs_s 331 | #endif 332 | 333 | /* Optimize dirent_set_errno() away on modern Microsoft compilers */ 334 | #if defined(_MSC_VER) && _MSC_VER >= 1400 335 | # define dirent_set_errno _set_errno 336 | #endif 337 | 338 | 339 | /* Internal utility functions */ 340 | static WIN32_FIND_DATAW *dirent_first(_WDIR *dirp); 341 | static WIN32_FIND_DATAW *dirent_next(_WDIR *dirp); 342 | 343 | #if !defined(_MSC_VER) || _MSC_VER < 1400 344 | static int dirent_mbstowcs_s( 345 | size_t *pReturnValue, wchar_t *wcstr, size_t sizeInWords, 346 | const char *mbstr, size_t count); 347 | #endif 348 | 349 | #if !defined(_MSC_VER) || _MSC_VER < 1400 350 | static int dirent_wcstombs_s( 351 | size_t *pReturnValue, char *mbstr, size_t sizeInBytes, 352 | const wchar_t *wcstr, size_t count); 353 | #endif 354 | 355 | #if !defined(_MSC_VER) || _MSC_VER < 1400 356 | static void dirent_set_errno(int error); 357 | #endif 358 | 359 | 360 | /* 361 | * Open directory stream DIRNAME for read and return a pointer to the 362 | * internal working area that is used to retrieve individual directory 363 | * entries. 364 | */ 365 | static _WDIR *_wopendir(const wchar_t *dirname) 366 | { 367 | wchar_t *p; 368 | 369 | /* Must have directory name */ 370 | if (dirname == NULL || dirname[0] == '\0') { 371 | dirent_set_errno(ENOENT); 372 | return NULL; 373 | } 374 | 375 | /* Allocate new _WDIR structure */ 376 | _WDIR *dirp = (_WDIR*) malloc(sizeof(struct _WDIR)); 377 | if (!dirp) 378 | return NULL; 379 | 380 | /* Reset _WDIR structure */ 381 | dirp->handle = INVALID_HANDLE_VALUE; 382 | dirp->patt = NULL; 383 | dirp->cached = 0; 384 | 385 | /* 386 | * Compute the length of full path plus zero terminator 387 | * 388 | * Note that on WinRT there's no way to convert relative paths 389 | * into absolute paths, so just assume it is an absolute path. 390 | */ 391 | #if WINAPI_FAMILY_PARTITION(WINAPI_PARTITION_DESKTOP) 392 | /* Desktop */ 393 | DWORD n = GetFullPathNameW(dirname, 0, NULL, NULL); 394 | #else 395 | /* WinRT */ 396 | size_t n = wcslen(dirname); 397 | #endif 398 | 399 | /* Allocate room for absolute directory name and search pattern */ 400 | dirp->patt = (wchar_t*) malloc(sizeof(wchar_t) * n + 16); 401 | if (dirp->patt == NULL) 402 | goto exit_closedir; 403 | 404 | /* 405 | * Convert relative directory name to an absolute one. This 406 | * allows rewinddir() to function correctly even when current 407 | * working directory is changed between opendir() and rewinddir(). 408 | * 409 | * Note that on WinRT there's no way to convert relative paths 410 | * into absolute paths, so just assume it is an absolute path. 411 | */ 412 | #if WINAPI_FAMILY_PARTITION(WINAPI_PARTITION_DESKTOP) 413 | /* Desktop */ 414 | n = GetFullPathNameW(dirname, n, dirp->patt, NULL); 415 | if (n <= 0) 416 | goto exit_closedir; 417 | #else 418 | /* WinRT */ 419 | wcsncpy_s(dirp->patt, n+1, dirname, n); 420 | #endif 421 | 422 | /* Append search pattern \* to the directory name */ 423 | p = dirp->patt + n; 424 | switch (p[-1]) { 425 | case '\\': 426 | case '/': 427 | case ':': 428 | /* Directory ends in path separator, e.g. c:\temp\ */ 429 | /*NOP*/; 430 | break; 431 | 432 | default: 433 | /* Directory name doesn't end in path separator */ 434 | *p++ = '\\'; 435 | } 436 | *p++ = '*'; 437 | *p = '\0'; 438 | 439 | /* Open directory stream and retrieve the first entry */ 440 | if (!dirent_first(dirp)) 441 | goto exit_closedir; 442 | 443 | /* Success */ 444 | return dirp; 445 | 446 | /* Failure */ 447 | exit_closedir: 448 | _wclosedir(dirp); 449 | return NULL; 450 | } 451 | 452 | /* 453 | * Read next directory entry. 454 | * 455 | * Returns pointer to static directory entry which may be overwritten by 456 | * subsequent calls to _wreaddir(). 457 | */ 458 | static struct _wdirent *_wreaddir(_WDIR *dirp) 459 | { 460 | /* 461 | * Read directory entry to buffer. We can safely ignore the return 462 | * value as entry will be set to NULL in case of error. 463 | */ 464 | struct _wdirent *entry; 465 | (void) _wreaddir_r(dirp, &dirp->ent, &entry); 466 | 467 | /* Return pointer to statically allocated directory entry */ 468 | return entry; 469 | } 470 | 471 | /* 472 | * Read next directory entry. 473 | * 474 | * Returns zero on success. If end of directory stream is reached, then sets 475 | * result to NULL and returns zero. 476 | */ 477 | static int _wreaddir_r( 478 | _WDIR *dirp, struct _wdirent *entry, struct _wdirent **result) 479 | { 480 | /* Read next directory entry */ 481 | WIN32_FIND_DATAW *datap = dirent_next(dirp); 482 | if (!datap) { 483 | /* Return NULL to indicate end of directory */ 484 | *result = NULL; 485 | return /*OK*/0; 486 | } 487 | 488 | /* 489 | * Copy file name as wide-character string. If the file name is too 490 | * long to fit in to the destination buffer, then truncate file name 491 | * to PATH_MAX characters and zero-terminate the buffer. 492 | */ 493 | size_t n = 0; 494 | while (n < PATH_MAX && datap->cFileName[n] != 0) { 495 | entry->d_name[n] = datap->cFileName[n]; 496 | n++; 497 | } 498 | entry->d_name[n] = 0; 499 | 500 | /* Length of file name excluding zero terminator */ 501 | entry->d_namlen = n; 502 | 503 | /* File type */ 504 | DWORD attr = datap->dwFileAttributes; 505 | if ((attr & FILE_ATTRIBUTE_DEVICE) != 0) 506 | entry->d_type = DT_CHR; 507 | else if ((attr & FILE_ATTRIBUTE_DIRECTORY) != 0) 508 | entry->d_type = DT_DIR; 509 | else 510 | entry->d_type = DT_REG; 511 | 512 | /* Reset dummy fields */ 513 | entry->d_ino = 0; 514 | entry->d_off = 0; 515 | entry->d_reclen = sizeof(struct _wdirent); 516 | 517 | /* Set result address */ 518 | *result = entry; 519 | return /*OK*/0; 520 | } 521 | 522 | /* 523 | * Close directory stream opened by opendir() function. This invalidates the 524 | * DIR structure as well as any directory entry read previously by 525 | * _wreaddir(). 526 | */ 527 | static int _wclosedir(_WDIR *dirp) 528 | { 529 | if (!dirp) { 530 | dirent_set_errno(EBADF); 531 | return /*failure*/-1; 532 | } 533 | 534 | /* Release search handle */ 535 | if (dirp->handle != INVALID_HANDLE_VALUE) 536 | FindClose(dirp->handle); 537 | 538 | /* Release search pattern */ 539 | free(dirp->patt); 540 | 541 | /* Release directory structure */ 542 | free(dirp); 543 | return /*success*/0; 544 | } 545 | 546 | /* 547 | * Rewind directory stream such that _wreaddir() returns the very first 548 | * file name again. 549 | */ 550 | static void _wrewinddir(_WDIR* dirp) 551 | { 552 | if (!dirp) 553 | return; 554 | 555 | /* Release existing search handle */ 556 | if (dirp->handle != INVALID_HANDLE_VALUE) 557 | FindClose(dirp->handle); 558 | 559 | /* Open new search handle */ 560 | dirent_first(dirp); 561 | } 562 | 563 | /* Get first directory entry */ 564 | static WIN32_FIND_DATAW *dirent_first(_WDIR *dirp) 565 | { 566 | if (!dirp) 567 | return NULL; 568 | 569 | /* Open directory and retrieve the first entry */ 570 | dirp->handle = FindFirstFileExW( 571 | dirp->patt, FindExInfoStandard, &dirp->data, 572 | FindExSearchNameMatch, NULL, 0); 573 | if (dirp->handle == INVALID_HANDLE_VALUE) 574 | goto error; 575 | 576 | /* A directory entry is now waiting in memory */ 577 | dirp->cached = 1; 578 | return &dirp->data; 579 | 580 | error: 581 | /* Failed to open directory: no directory entry in memory */ 582 | dirp->cached = 0; 583 | 584 | /* Set error code */ 585 | DWORD errorcode = GetLastError(); 586 | switch (errorcode) { 587 | case ERROR_ACCESS_DENIED: 588 | /* No read access to directory */ 589 | dirent_set_errno(EACCES); 590 | break; 591 | 592 | case ERROR_DIRECTORY: 593 | /* Directory name is invalid */ 594 | dirent_set_errno(ENOTDIR); 595 | break; 596 | 597 | case ERROR_PATH_NOT_FOUND: 598 | default: 599 | /* Cannot find the file */ 600 | dirent_set_errno(ENOENT); 601 | } 602 | return NULL; 603 | } 604 | 605 | /* Get next directory entry */ 606 | static WIN32_FIND_DATAW *dirent_next(_WDIR *dirp) 607 | { 608 | /* Is the next directory entry already in cache? */ 609 | if (dirp->cached) { 610 | /* Yes, a valid directory entry found in memory */ 611 | dirp->cached = 0; 612 | return &dirp->data; 613 | } 614 | 615 | /* No directory entry in cache */ 616 | if (dirp->handle == INVALID_HANDLE_VALUE) 617 | return NULL; 618 | 619 | /* Read the next directory entry from stream */ 620 | if (FindNextFileW(dirp->handle, &dirp->data) == FALSE) 621 | goto exit_close; 622 | 623 | /* Success */ 624 | return &dirp->data; 625 | 626 | /* Failure */ 627 | exit_close: 628 | FindClose(dirp->handle); 629 | dirp->handle = INVALID_HANDLE_VALUE; 630 | return NULL; 631 | } 632 | 633 | /* Open directory stream using plain old C-string */ 634 | static DIR *opendir(const char *dirname) 635 | { 636 | /* Must have directory name */ 637 | if (dirname == NULL || dirname[0] == '\0') { 638 | dirent_set_errno(ENOENT); 639 | return NULL; 640 | } 641 | 642 | /* Allocate memory for DIR structure */ 643 | struct DIR *dirp = (DIR*) malloc(sizeof(struct DIR)); 644 | if (!dirp) 645 | return NULL; 646 | 647 | /* Convert directory name to wide-character string */ 648 | wchar_t wname[PATH_MAX + 1]; 649 | size_t n; 650 | int error = mbstowcs_s(&n, wname, PATH_MAX + 1, dirname, PATH_MAX+1); 651 | if (error) 652 | goto exit_failure; 653 | 654 | /* Open directory stream using wide-character name */ 655 | dirp->wdirp = _wopendir(wname); 656 | if (!dirp->wdirp) 657 | goto exit_failure; 658 | 659 | /* Success */ 660 | return dirp; 661 | 662 | /* Failure */ 663 | exit_failure: 664 | free(dirp); 665 | return NULL; 666 | } 667 | 668 | /* Read next directory entry */ 669 | static struct dirent *readdir(DIR *dirp) 670 | { 671 | /* 672 | * Read directory entry to buffer. We can safely ignore the return 673 | * value as entry will be set to NULL in case of error. 674 | */ 675 | struct dirent *entry; 676 | (void) readdir_r(dirp, &dirp->ent, &entry); 677 | 678 | /* Return pointer to statically allocated directory entry */ 679 | return entry; 680 | } 681 | 682 | /* 683 | * Read next directory entry into called-allocated buffer. 684 | * 685 | * Returns zero on success. If the end of directory stream is reached, then 686 | * sets result to NULL and returns zero. 687 | */ 688 | static int readdir_r( 689 | DIR *dirp, struct dirent *entry, struct dirent **result) 690 | { 691 | /* Read next directory entry */ 692 | WIN32_FIND_DATAW *datap = dirent_next(dirp->wdirp); 693 | if (!datap) { 694 | /* No more directory entries */ 695 | *result = NULL; 696 | return /*OK*/0; 697 | } 698 | 699 | /* Attempt to convert file name to multi-byte string */ 700 | size_t n; 701 | int error = wcstombs_s( 702 | &n, entry->d_name, PATH_MAX + 1, 703 | datap->cFileName, PATH_MAX + 1); 704 | 705 | /* 706 | * If the file name cannot be represented by a multi-byte string, then 707 | * attempt to use old 8+3 file name. This allows the program to 708 | * access files although file names may seem unfamiliar to the user. 709 | * 710 | * Be ware that the code below cannot come up with a short file name 711 | * unless the file system provides one. At least VirtualBox shared 712 | * folders fail to do this. 713 | */ 714 | if (error && datap->cAlternateFileName[0] != '\0') { 715 | error = wcstombs_s( 716 | &n, entry->d_name, PATH_MAX + 1, 717 | datap->cAlternateFileName, PATH_MAX + 1); 718 | } 719 | 720 | if (!error) { 721 | /* Length of file name excluding zero terminator */ 722 | entry->d_namlen = n - 1; 723 | 724 | /* File attributes */ 725 | DWORD attr = datap->dwFileAttributes; 726 | if ((attr & FILE_ATTRIBUTE_DEVICE) != 0) 727 | entry->d_type = DT_CHR; 728 | else if ((attr & FILE_ATTRIBUTE_DIRECTORY) != 0) 729 | entry->d_type = DT_DIR; 730 | else 731 | entry->d_type = DT_REG; 732 | 733 | /* Reset dummy fields */ 734 | entry->d_ino = 0; 735 | entry->d_off = 0; 736 | entry->d_reclen = sizeof(struct dirent); 737 | } else { 738 | /* 739 | * Cannot convert file name to multi-byte string so construct 740 | * an erroneous directory entry and return that. Note that 741 | * we cannot return NULL as that would stop the processing 742 | * of directory entries completely. 743 | */ 744 | entry->d_name[0] = '?'; 745 | entry->d_name[1] = '\0'; 746 | entry->d_namlen = 1; 747 | entry->d_type = DT_UNKNOWN; 748 | entry->d_ino = 0; 749 | entry->d_off = -1; 750 | entry->d_reclen = 0; 751 | } 752 | 753 | /* Return pointer to directory entry */ 754 | *result = entry; 755 | return /*OK*/0; 756 | } 757 | 758 | /* Close directory stream */ 759 | static int closedir(DIR *dirp) 760 | { 761 | int ok; 762 | 763 | if (!dirp) 764 | goto exit_failure; 765 | 766 | /* Close wide-character directory stream */ 767 | ok = _wclosedir(dirp->wdirp); 768 | dirp->wdirp = NULL; 769 | 770 | /* Release multi-byte character version */ 771 | free(dirp); 772 | return ok; 773 | 774 | exit_failure: 775 | /* Invalid directory stream */ 776 | dirent_set_errno(EBADF); 777 | return /*failure*/-1; 778 | } 779 | 780 | /* Rewind directory stream to beginning */ 781 | static void rewinddir(DIR* dirp) 782 | { 783 | if (!dirp) 784 | return; 785 | 786 | /* Rewind wide-character string directory stream */ 787 | _wrewinddir(dirp->wdirp); 788 | } 789 | 790 | /* Scan directory for entries */ 791 | static int scandir( 792 | const char *dirname, struct dirent ***namelist, 793 | int (*filter)(const struct dirent*), 794 | int (*compare)(const struct dirent**, const struct dirent**)) 795 | { 796 | int result; 797 | 798 | /* Open directory stream */ 799 | DIR *dir = opendir(dirname); 800 | if (!dir) { 801 | /* Cannot open directory */ 802 | return /*Error*/ -1; 803 | } 804 | 805 | /* Read directory entries to memory */ 806 | struct dirent *tmp = NULL; 807 | struct dirent **files = NULL; 808 | size_t size = 0; 809 | size_t allocated = 0; 810 | while (1) { 811 | /* Allocate room for a temporary directory entry */ 812 | if (!tmp) { 813 | tmp = (struct dirent*) malloc(sizeof(struct dirent)); 814 | if (!tmp) 815 | goto exit_failure; 816 | } 817 | 818 | /* Read directory entry to temporary area */ 819 | struct dirent *entry; 820 | if (readdir_r(dir, tmp, &entry) != /*OK*/0) 821 | goto exit_failure; 822 | 823 | /* Stop if we already read the last directory entry */ 824 | if (entry == NULL) 825 | goto exit_success; 826 | 827 | /* Determine whether to include the entry in results */ 828 | if (filter && !filter(tmp)) 829 | continue; 830 | 831 | /* Enlarge pointer table to make room for another pointer */ 832 | if (size >= allocated) { 833 | /* Compute number of entries in the new table */ 834 | size_t num_entries = size * 2 + 16; 835 | 836 | /* Allocate new pointer table or enlarge existing */ 837 | void *p = realloc(files, sizeof(void*) * num_entries); 838 | if (!p) 839 | goto exit_failure; 840 | 841 | /* Got the memory */ 842 | files = (dirent**) p; 843 | allocated = num_entries; 844 | } 845 | 846 | /* Store the temporary entry to ptr table */ 847 | files[size++] = tmp; 848 | tmp = NULL; 849 | } 850 | 851 | exit_failure: 852 | /* Release allocated file entries */ 853 | for (size_t i = 0; i < size; i++) { 854 | free(files[i]); 855 | } 856 | 857 | /* Release the pointer table */ 858 | free(files); 859 | files = NULL; 860 | 861 | /* Exit with error code */ 862 | result = /*error*/ -1; 863 | goto exit_status; 864 | 865 | exit_success: 866 | /* Sort directory entries */ 867 | qsort(files, size, sizeof(void*), 868 | (int (*) (const void*, const void*)) compare); 869 | 870 | /* Pass pointer table to caller */ 871 | if (namelist) 872 | *namelist = files; 873 | 874 | /* Return the number of directory entries read */ 875 | result = (int) size; 876 | 877 | exit_status: 878 | /* Release temporary directory entry, if we had one */ 879 | free(tmp); 880 | 881 | /* Close directory stream */ 882 | closedir(dir); 883 | return result; 884 | } 885 | 886 | /* Alphabetical sorting */ 887 | static int alphasort(const struct dirent **a, const struct dirent **b) 888 | { 889 | return strcoll((*a)->d_name, (*b)->d_name); 890 | } 891 | 892 | /* Sort versions */ 893 | static int versionsort(const struct dirent **a, const struct dirent **b) 894 | { 895 | return strverscmp((*a)->d_name, (*b)->d_name); 896 | } 897 | 898 | /* Compare strings */ 899 | static int strverscmp(const char *a, const char *b) 900 | { 901 | size_t i = 0; 902 | size_t j; 903 | 904 | /* Find first difference */ 905 | while (a[i] == b[i]) { 906 | if (a[i] == '\0') { 907 | /* No difference */ 908 | return 0; 909 | } 910 | ++i; 911 | } 912 | 913 | /* Count backwards and find the leftmost digit */ 914 | j = i; 915 | while (j > 0 && isdigit(a[j-1])) { 916 | --j; 917 | } 918 | 919 | /* Determine mode of comparison */ 920 | if (a[j] == '0' || b[j] == '0') { 921 | /* Find the next non-zero digit */ 922 | while (a[j] == '0' && a[j] == b[j]) { 923 | j++; 924 | } 925 | 926 | /* String with more digits is smaller, e.g 002 < 01 */ 927 | if (isdigit(a[j])) { 928 | if (!isdigit(b[j])) { 929 | return -1; 930 | } 931 | } else if (isdigit(b[j])) { 932 | return 1; 933 | } 934 | } else if (isdigit(a[j]) && isdigit(b[j])) { 935 | /* Numeric comparison */ 936 | size_t k1 = j; 937 | size_t k2 = j; 938 | 939 | /* Compute number of digits in each string */ 940 | while (isdigit(a[k1])) { 941 | k1++; 942 | } 943 | while (isdigit(b[k2])) { 944 | k2++; 945 | } 946 | 947 | /* Number with more digits is bigger, e.g 999 < 1000 */ 948 | if (k1 < k2) 949 | return -1; 950 | else if (k1 > k2) 951 | return 1; 952 | } 953 | 954 | /* Alphabetical comparison */ 955 | return (int) ((unsigned char) a[i]) - ((unsigned char) b[i]); 956 | } 957 | 958 | /* Convert multi-byte string to wide character string */ 959 | #if !defined(_MSC_VER) || _MSC_VER < 1400 960 | static int dirent_mbstowcs_s( 961 | size_t *pReturnValue, wchar_t *wcstr, 962 | size_t sizeInWords, const char *mbstr, size_t count) 963 | { 964 | /* Older Visual Studio or non-Microsoft compiler */ 965 | size_t n = mbstowcs(wcstr, mbstr, sizeInWords); 966 | if (wcstr && n >= count) 967 | return /*error*/ 1; 968 | 969 | /* Zero-terminate output buffer */ 970 | if (wcstr && sizeInWords) { 971 | if (n >= sizeInWords) 972 | n = sizeInWords - 1; 973 | wcstr[n] = 0; 974 | } 975 | 976 | /* Length of multi-byte string with zero terminator */ 977 | if (pReturnValue) { 978 | *pReturnValue = n + 1; 979 | } 980 | 981 | /* Success */ 982 | return 0; 983 | } 984 | #endif 985 | 986 | /* Convert wide-character string to multi-byte string */ 987 | #if !defined(_MSC_VER) || _MSC_VER < 1400 988 | static int dirent_wcstombs_s( 989 | size_t *pReturnValue, char *mbstr, 990 | size_t sizeInBytes, const wchar_t *wcstr, size_t count) 991 | { 992 | /* Older Visual Studio or non-Microsoft compiler */ 993 | size_t n = wcstombs(mbstr, wcstr, sizeInBytes); 994 | if (mbstr && n >= count) 995 | return /*error*/1; 996 | 997 | /* Zero-terminate output buffer */ 998 | if (mbstr && sizeInBytes) { 999 | if (n >= sizeInBytes) { 1000 | n = sizeInBytes - 1; 1001 | } 1002 | mbstr[n] = '\0'; 1003 | } 1004 | 1005 | /* Length of resulting multi-bytes string WITH zero-terminator */ 1006 | if (pReturnValue) { 1007 | *pReturnValue = n + 1; 1008 | } 1009 | 1010 | /* Success */ 1011 | return 0; 1012 | } 1013 | #endif 1014 | 1015 | /* Set errno variable */ 1016 | #if !defined(_MSC_VER) || _MSC_VER < 1400 1017 | static void dirent_set_errno(int error) 1018 | { 1019 | /* Non-Microsoft compiler or older Microsoft compiler */ 1020 | errno = error; 1021 | } 1022 | #endif 1023 | 1024 | #ifdef __cplusplus 1025 | } 1026 | #endif 1027 | #endif /*DIRENT_H*/ 1028 | -------------------------------------------------------------------------------- /examples/batch.bat: -------------------------------------------------------------------------------- 1 | ..\bin\Release\noi.exe -b .\ .\test 16_1_1 2 | -------------------------------------------------------------------------------- /examples/lenna.bat: -------------------------------------------------------------------------------- 1 | @echo off 2 | 3 | ..\bin\Release\noi.exe -pc lenna.png lenna_16_1_1.noi 16_1_1 4 | ..\bin\Release\noi.exe -pd lenna_16_1_1.noi lenna_test_16_1_1.png 5 | lenna_test_16_1_1.png 6 | 7 | ..\bin\Release\noi.exe -pc lenna.png lenna_4_1_1.noi 4_1_1 8 | ..\bin\Release\noi.exe -pd lenna_4_1_1.noi lenna_test_4_1_1.png 9 | lenna_test_4_1_1.png 10 | 11 | ..\bin\Release\noi.exe -pc lenna.png lenna_2_1_1.noi 2_1_1 12 | ..\bin\Release\noi.exe -pd lenna_2_1_1.noi lenna_test_2_1_1.png 13 | lenna_test_2_1_1.png 14 | 15 | ..\bin\Release\noi.exe -pc lenna.png lenna_1_1_1.noi 1_1_1 16 | ..\bin\Release\noi.exe -pd lenna_1_1_1.noi lenna_test_1_1_1.png 17 | lenna_test_1_1_1.png 18 | 19 | ..\bin\Release\noi.exe -pc lenna.png lenna_1_0_0.noi 1_0_0 20 | ..\bin\Release\noi.exe -pd lenna_1_0_0.noi lenna_test_1_0_0.png 21 | lenna_test_1_0_0.png 22 | 23 | -------------------------------------------------------------------------------- /examples/lenna.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/borisbat/NotOkImageFormat/2eefc9f0e8c12486eac17df1cb1a56f9180e107d/examples/lenna.png -------------------------------------------------------------------------------- /examples/lenna_1_0_0.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/borisbat/NotOkImageFormat/2eefc9f0e8c12486eac17df1cb1a56f9180e107d/examples/lenna_1_0_0.png -------------------------------------------------------------------------------- /examples/lenna_osx.sh: -------------------------------------------------------------------------------- 1 | ../bin/noi -pc lenna.png lenna_16_1_1.noi 16_1_1 2 | ../bin/noi -pd lenna_16_1_1.noi lenna_test_16_1_1.png 3 | open lenna_test_16_1_1.png 4 | 5 | ../bin/noi -pc lenna.png lenna_4_1_1.noi 4_1_1 6 | ../bin/noi -pd lenna_4_1_1.noi lenna_test_4_1_1.png 7 | open lenna_test_4_1_1.png 8 | 9 | ../bin/noi -pc lenna.png lenna_2_1_1.noi 2_1_1 10 | ../bin/noi -pd lenna_2_1_1.noi lenna_test_2_1_1.png 11 | open lenna_test_2_1_1.png 12 | 13 | ../bin/noi -pc lenna.png lenna_1_1_1.noi 1_1_1 14 | ../bin/noi -pd lenna_1_1_1.noi lenna_test_1_1_1.png 15 | open lenna_test_1_1_1.png 16 | 17 | ../bin/noi -pc lenna.png lenna_1_0_0.noi 1_0_0 18 | ../bin/noi -pd lenna_1_0_0.noi lenna_1_0_0.png 19 | open lenna_1_0_0.png 20 | 21 | -------------------------------------------------------------------------------- /examples/lenna_test_16_1_1.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/borisbat/NotOkImageFormat/2eefc9f0e8c12486eac17df1cb1a56f9180e107d/examples/lenna_test_16_1_1.png -------------------------------------------------------------------------------- /examples/lenna_test_1_0_0.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/borisbat/NotOkImageFormat/2eefc9f0e8c12486eac17df1cb1a56f9180e107d/examples/lenna_test_1_0_0.png -------------------------------------------------------------------------------- /examples/lenna_test_1_1_1.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/borisbat/NotOkImageFormat/2eefc9f0e8c12486eac17df1cb1a56f9180e107d/examples/lenna_test_1_1_1.png -------------------------------------------------------------------------------- /examples/lenna_test_2_1_1.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/borisbat/NotOkImageFormat/2eefc9f0e8c12486eac17df1cb1a56f9180e107d/examples/lenna_test_2_1_1.png -------------------------------------------------------------------------------- /examples/lenna_test_4_1_1.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/borisbat/NotOkImageFormat/2eefc9f0e8c12486eac17df1cb1a56f9180e107d/examples/lenna_test_4_1_1.png -------------------------------------------------------------------------------- /noi_image.h: -------------------------------------------------------------------------------- 1 | // NotOkImageFormat revision 0 2 | // NOI_THREADS - when defined, compession will use posix (or windows) threads 3 | // NOI_IMAGE_IMPLEMENTATION - when defined, this also becomes implementation 4 | 5 | #ifndef NOI_IMAGE_H 6 | #define NOI_IMAGE_H 7 | 8 | #include 9 | #include 10 | #include 11 | #include 12 | #include 13 | 14 | #define NOI_KMEANS_NPASS 32 // this controls compression quality. 100 is excessive? 8 goes down, 32 is excellent 15 | #define NOI_MAX_THREADS 32 // when NOI_THREADS is defined, this is max number of threads 16 | #define NOI_MAGIC 0xBAD0 17 | 18 | #define NOI_YUV_16_1_1 18 // YUV 16:1:1 19 | #define NOI_YUV_4_1_1 6 // YUV 4:1:1 20 | #define NOI_YUV_2_1_1 4 // YUV 2:1:1 21 | #define NOI_RGB_1_1_1 3 // RGB 1:1:1 22 | #define NOI_RG_1_1_0 2 // RB 1:1:0 23 | #define NOI_Y_1_0_0 1 // Y 1:0:0 greyscale 24 | 25 | #define NOI_K3 1024 26 | #define NOI_K5 512 27 | #define NOI_K7 512 28 | 29 | #ifdef __cplusplus 30 | extern "C" { 31 | #endif 32 | 33 | typedef struct noi_header_s { // NOI file header 34 | uint16_t magic; 35 | uint16_t profile; 36 | uint16_t width; 37 | uint16_t height; 38 | } noi_header_t; 39 | 40 | void * noi_compress ( uint8_t * pixels, int w, int h, int * bytes, int profile ); // profile is NOI_YUV_16_1_1, etc 41 | uint8_t * noi_decompress ( void * bytes, int * W, int * H, uint8_t * pixels ); // if pixels are NULL, they are allocated 42 | void noi_image_size ( void * bytes, int * W, int * H ); 43 | const char * noi_profile_name ( int profile ); 44 | void noi_profile_block_size ( int profile, int * bx, int * by ); 45 | 46 | #ifdef __cplusplus 47 | } 48 | #endif 49 | 50 | #ifdef NOI_IMAGE_IMPLEMENTATION 51 | 52 | void noi_hdt2x2 ( int * a, int * b, int * c, int * d ) { 53 | int A = *a; int B = *b; int C = *c; int D = *d; 54 | int aab = A + B; int sab = A - B; int acd = C + D; int scd = C - D; 55 | *a = aab + acd; *b = sab + scd; *c = aab - acd; *d = sab - scd; 56 | } 57 | 58 | void noi_hdt4x4 ( int * m ) { 59 | noi_hdt2x2(m+ 0,m+ 1,m+ 4,m+ 5); noi_hdt2x2(m+ 2,m+ 3,m+ 6,m+ 7); noi_hdt2x2(m+ 8,m+ 9,m+12,m+13); noi_hdt2x2(m+10,m+11,m+14,m+15); 60 | noi_hdt2x2(m+ 0,m+ 2,m+ 8,m+10); noi_hdt2x2(m+ 1,m+ 3,m+ 9,m+11); noi_hdt2x2(m+ 4,m+ 6,m+12,m+14); noi_hdt2x2(m+ 5,m+ 7,m+13,m+15); 61 | } 62 | 63 | void noi_rgb2yuv ( int r, int g, int b, int * y, int * u, int * v) { 64 | *y = (( 66*r + 129*g + 25*b + 128)>>8) + 16; 65 | *u = ((-38*r - 74*g + 112*b + 128)>>8) + 128; 66 | *v = ((112*r - 94*g - 18*b + 128)>>8) + 128; 67 | } 68 | 69 | void noi_yuv2rgb ( int y, int u, int v, int * r, int * g, int * b ) { 70 | int c = y - 16; int d = u - 128; int e = v - 128; 71 | *r = (298*c + 409*e + 128) >> 8; 72 | *g = (298*c - 100*d - 209*e + 128) >> 8; 73 | *b = (298*c + 516*d + 128) >> 8; 74 | } 75 | 76 | uint8_t noi_saturate ( int t ) { 77 | return (t<0) ? 0 : (t>255 ? 255 : t); 78 | } 79 | 80 | // 0, 1, 2, 3 -> 0 2 1 3 81 | 82 | /* 83 | 0 2 1 3 84 | 8 10 9 11 85 | 4 6 5 7 86 | 12 14 13 15 87 | */ 88 | 89 | #define NOI_B3_0 2 90 | #define NOI_B3_1 8 91 | #define NOI_B3_2 10 92 | 93 | #define NOI_B5_0 1 94 | #define NOI_B5_1 4 95 | #define NOI_B5_2 5 96 | #define NOI_B5_3 6 97 | #define NOI_B5_4 9 98 | 99 | #define NOI_B7_0 3 100 | #define NOI_B7_1 7 101 | #define NOI_B7_2 11 102 | #define NOI_B7_3 12 103 | #define NOI_B7_4 13 104 | #define NOI_B7_5 14 105 | #define NOI_B7_6 15 106 | 107 | int noi_dist3 ( int * a, int * b ) { 108 | int d0 = a[NOI_B3_0] - b[NOI_B3_0]; int d1 = a[NOI_B3_1] - b[NOI_B3_1]; int d2 = a[NOI_B3_2] - b[NOI_B3_2]; 109 | return d0*d0 + d1*d1 + d2*d2; 110 | } 111 | 112 | int noi_dist5 ( int * a, int * b ) { 113 | int d0 = a[NOI_B5_0] - b[NOI_B5_0]; int d1 = a[NOI_B5_1] - b[NOI_B5_1]; int d2 = a[NOI_B5_2] - b[NOI_B5_2]; 114 | int d3 = a[NOI_B5_3] - b[NOI_B5_3]; int d4 = a[NOI_B5_4] - b[NOI_B5_4]; 115 | return d0*d0 + d1*d1 + d2*d2 + d3*d3 + d4*d4; 116 | } 117 | 118 | int noi_dist7 ( int * a, int * b ) { 119 | int d0 = a[NOI_B7_0] - b[NOI_B7_0]; int d1 = a[NOI_B7_1] - b[NOI_B7_1]; int d2 = a[NOI_B7_2] - b[NOI_B7_2]; 120 | int d3 = a[NOI_B7_3] - b[NOI_B7_3]; int d4 = a[NOI_B7_4] - b[NOI_B7_4]; int d5 = a[NOI_B7_5] - b[NOI_B7_5]; 121 | int d6 = a[NOI_B7_6] - b[NOI_B7_6]; 122 | return d0*d0 + d1*d1 + d2*d2 + d3*d3 + d4*d4 + d5*d5 + d6*d6; 123 | } 124 | 125 | typedef struct noi_kmeans_s { 126 | int * index; 127 | int * center; 128 | } noi_kmeans_t; 129 | 130 | void noi_free_means ( noi_kmeans_t * km ) { 131 | free(km->index); 132 | free(km->center); 133 | } 134 | 135 | typedef struct noi_kmeans_thread_s { 136 | int * blocks; 137 | int nblocks; 138 | int * index; 139 | int * center; 140 | int K; 141 | int ( *dist) ( int *, int * ); 142 | } noi_kmeans_thread_t; 143 | 144 | void noi_kmeans_single_thread ( noi_kmeans_thread_t * nkt ) { 145 | int * blocks = nkt->blocks; 146 | int nblocks = nkt->nblocks; 147 | int ( *dist ) ( int *, int * ) = nkt->dist; 148 | int * index = nkt->index; 149 | int * center = nkt->center; 150 | int K = nkt->K; 151 | for ( int o=0; o!=nblocks; ++o ) { 152 | int cur_dist = (*dist) ( blocks + o*16, center ); 153 | int cur_index = 0; 154 | for ( int k=1; k!=K; ++k ) { 155 | int d = (*dist) ( blocks + o*16, center + k*16 ); 156 | if ( d < cur_dist ) { 157 | cur_index = k; 158 | cur_dist = d; 159 | } 160 | } 161 | index[o] = cur_index; 162 | } 163 | } 164 | 165 | #ifdef NOI_THREADS 166 | #ifdef _WIN32 167 | #define WIN32_LEAN_AND_MEAN 168 | #include 169 | DWORD WINAPI noi_kmeans_thread ( void * _nkt ) { 170 | noi_kmeans_thread_t * nkt = (noi_kmeans_thread_t *) _nkt; 171 | noi_kmeans_single_thread(nkt); 172 | return 0; 173 | } 174 | #else 175 | #include 176 | #include 177 | void * noi_kmeans_thread ( void * _nkt ) { 178 | noi_kmeans_thread_t * nkt = (noi_kmeans_thread_t *) _nkt; 179 | noi_kmeans_single_thread(nkt); 180 | return NULL; 181 | } 182 | #endif 183 | #endif 184 | 185 | void noi_kmeans ( noi_kmeans_t * res, int * blocks, int nblocks, int K, int (* dist) (int * o, int * c) ) { 186 | int * index = (int *) malloc(nblocks*sizeof(int)); 187 | int * center = (int *) malloc(K*16*sizeof(int)); 188 | int * num = (int *) malloc(K*sizeof(int)); 189 | int dim_min[16], dim_max[16]; 190 | for ( int i=1; i!=16; ++i ) { 191 | dim_min[i] = INT_MAX; 192 | dim_max[i] = -INT_MAX; 193 | } 194 | for ( int o=0; o!=nblocks; ++o ) { 195 | for ( int i=1; i!=16; ++i ) { 196 | int d = blocks[o*16 + i]; 197 | if ( dim_min[i]>d ) dim_min[i] = d; 198 | if ( dim_max[i]>2)) % ( dim_max[i] - dim_min[i] ) + dim_min[i]; 207 | } else { 208 | center[k*16+i] = dim_min[i]; 209 | } 210 | } 211 | } 212 | memset(index, 0, nblocks*sizeof(int)); 213 | #ifdef NOI_THREADS 214 | #if _WIN32 215 | SYSTEM_INFO sysinfo; 216 | GetSystemInfo(&sysinfo); 217 | int nth = sysinfo.dwNumberOfProcessors - 2; 218 | #else 219 | int nth = sysconf(_SC_NPROCESSORS_ONLN) - 2; 220 | #endif 221 | nth = nth<=0 ? 1 : (nth>NOI_MAX_THREADS ? NOI_MAX_THREADS : nth); 222 | noi_kmeans_thread_t tinfo[NOI_MAX_THREADS]; 223 | int bpt = nblocks / nth; 224 | for ( int t=0; t!=nth; ++t ) { 225 | tinfo[t] = (noi_kmeans_thread_t ) { blocks+16*t*bpt,bpt,index+t*bpt,center,K,dist }; 226 | } 227 | tinfo[nth-1].nblocks = nblocks - (nth-1)*bpt; 228 | #endif 229 | for ( int iter=0; iter!=NOI_KMEANS_NPASS; ++iter ) { 230 | #ifdef NOI_THREADS 231 | #ifdef _WIN32 232 | HANDLE hth[NOI_MAX_THREADS]; 233 | for ( int t=1; t!=nth; ++t ) { 234 | hth[t] = CreateThread(NULL, 0, noi_kmeans_thread, tinfo + t, 0, NULL); 235 | } 236 | noi_kmeans_single_thread(tinfo); 237 | WaitForMultipleObjects(nth-1,hth+1,TRUE,INFINITE); 238 | #else 239 | pthread_t hth[NOI_MAX_THREADS]; 240 | for ( int t=1; t!=nth; ++t ) { 241 | pthread_create(&hth[t], NULL, noi_kmeans_thread, tinfo + t); 242 | } 243 | noi_kmeans_single_thread(tinfo); 244 | for ( int t=1; t!=nth; ++t ) { 245 | pthread_join(hth[t], NULL); 246 | } 247 | #endif 248 | #else 249 | noi_kmeans_thread_t nkt = (noi_kmeans_thread_t ) { blocks,nblocks,index,center,K,dist }; 250 | noi_kmeans_single_thread(&nkt); 251 | #endif 252 | memset(center, 0, K*16*sizeof(int) ); 253 | memset(num, 0, K*sizeof(int)); 254 | for ( int o=0; o!=nblocks; ++o ) { 255 | int k = index[o]; 256 | for ( int i=0; i!=16; ++i ) { 257 | center[k*16 + i] += blocks[o*16 + i]; 258 | } 259 | num[k] ++; 260 | } 261 | for ( int k=0; k!=K; ++k ) { 262 | int numK = num[k]; 263 | if ( numK ) { 264 | for ( int i=1; i!=16; ++i ) { 265 | center[k*16 + i] /= numK; 266 | } 267 | } 268 | } 269 | } 270 | res->index = index; 271 | res->center = center; 272 | free(num); 273 | } 274 | 275 | int * noi_compress_yuv ( int * nblocks, uint8_t * pixels, int w, int h, int profile ) { 276 | int psx, pxshift, psy, pyshift; 277 | switch ( profile ) { 278 | case NOI_YUV_16_1_1: psx = psy = 4; pxshift = pyshift = 2; break; 279 | case NOI_YUV_4_1_1: psx = psy = 2; pxshift = pyshift = 1; break; 280 | case NOI_YUV_2_1_1: psx = 2; psy = 1; pxshift = 1; pyshift = 0; break; 281 | } 282 | int bsx = psx * 4; 283 | int bsy = psy * 4; 284 | int bw = w / bsx; 285 | int bh = h / bsy; 286 | int stride = w * 4; 287 | int numBlocks = bw*bh*profile; *nblocks = numBlocks; 288 | int * blocks = (int *) malloc ( numBlocks*16*sizeof(int) ); 289 | int * block = blocks; 290 | for ( int by=0; by!=bh; ++by ) { 291 | for ( int bx=0; bx!=bw; ++bx ) { 292 | int * ublock = block; block += 16; 293 | int * vblock = block; block += 16; 294 | memset(ublock, 0, 16*sizeof(int)); 295 | memset(vblock, 0, 16*sizeof(int)); 296 | for ( int ty=0; ty!=psy; ty++ ) { 297 | for ( int tx=0; tx!=psx; tx++ ) { 298 | int * yblock = block; block += 16; 299 | for ( int y=0; y!=4; y++ ) { 300 | for ( int x=0; x!=4; x++ ) { 301 | int uvy = ty*4+y; 302 | int uvx = tx*4+x; 303 | int ofs = (by*bsy+uvy)*stride + (bx*bsx+uvx)*4; 304 | int r = pixels[ofs+0]; 305 | int g = pixels[ofs+1]; 306 | int b = pixels[ofs+2]; 307 | int Y, U, V; 308 | noi_rgb2yuv(r,g,b,&Y,&U,&V); 309 | yblock[y*4+x] = Y; 310 | int t = (uvy>>pyshift)*4 + (uvx>>pxshift); 311 | ublock[t] += U; 312 | vblock[t] += V; 313 | } 314 | } 315 | noi_hdt4x4(yblock); 316 | } 317 | } 318 | for ( int t=0; t!=16; ++t ) { 319 | ublock[t] /= (psy*psx); 320 | vblock[t] /= (psy*psx); 321 | } 322 | noi_hdt4x4(ublock); 323 | noi_hdt4x4(vblock); 324 | } 325 | } 326 | return blocks; 327 | } 328 | 329 | int * noi_compress_rgb ( int * nblocks, uint8_t * pixels, int w, int h, int profile, int hasB ) { 330 | int bw = w / 4; 331 | int bh = h / 4; 332 | int stride = w * 4; 333 | int numBlocks = bw*bh*profile; *nblocks = numBlocks; 334 | int * blocks = (int *) malloc ( numBlocks*16*sizeof(int) ); 335 | int * block = blocks; 336 | for ( int by=0; by!=bh; ++by ) { 337 | for ( int bx=0; bx!=bw; ++bx ) { 338 | int * rblock = block; block += 16; 339 | int * gblock = block; block += 16; 340 | int * bblock = block; if ( hasB ) block += 16; 341 | for ( int y=0; y!=4; y++ ) { 342 | for ( int x=0; x!=4; x++ ) { 343 | int ofs = (by*4+y)*stride + (bx*4+x)*4; 344 | int t = y*4 + x; 345 | rblock[t] = pixels[ofs+0]; 346 | gblock[t] = pixels[ofs+1]; 347 | if ( hasB ) bblock[t] = pixels[ofs+2]; 348 | } 349 | } 350 | noi_hdt4x4(rblock); 351 | noi_hdt4x4(gblock); 352 | if ( hasB ) noi_hdt4x4(bblock); 353 | } 354 | } 355 | return blocks; 356 | } 357 | 358 | int * noi_compress_greyscale ( int * nblocks, uint8_t * pixels, int w, int h, int profile ) { 359 | int bw = w / 4; 360 | int bh = h / 4; 361 | int stride = w * 4; 362 | int numBlocks = bw*bh*profile; *nblocks = numBlocks; 363 | int * blocks = (int *) malloc ( numBlocks*16*sizeof(int) ); 364 | int * block = blocks; 365 | for ( int by=0; by!=bh; ++by ) { 366 | for ( int bx=0; bx!=bw; ++bx ) { 367 | for ( int y=0; y!=4; y++ ) { 368 | for ( int x=0; x!=4; x++ ) { 369 | int ofs = (by*4+y)*stride + (bx*4+x)*4; 370 | int t = y*4 + x; 371 | int R = pixels[ofs+0]; 372 | int G = pixels[ofs+1]; 373 | int B = pixels[ofs+2]; 374 | block[t] = (int)(0.299*R + 0.587*G + 0.114*B); 375 | } 376 | } 377 | noi_hdt4x4(block); 378 | block += 16; 379 | } 380 | } 381 | return blocks; 382 | } 383 | 384 | void * noi_compress ( uint8_t * pixels, int w, int h, int * bytes, int profile ) { 385 | int numBlocks; int * blocks = NULL; 386 | switch ( profile ) { 387 | case NOI_YUV_16_1_1: case NOI_YUV_4_1_1: case NOI_YUV_2_1_1: 388 | blocks = noi_compress_yuv(&numBlocks, pixels, w, h, profile); break; 389 | case NOI_RGB_1_1_1: 390 | blocks = noi_compress_rgb(&numBlocks, pixels, w, h, profile, 1); break; 391 | case NOI_RG_1_1_0: 392 | blocks = noi_compress_rgb(&numBlocks, pixels, w, h, profile, 0); break; 393 | case NOI_Y_1_0_0: 394 | blocks = noi_compress_greyscale(&numBlocks, pixels, w, h, profile); break; 395 | default: 396 | return NULL; 397 | } 398 | noi_kmeans_t res3; 399 | noi_kmeans(&res3, blocks, numBlocks, NOI_K3, noi_dist3); 400 | noi_kmeans_t res5; 401 | noi_kmeans(&res5, blocks, numBlocks, NOI_K5, noi_dist5); 402 | noi_kmeans_t res7; 403 | noi_kmeans(&res7, blocks, numBlocks, NOI_K7, noi_dist7); 404 | int osize = numBlocks*5 + NOI_K3*3*sizeof(int16_t) + NOI_K5*5*sizeof(int16_t) + NOI_K7*7*sizeof(int16_t) + sizeof(noi_header_t); 405 | if ( bytes ) *bytes = osize; 406 | uint8_t * outbytes = (uint8_t *) malloc(osize); 407 | uint8_t * out = outbytes; 408 | noi_header_t * header = (noi_header_t *) out; out += sizeof(noi_header_t); 409 | *header = (noi_header_t) { NOI_MAGIC, profile, w, h }; 410 | for ( int o=0; o!=numBlocks; ++o ) { 411 | int i0 = blocks[o*16+0] & 4095; 412 | int i3 = res3.index[o]; 413 | int i5 = res5.index[o]; 414 | int i7 = res7.index[o]; 415 | *((uint16_t *)out) = i0 | ((i7&0x100)<<4) | ((i5&0x100)<<5) | ((i3&0x300)<<6); 416 | out[2] = i3; 417 | out[3] = i5; 418 | out[4] = i7; 419 | out += 5; 420 | } 421 | for ( int k=0; k!=NOI_K3; ++k ) { 422 | int16_t * B = (int16_t *) out; out += 3*sizeof(int16_t); 423 | int * c3 = res3.center + k*16; 424 | B[0] = c3[NOI_B3_0]; B[1] = c3[NOI_B3_1]; B[2] = c3[NOI_B3_2]; 425 | } 426 | for ( int k=0; k!=NOI_K5; ++k ) { 427 | int16_t * B = (int16_t *) out; out += 5*sizeof(int16_t); 428 | int * c5 = res5.center + k*16; 429 | B[0] = c5[NOI_B5_0]; B[1] = c5[NOI_B5_1]; B[2] = c5[NOI_B5_2]; B[3] = c5[NOI_B5_3]; B[4] = c5[NOI_B5_4]; 430 | } 431 | for ( int k=0; k!=NOI_K7; ++k ) { 432 | int16_t * B = (int16_t *) out; out += 7*sizeof(int16_t); 433 | int * c7 = res7.center + k*16; 434 | B[0] = c7[NOI_B7_0]; B[1] = c7[NOI_B7_1]; B[2] = c7[NOI_B7_2]; B[3] = c7[NOI_B7_3]; B[4] = c7[NOI_B7_4]; B[5] = c7[NOI_B7_5]; B[6] = c7[NOI_B7_6]; 435 | } 436 | free(blocks); 437 | noi_free_means(&res3); 438 | noi_free_means(&res5); 439 | noi_free_means(&res7); 440 | return outbytes; 441 | } 442 | 443 | #define NOI_HDT2X2(a,b,c,d) \ 444 | aab = a + b; sab = a - b; acd = c + d; scd = c - d; \ 445 | a = aab + acd; b = sab + scd; c = aab - acd; d = sab - scd; 446 | 447 | #define NOI_HDT2X2S(a,b,c,d) \ 448 | aab = a + b; sab = a - b; acd = c + d; scd = c - d; \ 449 | a = (aab + acd)>>4; b = (sab + scd)>>4; c = (aab - acd)>>4; d = (sab - scd)>>4; 450 | 451 | void noi_decompress_5_16 ( uint8_t * in, int * blocks, int16_t * center3, int profile ) { 452 | int16_t * center5 = center3 + NOI_K3*3; 453 | int16_t * center7 = center5 + NOI_K5*5; 454 | int * fb = blocks; 455 | for ( int o=0; o!=profile; ++o ) { 456 | int in0 = *((uint16_t *)in); 457 | int in3 = in[2] + ((in0 & 0xC000)>>6); 458 | int in5 = in[3] + ((in0 & 0x2000)>>5); 459 | int in7 = in[4] + ((in0 & 0x1000)>>4); 460 | int16_t * c3 = center3 + in3*3; 461 | int16_t * c5 = center5 + in5*5; 462 | int16_t * c7 = center7 + in7*7; 463 | int a00 = in0 & 4095; 464 | int a02 = c3[0]; int a20 = c3[1]; int a22 = c3[2]; 465 | int a01 = c5[0]; int a10 = c5[1]; int a11 = c5[2]; int a12 = c5[3]; int a21 = c5[4]; 466 | int a03 = c7[0]; int a13 = c7[1]; int a23 = c7[2]; int a30 = c7[3]; int a31 = c7[4]; int a32 = c7[5]; int a33 = c7[6]; 467 | int aab, sab, acd, scd; 468 | NOI_HDT2X2(a00,a01,a10,a11); NOI_HDT2X2(a02,a03,a12,a13); NOI_HDT2X2(a20,a21,a30,a31); NOI_HDT2X2(a22,a23,a32,a33); 469 | NOI_HDT2X2S(a00,a02,a20,a22); NOI_HDT2X2S(a01,a03,a21,a23); NOI_HDT2X2S(a10,a12,a30,a32); NOI_HDT2X2S(a11,a13,a31,a33); 470 | fb[0] = a00; fb[1] = a01; fb[2] = a02; fb[3] = a03; 471 | fb[4] = a10; fb[5] = a11; fb[6] = a12; fb[7] = a13; 472 | fb[8] = a20; fb[9] = a21; fb[10] = a22; fb[11] = a23; 473 | fb[12] = a30; fb[13] = a31; fb[14] = a32; fb[15] = a33; 474 | fb += 16; 475 | in += 5; 476 | } 477 | } 478 | 479 | #define NOI_YUV2RGB(ofsx,AA) \ 480 | C = 298*(AA-16); \ 481 | bpixels[ofsx*4+0] = noi_saturate((C+DR) >> 8); \ 482 | bpixels[ofsx*4+1] = noi_saturate((C+DG) >> 8); \ 483 | bpixels[ofsx*4+2] = noi_saturate((C+DB) >> 8); \ 484 | bpixels[ofsx*4+3] = 255; 485 | 486 | void noi_convert_colors_YUV_16_1_1 ( uint8_t * in, int16_t * center3, uint8_t * pixels, int stride ) { 487 | int blocks[16*2]; 488 | noi_decompress_5_16(in, blocks, center3, 2); 489 | int * ublock = blocks; 490 | int * vblock = blocks + 16; 491 | int16_t * center5 = center3 + NOI_K3*3; 492 | int16_t * center7 = center5 + NOI_K5*5; 493 | in += 10; 494 | for ( int ty=0; ty!=4; ty++ ) { 495 | for ( int tx=0; tx!=4; tx++ ) { 496 | int in0 = *((uint16_t *)in); 497 | int in3 = in[2] + ((in0 & 0xC000)>>6); 498 | int in5 = in[3] + ((in0 & 0x2000)>>5); 499 | int in7 = in[4] + ((in0 & 0x1000)>>4); 500 | int16_t * c3 = center3 + in3*3; 501 | int16_t * c5 = center5 + in5*5; 502 | int16_t * c7 = center7 + in7*7; 503 | int a00 = in0 & 4095; 504 | int a02 = c3[0]; int a20 = c3[1]; int a22 = c3[2]; 505 | int a01 = c5[0]; int a10 = c5[1]; int a11 = c5[2]; int a12 = c5[3]; int a21 = c5[4]; 506 | int a03 = c7[0]; int a13 = c7[1]; int a23 = c7[2]; int a30 = c7[3]; int a31 = c7[4]; int a32 = c7[5]; int a33 = c7[6]; 507 | int aab, sab, acd, scd; 508 | NOI_HDT2X2(a00,a01,a10,a11); NOI_HDT2X2(a02,a03,a12,a13); NOI_HDT2X2(a20,a21,a30,a31); NOI_HDT2X2(a22,a23,a32,a33); 509 | NOI_HDT2X2S(a00,a02,a20,a22); NOI_HDT2X2S(a01,a03,a21,a23); NOI_HDT2X2S(a10,a12,a30,a32); NOI_HDT2X2S(a11,a13,a31,a33); 510 | // convert Y-block directly from HDT output 511 | uint8_t * bpixels = pixels + ty*4*stride + tx*4*4; int C; 512 | int d = *ublock++ - 128; int e = *vblock++ - 128; 513 | int DR = + 409*e + 128; int DG = - 100*d - 209*e + 128; int DB = + 516*d + 128; 514 | NOI_YUV2RGB(0,a00); NOI_YUV2RGB(1,a01); NOI_YUV2RGB(2,a02); NOI_YUV2RGB(3,a03); bpixels += stride; 515 | NOI_YUV2RGB(0,a10); NOI_YUV2RGB(1,a11); NOI_YUV2RGB(2,a12); NOI_YUV2RGB(3,a13); bpixels += stride; 516 | NOI_YUV2RGB(0,a20); NOI_YUV2RGB(1,a21); NOI_YUV2RGB(2,a22); NOI_YUV2RGB(3,a23); bpixels += stride; 517 | NOI_YUV2RGB(0,a30); NOI_YUV2RGB(1,a31); NOI_YUV2RGB(2,a32); NOI_YUV2RGB(3,a33); 518 | in += 5; 519 | } 520 | } 521 | } 522 | 523 | #undef NOI_YUV2RGB 524 | 525 | #define NOI_Y2RGB(ofsx,AA) \ 526 | bpixels[ofsx*4+0] = bpixels[ofsx*4+1] = bpixels[ofsx*4+2] = noi_saturate(AA); bpixels[ofsx*4+3] = 255; 527 | 528 | void noi_convert_colors_greyscale ( uint8_t * in, int16_t * center3, uint8_t * pixels, int stride ) { 529 | int16_t * center5 = center3 + NOI_K3*3; 530 | int16_t * center7 = center5 + NOI_K5*5; 531 | int in0 = *((uint16_t *)in); 532 | int in3 = in[2] + ((in0 & 0xC000)>>6); 533 | int in5 = in[3] + ((in0 & 0x2000)>>5); 534 | int in7 = in[4] + ((in0 & 0x1000)>>4); 535 | int16_t * c3 = center3 + in3*3; 536 | int16_t * c5 = center5 + in5*5; 537 | int16_t * c7 = center7 + in7*7; 538 | int a00 = in0 & 4095; 539 | int a02 = c3[0]; int a20 = c3[1]; int a22 = c3[2]; 540 | int a01 = c5[0]; int a10 = c5[1]; int a11 = c5[2]; int a12 = c5[3]; int a21 = c5[4]; 541 | int a03 = c7[0]; int a13 = c7[1]; int a23 = c7[2]; int a30 = c7[3]; int a31 = c7[4]; int a32 = c7[5]; int a33 = c7[6]; 542 | int aab, sab, acd, scd; 543 | NOI_HDT2X2(a00,a01,a10,a11); NOI_HDT2X2(a02,a03,a12,a13); NOI_HDT2X2(a20,a21,a30,a31); NOI_HDT2X2(a22,a23,a32,a33); 544 | NOI_HDT2X2S(a00,a02,a20,a22); NOI_HDT2X2S(a01,a03,a21,a23); NOI_HDT2X2S(a10,a12,a30,a32); NOI_HDT2X2S(a11,a13,a31,a33); 545 | // convert Y-block directly from HDT output 546 | uint8_t * bpixels = pixels; 547 | NOI_Y2RGB(0,a00); NOI_Y2RGB(1,a01); NOI_Y2RGB(2,a02); NOI_Y2RGB(3,a03); bpixels += stride; 548 | NOI_Y2RGB(0,a10); NOI_Y2RGB(1,a11); NOI_Y2RGB(2,a12); NOI_Y2RGB(3,a13); bpixels += stride; 549 | NOI_Y2RGB(0,a20); NOI_Y2RGB(1,a21); NOI_Y2RGB(2,a22); NOI_Y2RGB(3,a23); bpixels += stride; 550 | NOI_Y2RGB(0,a30); NOI_Y2RGB(1,a31); NOI_Y2RGB(2,a32); NOI_Y2RGB(3,a33); 551 | in += 5; 552 | } 553 | 554 | #undef NOI_Y2RGB 555 | 556 | #undef NOI_HDT2X2 557 | #undef NOI_HDT2X2S 558 | 559 | void noi_convert_colors_YUV_4_1_1 ( uint8_t * in, int16_t * c3, uint8_t * pixels, int stride ) { 560 | int blocks[16*6]; 561 | noi_decompress_5_16(in, blocks, c3, 6 ); 562 | int * ublock = blocks; 563 | int * vblock = blocks + 16; 564 | int * yblock = blocks + 32; 565 | for ( int ty=0; ty!=2; ty++ ) { 566 | for ( int tx=0; tx!=2; tx++ ) { 567 | uint8_t * bpixels = pixels + ty*4*stride + tx*4*4; 568 | for ( int y=0; y!=4; y++ ) { 569 | for ( int x=0; x!=4; x++ ) { 570 | int uvy = ty*4+y; 571 | int uvx = tx*4+x; 572 | int t = (uvy>>1)*4 + (uvx>>1); 573 | int U = ublock[t]; int V = vblock[t]; int Y = yblock[y*4+x]; 574 | int R, G, B; noi_yuv2rgb(Y, U, V, &R, &G, &B); 575 | bpixels[x*4+0] = noi_saturate(R); bpixels[x*4+1] = noi_saturate(G); bpixels[x*4+2] = noi_saturate(B); bpixels[x*4+3] = 255; 576 | } 577 | bpixels += stride; 578 | } 579 | yblock += 16; 580 | } 581 | } 582 | } 583 | 584 | void noi_convert_colors_YUV_2_1_1 ( uint8_t * in, int16_t * c3, uint8_t * pixels, int stride ) { 585 | int blocks[16*4]; 586 | noi_decompress_5_16(in, blocks, c3, 4 ); 587 | int * ublock = blocks; 588 | int * vblock = blocks + 16; 589 | int * yblock = blocks + 32; 590 | for ( int tx=0; tx!=2; tx++ ) { 591 | uint8_t * bpixels = pixels + tx*4*4; 592 | for ( int y=0; y!=4; y++ ) { 593 | for ( int x=0; x!=4; x++ ) { 594 | int uvy = y; 595 | int uvx = tx*4+x; 596 | int t = uvy*4 + (uvx>>1); 597 | int U = ublock[t]; int V = vblock[t]; int Y = yblock[y*4+x]; 598 | int R, G, B; noi_yuv2rgb(Y, U, V, &R, &G, &B); 599 | bpixels[x*4+0] = noi_saturate(R); bpixels[x*4+1] = noi_saturate(G); bpixels[x*4+2] = noi_saturate(B); bpixels[x*4+3] = 255; 600 | } 601 | bpixels += stride; 602 | } 603 | yblock += 16; 604 | } 605 | } 606 | 607 | void noi_convert_colors_RGB_1_1_1 ( uint8_t * in, int16_t * c3, uint8_t * pixels, int stride ) { 608 | int blocks[16*3]; 609 | noi_decompress_5_16(in, blocks, c3, 3 ); 610 | int * rblock = blocks; 611 | int * gblock = blocks + 16; 612 | int * bblock = blocks + 32; 613 | uint8_t * bpixels = pixels; 614 | for ( int y=0; y!=4; y++ ) { 615 | for ( int x=0; x!=4; x++ ) { 616 | int t = y*4 + x; 617 | int R = rblock[t]; int G = gblock[t]; int B = bblock[t]; 618 | bpixels[x*4+0] = noi_saturate(R); bpixels[x*4+1] = noi_saturate(G); bpixels[x*4+2] = noi_saturate(B); bpixels[x*4+3] = 255; 619 | } 620 | bpixels += stride; 621 | } 622 | } 623 | 624 | void noi_convert_colors_RG_1_1_0 ( uint8_t * in, int16_t * c3, uint8_t * pixels, int stride ) { 625 | int blocks[16*2]; 626 | noi_decompress_5_16(in, blocks, c3, 2 ); 627 | int * rblock = blocks; 628 | int * gblock = blocks + 16; 629 | uint8_t * bpixels = pixels; 630 | for ( int y=0; y!=4; y++ ) { 631 | for ( int x=0; x!=4; x++ ) { 632 | int t = y*4 + x; 633 | int R = rblock[t]; int G = gblock[t]; 634 | bpixels[x*4+0] = noi_saturate(R); bpixels[x*4+1] = noi_saturate(G); bpixels[x*4+2] = bpixels[x*4+3] = 255; 635 | } 636 | bpixels += stride; 637 | } 638 | } 639 | 640 | uint8_t * noi_decompress ( void * bytes, int * W, int * H, uint8_t * pixels ) { 641 | uint8_t * in = (uint8_t *) bytes; 642 | noi_header_t * header = (noi_header_t *) in; in += sizeof(noi_header_t); 643 | if ( header->magic != NOI_MAGIC ) return NULL; 644 | int w = header->width; if ( W ) *W = w; 645 | int h = header->height; if ( H ) *H = h; 646 | int profile = header->profile; 647 | int psx, psy; 648 | void ( *decmpress_block ) ( uint8_t * in, int16_t * c3, uint8_t * pixels, int stride ); 649 | switch ( profile ) { 650 | case NOI_YUV_16_1_1: psx = psy = 4; decmpress_block = noi_convert_colors_YUV_16_1_1; break; 651 | case NOI_YUV_4_1_1: psx = psy = 2; decmpress_block = noi_convert_colors_YUV_4_1_1; break; 652 | case NOI_YUV_2_1_1: psx = 2; psy = 1; decmpress_block = noi_convert_colors_YUV_2_1_1; break; 653 | case NOI_RGB_1_1_1: psx = psy = 1; decmpress_block = noi_convert_colors_RGB_1_1_1; break; 654 | case NOI_RG_1_1_0: psx = psy = 1; decmpress_block = noi_convert_colors_RG_1_1_0; break; 655 | case NOI_Y_1_0_0: psx = psy = 1; decmpress_block = noi_convert_colors_greyscale; break; 656 | } 657 | int bsx = psx * 4; 658 | int bsy = psy * 4; 659 | int bw = w / bsx; 660 | int bh = h / bsy; 661 | int numBlocks = bw*bh*profile; 662 | if ( !pixels ) pixels = (uint8_t *) malloc(w*h*4); 663 | int stride = w * 4; 664 | int16_t * c3 = (int16_t *)(in + numBlocks * 5); 665 | for ( int by=0; by!=bh; ++by ) { 666 | for ( int bx=0; bx!=bw; ++bx ) { 667 | (*decmpress_block) (in, c3, pixels + (bx*bsx*4) + (by*bsy)*stride, stride ); 668 | in += profile*5; 669 | } 670 | } 671 | return pixels; 672 | } 673 | 674 | void noi_image_size ( void * bytes, int * W, int * H ) { 675 | noi_header_t * header = (noi_header_t *) bytes; 676 | if ( header->magic != NOI_MAGIC ) return; 677 | int w = header->width; if ( W ) *W = w; 678 | int h = header->height; if ( H ) *H = h; 679 | } 680 | 681 | const char * noi_profile_name ( int profile ) { 682 | switch ( profile ) { 683 | case NOI_YUV_16_1_1: return "YUV_16_1_1"; 684 | case NOI_YUV_4_1_1: return "YUV_4_1_1"; 685 | case NOI_YUV_2_1_1: return "YUV_2_1_1"; 686 | case NOI_RGB_1_1_1: return "RGB_1_1_1"; 687 | case NOI_RG_1_1_0: return "RG_1_1_0"; 688 | case NOI_Y_1_0_0: return "Y_1_0_0"; 689 | default: return NULL; 690 | } 691 | } 692 | 693 | void noi_profile_block_size ( int profile, int * bx, int * by ) { 694 | switch ( profile ) { 695 | case NOI_YUV_16_1_1: *bx=16; *by=16; break; 696 | case NOI_YUV_4_1_1: *bx=8; *by=8; break; 697 | case NOI_YUV_2_1_1: *bx=8; *by=4; break; 698 | case NOI_RGB_1_1_1: *bx=4; *by=4; break; 699 | case NOI_RG_1_1_0: *bx=4; *by=4; break; 700 | case NOI_Y_1_0_0: *bx=4; *by=4; break; 701 | } 702 | } 703 | 704 | #endif 705 | #endif 706 | -------------------------------------------------------------------------------- /stb/stb_image_write.h: -------------------------------------------------------------------------------- 1 | /* stb_image_write - v1.16 - public domain - http://nothings.org/stb 2 | writes out PNG/BMP/TGA/JPEG/HDR images to C stdio - Sean Barrett 2010-2015 3 | no warranty implied; use at your own risk 4 | 5 | Before #including, 6 | 7 | #define STB_IMAGE_WRITE_IMPLEMENTATION 8 | 9 | in the file that you want to have the implementation. 10 | 11 | Will probably not work correctly with strict-aliasing optimizations. 12 | 13 | ABOUT: 14 | 15 | This header file is a library for writing images to C stdio or a callback. 16 | 17 | The PNG output is not optimal; it is 20-50% larger than the file 18 | written by a decent optimizing implementation; though providing a custom 19 | zlib compress function (see STBIW_ZLIB_COMPRESS) can mitigate that. 20 | This library is designed for source code compactness and simplicity, 21 | not optimal image file size or run-time performance. 22 | 23 | BUILDING: 24 | 25 | You can #define STBIW_ASSERT(x) before the #include to avoid using assert.h. 26 | You can #define STBIW_MALLOC(), STBIW_REALLOC(), and STBIW_FREE() to replace 27 | malloc,realloc,free. 28 | You can #define STBIW_MEMMOVE() to replace memmove() 29 | You can #define STBIW_ZLIB_COMPRESS to use a custom zlib-style compress function 30 | for PNG compression (instead of the builtin one), it must have the following signature: 31 | unsigned char * my_compress(unsigned char *data, int data_len, int *out_len, int quality); 32 | The returned data will be freed with STBIW_FREE() (free() by default), 33 | so it must be heap allocated with STBIW_MALLOC() (malloc() by default), 34 | 35 | UNICODE: 36 | 37 | If compiling for Windows and you wish to use Unicode filenames, compile 38 | with 39 | #define STBIW_WINDOWS_UTF8 40 | and pass utf8-encoded filenames. Call stbiw_convert_wchar_to_utf8 to convert 41 | Windows wchar_t filenames to utf8. 42 | 43 | USAGE: 44 | 45 | There are five functions, one for each image file format: 46 | 47 | int stbi_write_png(char const *filename, int w, int h, int comp, const void *data, int stride_in_bytes); 48 | int stbi_write_bmp(char const *filename, int w, int h, int comp, const void *data); 49 | int stbi_write_tga(char const *filename, int w, int h, int comp, const void *data); 50 | int stbi_write_jpg(char const *filename, int w, int h, int comp, const void *data, int quality); 51 | int stbi_write_hdr(char const *filename, int w, int h, int comp, const float *data); 52 | 53 | void stbi_flip_vertically_on_write(int flag); // flag is non-zero to flip data vertically 54 | 55 | There are also five equivalent functions that use an arbitrary write function. You are 56 | expected to open/close your file-equivalent before and after calling these: 57 | 58 | int stbi_write_png_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const void *data, int stride_in_bytes); 59 | int stbi_write_bmp_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const void *data); 60 | int stbi_write_tga_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const void *data); 61 | int stbi_write_hdr_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const float *data); 62 | int stbi_write_jpg_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data, int quality); 63 | 64 | where the callback is: 65 | void stbi_write_func(void *context, void *data, int size); 66 | 67 | You can configure it with these global variables: 68 | int stbi_write_tga_with_rle; // defaults to true; set to 0 to disable RLE 69 | int stbi_write_png_compression_level; // defaults to 8; set to higher for more compression 70 | int stbi_write_force_png_filter; // defaults to -1; set to 0..5 to force a filter mode 71 | 72 | 73 | You can define STBI_WRITE_NO_STDIO to disable the file variant of these 74 | functions, so the library will not use stdio.h at all. However, this will 75 | also disable HDR writing, because it requires stdio for formatted output. 76 | 77 | Each function returns 0 on failure and non-0 on success. 78 | 79 | The functions create an image file defined by the parameters. The image 80 | is a rectangle of pixels stored from left-to-right, top-to-bottom. 81 | Each pixel contains 'comp' channels of data stored interleaved with 8-bits 82 | per channel, in the following order: 1=Y, 2=YA, 3=RGB, 4=RGBA. (Y is 83 | monochrome color.) The rectangle is 'w' pixels wide and 'h' pixels tall. 84 | The *data pointer points to the first byte of the top-left-most pixel. 85 | For PNG, "stride_in_bytes" is the distance in bytes from the first byte of 86 | a row of pixels to the first byte of the next row of pixels. 87 | 88 | PNG creates output files with the same number of components as the input. 89 | The BMP format expands Y to RGB in the file format and does not 90 | output alpha. 91 | 92 | PNG supports writing rectangles of data even when the bytes storing rows of 93 | data are not consecutive in memory (e.g. sub-rectangles of a larger image), 94 | by supplying the stride between the beginning of adjacent rows. The other 95 | formats do not. (Thus you cannot write a native-format BMP through the BMP 96 | writer, both because it is in BGR order and because it may have padding 97 | at the end of the line.) 98 | 99 | PNG allows you to set the deflate compression level by setting the global 100 | variable 'stbi_write_png_compression_level' (it defaults to 8). 101 | 102 | HDR expects linear float data. Since the format is always 32-bit rgb(e) 103 | data, alpha (if provided) is discarded, and for monochrome data it is 104 | replicated across all three channels. 105 | 106 | TGA supports RLE or non-RLE compressed data. To use non-RLE-compressed 107 | data, set the global variable 'stbi_write_tga_with_rle' to 0. 108 | 109 | JPEG does ignore alpha channels in input data; quality is between 1 and 100. 110 | Higher quality looks better but results in a bigger image. 111 | JPEG baseline (no JPEG progressive). 112 | 113 | CREDITS: 114 | 115 | 116 | Sean Barrett - PNG/BMP/TGA 117 | Baldur Karlsson - HDR 118 | Jean-Sebastien Guay - TGA monochrome 119 | Tim Kelsey - misc enhancements 120 | Alan Hickman - TGA RLE 121 | Emmanuel Julien - initial file IO callback implementation 122 | Jon Olick - original jo_jpeg.cpp code 123 | Daniel Gibson - integrate JPEG, allow external zlib 124 | Aarni Koskela - allow choosing PNG filter 125 | 126 | bugfixes: 127 | github:Chribba 128 | Guillaume Chereau 129 | github:jry2 130 | github:romigrou 131 | Sergio Gonzalez 132 | Jonas Karlsson 133 | Filip Wasil 134 | Thatcher Ulrich 135 | github:poppolopoppo 136 | Patrick Boettcher 137 | github:xeekworx 138 | Cap Petschulat 139 | Simon Rodriguez 140 | Ivan Tikhonov 141 | github:ignotion 142 | Adam Schackart 143 | Andrew Kensler 144 | 145 | LICENSE 146 | 147 | See end of file for license information. 148 | 149 | */ 150 | 151 | #ifndef INCLUDE_STB_IMAGE_WRITE_H 152 | #define INCLUDE_STB_IMAGE_WRITE_H 153 | 154 | #include 155 | 156 | // if STB_IMAGE_WRITE_STATIC causes problems, try defining STBIWDEF to 'inline' or 'static inline' 157 | #ifndef STBIWDEF 158 | #ifdef STB_IMAGE_WRITE_STATIC 159 | #define STBIWDEF static 160 | #else 161 | #ifdef __cplusplus 162 | #define STBIWDEF extern "C" 163 | #else 164 | #define STBIWDEF extern 165 | #endif 166 | #endif 167 | #endif 168 | 169 | #ifndef STB_IMAGE_WRITE_STATIC // C++ forbids static forward declarations 170 | STBIWDEF int stbi_write_tga_with_rle; 171 | STBIWDEF int stbi_write_png_compression_level; 172 | STBIWDEF int stbi_write_force_png_filter; 173 | #endif 174 | 175 | #ifndef STBI_WRITE_NO_STDIO 176 | STBIWDEF int stbi_write_png(char const *filename, int w, int h, int comp, const void *data, int stride_in_bytes); 177 | STBIWDEF int stbi_write_bmp(char const *filename, int w, int h, int comp, const void *data); 178 | STBIWDEF int stbi_write_tga(char const *filename, int w, int h, int comp, const void *data); 179 | STBIWDEF int stbi_write_hdr(char const *filename, int w, int h, int comp, const float *data); 180 | STBIWDEF int stbi_write_jpg(char const *filename, int x, int y, int comp, const void *data, int quality); 181 | 182 | #ifdef STBIW_WINDOWS_UTF8 183 | STBIWDEF int stbiw_convert_wchar_to_utf8(char *buffer, size_t bufferlen, const wchar_t* input); 184 | #endif 185 | #endif 186 | 187 | typedef void stbi_write_func(void *context, void *data, int size); 188 | 189 | STBIWDEF int stbi_write_png_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const void *data, int stride_in_bytes); 190 | STBIWDEF int stbi_write_bmp_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const void *data); 191 | STBIWDEF int stbi_write_tga_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const void *data); 192 | STBIWDEF int stbi_write_hdr_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const float *data); 193 | STBIWDEF int stbi_write_jpg_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data, int quality); 194 | 195 | STBIWDEF void stbi_flip_vertically_on_write(int flip_boolean); 196 | 197 | #endif//INCLUDE_STB_IMAGE_WRITE_H 198 | 199 | #ifdef STB_IMAGE_WRITE_IMPLEMENTATION 200 | 201 | #ifdef _WIN32 202 | #ifndef _CRT_SECURE_NO_WARNINGS 203 | #define _CRT_SECURE_NO_WARNINGS 204 | #endif 205 | #ifndef _CRT_NONSTDC_NO_DEPRECATE 206 | #define _CRT_NONSTDC_NO_DEPRECATE 207 | #endif 208 | #endif 209 | 210 | #ifndef STBI_WRITE_NO_STDIO 211 | #include 212 | #endif // STBI_WRITE_NO_STDIO 213 | 214 | #include 215 | #include 216 | #include 217 | #include 218 | 219 | #if defined(STBIW_MALLOC) && defined(STBIW_FREE) && (defined(STBIW_REALLOC) || defined(STBIW_REALLOC_SIZED)) 220 | // ok 221 | #elif !defined(STBIW_MALLOC) && !defined(STBIW_FREE) && !defined(STBIW_REALLOC) && !defined(STBIW_REALLOC_SIZED) 222 | // ok 223 | #else 224 | #error "Must define all or none of STBIW_MALLOC, STBIW_FREE, and STBIW_REALLOC (or STBIW_REALLOC_SIZED)." 225 | #endif 226 | 227 | #ifndef STBIW_MALLOC 228 | #define STBIW_MALLOC(sz) malloc(sz) 229 | #define STBIW_REALLOC(p,newsz) realloc(p,newsz) 230 | #define STBIW_FREE(p) free(p) 231 | #endif 232 | 233 | #ifndef STBIW_REALLOC_SIZED 234 | #define STBIW_REALLOC_SIZED(p,oldsz,newsz) STBIW_REALLOC(p,newsz) 235 | #endif 236 | 237 | 238 | #ifndef STBIW_MEMMOVE 239 | #define STBIW_MEMMOVE(a,b,sz) memmove(a,b,sz) 240 | #endif 241 | 242 | 243 | #ifndef STBIW_ASSERT 244 | #include 245 | #define STBIW_ASSERT(x) assert(x) 246 | #endif 247 | 248 | #define STBIW_UCHAR(x) (unsigned char) ((x) & 0xff) 249 | 250 | #ifdef STB_IMAGE_WRITE_STATIC 251 | static int stbi_write_png_compression_level = 8; 252 | static int stbi_write_tga_with_rle = 1; 253 | static int stbi_write_force_png_filter = -1; 254 | #else 255 | int stbi_write_png_compression_level = 8; 256 | int stbi_write_tga_with_rle = 1; 257 | int stbi_write_force_png_filter = -1; 258 | #endif 259 | 260 | static int stbi__flip_vertically_on_write = 0; 261 | 262 | STBIWDEF void stbi_flip_vertically_on_write(int flag) 263 | { 264 | stbi__flip_vertically_on_write = flag; 265 | } 266 | 267 | typedef struct 268 | { 269 | stbi_write_func *func; 270 | void *context; 271 | unsigned char buffer[64]; 272 | int buf_used; 273 | } stbi__write_context; 274 | 275 | // initialize a callback-based context 276 | static void stbi__start_write_callbacks(stbi__write_context *s, stbi_write_func *c, void *context) 277 | { 278 | s->func = c; 279 | s->context = context; 280 | } 281 | 282 | #ifndef STBI_WRITE_NO_STDIO 283 | 284 | static void stbi__stdio_write(void *context, void *data, int size) 285 | { 286 | fwrite(data,1,size,(FILE*) context); 287 | } 288 | 289 | #if defined(_WIN32) && defined(STBIW_WINDOWS_UTF8) 290 | #ifdef __cplusplus 291 | #define STBIW_EXTERN extern "C" 292 | #else 293 | #define STBIW_EXTERN extern 294 | #endif 295 | STBIW_EXTERN __declspec(dllimport) int __stdcall MultiByteToWideChar(unsigned int cp, unsigned long flags, const char *str, int cbmb, wchar_t *widestr, int cchwide); 296 | STBIW_EXTERN __declspec(dllimport) int __stdcall WideCharToMultiByte(unsigned int cp, unsigned long flags, const wchar_t *widestr, int cchwide, char *str, int cbmb, const char *defchar, int *used_default); 297 | 298 | STBIWDEF int stbiw_convert_wchar_to_utf8(char *buffer, size_t bufferlen, const wchar_t* input) 299 | { 300 | return WideCharToMultiByte(65001 /* UTF8 */, 0, input, -1, buffer, (int) bufferlen, NULL, NULL); 301 | } 302 | #endif 303 | 304 | static FILE *stbiw__fopen(char const *filename, char const *mode) 305 | { 306 | FILE *f; 307 | #if defined(_WIN32) && defined(STBIW_WINDOWS_UTF8) 308 | wchar_t wMode[64]; 309 | wchar_t wFilename[1024]; 310 | if (0 == MultiByteToWideChar(65001 /* UTF8 */, 0, filename, -1, wFilename, sizeof(wFilename)/sizeof(*wFilename))) 311 | return 0; 312 | 313 | if (0 == MultiByteToWideChar(65001 /* UTF8 */, 0, mode, -1, wMode, sizeof(wMode)/sizeof(*wMode))) 314 | return 0; 315 | 316 | #if defined(_MSC_VER) && _MSC_VER >= 1400 317 | if (0 != _wfopen_s(&f, wFilename, wMode)) 318 | f = 0; 319 | #else 320 | f = _wfopen(wFilename, wMode); 321 | #endif 322 | 323 | #elif defined(_MSC_VER) && _MSC_VER >= 1400 324 | if (0 != fopen_s(&f, filename, mode)) 325 | f=0; 326 | #else 327 | f = fopen(filename, mode); 328 | #endif 329 | return f; 330 | } 331 | 332 | static int stbi__start_write_file(stbi__write_context *s, const char *filename) 333 | { 334 | FILE *f = stbiw__fopen(filename, "wb"); 335 | stbi__start_write_callbacks(s, stbi__stdio_write, (void *) f); 336 | return f != NULL; 337 | } 338 | 339 | static void stbi__end_write_file(stbi__write_context *s) 340 | { 341 | fclose((FILE *)s->context); 342 | } 343 | 344 | #endif // !STBI_WRITE_NO_STDIO 345 | 346 | typedef unsigned int stbiw_uint32; 347 | typedef int stb_image_write_test[sizeof(stbiw_uint32)==4 ? 1 : -1]; 348 | 349 | static void stbiw__writefv(stbi__write_context *s, const char *fmt, va_list v) 350 | { 351 | while (*fmt) { 352 | switch (*fmt++) { 353 | case ' ': break; 354 | case '1': { unsigned char x = STBIW_UCHAR(va_arg(v, int)); 355 | s->func(s->context,&x,1); 356 | break; } 357 | case '2': { int x = va_arg(v,int); 358 | unsigned char b[2]; 359 | b[0] = STBIW_UCHAR(x); 360 | b[1] = STBIW_UCHAR(x>>8); 361 | s->func(s->context,b,2); 362 | break; } 363 | case '4': { stbiw_uint32 x = va_arg(v,int); 364 | unsigned char b[4]; 365 | b[0]=STBIW_UCHAR(x); 366 | b[1]=STBIW_UCHAR(x>>8); 367 | b[2]=STBIW_UCHAR(x>>16); 368 | b[3]=STBIW_UCHAR(x>>24); 369 | s->func(s->context,b,4); 370 | break; } 371 | default: 372 | STBIW_ASSERT(0); 373 | return; 374 | } 375 | } 376 | } 377 | 378 | static void stbiw__writef(stbi__write_context *s, const char *fmt, ...) 379 | { 380 | va_list v; 381 | va_start(v, fmt); 382 | stbiw__writefv(s, fmt, v); 383 | va_end(v); 384 | } 385 | 386 | static void stbiw__write_flush(stbi__write_context *s) 387 | { 388 | if (s->buf_used) { 389 | s->func(s->context, &s->buffer, s->buf_used); 390 | s->buf_used = 0; 391 | } 392 | } 393 | 394 | static void stbiw__putc(stbi__write_context *s, unsigned char c) 395 | { 396 | s->func(s->context, &c, 1); 397 | } 398 | 399 | static void stbiw__write1(stbi__write_context *s, unsigned char a) 400 | { 401 | if ((size_t)s->buf_used + 1 > sizeof(s->buffer)) 402 | stbiw__write_flush(s); 403 | s->buffer[s->buf_used++] = a; 404 | } 405 | 406 | static void stbiw__write3(stbi__write_context *s, unsigned char a, unsigned char b, unsigned char c) 407 | { 408 | int n; 409 | if ((size_t)s->buf_used + 3 > sizeof(s->buffer)) 410 | stbiw__write_flush(s); 411 | n = s->buf_used; 412 | s->buf_used = n+3; 413 | s->buffer[n+0] = a; 414 | s->buffer[n+1] = b; 415 | s->buffer[n+2] = c; 416 | } 417 | 418 | static void stbiw__write_pixel(stbi__write_context *s, int rgb_dir, int comp, int write_alpha, int expand_mono, unsigned char *d) 419 | { 420 | unsigned char bg[3] = { 255, 0, 255}, px[3]; 421 | int k; 422 | 423 | if (write_alpha < 0) 424 | stbiw__write1(s, d[comp - 1]); 425 | 426 | switch (comp) { 427 | case 2: // 2 pixels = mono + alpha, alpha is written separately, so same as 1-channel case 428 | case 1: 429 | if (expand_mono) 430 | stbiw__write3(s, d[0], d[0], d[0]); // monochrome bmp 431 | else 432 | stbiw__write1(s, d[0]); // monochrome TGA 433 | break; 434 | case 4: 435 | if (!write_alpha) { 436 | // composite against pink background 437 | for (k = 0; k < 3; ++k) 438 | px[k] = bg[k] + ((d[k] - bg[k]) * d[3]) / 255; 439 | stbiw__write3(s, px[1 - rgb_dir], px[1], px[1 + rgb_dir]); 440 | break; 441 | } 442 | /* FALLTHROUGH */ 443 | case 3: 444 | stbiw__write3(s, d[1 - rgb_dir], d[1], d[1 + rgb_dir]); 445 | break; 446 | } 447 | if (write_alpha > 0) 448 | stbiw__write1(s, d[comp - 1]); 449 | } 450 | 451 | static void stbiw__write_pixels(stbi__write_context *s, int rgb_dir, int vdir, int x, int y, int comp, void *data, int write_alpha, int scanline_pad, int expand_mono) 452 | { 453 | stbiw_uint32 zero = 0; 454 | int i,j, j_end; 455 | 456 | if (y <= 0) 457 | return; 458 | 459 | if (stbi__flip_vertically_on_write) 460 | vdir *= -1; 461 | 462 | if (vdir < 0) { 463 | j_end = -1; j = y-1; 464 | } else { 465 | j_end = y; j = 0; 466 | } 467 | 468 | for (; j != j_end; j += vdir) { 469 | for (i=0; i < x; ++i) { 470 | unsigned char *d = (unsigned char *) data + (j*x+i)*comp; 471 | stbiw__write_pixel(s, rgb_dir, comp, write_alpha, expand_mono, d); 472 | } 473 | stbiw__write_flush(s); 474 | s->func(s->context, &zero, scanline_pad); 475 | } 476 | } 477 | 478 | static int stbiw__outfile(stbi__write_context *s, int rgb_dir, int vdir, int x, int y, int comp, int expand_mono, void *data, int alpha, int pad, const char *fmt, ...) 479 | { 480 | if (y < 0 || x < 0) { 481 | return 0; 482 | } else { 483 | va_list v; 484 | va_start(v, fmt); 485 | stbiw__writefv(s, fmt, v); 486 | va_end(v); 487 | stbiw__write_pixels(s,rgb_dir,vdir,x,y,comp,data,alpha,pad, expand_mono); 488 | return 1; 489 | } 490 | } 491 | 492 | static int stbi_write_bmp_core(stbi__write_context *s, int x, int y, int comp, const void *data) 493 | { 494 | if (comp != 4) { 495 | // write RGB bitmap 496 | int pad = (-x*3) & 3; 497 | return stbiw__outfile(s,-1,-1,x,y,comp,1,(void *) data,0,pad, 498 | "11 4 22 4" "4 44 22 444444", 499 | 'B', 'M', 14+40+(x*3+pad)*y, 0,0, 14+40, // file header 500 | 40, x,y, 1,24, 0,0,0,0,0,0); // bitmap header 501 | } else { 502 | // RGBA bitmaps need a v4 header 503 | // use BI_BITFIELDS mode with 32bpp and alpha mask 504 | // (straight BI_RGB with alpha mask doesn't work in most readers) 505 | return stbiw__outfile(s,-1,-1,x,y,comp,1,(void *)data,1,0, 506 | "11 4 22 4" "4 44 22 444444 4444 4 444 444 444 444", 507 | 'B', 'M', 14+108+x*y*4, 0, 0, 14+108, // file header 508 | 108, x,y, 1,32, 3,0,0,0,0,0, 0xff0000,0xff00,0xff,0xff000000u, 0, 0,0,0, 0,0,0, 0,0,0, 0,0,0); // bitmap V4 header 509 | } 510 | } 511 | 512 | STBIWDEF int stbi_write_bmp_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data) 513 | { 514 | stbi__write_context s = { 0 }; 515 | stbi__start_write_callbacks(&s, func, context); 516 | return stbi_write_bmp_core(&s, x, y, comp, data); 517 | } 518 | 519 | #ifndef STBI_WRITE_NO_STDIO 520 | STBIWDEF int stbi_write_bmp(char const *filename, int x, int y, int comp, const void *data) 521 | { 522 | stbi__write_context s = { 0 }; 523 | if (stbi__start_write_file(&s,filename)) { 524 | int r = stbi_write_bmp_core(&s, x, y, comp, data); 525 | stbi__end_write_file(&s); 526 | return r; 527 | } else 528 | return 0; 529 | } 530 | #endif //!STBI_WRITE_NO_STDIO 531 | 532 | static int stbi_write_tga_core(stbi__write_context *s, int x, int y, int comp, void *data) 533 | { 534 | int has_alpha = (comp == 2 || comp == 4); 535 | int colorbytes = has_alpha ? comp-1 : comp; 536 | int format = colorbytes < 2 ? 3 : 2; // 3 color channels (RGB/RGBA) = 2, 1 color channel (Y/YA) = 3 537 | 538 | if (y < 0 || x < 0) 539 | return 0; 540 | 541 | if (!stbi_write_tga_with_rle) { 542 | return stbiw__outfile(s, -1, -1, x, y, comp, 0, (void *) data, has_alpha, 0, 543 | "111 221 2222 11", 0, 0, format, 0, 0, 0, 0, 0, x, y, (colorbytes + has_alpha) * 8, has_alpha * 8); 544 | } else { 545 | int i,j,k; 546 | int jend, jdir; 547 | 548 | stbiw__writef(s, "111 221 2222 11", 0,0,format+8, 0,0,0, 0,0,x,y, (colorbytes + has_alpha) * 8, has_alpha * 8); 549 | 550 | if (stbi__flip_vertically_on_write) { 551 | j = 0; 552 | jend = y; 553 | jdir = 1; 554 | } else { 555 | j = y-1; 556 | jend = -1; 557 | jdir = -1; 558 | } 559 | for (; j != jend; j += jdir) { 560 | unsigned char *row = (unsigned char *) data + j * x * comp; 561 | int len; 562 | 563 | for (i = 0; i < x; i += len) { 564 | unsigned char *begin = row + i * comp; 565 | int diff = 1; 566 | len = 1; 567 | 568 | if (i < x - 1) { 569 | ++len; 570 | diff = memcmp(begin, row + (i + 1) * comp, comp); 571 | if (diff) { 572 | const unsigned char *prev = begin; 573 | for (k = i + 2; k < x && len < 128; ++k) { 574 | if (memcmp(prev, row + k * comp, comp)) { 575 | prev += comp; 576 | ++len; 577 | } else { 578 | --len; 579 | break; 580 | } 581 | } 582 | } else { 583 | for (k = i + 2; k < x && len < 128; ++k) { 584 | if (!memcmp(begin, row + k * comp, comp)) { 585 | ++len; 586 | } else { 587 | break; 588 | } 589 | } 590 | } 591 | } 592 | 593 | if (diff) { 594 | unsigned char header = STBIW_UCHAR(len - 1); 595 | stbiw__write1(s, header); 596 | for (k = 0; k < len; ++k) { 597 | stbiw__write_pixel(s, -1, comp, has_alpha, 0, begin + k * comp); 598 | } 599 | } else { 600 | unsigned char header = STBIW_UCHAR(len - 129); 601 | stbiw__write1(s, header); 602 | stbiw__write_pixel(s, -1, comp, has_alpha, 0, begin); 603 | } 604 | } 605 | } 606 | stbiw__write_flush(s); 607 | } 608 | return 1; 609 | } 610 | 611 | STBIWDEF int stbi_write_tga_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data) 612 | { 613 | stbi__write_context s = { 0 }; 614 | stbi__start_write_callbacks(&s, func, context); 615 | return stbi_write_tga_core(&s, x, y, comp, (void *) data); 616 | } 617 | 618 | #ifndef STBI_WRITE_NO_STDIO 619 | STBIWDEF int stbi_write_tga(char const *filename, int x, int y, int comp, const void *data) 620 | { 621 | stbi__write_context s = { 0 }; 622 | if (stbi__start_write_file(&s,filename)) { 623 | int r = stbi_write_tga_core(&s, x, y, comp, (void *) data); 624 | stbi__end_write_file(&s); 625 | return r; 626 | } else 627 | return 0; 628 | } 629 | #endif 630 | 631 | // ************************************************************************************************* 632 | // Radiance RGBE HDR writer 633 | // by Baldur Karlsson 634 | 635 | #define stbiw__max(a, b) ((a) > (b) ? (a) : (b)) 636 | 637 | #ifndef STBI_WRITE_NO_STDIO 638 | 639 | static void stbiw__linear_to_rgbe(unsigned char *rgbe, float *linear) 640 | { 641 | int exponent; 642 | float maxcomp = stbiw__max(linear[0], stbiw__max(linear[1], linear[2])); 643 | 644 | if (maxcomp < 1e-32f) { 645 | rgbe[0] = rgbe[1] = rgbe[2] = rgbe[3] = 0; 646 | } else { 647 | float normalize = (float) frexp(maxcomp, &exponent) * 256.0f/maxcomp; 648 | 649 | rgbe[0] = (unsigned char)(linear[0] * normalize); 650 | rgbe[1] = (unsigned char)(linear[1] * normalize); 651 | rgbe[2] = (unsigned char)(linear[2] * normalize); 652 | rgbe[3] = (unsigned char)(exponent + 128); 653 | } 654 | } 655 | 656 | static void stbiw__write_run_data(stbi__write_context *s, int length, unsigned char databyte) 657 | { 658 | unsigned char lengthbyte = STBIW_UCHAR(length+128); 659 | STBIW_ASSERT(length+128 <= 255); 660 | s->func(s->context, &lengthbyte, 1); 661 | s->func(s->context, &databyte, 1); 662 | } 663 | 664 | static void stbiw__write_dump_data(stbi__write_context *s, int length, unsigned char *data) 665 | { 666 | unsigned char lengthbyte = STBIW_UCHAR(length); 667 | STBIW_ASSERT(length <= 128); // inconsistent with spec but consistent with official code 668 | s->func(s->context, &lengthbyte, 1); 669 | s->func(s->context, data, length); 670 | } 671 | 672 | static void stbiw__write_hdr_scanline(stbi__write_context *s, int width, int ncomp, unsigned char *scratch, float *scanline) 673 | { 674 | unsigned char scanlineheader[4] = { 2, 2, 0, 0 }; 675 | unsigned char rgbe[4]; 676 | float linear[3]; 677 | int x; 678 | 679 | scanlineheader[2] = (width&0xff00)>>8; 680 | scanlineheader[3] = (width&0x00ff); 681 | 682 | /* skip RLE for images too small or large */ 683 | if (width < 8 || width >= 32768) { 684 | for (x=0; x < width; x++) { 685 | switch (ncomp) { 686 | case 4: /* fallthrough */ 687 | case 3: linear[2] = scanline[x*ncomp + 2]; 688 | linear[1] = scanline[x*ncomp + 1]; 689 | linear[0] = scanline[x*ncomp + 0]; 690 | break; 691 | default: 692 | linear[0] = linear[1] = linear[2] = scanline[x*ncomp + 0]; 693 | break; 694 | } 695 | stbiw__linear_to_rgbe(rgbe, linear); 696 | s->func(s->context, rgbe, 4); 697 | } 698 | } else { 699 | int c,r; 700 | /* encode into scratch buffer */ 701 | for (x=0; x < width; x++) { 702 | switch(ncomp) { 703 | case 4: /* fallthrough */ 704 | case 3: linear[2] = scanline[x*ncomp + 2]; 705 | linear[1] = scanline[x*ncomp + 1]; 706 | linear[0] = scanline[x*ncomp + 0]; 707 | break; 708 | default: 709 | linear[0] = linear[1] = linear[2] = scanline[x*ncomp + 0]; 710 | break; 711 | } 712 | stbiw__linear_to_rgbe(rgbe, linear); 713 | scratch[x + width*0] = rgbe[0]; 714 | scratch[x + width*1] = rgbe[1]; 715 | scratch[x + width*2] = rgbe[2]; 716 | scratch[x + width*3] = rgbe[3]; 717 | } 718 | 719 | s->func(s->context, scanlineheader, 4); 720 | 721 | /* RLE each component separately */ 722 | for (c=0; c < 4; c++) { 723 | unsigned char *comp = &scratch[width*c]; 724 | 725 | x = 0; 726 | while (x < width) { 727 | // find first run 728 | r = x; 729 | while (r+2 < width) { 730 | if (comp[r] == comp[r+1] && comp[r] == comp[r+2]) 731 | break; 732 | ++r; 733 | } 734 | if (r+2 >= width) 735 | r = width; 736 | // dump up to first run 737 | while (x < r) { 738 | int len = r-x; 739 | if (len > 128) len = 128; 740 | stbiw__write_dump_data(s, len, &comp[x]); 741 | x += len; 742 | } 743 | // if there's a run, output it 744 | if (r+2 < width) { // same test as what we break out of in search loop, so only true if we break'd 745 | // find next byte after run 746 | while (r < width && comp[r] == comp[x]) 747 | ++r; 748 | // output run up to r 749 | while (x < r) { 750 | int len = r-x; 751 | if (len > 127) len = 127; 752 | stbiw__write_run_data(s, len, comp[x]); 753 | x += len; 754 | } 755 | } 756 | } 757 | } 758 | } 759 | } 760 | 761 | static int stbi_write_hdr_core(stbi__write_context *s, int x, int y, int comp, float *data) 762 | { 763 | if (y <= 0 || x <= 0 || data == NULL) 764 | return 0; 765 | else { 766 | // Each component is stored separately. Allocate scratch space for full output scanline. 767 | unsigned char *scratch = (unsigned char *) STBIW_MALLOC(x*4); 768 | int i, len; 769 | char buffer[128]; 770 | char header[] = "#?RADIANCE\n# Written by stb_image_write.h\nFORMAT=32-bit_rle_rgbe\n"; 771 | s->func(s->context, header, sizeof(header)-1); 772 | 773 | #ifdef __STDC_LIB_EXT1__ 774 | len = sprintf_s(buffer, sizeof(buffer), "EXPOSURE= 1.0000000000000\n\n-Y %d +X %d\n", y, x); 775 | #else 776 | len = sprintf(buffer, "EXPOSURE= 1.0000000000000\n\n-Y %d +X %d\n", y, x); 777 | #endif 778 | s->func(s->context, buffer, len); 779 | 780 | for(i=0; i < y; i++) 781 | stbiw__write_hdr_scanline(s, x, comp, scratch, data + comp*x*(stbi__flip_vertically_on_write ? y-1-i : i)); 782 | STBIW_FREE(scratch); 783 | return 1; 784 | } 785 | } 786 | 787 | STBIWDEF int stbi_write_hdr_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const float *data) 788 | { 789 | stbi__write_context s = { 0 }; 790 | stbi__start_write_callbacks(&s, func, context); 791 | return stbi_write_hdr_core(&s, x, y, comp, (float *) data); 792 | } 793 | 794 | STBIWDEF int stbi_write_hdr(char const *filename, int x, int y, int comp, const float *data) 795 | { 796 | stbi__write_context s = { 0 }; 797 | if (stbi__start_write_file(&s,filename)) { 798 | int r = stbi_write_hdr_core(&s, x, y, comp, (float *) data); 799 | stbi__end_write_file(&s); 800 | return r; 801 | } else 802 | return 0; 803 | } 804 | #endif // STBI_WRITE_NO_STDIO 805 | 806 | 807 | ////////////////////////////////////////////////////////////////////////////// 808 | // 809 | // PNG writer 810 | // 811 | 812 | #ifndef STBIW_ZLIB_COMPRESS 813 | // stretchy buffer; stbiw__sbpush() == vector<>::push_back() -- stbiw__sbcount() == vector<>::size() 814 | #define stbiw__sbraw(a) ((int *) (void *) (a) - 2) 815 | #define stbiw__sbm(a) stbiw__sbraw(a)[0] 816 | #define stbiw__sbn(a) stbiw__sbraw(a)[1] 817 | 818 | #define stbiw__sbneedgrow(a,n) ((a)==0 || stbiw__sbn(a)+n >= stbiw__sbm(a)) 819 | #define stbiw__sbmaybegrow(a,n) (stbiw__sbneedgrow(a,(n)) ? stbiw__sbgrow(a,n) : 0) 820 | #define stbiw__sbgrow(a,n) stbiw__sbgrowf((void **) &(a), (n), sizeof(*(a))) 821 | 822 | #define stbiw__sbpush(a, v) (stbiw__sbmaybegrow(a,1), (a)[stbiw__sbn(a)++] = (v)) 823 | #define stbiw__sbcount(a) ((a) ? stbiw__sbn(a) : 0) 824 | #define stbiw__sbfree(a) ((a) ? STBIW_FREE(stbiw__sbraw(a)),0 : 0) 825 | 826 | static void *stbiw__sbgrowf(void **arr, int increment, int itemsize) 827 | { 828 | int m = *arr ? 2*stbiw__sbm(*arr)+increment : increment+1; 829 | void *p = STBIW_REALLOC_SIZED(*arr ? stbiw__sbraw(*arr) : 0, *arr ? (stbiw__sbm(*arr)*itemsize + sizeof(int)*2) : 0, itemsize * m + sizeof(int)*2); 830 | STBIW_ASSERT(p); 831 | if (p) { 832 | if (!*arr) ((int *) p)[1] = 0; 833 | *arr = (void *) ((int *) p + 2); 834 | stbiw__sbm(*arr) = m; 835 | } 836 | return *arr; 837 | } 838 | 839 | static unsigned char *stbiw__zlib_flushf(unsigned char *data, unsigned int *bitbuffer, int *bitcount) 840 | { 841 | while (*bitcount >= 8) { 842 | stbiw__sbpush(data, STBIW_UCHAR(*bitbuffer)); 843 | *bitbuffer >>= 8; 844 | *bitcount -= 8; 845 | } 846 | return data; 847 | } 848 | 849 | static int stbiw__zlib_bitrev(int code, int codebits) 850 | { 851 | int res=0; 852 | while (codebits--) { 853 | res = (res << 1) | (code & 1); 854 | code >>= 1; 855 | } 856 | return res; 857 | } 858 | 859 | static unsigned int stbiw__zlib_countm(unsigned char *a, unsigned char *b, int limit) 860 | { 861 | int i; 862 | for (i=0; i < limit && i < 258; ++i) 863 | if (a[i] != b[i]) break; 864 | return i; 865 | } 866 | 867 | static unsigned int stbiw__zhash(unsigned char *data) 868 | { 869 | stbiw_uint32 hash = data[0] + (data[1] << 8) + (data[2] << 16); 870 | hash ^= hash << 3; 871 | hash += hash >> 5; 872 | hash ^= hash << 4; 873 | hash += hash >> 17; 874 | hash ^= hash << 25; 875 | hash += hash >> 6; 876 | return hash; 877 | } 878 | 879 | #define stbiw__zlib_flush() (out = stbiw__zlib_flushf(out, &bitbuf, &bitcount)) 880 | #define stbiw__zlib_add(code,codebits) \ 881 | (bitbuf |= (code) << bitcount, bitcount += (codebits), stbiw__zlib_flush()) 882 | #define stbiw__zlib_huffa(b,c) stbiw__zlib_add(stbiw__zlib_bitrev(b,c),c) 883 | // default huffman tables 884 | #define stbiw__zlib_huff1(n) stbiw__zlib_huffa(0x30 + (n), 8) 885 | #define stbiw__zlib_huff2(n) stbiw__zlib_huffa(0x190 + (n)-144, 9) 886 | #define stbiw__zlib_huff3(n) stbiw__zlib_huffa(0 + (n)-256,7) 887 | #define stbiw__zlib_huff4(n) stbiw__zlib_huffa(0xc0 + (n)-280,8) 888 | #define stbiw__zlib_huff(n) ((n) <= 143 ? stbiw__zlib_huff1(n) : (n) <= 255 ? stbiw__zlib_huff2(n) : (n) <= 279 ? stbiw__zlib_huff3(n) : stbiw__zlib_huff4(n)) 889 | #define stbiw__zlib_huffb(n) ((n) <= 143 ? stbiw__zlib_huff1(n) : stbiw__zlib_huff2(n)) 890 | 891 | #define stbiw__ZHASH 16384 892 | 893 | #endif // STBIW_ZLIB_COMPRESS 894 | 895 | STBIWDEF unsigned char * stbi_zlib_compress(unsigned char *data, int data_len, int *out_len, int quality) 896 | { 897 | #ifdef STBIW_ZLIB_COMPRESS 898 | // user provided a zlib compress implementation, use that 899 | return STBIW_ZLIB_COMPRESS(data, data_len, out_len, quality); 900 | #else // use builtin 901 | static unsigned short lengthc[] = { 3,4,5,6,7,8,9,10,11,13,15,17,19,23,27,31,35,43,51,59,67,83,99,115,131,163,195,227,258, 259 }; 902 | static unsigned char lengtheb[]= { 0,0,0,0,0,0,0, 0, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 5, 0 }; 903 | static unsigned short distc[] = { 1,2,3,4,5,7,9,13,17,25,33,49,65,97,129,193,257,385,513,769,1025,1537,2049,3073,4097,6145,8193,12289,16385,24577, 32768 }; 904 | static unsigned char disteb[] = { 0,0,0,0,1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,9,9,10,10,11,11,12,12,13,13 }; 905 | unsigned int bitbuf=0; 906 | int i,j, bitcount=0; 907 | unsigned char *out = NULL; 908 | unsigned char ***hash_table = (unsigned char***) STBIW_MALLOC(stbiw__ZHASH * sizeof(unsigned char**)); 909 | if (hash_table == NULL) 910 | return NULL; 911 | if (quality < 5) quality = 5; 912 | 913 | stbiw__sbpush(out, 0x78); // DEFLATE 32K window 914 | stbiw__sbpush(out, 0x5e); // FLEVEL = 1 915 | stbiw__zlib_add(1,1); // BFINAL = 1 916 | stbiw__zlib_add(1,2); // BTYPE = 1 -- fixed huffman 917 | 918 | for (i=0; i < stbiw__ZHASH; ++i) 919 | hash_table[i] = NULL; 920 | 921 | i=0; 922 | while (i < data_len-3) { 923 | // hash next 3 bytes of data to be compressed 924 | int h = stbiw__zhash(data+i)&(stbiw__ZHASH-1), best=3; 925 | unsigned char *bestloc = 0; 926 | unsigned char **hlist = hash_table[h]; 927 | int n = stbiw__sbcount(hlist); 928 | for (j=0; j < n; ++j) { 929 | if (hlist[j]-data > i-32768) { // if entry lies within window 930 | int d = stbiw__zlib_countm(hlist[j], data+i, data_len-i); 931 | if (d >= best) { best=d; bestloc=hlist[j]; } 932 | } 933 | } 934 | // when hash table entry is too long, delete half the entries 935 | if (hash_table[h] && stbiw__sbn(hash_table[h]) == 2*quality) { 936 | STBIW_MEMMOVE(hash_table[h], hash_table[h]+quality, sizeof(hash_table[h][0])*quality); 937 | stbiw__sbn(hash_table[h]) = quality; 938 | } 939 | stbiw__sbpush(hash_table[h],data+i); 940 | 941 | if (bestloc) { 942 | // "lazy matching" - check match at *next* byte, and if it's better, do cur byte as literal 943 | h = stbiw__zhash(data+i+1)&(stbiw__ZHASH-1); 944 | hlist = hash_table[h]; 945 | n = stbiw__sbcount(hlist); 946 | for (j=0; j < n; ++j) { 947 | if (hlist[j]-data > i-32767) { 948 | int e = stbiw__zlib_countm(hlist[j], data+i+1, data_len-i-1); 949 | if (e > best) { // if next match is better, bail on current match 950 | bestloc = NULL; 951 | break; 952 | } 953 | } 954 | } 955 | } 956 | 957 | if (bestloc) { 958 | int d = (int) (data+i - bestloc); // distance back 959 | STBIW_ASSERT(d <= 32767 && best <= 258); 960 | for (j=0; best > lengthc[j+1]-1; ++j); 961 | stbiw__zlib_huff(j+257); 962 | if (lengtheb[j]) stbiw__zlib_add(best - lengthc[j], lengtheb[j]); 963 | for (j=0; d > distc[j+1]-1; ++j); 964 | stbiw__zlib_add(stbiw__zlib_bitrev(j,5),5); 965 | if (disteb[j]) stbiw__zlib_add(d - distc[j], disteb[j]); 966 | i += best; 967 | } else { 968 | stbiw__zlib_huffb(data[i]); 969 | ++i; 970 | } 971 | } 972 | // write out final bytes 973 | for (;i < data_len; ++i) 974 | stbiw__zlib_huffb(data[i]); 975 | stbiw__zlib_huff(256); // end of block 976 | // pad with 0 bits to byte boundary 977 | while (bitcount) 978 | stbiw__zlib_add(0,1); 979 | 980 | for (i=0; i < stbiw__ZHASH; ++i) 981 | (void) stbiw__sbfree(hash_table[i]); 982 | STBIW_FREE(hash_table); 983 | 984 | // store uncompressed instead if compression was worse 985 | if (stbiw__sbn(out) > data_len + 2 + ((data_len+32766)/32767)*5) { 986 | stbiw__sbn(out) = 2; // truncate to DEFLATE 32K window and FLEVEL = 1 987 | for (j = 0; j < data_len;) { 988 | int blocklen = data_len - j; 989 | if (blocklen > 32767) blocklen = 32767; 990 | stbiw__sbpush(out, data_len - j == blocklen); // BFINAL = ?, BTYPE = 0 -- no compression 991 | stbiw__sbpush(out, STBIW_UCHAR(blocklen)); // LEN 992 | stbiw__sbpush(out, STBIW_UCHAR(blocklen >> 8)); 993 | stbiw__sbpush(out, STBIW_UCHAR(~blocklen)); // NLEN 994 | stbiw__sbpush(out, STBIW_UCHAR(~blocklen >> 8)); 995 | memcpy(out+stbiw__sbn(out), data+j, blocklen); 996 | stbiw__sbn(out) += blocklen; 997 | j += blocklen; 998 | } 999 | } 1000 | 1001 | { 1002 | // compute adler32 on input 1003 | unsigned int s1=1, s2=0; 1004 | int blocklen = (int) (data_len % 5552); 1005 | j=0; 1006 | while (j < data_len) { 1007 | for (i=0; i < blocklen; ++i) { s1 += data[j+i]; s2 += s1; } 1008 | s1 %= 65521; s2 %= 65521; 1009 | j += blocklen; 1010 | blocklen = 5552; 1011 | } 1012 | stbiw__sbpush(out, STBIW_UCHAR(s2 >> 8)); 1013 | stbiw__sbpush(out, STBIW_UCHAR(s2)); 1014 | stbiw__sbpush(out, STBIW_UCHAR(s1 >> 8)); 1015 | stbiw__sbpush(out, STBIW_UCHAR(s1)); 1016 | } 1017 | *out_len = stbiw__sbn(out); 1018 | // make returned pointer freeable 1019 | STBIW_MEMMOVE(stbiw__sbraw(out), out, *out_len); 1020 | return (unsigned char *) stbiw__sbraw(out); 1021 | #endif // STBIW_ZLIB_COMPRESS 1022 | } 1023 | 1024 | static unsigned int stbiw__crc32(unsigned char *buffer, int len) 1025 | { 1026 | #ifdef STBIW_CRC32 1027 | return STBIW_CRC32(buffer, len); 1028 | #else 1029 | static unsigned int crc_table[256] = 1030 | { 1031 | 0x00000000, 0x77073096, 0xEE0E612C, 0x990951BA, 0x076DC419, 0x706AF48F, 0xE963A535, 0x9E6495A3, 1032 | 0x0eDB8832, 0x79DCB8A4, 0xE0D5E91E, 0x97D2D988, 0x09B64C2B, 0x7EB17CBD, 0xE7B82D07, 0x90BF1D91, 1033 | 0x1DB71064, 0x6AB020F2, 0xF3B97148, 0x84BE41DE, 0x1ADAD47D, 0x6DDDE4EB, 0xF4D4B551, 0x83D385C7, 1034 | 0x136C9856, 0x646BA8C0, 0xFD62F97A, 0x8A65C9EC, 0x14015C4F, 0x63066CD9, 0xFA0F3D63, 0x8D080DF5, 1035 | 0x3B6E20C8, 0x4C69105E, 0xD56041E4, 0xA2677172, 0x3C03E4D1, 0x4B04D447, 0xD20D85FD, 0xA50AB56B, 1036 | 0x35B5A8FA, 0x42B2986C, 0xDBBBC9D6, 0xACBCF940, 0x32D86CE3, 0x45DF5C75, 0xDCD60DCF, 0xABD13D59, 1037 | 0x26D930AC, 0x51DE003A, 0xC8D75180, 0xBFD06116, 0x21B4F4B5, 0x56B3C423, 0xCFBA9599, 0xB8BDA50F, 1038 | 0x2802B89E, 0x5F058808, 0xC60CD9B2, 0xB10BE924, 0x2F6F7C87, 0x58684C11, 0xC1611DAB, 0xB6662D3D, 1039 | 0x76DC4190, 0x01DB7106, 0x98D220BC, 0xEFD5102A, 0x71B18589, 0x06B6B51F, 0x9FBFE4A5, 0xE8B8D433, 1040 | 0x7807C9A2, 0x0F00F934, 0x9609A88E, 0xE10E9818, 0x7F6A0DBB, 0x086D3D2D, 0x91646C97, 0xE6635C01, 1041 | 0x6B6B51F4, 0x1C6C6162, 0x856530D8, 0xF262004E, 0x6C0695ED, 0x1B01A57B, 0x8208F4C1, 0xF50FC457, 1042 | 0x65B0D9C6, 0x12B7E950, 0x8BBEB8EA, 0xFCB9887C, 0x62DD1DDF, 0x15DA2D49, 0x8CD37CF3, 0xFBD44C65, 1043 | 0x4DB26158, 0x3AB551CE, 0xA3BC0074, 0xD4BB30E2, 0x4ADFA541, 0x3DD895D7, 0xA4D1C46D, 0xD3D6F4FB, 1044 | 0x4369E96A, 0x346ED9FC, 0xAD678846, 0xDA60B8D0, 0x44042D73, 0x33031DE5, 0xAA0A4C5F, 0xDD0D7CC9, 1045 | 0x5005713C, 0x270241AA, 0xBE0B1010, 0xC90C2086, 0x5768B525, 0x206F85B3, 0xB966D409, 0xCE61E49F, 1046 | 0x5EDEF90E, 0x29D9C998, 0xB0D09822, 0xC7D7A8B4, 0x59B33D17, 0x2EB40D81, 0xB7BD5C3B, 0xC0BA6CAD, 1047 | 0xEDB88320, 0x9ABFB3B6, 0x03B6E20C, 0x74B1D29A, 0xEAD54739, 0x9DD277AF, 0x04DB2615, 0x73DC1683, 1048 | 0xE3630B12, 0x94643B84, 0x0D6D6A3E, 0x7A6A5AA8, 0xE40ECF0B, 0x9309FF9D, 0x0A00AE27, 0x7D079EB1, 1049 | 0xF00F9344, 0x8708A3D2, 0x1E01F268, 0x6906C2FE, 0xF762575D, 0x806567CB, 0x196C3671, 0x6E6B06E7, 1050 | 0xFED41B76, 0x89D32BE0, 0x10DA7A5A, 0x67DD4ACC, 0xF9B9DF6F, 0x8EBEEFF9, 0x17B7BE43, 0x60B08ED5, 1051 | 0xD6D6A3E8, 0xA1D1937E, 0x38D8C2C4, 0x4FDFF252, 0xD1BB67F1, 0xA6BC5767, 0x3FB506DD, 0x48B2364B, 1052 | 0xD80D2BDA, 0xAF0A1B4C, 0x36034AF6, 0x41047A60, 0xDF60EFC3, 0xA867DF55, 0x316E8EEF, 0x4669BE79, 1053 | 0xCB61B38C, 0xBC66831A, 0x256FD2A0, 0x5268E236, 0xCC0C7795, 0xBB0B4703, 0x220216B9, 0x5505262F, 1054 | 0xC5BA3BBE, 0xB2BD0B28, 0x2BB45A92, 0x5CB36A04, 0xC2D7FFA7, 0xB5D0CF31, 0x2CD99E8B, 0x5BDEAE1D, 1055 | 0x9B64C2B0, 0xEC63F226, 0x756AA39C, 0x026D930A, 0x9C0906A9, 0xEB0E363F, 0x72076785, 0x05005713, 1056 | 0x95BF4A82, 0xE2B87A14, 0x7BB12BAE, 0x0CB61B38, 0x92D28E9B, 0xE5D5BE0D, 0x7CDCEFB7, 0x0BDBDF21, 1057 | 0x86D3D2D4, 0xF1D4E242, 0x68DDB3F8, 0x1FDA836E, 0x81BE16CD, 0xF6B9265B, 0x6FB077E1, 0x18B74777, 1058 | 0x88085AE6, 0xFF0F6A70, 0x66063BCA, 0x11010B5C, 0x8F659EFF, 0xF862AE69, 0x616BFFD3, 0x166CCF45, 1059 | 0xA00AE278, 0xD70DD2EE, 0x4E048354, 0x3903B3C2, 0xA7672661, 0xD06016F7, 0x4969474D, 0x3E6E77DB, 1060 | 0xAED16A4A, 0xD9D65ADC, 0x40DF0B66, 0x37D83BF0, 0xA9BCAE53, 0xDEBB9EC5, 0x47B2CF7F, 0x30B5FFE9, 1061 | 0xBDBDF21C, 0xCABAC28A, 0x53B39330, 0x24B4A3A6, 0xBAD03605, 0xCDD70693, 0x54DE5729, 0x23D967BF, 1062 | 0xB3667A2E, 0xC4614AB8, 0x5D681B02, 0x2A6F2B94, 0xB40BBE37, 0xC30C8EA1, 0x5A05DF1B, 0x2D02EF8D 1063 | }; 1064 | 1065 | unsigned int crc = ~0u; 1066 | int i; 1067 | for (i=0; i < len; ++i) 1068 | crc = (crc >> 8) ^ crc_table[buffer[i] ^ (crc & 0xff)]; 1069 | return ~crc; 1070 | #endif 1071 | } 1072 | 1073 | #define stbiw__wpng4(o,a,b,c,d) ((o)[0]=STBIW_UCHAR(a),(o)[1]=STBIW_UCHAR(b),(o)[2]=STBIW_UCHAR(c),(o)[3]=STBIW_UCHAR(d),(o)+=4) 1074 | #define stbiw__wp32(data,v) stbiw__wpng4(data, (v)>>24,(v)>>16,(v)>>8,(v)); 1075 | #define stbiw__wptag(data,s) stbiw__wpng4(data, s[0],s[1],s[2],s[3]) 1076 | 1077 | static void stbiw__wpcrc(unsigned char **data, int len) 1078 | { 1079 | unsigned int crc = stbiw__crc32(*data - len - 4, len+4); 1080 | stbiw__wp32(*data, crc); 1081 | } 1082 | 1083 | static unsigned char stbiw__paeth(int a, int b, int c) 1084 | { 1085 | int p = a + b - c, pa = abs(p-a), pb = abs(p-b), pc = abs(p-c); 1086 | if (pa <= pb && pa <= pc) return STBIW_UCHAR(a); 1087 | if (pb <= pc) return STBIW_UCHAR(b); 1088 | return STBIW_UCHAR(c); 1089 | } 1090 | 1091 | // @OPTIMIZE: provide an option that always forces left-predict or paeth predict 1092 | static void stbiw__encode_png_line(unsigned char *pixels, int stride_bytes, int width, int height, int y, int n, int filter_type, signed char *line_buffer) 1093 | { 1094 | static int mapping[] = { 0,1,2,3,4 }; 1095 | static int firstmap[] = { 0,1,0,5,6 }; 1096 | int *mymap = (y != 0) ? mapping : firstmap; 1097 | int i; 1098 | int type = mymap[filter_type]; 1099 | unsigned char *z = pixels + stride_bytes * (stbi__flip_vertically_on_write ? height-1-y : y); 1100 | int signed_stride = stbi__flip_vertically_on_write ? -stride_bytes : stride_bytes; 1101 | 1102 | if (type==0) { 1103 | memcpy(line_buffer, z, width*n); 1104 | return; 1105 | } 1106 | 1107 | // first loop isn't optimized since it's just one pixel 1108 | for (i = 0; i < n; ++i) { 1109 | switch (type) { 1110 | case 1: line_buffer[i] = z[i]; break; 1111 | case 2: line_buffer[i] = z[i] - z[i-signed_stride]; break; 1112 | case 3: line_buffer[i] = z[i] - (z[i-signed_stride]>>1); break; 1113 | case 4: line_buffer[i] = (signed char) (z[i] - stbiw__paeth(0,z[i-signed_stride],0)); break; 1114 | case 5: line_buffer[i] = z[i]; break; 1115 | case 6: line_buffer[i] = z[i]; break; 1116 | } 1117 | } 1118 | switch (type) { 1119 | case 1: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - z[i-n]; break; 1120 | case 2: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - z[i-signed_stride]; break; 1121 | case 3: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - ((z[i-n] + z[i-signed_stride])>>1); break; 1122 | case 4: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - stbiw__paeth(z[i-n], z[i-signed_stride], z[i-signed_stride-n]); break; 1123 | case 5: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - (z[i-n]>>1); break; 1124 | case 6: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - stbiw__paeth(z[i-n], 0,0); break; 1125 | } 1126 | } 1127 | 1128 | STBIWDEF unsigned char *stbi_write_png_to_mem(const unsigned char *pixels, int stride_bytes, int x, int y, int n, int *out_len) 1129 | { 1130 | int force_filter = stbi_write_force_png_filter; 1131 | int ctype[5] = { -1, 0, 4, 2, 6 }; 1132 | unsigned char sig[8] = { 137,80,78,71,13,10,26,10 }; 1133 | unsigned char *out,*o, *filt, *zlib; 1134 | signed char *line_buffer; 1135 | int j,zlen; 1136 | 1137 | if (stride_bytes == 0) 1138 | stride_bytes = x * n; 1139 | 1140 | if (force_filter >= 5) { 1141 | force_filter = -1; 1142 | } 1143 | 1144 | filt = (unsigned char *) STBIW_MALLOC((x*n+1) * y); if (!filt) return 0; 1145 | line_buffer = (signed char *) STBIW_MALLOC(x * n); if (!line_buffer) { STBIW_FREE(filt); return 0; } 1146 | for (j=0; j < y; ++j) { 1147 | int filter_type; 1148 | if (force_filter > -1) { 1149 | filter_type = force_filter; 1150 | stbiw__encode_png_line((unsigned char*)(pixels), stride_bytes, x, y, j, n, force_filter, line_buffer); 1151 | } else { // Estimate the best filter by running through all of them: 1152 | int best_filter = 0, best_filter_val = 0x7fffffff, est, i; 1153 | for (filter_type = 0; filter_type < 5; filter_type++) { 1154 | stbiw__encode_png_line((unsigned char*)(pixels), stride_bytes, x, y, j, n, filter_type, line_buffer); 1155 | 1156 | // Estimate the entropy of the line using this filter; the less, the better. 1157 | est = 0; 1158 | for (i = 0; i < x*n; ++i) { 1159 | est += abs((signed char) line_buffer[i]); 1160 | } 1161 | if (est < best_filter_val) { 1162 | best_filter_val = est; 1163 | best_filter = filter_type; 1164 | } 1165 | } 1166 | if (filter_type != best_filter) { // If the last iteration already got us the best filter, don't redo it 1167 | stbiw__encode_png_line((unsigned char*)(pixels), stride_bytes, x, y, j, n, best_filter, line_buffer); 1168 | filter_type = best_filter; 1169 | } 1170 | } 1171 | // when we get here, filter_type contains the filter type, and line_buffer contains the data 1172 | filt[j*(x*n+1)] = (unsigned char) filter_type; 1173 | STBIW_MEMMOVE(filt+j*(x*n+1)+1, line_buffer, x*n); 1174 | } 1175 | STBIW_FREE(line_buffer); 1176 | zlib = stbi_zlib_compress(filt, y*( x*n+1), &zlen, stbi_write_png_compression_level); 1177 | STBIW_FREE(filt); 1178 | if (!zlib) return 0; 1179 | 1180 | // each tag requires 12 bytes of overhead 1181 | out = (unsigned char *) STBIW_MALLOC(8 + 12+13 + 12+zlen + 12); 1182 | if (!out) return 0; 1183 | *out_len = 8 + 12+13 + 12+zlen + 12; 1184 | 1185 | o=out; 1186 | STBIW_MEMMOVE(o,sig,8); o+= 8; 1187 | stbiw__wp32(o, 13); // header length 1188 | stbiw__wptag(o, "IHDR"); 1189 | stbiw__wp32(o, x); 1190 | stbiw__wp32(o, y); 1191 | *o++ = 8; 1192 | *o++ = STBIW_UCHAR(ctype[n]); 1193 | *o++ = 0; 1194 | *o++ = 0; 1195 | *o++ = 0; 1196 | stbiw__wpcrc(&o,13); 1197 | 1198 | stbiw__wp32(o, zlen); 1199 | stbiw__wptag(o, "IDAT"); 1200 | STBIW_MEMMOVE(o, zlib, zlen); 1201 | o += zlen; 1202 | STBIW_FREE(zlib); 1203 | stbiw__wpcrc(&o, zlen); 1204 | 1205 | stbiw__wp32(o,0); 1206 | stbiw__wptag(o, "IEND"); 1207 | stbiw__wpcrc(&o,0); 1208 | 1209 | STBIW_ASSERT(o == out + *out_len); 1210 | 1211 | return out; 1212 | } 1213 | 1214 | #ifndef STBI_WRITE_NO_STDIO 1215 | STBIWDEF int stbi_write_png(char const *filename, int x, int y, int comp, const void *data, int stride_bytes) 1216 | { 1217 | FILE *f; 1218 | int len; 1219 | unsigned char *png = stbi_write_png_to_mem((const unsigned char *) data, stride_bytes, x, y, comp, &len); 1220 | if (png == NULL) return 0; 1221 | 1222 | f = stbiw__fopen(filename, "wb"); 1223 | if (!f) { STBIW_FREE(png); return 0; } 1224 | fwrite(png, 1, len, f); 1225 | fclose(f); 1226 | STBIW_FREE(png); 1227 | return 1; 1228 | } 1229 | #endif 1230 | 1231 | STBIWDEF int stbi_write_png_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data, int stride_bytes) 1232 | { 1233 | int len; 1234 | unsigned char *png = stbi_write_png_to_mem((const unsigned char *) data, stride_bytes, x, y, comp, &len); 1235 | if (png == NULL) return 0; 1236 | func(context, png, len); 1237 | STBIW_FREE(png); 1238 | return 1; 1239 | } 1240 | 1241 | 1242 | /* *************************************************************************** 1243 | * 1244 | * JPEG writer 1245 | * 1246 | * This is based on Jon Olick's jo_jpeg.cpp: 1247 | * public domain Simple, Minimalistic JPEG writer - http://www.jonolick.com/code.html 1248 | */ 1249 | 1250 | static const unsigned char stbiw__jpg_ZigZag[] = { 0,1,5,6,14,15,27,28,2,4,7,13,16,26,29,42,3,8,12,17,25,30,41,43,9,11,18, 1251 | 24,31,40,44,53,10,19,23,32,39,45,52,54,20,22,33,38,46,51,55,60,21,34,37,47,50,56,59,61,35,36,48,49,57,58,62,63 }; 1252 | 1253 | static void stbiw__jpg_writeBits(stbi__write_context *s, int *bitBufP, int *bitCntP, const unsigned short *bs) { 1254 | int bitBuf = *bitBufP, bitCnt = *bitCntP; 1255 | bitCnt += bs[1]; 1256 | bitBuf |= bs[0] << (24 - bitCnt); 1257 | while(bitCnt >= 8) { 1258 | unsigned char c = (bitBuf >> 16) & 255; 1259 | stbiw__putc(s, c); 1260 | if(c == 255) { 1261 | stbiw__putc(s, 0); 1262 | } 1263 | bitBuf <<= 8; 1264 | bitCnt -= 8; 1265 | } 1266 | *bitBufP = bitBuf; 1267 | *bitCntP = bitCnt; 1268 | } 1269 | 1270 | static void stbiw__jpg_DCT(float *d0p, float *d1p, float *d2p, float *d3p, float *d4p, float *d5p, float *d6p, float *d7p) { 1271 | float d0 = *d0p, d1 = *d1p, d2 = *d2p, d3 = *d3p, d4 = *d4p, d5 = *d5p, d6 = *d6p, d7 = *d7p; 1272 | float z1, z2, z3, z4, z5, z11, z13; 1273 | 1274 | float tmp0 = d0 + d7; 1275 | float tmp7 = d0 - d7; 1276 | float tmp1 = d1 + d6; 1277 | float tmp6 = d1 - d6; 1278 | float tmp2 = d2 + d5; 1279 | float tmp5 = d2 - d5; 1280 | float tmp3 = d3 + d4; 1281 | float tmp4 = d3 - d4; 1282 | 1283 | // Even part 1284 | float tmp10 = tmp0 + tmp3; // phase 2 1285 | float tmp13 = tmp0 - tmp3; 1286 | float tmp11 = tmp1 + tmp2; 1287 | float tmp12 = tmp1 - tmp2; 1288 | 1289 | d0 = tmp10 + tmp11; // phase 3 1290 | d4 = tmp10 - tmp11; 1291 | 1292 | z1 = (tmp12 + tmp13) * 0.707106781f; // c4 1293 | d2 = tmp13 + z1; // phase 5 1294 | d6 = tmp13 - z1; 1295 | 1296 | // Odd part 1297 | tmp10 = tmp4 + tmp5; // phase 2 1298 | tmp11 = tmp5 + tmp6; 1299 | tmp12 = tmp6 + tmp7; 1300 | 1301 | // The rotator is modified from fig 4-8 to avoid extra negations. 1302 | z5 = (tmp10 - tmp12) * 0.382683433f; // c6 1303 | z2 = tmp10 * 0.541196100f + z5; // c2-c6 1304 | z4 = tmp12 * 1.306562965f + z5; // c2+c6 1305 | z3 = tmp11 * 0.707106781f; // c4 1306 | 1307 | z11 = tmp7 + z3; // phase 5 1308 | z13 = tmp7 - z3; 1309 | 1310 | *d5p = z13 + z2; // phase 6 1311 | *d3p = z13 - z2; 1312 | *d1p = z11 + z4; 1313 | *d7p = z11 - z4; 1314 | 1315 | *d0p = d0; *d2p = d2; *d4p = d4; *d6p = d6; 1316 | } 1317 | 1318 | static void stbiw__jpg_calcBits(int val, unsigned short bits[2]) { 1319 | int tmp1 = val < 0 ? -val : val; 1320 | val = val < 0 ? val-1 : val; 1321 | bits[1] = 1; 1322 | while(tmp1 >>= 1) { 1323 | ++bits[1]; 1324 | } 1325 | bits[0] = val & ((1<0)&&(DU[end0pos]==0); --end0pos) { 1368 | } 1369 | // end0pos = first element in reverse order !=0 1370 | if(end0pos == 0) { 1371 | stbiw__jpg_writeBits(s, bitBuf, bitCnt, EOB); 1372 | return DU[0]; 1373 | } 1374 | for(i = 1; i <= end0pos; ++i) { 1375 | int startpos = i; 1376 | int nrzeroes; 1377 | unsigned short bits[2]; 1378 | for (; DU[i]==0 && i<=end0pos; ++i) { 1379 | } 1380 | nrzeroes = i-startpos; 1381 | if ( nrzeroes >= 16 ) { 1382 | int lng = nrzeroes>>4; 1383 | int nrmarker; 1384 | for (nrmarker=1; nrmarker <= lng; ++nrmarker) 1385 | stbiw__jpg_writeBits(s, bitBuf, bitCnt, M16zeroes); 1386 | nrzeroes &= 15; 1387 | } 1388 | stbiw__jpg_calcBits(DU[i], bits); 1389 | stbiw__jpg_writeBits(s, bitBuf, bitCnt, HTAC[(nrzeroes<<4)+bits[1]]); 1390 | stbiw__jpg_writeBits(s, bitBuf, bitCnt, bits); 1391 | } 1392 | if(end0pos != 63) { 1393 | stbiw__jpg_writeBits(s, bitBuf, bitCnt, EOB); 1394 | } 1395 | return DU[0]; 1396 | } 1397 | 1398 | static int stbi_write_jpg_core(stbi__write_context *s, int width, int height, int comp, const void* data, int quality) { 1399 | // Constants that don't pollute global namespace 1400 | static const unsigned char std_dc_luminance_nrcodes[] = {0,0,1,5,1,1,1,1,1,1,0,0,0,0,0,0,0}; 1401 | static const unsigned char std_dc_luminance_values[] = {0,1,2,3,4,5,6,7,8,9,10,11}; 1402 | static const unsigned char std_ac_luminance_nrcodes[] = {0,0,2,1,3,3,2,4,3,5,5,4,4,0,0,1,0x7d}; 1403 | static const unsigned char std_ac_luminance_values[] = { 1404 | 0x01,0x02,0x03,0x00,0x04,0x11,0x05,0x12,0x21,0x31,0x41,0x06,0x13,0x51,0x61,0x07,0x22,0x71,0x14,0x32,0x81,0x91,0xa1,0x08, 1405 | 0x23,0x42,0xb1,0xc1,0x15,0x52,0xd1,0xf0,0x24,0x33,0x62,0x72,0x82,0x09,0x0a,0x16,0x17,0x18,0x19,0x1a,0x25,0x26,0x27,0x28, 1406 | 0x29,0x2a,0x34,0x35,0x36,0x37,0x38,0x39,0x3a,0x43,0x44,0x45,0x46,0x47,0x48,0x49,0x4a,0x53,0x54,0x55,0x56,0x57,0x58,0x59, 1407 | 0x5a,0x63,0x64,0x65,0x66,0x67,0x68,0x69,0x6a,0x73,0x74,0x75,0x76,0x77,0x78,0x79,0x7a,0x83,0x84,0x85,0x86,0x87,0x88,0x89, 1408 | 0x8a,0x92,0x93,0x94,0x95,0x96,0x97,0x98,0x99,0x9a,0xa2,0xa3,0xa4,0xa5,0xa6,0xa7,0xa8,0xa9,0xaa,0xb2,0xb3,0xb4,0xb5,0xb6, 1409 | 0xb7,0xb8,0xb9,0xba,0xc2,0xc3,0xc4,0xc5,0xc6,0xc7,0xc8,0xc9,0xca,0xd2,0xd3,0xd4,0xd5,0xd6,0xd7,0xd8,0xd9,0xda,0xe1,0xe2, 1410 | 0xe3,0xe4,0xe5,0xe6,0xe7,0xe8,0xe9,0xea,0xf1,0xf2,0xf3,0xf4,0xf5,0xf6,0xf7,0xf8,0xf9,0xfa 1411 | }; 1412 | static const unsigned char std_dc_chrominance_nrcodes[] = {0,0,3,1,1,1,1,1,1,1,1,1,0,0,0,0,0}; 1413 | static const unsigned char std_dc_chrominance_values[] = {0,1,2,3,4,5,6,7,8,9,10,11}; 1414 | static const unsigned char std_ac_chrominance_nrcodes[] = {0,0,2,1,2,4,4,3,4,7,5,4,4,0,1,2,0x77}; 1415 | static const unsigned char std_ac_chrominance_values[] = { 1416 | 0x00,0x01,0x02,0x03,0x11,0x04,0x05,0x21,0x31,0x06,0x12,0x41,0x51,0x07,0x61,0x71,0x13,0x22,0x32,0x81,0x08,0x14,0x42,0x91, 1417 | 0xa1,0xb1,0xc1,0x09,0x23,0x33,0x52,0xf0,0x15,0x62,0x72,0xd1,0x0a,0x16,0x24,0x34,0xe1,0x25,0xf1,0x17,0x18,0x19,0x1a,0x26, 1418 | 0x27,0x28,0x29,0x2a,0x35,0x36,0x37,0x38,0x39,0x3a,0x43,0x44,0x45,0x46,0x47,0x48,0x49,0x4a,0x53,0x54,0x55,0x56,0x57,0x58, 1419 | 0x59,0x5a,0x63,0x64,0x65,0x66,0x67,0x68,0x69,0x6a,0x73,0x74,0x75,0x76,0x77,0x78,0x79,0x7a,0x82,0x83,0x84,0x85,0x86,0x87, 1420 | 0x88,0x89,0x8a,0x92,0x93,0x94,0x95,0x96,0x97,0x98,0x99,0x9a,0xa2,0xa3,0xa4,0xa5,0xa6,0xa7,0xa8,0xa9,0xaa,0xb2,0xb3,0xb4, 1421 | 0xb5,0xb6,0xb7,0xb8,0xb9,0xba,0xc2,0xc3,0xc4,0xc5,0xc6,0xc7,0xc8,0xc9,0xca,0xd2,0xd3,0xd4,0xd5,0xd6,0xd7,0xd8,0xd9,0xda, 1422 | 0xe2,0xe3,0xe4,0xe5,0xe6,0xe7,0xe8,0xe9,0xea,0xf2,0xf3,0xf4,0xf5,0xf6,0xf7,0xf8,0xf9,0xfa 1423 | }; 1424 | // Huffman tables 1425 | static const unsigned short YDC_HT[256][2] = { {0,2},{2,3},{3,3},{4,3},{5,3},{6,3},{14,4},{30,5},{62,6},{126,7},{254,8},{510,9}}; 1426 | static const unsigned short UVDC_HT[256][2] = { {0,2},{1,2},{2,2},{6,3},{14,4},{30,5},{62,6},{126,7},{254,8},{510,9},{1022,10},{2046,11}}; 1427 | static const unsigned short YAC_HT[256][2] = { 1428 | {10,4},{0,2},{1,2},{4,3},{11,4},{26,5},{120,7},{248,8},{1014,10},{65410,16},{65411,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1429 | {12,4},{27,5},{121,7},{502,9},{2038,11},{65412,16},{65413,16},{65414,16},{65415,16},{65416,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1430 | {28,5},{249,8},{1015,10},{4084,12},{65417,16},{65418,16},{65419,16},{65420,16},{65421,16},{65422,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1431 | {58,6},{503,9},{4085,12},{65423,16},{65424,16},{65425,16},{65426,16},{65427,16},{65428,16},{65429,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1432 | {59,6},{1016,10},{65430,16},{65431,16},{65432,16},{65433,16},{65434,16},{65435,16},{65436,16},{65437,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1433 | {122,7},{2039,11},{65438,16},{65439,16},{65440,16},{65441,16},{65442,16},{65443,16},{65444,16},{65445,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1434 | {123,7},{4086,12},{65446,16},{65447,16},{65448,16},{65449,16},{65450,16},{65451,16},{65452,16},{65453,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1435 | {250,8},{4087,12},{65454,16},{65455,16},{65456,16},{65457,16},{65458,16},{65459,16},{65460,16},{65461,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1436 | {504,9},{32704,15},{65462,16},{65463,16},{65464,16},{65465,16},{65466,16},{65467,16},{65468,16},{65469,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1437 | {505,9},{65470,16},{65471,16},{65472,16},{65473,16},{65474,16},{65475,16},{65476,16},{65477,16},{65478,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1438 | {506,9},{65479,16},{65480,16},{65481,16},{65482,16},{65483,16},{65484,16},{65485,16},{65486,16},{65487,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1439 | {1017,10},{65488,16},{65489,16},{65490,16},{65491,16},{65492,16},{65493,16},{65494,16},{65495,16},{65496,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1440 | {1018,10},{65497,16},{65498,16},{65499,16},{65500,16},{65501,16},{65502,16},{65503,16},{65504,16},{65505,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1441 | {2040,11},{65506,16},{65507,16},{65508,16},{65509,16},{65510,16},{65511,16},{65512,16},{65513,16},{65514,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1442 | {65515,16},{65516,16},{65517,16},{65518,16},{65519,16},{65520,16},{65521,16},{65522,16},{65523,16},{65524,16},{0,0},{0,0},{0,0},{0,0},{0,0}, 1443 | {2041,11},{65525,16},{65526,16},{65527,16},{65528,16},{65529,16},{65530,16},{65531,16},{65532,16},{65533,16},{65534,16},{0,0},{0,0},{0,0},{0,0},{0,0} 1444 | }; 1445 | static const unsigned short UVAC_HT[256][2] = { 1446 | {0,2},{1,2},{4,3},{10,4},{24,5},{25,5},{56,6},{120,7},{500,9},{1014,10},{4084,12},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1447 | {11,4},{57,6},{246,8},{501,9},{2038,11},{4085,12},{65416,16},{65417,16},{65418,16},{65419,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1448 | {26,5},{247,8},{1015,10},{4086,12},{32706,15},{65420,16},{65421,16},{65422,16},{65423,16},{65424,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1449 | {27,5},{248,8},{1016,10},{4087,12},{65425,16},{65426,16},{65427,16},{65428,16},{65429,16},{65430,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1450 | {58,6},{502,9},{65431,16},{65432,16},{65433,16},{65434,16},{65435,16},{65436,16},{65437,16},{65438,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1451 | {59,6},{1017,10},{65439,16},{65440,16},{65441,16},{65442,16},{65443,16},{65444,16},{65445,16},{65446,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1452 | {121,7},{2039,11},{65447,16},{65448,16},{65449,16},{65450,16},{65451,16},{65452,16},{65453,16},{65454,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1453 | {122,7},{2040,11},{65455,16},{65456,16},{65457,16},{65458,16},{65459,16},{65460,16},{65461,16},{65462,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1454 | {249,8},{65463,16},{65464,16},{65465,16},{65466,16},{65467,16},{65468,16},{65469,16},{65470,16},{65471,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1455 | {503,9},{65472,16},{65473,16},{65474,16},{65475,16},{65476,16},{65477,16},{65478,16},{65479,16},{65480,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1456 | {504,9},{65481,16},{65482,16},{65483,16},{65484,16},{65485,16},{65486,16},{65487,16},{65488,16},{65489,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1457 | {505,9},{65490,16},{65491,16},{65492,16},{65493,16},{65494,16},{65495,16},{65496,16},{65497,16},{65498,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1458 | {506,9},{65499,16},{65500,16},{65501,16},{65502,16},{65503,16},{65504,16},{65505,16},{65506,16},{65507,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1459 | {2041,11},{65508,16},{65509,16},{65510,16},{65511,16},{65512,16},{65513,16},{65514,16},{65515,16},{65516,16},{0,0},{0,0},{0,0},{0,0},{0,0},{0,0}, 1460 | {16352,14},{65517,16},{65518,16},{65519,16},{65520,16},{65521,16},{65522,16},{65523,16},{65524,16},{65525,16},{0,0},{0,0},{0,0},{0,0},{0,0}, 1461 | {1018,10},{32707,15},{65526,16},{65527,16},{65528,16},{65529,16},{65530,16},{65531,16},{65532,16},{65533,16},{65534,16},{0,0},{0,0},{0,0},{0,0},{0,0} 1462 | }; 1463 | static const int YQT[] = {16,11,10,16,24,40,51,61,12,12,14,19,26,58,60,55,14,13,16,24,40,57,69,56,14,17,22,29,51,87,80,62,18,22, 1464 | 37,56,68,109,103,77,24,35,55,64,81,104,113,92,49,64,78,87,103,121,120,101,72,92,95,98,112,100,103,99}; 1465 | static const int UVQT[] = {17,18,24,47,99,99,99,99,18,21,26,66,99,99,99,99,24,26,56,99,99,99,99,99,47,66,99,99,99,99,99,99, 1466 | 99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99,99}; 1467 | static const float aasf[] = { 1.0f * 2.828427125f, 1.387039845f * 2.828427125f, 1.306562965f * 2.828427125f, 1.175875602f * 2.828427125f, 1468 | 1.0f * 2.828427125f, 0.785694958f * 2.828427125f, 0.541196100f * 2.828427125f, 0.275899379f * 2.828427125f }; 1469 | 1470 | int row, col, i, k, subsample; 1471 | float fdtbl_Y[64], fdtbl_UV[64]; 1472 | unsigned char YTable[64], UVTable[64]; 1473 | 1474 | if(!data || !width || !height || comp > 4 || comp < 1) { 1475 | return 0; 1476 | } 1477 | 1478 | quality = quality ? quality : 90; 1479 | subsample = quality <= 90 ? 1 : 0; 1480 | quality = quality < 1 ? 1 : quality > 100 ? 100 : quality; 1481 | quality = quality < 50 ? 5000 / quality : 200 - quality * 2; 1482 | 1483 | for(i = 0; i < 64; ++i) { 1484 | int uvti, yti = (YQT[i]*quality+50)/100; 1485 | YTable[stbiw__jpg_ZigZag[i]] = (unsigned char) (yti < 1 ? 1 : yti > 255 ? 255 : yti); 1486 | uvti = (UVQT[i]*quality+50)/100; 1487 | UVTable[stbiw__jpg_ZigZag[i]] = (unsigned char) (uvti < 1 ? 1 : uvti > 255 ? 255 : uvti); 1488 | } 1489 | 1490 | for(row = 0, k = 0; row < 8; ++row) { 1491 | for(col = 0; col < 8; ++col, ++k) { 1492 | fdtbl_Y[k] = 1 / (YTable [stbiw__jpg_ZigZag[k]] * aasf[row] * aasf[col]); 1493 | fdtbl_UV[k] = 1 / (UVTable[stbiw__jpg_ZigZag[k]] * aasf[row] * aasf[col]); 1494 | } 1495 | } 1496 | 1497 | // Write Headers 1498 | { 1499 | static const unsigned char head0[] = { 0xFF,0xD8,0xFF,0xE0,0,0x10,'J','F','I','F',0,1,1,0,0,1,0,1,0,0,0xFF,0xDB,0,0x84,0 }; 1500 | static const unsigned char head2[] = { 0xFF,0xDA,0,0xC,3,1,0,2,0x11,3,0x11,0,0x3F,0 }; 1501 | const unsigned char head1[] = { 0xFF,0xC0,0,0x11,8,(unsigned char)(height>>8),STBIW_UCHAR(height),(unsigned char)(width>>8),STBIW_UCHAR(width), 1502 | 3,1,(unsigned char)(subsample?0x22:0x11),0,2,0x11,1,3,0x11,1,0xFF,0xC4,0x01,0xA2,0 }; 1503 | s->func(s->context, (void*)head0, sizeof(head0)); 1504 | s->func(s->context, (void*)YTable, sizeof(YTable)); 1505 | stbiw__putc(s, 1); 1506 | s->func(s->context, UVTable, sizeof(UVTable)); 1507 | s->func(s->context, (void*)head1, sizeof(head1)); 1508 | s->func(s->context, (void*)(std_dc_luminance_nrcodes+1), sizeof(std_dc_luminance_nrcodes)-1); 1509 | s->func(s->context, (void*)std_dc_luminance_values, sizeof(std_dc_luminance_values)); 1510 | stbiw__putc(s, 0x10); // HTYACinfo 1511 | s->func(s->context, (void*)(std_ac_luminance_nrcodes+1), sizeof(std_ac_luminance_nrcodes)-1); 1512 | s->func(s->context, (void*)std_ac_luminance_values, sizeof(std_ac_luminance_values)); 1513 | stbiw__putc(s, 1); // HTUDCinfo 1514 | s->func(s->context, (void*)(std_dc_chrominance_nrcodes+1), sizeof(std_dc_chrominance_nrcodes)-1); 1515 | s->func(s->context, (void*)std_dc_chrominance_values, sizeof(std_dc_chrominance_values)); 1516 | stbiw__putc(s, 0x11); // HTUACinfo 1517 | s->func(s->context, (void*)(std_ac_chrominance_nrcodes+1), sizeof(std_ac_chrominance_nrcodes)-1); 1518 | s->func(s->context, (void*)std_ac_chrominance_values, sizeof(std_ac_chrominance_values)); 1519 | s->func(s->context, (void*)head2, sizeof(head2)); 1520 | } 1521 | 1522 | // Encode 8x8 macroblocks 1523 | { 1524 | static const unsigned short fillBits[] = {0x7F, 7}; 1525 | int DCY=0, DCU=0, DCV=0; 1526 | int bitBuf=0, bitCnt=0; 1527 | // comp == 2 is grey+alpha (alpha is ignored) 1528 | int ofsG = comp > 2 ? 1 : 0, ofsB = comp > 2 ? 2 : 0; 1529 | const unsigned char *dataR = (const unsigned char *)data; 1530 | const unsigned char *dataG = dataR + ofsG; 1531 | const unsigned char *dataB = dataR + ofsB; 1532 | int x, y, pos; 1533 | if(subsample) { 1534 | for(y = 0; y < height; y += 16) { 1535 | for(x = 0; x < width; x += 16) { 1536 | float Y[256], U[256], V[256]; 1537 | for(row = y, pos = 0; row < y+16; ++row) { 1538 | // row >= height => use last input row 1539 | int clamped_row = (row < height) ? row : height - 1; 1540 | int base_p = (stbi__flip_vertically_on_write ? (height-1-clamped_row) : clamped_row)*width*comp; 1541 | for(col = x; col < x+16; ++col, ++pos) { 1542 | // if col >= width => use pixel from last input column 1543 | int p = base_p + ((col < width) ? col : (width-1))*comp; 1544 | float r = dataR[p], g = dataG[p], b = dataB[p]; 1545 | Y[pos]= +0.29900f*r + 0.58700f*g + 0.11400f*b - 128; 1546 | U[pos]= -0.16874f*r - 0.33126f*g + 0.50000f*b; 1547 | V[pos]= +0.50000f*r - 0.41869f*g - 0.08131f*b; 1548 | } 1549 | } 1550 | DCY = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, Y+0, 16, fdtbl_Y, DCY, YDC_HT, YAC_HT); 1551 | DCY = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, Y+8, 16, fdtbl_Y, DCY, YDC_HT, YAC_HT); 1552 | DCY = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, Y+128, 16, fdtbl_Y, DCY, YDC_HT, YAC_HT); 1553 | DCY = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, Y+136, 16, fdtbl_Y, DCY, YDC_HT, YAC_HT); 1554 | 1555 | // subsample U,V 1556 | { 1557 | float subU[64], subV[64]; 1558 | int yy, xx; 1559 | for(yy = 0, pos = 0; yy < 8; ++yy) { 1560 | for(xx = 0; xx < 8; ++xx, ++pos) { 1561 | int j = yy*32+xx*2; 1562 | subU[pos] = (U[j+0] + U[j+1] + U[j+16] + U[j+17]) * 0.25f; 1563 | subV[pos] = (V[j+0] + V[j+1] + V[j+16] + V[j+17]) * 0.25f; 1564 | } 1565 | } 1566 | DCU = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, subU, 8, fdtbl_UV, DCU, UVDC_HT, UVAC_HT); 1567 | DCV = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, subV, 8, fdtbl_UV, DCV, UVDC_HT, UVAC_HT); 1568 | } 1569 | } 1570 | } 1571 | } else { 1572 | for(y = 0; y < height; y += 8) { 1573 | for(x = 0; x < width; x += 8) { 1574 | float Y[64], U[64], V[64]; 1575 | for(row = y, pos = 0; row < y+8; ++row) { 1576 | // row >= height => use last input row 1577 | int clamped_row = (row < height) ? row : height - 1; 1578 | int base_p = (stbi__flip_vertically_on_write ? (height-1-clamped_row) : clamped_row)*width*comp; 1579 | for(col = x; col < x+8; ++col, ++pos) { 1580 | // if col >= width => use pixel from last input column 1581 | int p = base_p + ((col < width) ? col : (width-1))*comp; 1582 | float r = dataR[p], g = dataG[p], b = dataB[p]; 1583 | Y[pos]= +0.29900f*r + 0.58700f*g + 0.11400f*b - 128; 1584 | U[pos]= -0.16874f*r - 0.33126f*g + 0.50000f*b; 1585 | V[pos]= +0.50000f*r - 0.41869f*g - 0.08131f*b; 1586 | } 1587 | } 1588 | 1589 | DCY = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, Y, 8, fdtbl_Y, DCY, YDC_HT, YAC_HT); 1590 | DCU = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, U, 8, fdtbl_UV, DCU, UVDC_HT, UVAC_HT); 1591 | DCV = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, V, 8, fdtbl_UV, DCV, UVDC_HT, UVAC_HT); 1592 | } 1593 | } 1594 | } 1595 | 1596 | // Do the bit alignment of the EOI marker 1597 | stbiw__jpg_writeBits(s, &bitBuf, &bitCnt, fillBits); 1598 | } 1599 | 1600 | // EOI 1601 | stbiw__putc(s, 0xFF); 1602 | stbiw__putc(s, 0xD9); 1603 | 1604 | return 1; 1605 | } 1606 | 1607 | STBIWDEF int stbi_write_jpg_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data, int quality) 1608 | { 1609 | stbi__write_context s = { 0 }; 1610 | stbi__start_write_callbacks(&s, func, context); 1611 | return stbi_write_jpg_core(&s, x, y, comp, (void *) data, quality); 1612 | } 1613 | 1614 | 1615 | #ifndef STBI_WRITE_NO_STDIO 1616 | STBIWDEF int stbi_write_jpg(char const *filename, int x, int y, int comp, const void *data, int quality) 1617 | { 1618 | stbi__write_context s = { 0 }; 1619 | if (stbi__start_write_file(&s,filename)) { 1620 | int r = stbi_write_jpg_core(&s, x, y, comp, data, quality); 1621 | stbi__end_write_file(&s); 1622 | return r; 1623 | } else 1624 | return 0; 1625 | } 1626 | #endif 1627 | 1628 | #endif // STB_IMAGE_WRITE_IMPLEMENTATION 1629 | 1630 | /* Revision history 1631 | 1.16 (2021-07-11) 1632 | make Deflate code emit uncompressed blocks when it would otherwise expand 1633 | support writing BMPs with alpha channel 1634 | 1.15 (2020-07-13) unknown 1635 | 1.14 (2020-02-02) updated JPEG writer to downsample chroma channels 1636 | 1.13 1637 | 1.12 1638 | 1.11 (2019-08-11) 1639 | 1640 | 1.10 (2019-02-07) 1641 | support utf8 filenames in Windows; fix warnings and platform ifdefs 1642 | 1.09 (2018-02-11) 1643 | fix typo in zlib quality API, improve STB_I_W_STATIC in C++ 1644 | 1.08 (2018-01-29) 1645 | add stbi__flip_vertically_on_write, external zlib, zlib quality, choose PNG filter 1646 | 1.07 (2017-07-24) 1647 | doc fix 1648 | 1.06 (2017-07-23) 1649 | writing JPEG (using Jon Olick's code) 1650 | 1.05 ??? 1651 | 1.04 (2017-03-03) 1652 | monochrome BMP expansion 1653 | 1.03 ??? 1654 | 1.02 (2016-04-02) 1655 | avoid allocating large structures on the stack 1656 | 1.01 (2016-01-16) 1657 | STBIW_REALLOC_SIZED: support allocators with no realloc support 1658 | avoid race-condition in crc initialization 1659 | minor compile issues 1660 | 1.00 (2015-09-14) 1661 | installable file IO function 1662 | 0.99 (2015-09-13) 1663 | warning fixes; TGA rle support 1664 | 0.98 (2015-04-08) 1665 | added STBIW_MALLOC, STBIW_ASSERT etc 1666 | 0.97 (2015-01-18) 1667 | fixed HDR asserts, rewrote HDR rle logic 1668 | 0.96 (2015-01-17) 1669 | add HDR output 1670 | fix monochrome BMP 1671 | 0.95 (2014-08-17) 1672 | add monochrome TGA output 1673 | 0.94 (2014-05-31) 1674 | rename private functions to avoid conflicts with stb_image.h 1675 | 0.93 (2014-05-27) 1676 | warning fixes 1677 | 0.92 (2010-08-01) 1678 | casts to unsigned char to fix warnings 1679 | 0.91 (2010-07-17) 1680 | first public release 1681 | 0.90 first internal release 1682 | */ 1683 | 1684 | /* 1685 | ------------------------------------------------------------------------------ 1686 | This software is available under 2 licenses -- choose whichever you prefer. 1687 | ------------------------------------------------------------------------------ 1688 | ALTERNATIVE A - MIT License 1689 | Copyright (c) 2017 Sean Barrett 1690 | Permission is hereby granted, free of charge, to any person obtaining a copy of 1691 | this software and associated documentation files (the "Software"), to deal in 1692 | the Software without restriction, including without limitation the rights to 1693 | use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies 1694 | of the Software, and to permit persons to whom the Software is furnished to do 1695 | so, subject to the following conditions: 1696 | The above copyright notice and this permission notice shall be included in all 1697 | copies or substantial portions of the Software. 1698 | THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 1699 | IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 1700 | FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE 1701 | AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER 1702 | LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, 1703 | OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE 1704 | SOFTWARE. 1705 | ------------------------------------------------------------------------------ 1706 | ALTERNATIVE B - Public Domain (www.unlicense.org) 1707 | This is free and unencumbered software released into the public domain. 1708 | Anyone is free to copy, modify, publish, use, compile, sell, or distribute this 1709 | software, either in source code form or as a compiled binary, for any purpose, 1710 | commercial or non-commercial, and by any means. 1711 | In jurisdictions that recognize copyright laws, the author or authors of this 1712 | software dedicate any and all copyright interest in the software to the public 1713 | domain. We make this dedication for the benefit of the public at large and to 1714 | the detriment of our heirs and successors. We intend this dedication to be an 1715 | overt act of relinquishment in perpetuity of all present and future rights to 1716 | this software under copyright law. 1717 | THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR 1718 | IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, 1719 | FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE 1720 | AUTHORS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN 1721 | ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION 1722 | WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. 1723 | ------------------------------------------------------------------------------ 1724 | */ 1725 | -------------------------------------------------------------------------------- /test/impl.c: -------------------------------------------------------------------------------- 1 | 2 | #define NOI_IMAGE_IMPLEMENTATION 3 | #include 4 | 5 | #define STB_IMAGE_IMPLEMENTATION 6 | #include 7 | 8 | #define STB_IMAGE_WRITE_IMPLEMENTATION 9 | #include 10 | -------------------------------------------------------------------------------- /test/main.c: -------------------------------------------------------------------------------- 1 | #ifdef _MSC_VER 2 | #define _CRT_SECURE_NO_WARNINGS 3 | #endif 4 | 5 | #include 6 | #include 7 | #include 8 | 9 | #include 10 | 11 | #ifdef _MSC_VER 12 | #define WIN32_LEAN_AND_MEAN 13 | #include 14 | int64_t ref_time_ticks () { 15 | LARGE_INTEGER t0; 16 | QueryPerformanceCounter(&t0); 17 | return t0.QuadPart; 18 | } 19 | int get_time_usec ( int64_t reft ) { 20 | int64_t t0 = ref_time_ticks(); 21 | LARGE_INTEGER freq; 22 | QueryPerformanceFrequency(&freq); 23 | return (int)(((t0-reft)*1000000) / freq.QuadPart); 24 | } 25 | #elif __linux__ 26 | #include 27 | const uint64_t NSEC_IN_SEC = 1000000000LL; 28 | int64_t ref_time_ticks () { 29 | struct timespec ts; 30 | clock_gettime(CLOCK_MONOTONIC, &ts); 31 | return ts.tv_sec * NSEC_IN_SEC + ts.tv_nsec; 32 | } 33 | int get_time_usec ( int64_t reft ) { 34 | return (int) ((ref_time_ticks() - reft) / (NSEC_IN_SEC/1000000)); 35 | } 36 | #else // osx 37 | #include 38 | #include 39 | int64_t ref_time_ticks() { 40 | return mach_absolute_time(); 41 | } 42 | int get_time_usec ( int64_t reft ) { 43 | int64_t relt = ref_time_ticks() - reft; 44 | mach_timebase_info_data_t s_timebase_info; 45 | mach_timebase_info(&s_timebase_info); 46 | return relt * s_timebase_info.numer/s_timebase_info.denom/1000; 47 | } 48 | #endif 49 | 50 | void print_use ( void ) { 51 | printf( 52 | "NotOk image compression\n" 53 | " noi -b path_to_images/* dest_folder batch compress and output PSNR\n" 54 | " noi -c source_image dest_noi {profile} compress\n" 55 | " noi -pc source_image dest_noi {profile} compress and output stats\n" 56 | " noi -d source_noi dest_png decompress\n" 57 | " noi -pd source_noi dest_png profile decompression\n" 58 | " noi -stbjpg source_image dest_jpg profile stbimage jpeg decoder\n" 59 | ); 60 | } 61 | 62 | void rgb2yuv ( double R, double G, double B, double * Y, double * U, double * V ) { 63 | *Y = 0.257 * R + 0.504 * G + 0.098 * B + 16; 64 | *U = -0.148 * R - 0.291 * G + 0.439 * B + 128; 65 | *V = 0.439 * R - 0.368 * G - 0.071 * B + 128; 66 | } 67 | 68 | void psnr_greyscale ( uint8_t * oi, uint8_t * ci, int npixels, double * PSNR ) { 69 | double mse = 0.; 70 | double mse_yuv = 0.; 71 | double npix = npixels; 72 | while ( npixels-- ) { 73 | int R = oi[0]; 74 | int G = oi[1]; 75 | int B = oi[2]; 76 | int Y = (int)(0.299*R + 0.587*G + 0.114*B); 77 | double dg = ((double)Y) - ((double)ci[1]); 78 | mse += dg*dg; 79 | oi += 4; 80 | ci += 4; 81 | } 82 | mse /= npix; 83 | double D = 255.; 84 | *PSNR = (10 * log10((D*D) / mse)); 85 | } 86 | 87 | void psnr ( uint8_t * oi, uint8_t * ci, int npixels, double * PSNR, double * PSNR_YUV ) { 88 | double mse = 0.; 89 | double mse_yuv = 0.; 90 | double npix = npixels; 91 | while ( npixels-- ) { 92 | double dr = ((double)oi[0]) - ((double)ci[0]); 93 | double dg = ((double)oi[1]) - ((double)ci[1]); 94 | double db = ((double)oi[2]) - ((double)ci[2]); 95 | mse += dr*dr + dg*dg + db*db; 96 | double Y1,U1,V1,Y2,U2,V2; 97 | rgb2yuv(oi[0],oi[1],oi[2],&Y1,&U1,&V1); 98 | rgb2yuv(ci[0],ci[1],ci[2],&Y2,&U2,&V2); 99 | double dy = Y1 - Y2; 100 | double du = U1 - U2; 101 | double dv = V1 - V2; 102 | mse_yuv += dy*dy + du*du + dv*dv; 103 | oi += 4; 104 | ci += 4; 105 | } 106 | mse /= npix; 107 | mse_yuv /= npix; 108 | double D = 255.; 109 | *PSNR = (10 * log10((3*D*D) / mse)); 110 | *PSNR_YUV = (10 * log10((3*D*D) / mse_yuv)); 111 | } 112 | 113 | void * load_file ( const char * fileName, int * fsize ) { 114 | FILE * f = fopen(fileName,"rb"); 115 | if ( !f ) return NULL; 116 | fseek(f,0,SEEK_END); 117 | int csize = ftell(f); 118 | fseek(f,0,SEEK_SET); 119 | void * cbytes = malloc ( csize ); 120 | int cdsize = fread(cbytes, 1, csize, f); 121 | if ( csize != cdsize ) { 122 | free(cbytes); 123 | return NULL; 124 | } 125 | fclose(f); 126 | if ( fsize ) *fsize = csize; 127 | return cbytes; 128 | } 129 | 130 | int get_profile ( const char * arg ) { 131 | int profile = NOI_YUV_16_1_1; 132 | if ( strcmp(arg,"4_1_1")==0 ) profile = NOI_YUV_4_1_1; 133 | else if ( strcmp(arg,"2_1_1")==0 ) profile = NOI_YUV_2_1_1; 134 | else if ( strcmp(arg,"16_1_1")==0 ) profile = NOI_YUV_16_1_1; 135 | else if ( strcmp(arg,"1_1_1")==0 ) profile = NOI_RGB_1_1_1; 136 | else if ( strcmp(arg,"1_1_0")==0 ) profile = NOI_RG_1_1_0; 137 | else if ( strcmp(arg,"1_0_0")==0 ) profile = NOI_Y_1_0_0; 138 | else { 139 | printf("unsupported profile %s\n", arg); 140 | exit(-9); 141 | } 142 | return profile; 143 | } 144 | 145 | #ifdef _WIN32 146 | #include "dirent/dirent.h" 147 | #else 148 | #include 149 | #endif 150 | 151 | int main(int argc, char** argv) { 152 | if ( !(argc==4 || argc==5) ) { 153 | print_use(); 154 | return -1; 155 | } 156 | if ( strcmp(argv[1],"-b")==0 ) { 157 | int profile = argc==5 ? get_profile(argv[4]) : NOI_YUV_16_1_1; 158 | int bx, by; 159 | noi_profile_block_size(profile, &bx, &by); 160 | struct dirent **namelist; 161 | int n = scandir(argv[2], &namelist, 0, alphasort); 162 | if ( n>=0 ) { 163 | printf("IMAGE, PSNR, PSNR_YUV, PROFILE, W, H\n"); 164 | while( n-- ) { 165 | size_t sl = strlen(namelist[n]->d_name); 166 | if ( sl>4 && strcmp(namelist[n]->d_name+sl-4,".png")==0 ) { 167 | char img_name[1024], test_name[1024]; 168 | snprintf(img_name, sizeof(img_name), "%s\\%s", argv[2], namelist[n]->d_name); 169 | snprintf(test_name, sizeof(test_name), "%s\\%s", argv[3], namelist[n]->d_name); 170 | int w, h; 171 | uint8_t * pixels = stbi_load(img_name, &w, &h, NULL, 4); 172 | if ( pixels!=NULL ) { 173 | if ( (w%bx) || (h%by) ) { 174 | printf("\"%s\",\"skipped\",\"skipped\",%i,%i\n", namelist[n]->d_name, w, h ); 175 | } else { 176 | int csize = 0; 177 | void * cbytes = noi_compress(pixels, w, h, &csize, profile); 178 | uint8_t * cpixels = noi_decompress(cbytes, NULL, NULL, NULL); 179 | double PSNR = 0., PSNR_YUV = 0.; 180 | if ( profile==NOI_Y_1_0_0 ) { 181 | psnr_greyscale(pixels, cpixels, w*h, &PSNR); 182 | PSNR_YUV = PSNR; 183 | } else { 184 | psnr(pixels, cpixels, w*h, &PSNR, &PSNR_YUV); 185 | } 186 | printf("\"%s\", -%.1f, -%.1f, \"%s\", %i, %i\n", namelist[n]->d_name, PSNR, PSNR_YUV, noi_profile_name(profile), w, h ); 187 | stbi_write_png(test_name, w, h, 4, cpixels, w*4); 188 | free(cpixels); 189 | free(cbytes); 190 | } 191 | stbi_image_free(pixels); 192 | } 193 | } 194 | free(namelist[n]); 195 | } 196 | } 197 | free(namelist); 198 | return 0; 199 | } else if ( strcmp(argv[1],"-c")==0 || strcmp(argv[1],"-pc")==0 ) { 200 | int profile = argc==5 ? get_profile(argv[4]) : NOI_YUV_16_1_1; 201 | int bx, by; 202 | noi_profile_block_size(profile, &bx, &by); 203 | int w, h; 204 | uint8_t * pixels = stbi_load(argv[2], &w, &h, NULL, 4); 205 | if ( !pixels ) { 206 | printf("can't load image from %s\n", argv[2]); 207 | return -2; 208 | } 209 | if ( (w%bx) || (h%by) ) { 210 | printf("image dimensions need to be proportional to %ix%i block, and not %ix%i\n", bx, by, w, h ); 211 | return -3; 212 | } 213 | printf("noi_compress %i x %i profile %s\n", w, h, noi_profile_name(profile)); 214 | int csize; 215 | uint64_t t0 = ref_time_ticks(); 216 | void * cbytes = noi_compress(pixels, w, h, &csize, profile); 217 | double sec = get_time_usec(t0) / 1000000.0; 218 | if ( !cbytes ) { 219 | printf("compression failed\n"); 220 | return -4; 221 | } 222 | double mb = ((double)(w*h*3))/1024./1024.; 223 | printf("%i mb in %.2f sec, %.3fmb/sec\n", ((int)mb), sec, mb/sec ); 224 | if ( strcmp(argv[1],"-pc")==0 ) { 225 | uint8_t * cpixels = noi_decompress(cbytes, NULL, NULL, NULL); 226 | double PSNR, PSNR_YUV; 227 | if ( profile==NOI_Y_1_0_0 ) { 228 | psnr_greyscale(pixels, cpixels, w*h, &PSNR); 229 | PSNR_YUV = PSNR; 230 | } else { 231 | psnr(pixels, cpixels, w*h, &PSNR, &PSNR_YUV); 232 | } 233 | printf("PSNR = -%.1f PSNR(YUV) = -%.1f\n", PSNR, PSNR_YUV); 234 | free(cpixels); 235 | } 236 | stbi_image_free(pixels); 237 | FILE * f = fopen(argv[3], "wb"); 238 | if ( !f ) { 239 | printf("can't save noi to %s\n", argv[2]); 240 | return -5; 241 | } 242 | fwrite(cbytes, 1, csize, f); 243 | fclose(f); 244 | free(cbytes); 245 | return 0; 246 | } else if ( strcmp(argv[1],"-d")==0 || strcmp(argv[1],"-pd")==0 ) { 247 | int csize; 248 | void * cbytes = load_file(argv[2], &csize); 249 | if ( !cbytes ) { 250 | printf("can't load NOI %s\n", argv[2]); 251 | return -6; 252 | } 253 | int w, h; 254 | if ( strcmp(argv[1],"-pd")==0 ) { 255 | int nTimes = 100; 256 | printf("running noi_decompressing %i times, %i bytes\n", nTimes, csize); 257 | noi_image_size(cbytes, &w, &h); 258 | uint8_t * pixels = (uint8_t *) malloc(w*h*4); 259 | double minsec = 100500.; 260 | for ( int t=0; t!=nTimes; ++t ) { 261 | uint64_t t0 = ref_time_ticks(); 262 | noi_decompress(cbytes, &w, &h, pixels); 263 | double sec = get_time_usec(t0) / 1000000.0; 264 | if ( sec