g30rv17ys commited on
Commit
eec42bd
·
verified ·
1 Parent(s): fd4b932

Add files using upload-large-folder tool

Browse files
This view is limited to 50 files because it contains too many changes.   See raw diff
Files changed (50) hide show
  1. Dockerfile +6 -0
  2. LICENSE +674 -0
  3. README.md +92 -0
  4. convert_to_onnx.sh +8 -0
  5. convert_to_tfjs.py +58 -0
  6. data/NN_human_mouse_eyes/annotation/annotations.csv +0 -0
  7. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0001.png +0 -0
  8. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0002.png +0 -0
  9. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0003.png +0 -0
  10. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0004.png +0 -0
  11. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0005.png +0 -0
  12. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0006.png +0 -0
  13. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0007.png +0 -0
  14. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0008.png +0 -0
  15. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0009.png +0 -0
  16. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0010.png +0 -0
  17. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0011.png +0 -0
  18. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0012.png +0 -0
  19. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0013.png +0 -0
  20. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0014.png +0 -0
  21. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0015.png +0 -0
  22. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0016.png +0 -0
  23. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0017.png +0 -0
  24. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0018.png +0 -0
  25. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0019.png +0 -0
  26. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0020.png +0 -0
  27. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0021.png +0 -0
  28. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0022.png +0 -0
  29. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0023.png +0 -0
  30. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0024.png +0 -0
  31. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0025.png +0 -0
  32. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0026.png +0 -0
  33. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0027.png +0 -0
  34. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0028.png +0 -0
  35. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0029.png +0 -0
  36. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0040.png +0 -0
  37. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0041.png +0 -0
  38. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0042.png +0 -0
  39. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0043.png +0 -0
  40. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0044.png +0 -0
  41. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0045.png +0 -0
  42. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0046.png +0 -0
  43. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0047.png +0 -0
  44. data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0048.png +0 -0
  45. dataloader.py +240 -0
  46. evaluate.py +339 -0
  47. models/__pycache__/unet.cpython-311.pyc +0 -0
  48. models/deeplab/inference/tflite/stb_image.h +0 -0
  49. models/deeplab/inference/tflite/stb_image_resize.h +0 -0
  50. models/deeplab/inference/tflite/stb_image_write.h +1622 -0
Dockerfile ADDED
@@ -0,0 +1,6 @@
 
 
 
 
 
 
 
1
+ FROM tensorflow/tensorflow:2.4.1-gpu-jupyter
2
+
3
+ RUN apt update && apt install -y nvidia-modprobe
4
+ RUN pip install adabelief-tf pandas imageio imageio-ffmpeg seaborn sklearn tensorflow_addons tqdm
5
+ RUN pip install keras_applications tensorflowjs
6
+ RUN pip install tf2onnx
LICENSE ADDED
@@ -0,0 +1,674 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ GNU GENERAL PUBLIC LICENSE
2
+ Version 3, 29 June 2007
3
+
4
+ Copyright (C) 2007 Free Software Foundation, Inc. <https://fsf.org/>
5
+ Everyone is permitted to copy and distribute verbatim copies
6
+ of this license document, but changing it is not allowed.
7
+
8
+ Preamble
9
+
10
+ The GNU General Public License is a free, copyleft license for
11
+ software and other kinds of works.
12
+
13
+ The licenses for most software and other practical works are designed
14
+ to take away your freedom to share and change the works. By contrast,
15
+ the GNU General Public License is intended to guarantee your freedom to
16
+ share and change all versions of a program--to make sure it remains free
17
+ software for all its users. We, the Free Software Foundation, use the
18
+ GNU General Public License for most of our software; it applies also to
19
+ any other work released this way by its authors. You can apply it to
20
+ your programs, too.
21
+
22
+ When we speak of free software, we are referring to freedom, not
23
+ price. Our General Public Licenses are designed to make sure that you
24
+ have the freedom to distribute copies of free software (and charge for
25
+ them if you wish), that you receive source code or can get it if you
26
+ want it, that you can change the software or use pieces of it in new
27
+ free programs, and that you know you can do these things.
28
+
29
+ To protect your rights, we need to prevent others from denying you
30
+ these rights or asking you to surrender the rights. Therefore, you have
31
+ certain responsibilities if you distribute copies of the software, or if
32
+ you modify it: responsibilities to respect the freedom of others.
33
+
34
+ For example, if you distribute copies of such a program, whether
35
+ gratis or for a fee, you must pass on to the recipients the same
36
+ freedoms that you received. You must make sure that they, too, receive
37
+ or can get the source code. And you must show them these terms so they
38
+ know their rights.
39
+
40
+ Developers that use the GNU GPL protect your rights with two steps:
41
+ (1) assert copyright on the software, and (2) offer you this License
42
+ giving you legal permission to copy, distribute and/or modify it.
43
+
44
+ For the developers' and authors' protection, the GPL clearly explains
45
+ that there is no warranty for this free software. For both users' and
46
+ authors' sake, the GPL requires that modified versions be marked as
47
+ changed, so that their problems will not be attributed erroneously to
48
+ authors of previous versions.
49
+
50
+ Some devices are designed to deny users access to install or run
51
+ modified versions of the software inside them, although the manufacturer
52
+ can do so. This is fundamentally incompatible with the aim of
53
+ protecting users' freedom to change the software. The systematic
54
+ pattern of such abuse occurs in the area of products for individuals to
55
+ use, which is precisely where it is most unacceptable. Therefore, we
56
+ have designed this version of the GPL to prohibit the practice for those
57
+ products. If such problems arise substantially in other domains, we
58
+ stand ready to extend this provision to those domains in future versions
59
+ of the GPL, as needed to protect the freedom of users.
60
+
61
+ Finally, every program is threatened constantly by software patents.
62
+ States should not allow patents to restrict development and use of
63
+ software on general-purpose computers, but in those that do, we wish to
64
+ avoid the special danger that patents applied to a free program could
65
+ make it effectively proprietary. To prevent this, the GPL assures that
66
+ patents cannot be used to render the program non-free.
67
+
68
+ The precise terms and conditions for copying, distribution and
69
+ modification follow.
70
+
71
+ TERMS AND CONDITIONS
72
+
73
+ 0. Definitions.
74
+
75
+ "This License" refers to version 3 of the GNU General Public License.
76
+
77
+ "Copyright" also means copyright-like laws that apply to other kinds of
78
+ works, such as semiconductor masks.
79
+
80
+ "The Program" refers to any copyrightable work licensed under this
81
+ License. Each licensee is addressed as "you". "Licensees" and
82
+ "recipients" may be individuals or organizations.
83
+
84
+ To "modify" a work means to copy from or adapt all or part of the work
85
+ in a fashion requiring copyright permission, other than the making of an
86
+ exact copy. The resulting work is called a "modified version" of the
87
+ earlier work or a work "based on" the earlier work.
88
+
89
+ A "covered work" means either the unmodified Program or a work based
90
+ on the Program.
91
+
92
+ To "propagate" a work means to do anything with it that, without
93
+ permission, would make you directly or secondarily liable for
94
+ infringement under applicable copyright law, except executing it on a
95
+ computer or modifying a private copy. Propagation includes copying,
96
+ distribution (with or without modification), making available to the
97
+ public, and in some countries other activities as well.
98
+
99
+ To "convey" a work means any kind of propagation that enables other
100
+ parties to make or receive copies. Mere interaction with a user through
101
+ a computer network, with no transfer of a copy, is not conveying.
102
+
103
+ An interactive user interface displays "Appropriate Legal Notices"
104
+ to the extent that it includes a convenient and prominently visible
105
+ feature that (1) displays an appropriate copyright notice, and (2)
106
+ tells the user that there is no warranty for the work (except to the
107
+ extent that warranties are provided), that licensees may convey the
108
+ work under this License, and how to view a copy of this License. If
109
+ the interface presents a list of user commands or options, such as a
110
+ menu, a prominent item in the list meets this criterion.
111
+
112
+ 1. Source Code.
113
+
114
+ The "source code" for a work means the preferred form of the work
115
+ for making modifications to it. "Object code" means any non-source
116
+ form of a work.
117
+
118
+ A "Standard Interface" means an interface that either is an official
119
+ standard defined by a recognized standards body, or, in the case of
120
+ interfaces specified for a particular programming language, one that
121
+ is widely used among developers working in that language.
122
+
123
+ The "System Libraries" of an executable work include anything, other
124
+ than the work as a whole, that (a) is included in the normal form of
125
+ packaging a Major Component, but which is not part of that Major
126
+ Component, and (b) serves only to enable use of the work with that
127
+ Major Component, or to implement a Standard Interface for which an
128
+ implementation is available to the public in source code form. A
129
+ "Major Component", in this context, means a major essential component
130
+ (kernel, window system, and so on) of the specific operating system
131
+ (if any) on which the executable work runs, or a compiler used to
132
+ produce the work, or an object code interpreter used to run it.
133
+
134
+ The "Corresponding Source" for a work in object code form means all
135
+ the source code needed to generate, install, and (for an executable
136
+ work) run the object code and to modify the work, including scripts to
137
+ control those activities. However, it does not include the work's
138
+ System Libraries, or general-purpose tools or generally available free
139
+ programs which are used unmodified in performing those activities but
140
+ which are not part of the work. For example, Corresponding Source
141
+ includes interface definition files associated with source files for
142
+ the work, and the source code for shared libraries and dynamically
143
+ linked subprograms that the work is specifically designed to require,
144
+ such as by intimate data communication or control flow between those
145
+ subprograms and other parts of the work.
146
+
147
+ The Corresponding Source need not include anything that users
148
+ can regenerate automatically from other parts of the Corresponding
149
+ Source.
150
+
151
+ The Corresponding Source for a work in source code form is that
152
+ same work.
153
+
154
+ 2. Basic Permissions.
155
+
156
+ All rights granted under this License are granted for the term of
157
+ copyright on the Program, and are irrevocable provided the stated
158
+ conditions are met. This License explicitly affirms your unlimited
159
+ permission to run the unmodified Program. The output from running a
160
+ covered work is covered by this License only if the output, given its
161
+ content, constitutes a covered work. This License acknowledges your
162
+ rights of fair use or other equivalent, as provided by copyright law.
163
+
164
+ You may make, run and propagate covered works that you do not
165
+ convey, without conditions so long as your license otherwise remains
166
+ in force. You may convey covered works to others for the sole purpose
167
+ of having them make modifications exclusively for you, or provide you
168
+ with facilities for running those works, provided that you comply with
169
+ the terms of this License in conveying all material for which you do
170
+ not control copyright. Those thus making or running the covered works
171
+ for you must do so exclusively on your behalf, under your direction
172
+ and control, on terms that prohibit them from making any copies of
173
+ your copyrighted material outside their relationship with you.
174
+
175
+ Conveying under any other circumstances is permitted solely under
176
+ the conditions stated below. Sublicensing is not allowed; section 10
177
+ makes it unnecessary.
178
+
179
+ 3. Protecting Users' Legal Rights From Anti-Circumvention Law.
180
+
181
+ No covered work shall be deemed part of an effective technological
182
+ measure under any applicable law fulfilling obligations under article
183
+ 11 of the WIPO copyright treaty adopted on 20 December 1996, or
184
+ similar laws prohibiting or restricting circumvention of such
185
+ measures.
186
+
187
+ When you convey a covered work, you waive any legal power to forbid
188
+ circumvention of technological measures to the extent such circumvention
189
+ is effected by exercising rights under this License with respect to
190
+ the covered work, and you disclaim any intention to limit operation or
191
+ modification of the work as a means of enforcing, against the work's
192
+ users, your or third parties' legal rights to forbid circumvention of
193
+ technological measures.
194
+
195
+ 4. Conveying Verbatim Copies.
196
+
197
+ You may convey verbatim copies of the Program's source code as you
198
+ receive it, in any medium, provided that you conspicuously and
199
+ appropriately publish on each copy an appropriate copyright notice;
200
+ keep intact all notices stating that this License and any
201
+ non-permissive terms added in accord with section 7 apply to the code;
202
+ keep intact all notices of the absence of any warranty; and give all
203
+ recipients a copy of this License along with the Program.
204
+
205
+ You may charge any price or no price for each copy that you convey,
206
+ and you may offer support or warranty protection for a fee.
207
+
208
+ 5. Conveying Modified Source Versions.
209
+
210
+ You may convey a work based on the Program, or the modifications to
211
+ produce it from the Program, in the form of source code under the
212
+ terms of section 4, provided that you also meet all of these conditions:
213
+
214
+ a) The work must carry prominent notices stating that you modified
215
+ it, and giving a relevant date.
216
+
217
+ b) The work must carry prominent notices stating that it is
218
+ released under this License and any conditions added under section
219
+ 7. This requirement modifies the requirement in section 4 to
220
+ "keep intact all notices".
221
+
222
+ c) You must license the entire work, as a whole, under this
223
+ License to anyone who comes into possession of a copy. This
224
+ License will therefore apply, along with any applicable section 7
225
+ additional terms, to the whole of the work, and all its parts,
226
+ regardless of how they are packaged. This License gives no
227
+ permission to license the work in any other way, but it does not
228
+ invalidate such permission if you have separately received it.
229
+
230
+ d) If the work has interactive user interfaces, each must display
231
+ Appropriate Legal Notices; however, if the Program has interactive
232
+ interfaces that do not display Appropriate Legal Notices, your
233
+ work need not make them do so.
234
+
235
+ A compilation of a covered work with other separate and independent
236
+ works, which are not by their nature extensions of the covered work,
237
+ and which are not combined with it such as to form a larger program,
238
+ in or on a volume of a storage or distribution medium, is called an
239
+ "aggregate" if the compilation and its resulting copyright are not
240
+ used to limit the access or legal rights of the compilation's users
241
+ beyond what the individual works permit. Inclusion of a covered work
242
+ in an aggregate does not cause this License to apply to the other
243
+ parts of the aggregate.
244
+
245
+ 6. Conveying Non-Source Forms.
246
+
247
+ You may convey a covered work in object code form under the terms
248
+ of sections 4 and 5, provided that you also convey the
249
+ machine-readable Corresponding Source under the terms of this License,
250
+ in one of these ways:
251
+
252
+ a) Convey the object code in, or embodied in, a physical product
253
+ (including a physical distribution medium), accompanied by the
254
+ Corresponding Source fixed on a durable physical medium
255
+ customarily used for software interchange.
256
+
257
+ b) Convey the object code in, or embodied in, a physical product
258
+ (including a physical distribution medium), accompanied by a
259
+ written offer, valid for at least three years and valid for as
260
+ long as you offer spare parts or customer support for that product
261
+ model, to give anyone who possesses the object code either (1) a
262
+ copy of the Corresponding Source for all the software in the
263
+ product that is covered by this License, on a durable physical
264
+ medium customarily used for software interchange, for a price no
265
+ more than your reasonable cost of physically performing this
266
+ conveying of source, or (2) access to copy the
267
+ Corresponding Source from a network server at no charge.
268
+
269
+ c) Convey individual copies of the object code with a copy of the
270
+ written offer to provide the Corresponding Source. This
271
+ alternative is allowed only occasionally and noncommercially, and
272
+ only if you received the object code with such an offer, in accord
273
+ with subsection 6b.
274
+
275
+ d) Convey the object code by offering access from a designated
276
+ place (gratis or for a charge), and offer equivalent access to the
277
+ Corresponding Source in the same way through the same place at no
278
+ further charge. You need not require recipients to copy the
279
+ Corresponding Source along with the object code. If the place to
280
+ copy the object code is a network server, the Corresponding Source
281
+ may be on a different server (operated by you or a third party)
282
+ that supports equivalent copying facilities, provided you maintain
283
+ clear directions next to the object code saying where to find the
284
+ Corresponding Source. Regardless of what server hosts the
285
+ Corresponding Source, you remain obligated to ensure that it is
286
+ available for as long as needed to satisfy these requirements.
287
+
288
+ e) Convey the object code using peer-to-peer transmission, provided
289
+ you inform other peers where the object code and Corresponding
290
+ Source of the work are being offered to the general public at no
291
+ charge under subsection 6d.
292
+
293
+ A separable portion of the object code, whose source code is excluded
294
+ from the Corresponding Source as a System Library, need not be
295
+ included in conveying the object code work.
296
+
297
+ A "User Product" is either (1) a "consumer product", which means any
298
+ tangible personal property which is normally used for personal, family,
299
+ or household purposes, or (2) anything designed or sold for incorporation
300
+ into a dwelling. In determining whether a product is a consumer product,
301
+ doubtful cases shall be resolved in favor of coverage. For a particular
302
+ product received by a particular user, "normally used" refers to a
303
+ typical or common use of that class of product, regardless of the status
304
+ of the particular user or of the way in which the particular user
305
+ actually uses, or expects or is expected to use, the product. A product
306
+ is a consumer product regardless of whether the product has substantial
307
+ commercial, industrial or non-consumer uses, unless such uses represent
308
+ the only significant mode of use of the product.
309
+
310
+ "Installation Information" for a User Product means any methods,
311
+ procedures, authorization keys, or other information required to install
312
+ and execute modified versions of a covered work in that User Product from
313
+ a modified version of its Corresponding Source. The information must
314
+ suffice to ensure that the continued functioning of the modified object
315
+ code is in no case prevented or interfered with solely because
316
+ modification has been made.
317
+
318
+ If you convey an object code work under this section in, or with, or
319
+ specifically for use in, a User Product, and the conveying occurs as
320
+ part of a transaction in which the right of possession and use of the
321
+ User Product is transferred to the recipient in perpetuity or for a
322
+ fixed term (regardless of how the transaction is characterized), the
323
+ Corresponding Source conveyed under this section must be accompanied
324
+ by the Installation Information. But this requirement does not apply
325
+ if neither you nor any third party retains the ability to install
326
+ modified object code on the User Product (for example, the work has
327
+ been installed in ROM).
328
+
329
+ The requirement to provide Installation Information does not include a
330
+ requirement to continue to provide support service, warranty, or updates
331
+ for a work that has been modified or installed by the recipient, or for
332
+ the User Product in which it has been modified or installed. Access to a
333
+ network may be denied when the modification itself materially and
334
+ adversely affects the operation of the network or violates the rules and
335
+ protocols for communication across the network.
336
+
337
+ Corresponding Source conveyed, and Installation Information provided,
338
+ in accord with this section must be in a format that is publicly
339
+ documented (and with an implementation available to the public in
340
+ source code form), and must require no special password or key for
341
+ unpacking, reading or copying.
342
+
343
+ 7. Additional Terms.
344
+
345
+ "Additional permissions" are terms that supplement the terms of this
346
+ License by making exceptions from one or more of its conditions.
347
+ Additional permissions that are applicable to the entire Program shall
348
+ be treated as though they were included in this License, to the extent
349
+ that they are valid under applicable law. If additional permissions
350
+ apply only to part of the Program, that part may be used separately
351
+ under those permissions, but the entire Program remains governed by
352
+ this License without regard to the additional permissions.
353
+
354
+ When you convey a copy of a covered work, you may at your option
355
+ remove any additional permissions from that copy, or from any part of
356
+ it. (Additional permissions may be written to require their own
357
+ removal in certain cases when you modify the work.) You may place
358
+ additional permissions on material, added by you to a covered work,
359
+ for which you have or can give appropriate copyright permission.
360
+
361
+ Notwithstanding any other provision of this License, for material you
362
+ add to a covered work, you may (if authorized by the copyright holders of
363
+ that material) supplement the terms of this License with terms:
364
+
365
+ a) Disclaiming warranty or limiting liability differently from the
366
+ terms of sections 15 and 16 of this License; or
367
+
368
+ b) Requiring preservation of specified reasonable legal notices or
369
+ author attributions in that material or in the Appropriate Legal
370
+ Notices displayed by works containing it; or
371
+
372
+ c) Prohibiting misrepresentation of the origin of that material, or
373
+ requiring that modified versions of such material be marked in
374
+ reasonable ways as different from the original version; or
375
+
376
+ d) Limiting the use for publicity purposes of names of licensors or
377
+ authors of the material; or
378
+
379
+ e) Declining to grant rights under trademark law for use of some
380
+ trade names, trademarks, or service marks; or
381
+
382
+ f) Requiring indemnification of licensors and authors of that
383
+ material by anyone who conveys the material (or modified versions of
384
+ it) with contractual assumptions of liability to the recipient, for
385
+ any liability that these contractual assumptions directly impose on
386
+ those licensors and authors.
387
+
388
+ All other non-permissive additional terms are considered "further
389
+ restrictions" within the meaning of section 10. If the Program as you
390
+ received it, or any part of it, contains a notice stating that it is
391
+ governed by this License along with a term that is a further
392
+ restriction, you may remove that term. If a license document contains
393
+ a further restriction but permits relicensing or conveying under this
394
+ License, you may add to a covered work material governed by the terms
395
+ of that license document, provided that the further restriction does
396
+ not survive such relicensing or conveying.
397
+
398
+ If you add terms to a covered work in accord with this section, you
399
+ must place, in the relevant source files, a statement of the
400
+ additional terms that apply to those files, or a notice indicating
401
+ where to find the applicable terms.
402
+
403
+ Additional terms, permissive or non-permissive, may be stated in the
404
+ form of a separately written license, or stated as exceptions;
405
+ the above requirements apply either way.
406
+
407
+ 8. Termination.
408
+
409
+ You may not propagate or modify a covered work except as expressly
410
+ provided under this License. Any attempt otherwise to propagate or
411
+ modify it is void, and will automatically terminate your rights under
412
+ this License (including any patent licenses granted under the third
413
+ paragraph of section 11).
414
+
415
+ However, if you cease all violation of this License, then your
416
+ license from a particular copyright holder is reinstated (a)
417
+ provisionally, unless and until the copyright holder explicitly and
418
+ finally terminates your license, and (b) permanently, if the copyright
419
+ holder fails to notify you of the violation by some reasonable means
420
+ prior to 60 days after the cessation.
421
+
422
+ Moreover, your license from a particular copyright holder is
423
+ reinstated permanently if the copyright holder notifies you of the
424
+ violation by some reasonable means, this is the first time you have
425
+ received notice of violation of this License (for any work) from that
426
+ copyright holder, and you cure the violation prior to 30 days after
427
+ your receipt of the notice.
428
+
429
+ Termination of your rights under this section does not terminate the
430
+ licenses of parties who have received copies or rights from you under
431
+ this License. If your rights have been terminated and not permanently
432
+ reinstated, you do not qualify to receive new licenses for the same
433
+ material under section 10.
434
+
435
+ 9. Acceptance Not Required for Having Copies.
436
+
437
+ You are not required to accept this License in order to receive or
438
+ run a copy of the Program. Ancillary propagation of a covered work
439
+ occurring solely as a consequence of using peer-to-peer transmission
440
+ to receive a copy likewise does not require acceptance. However,
441
+ nothing other than this License grants you permission to propagate or
442
+ modify any covered work. These actions infringe copyright if you do
443
+ not accept this License. Therefore, by modifying or propagating a
444
+ covered work, you indicate your acceptance of this License to do so.
445
+
446
+ 10. Automatic Licensing of Downstream Recipients.
447
+
448
+ Each time you convey a covered work, the recipient automatically
449
+ receives a license from the original licensors, to run, modify and
450
+ propagate that work, subject to this License. You are not responsible
451
+ for enforcing compliance by third parties with this License.
452
+
453
+ An "entity transaction" is a transaction transferring control of an
454
+ organization, or substantially all assets of one, or subdividing an
455
+ organization, or merging organizations. If propagation of a covered
456
+ work results from an entity transaction, each party to that
457
+ transaction who receives a copy of the work also receives whatever
458
+ licenses to the work the party's predecessor in interest had or could
459
+ give under the previous paragraph, plus a right to possession of the
460
+ Corresponding Source of the work from the predecessor in interest, if
461
+ the predecessor has it or can get it with reasonable efforts.
462
+
463
+ You may not impose any further restrictions on the exercise of the
464
+ rights granted or affirmed under this License. For example, you may
465
+ not impose a license fee, royalty, or other charge for exercise of
466
+ rights granted under this License, and you may not initiate litigation
467
+ (including a cross-claim or counterclaim in a lawsuit) alleging that
468
+ any patent claim is infringed by making, using, selling, offering for
469
+ sale, or importing the Program or any portion of it.
470
+
471
+ 11. Patents.
472
+
473
+ A "contributor" is a copyright holder who authorizes use under this
474
+ License of the Program or a work on which the Program is based. The
475
+ work thus licensed is called the contributor's "contributor version".
476
+
477
+ A contributor's "essential patent claims" are all patent claims
478
+ owned or controlled by the contributor, whether already acquired or
479
+ hereafter acquired, that would be infringed by some manner, permitted
480
+ by this License, of making, using, or selling its contributor version,
481
+ but do not include claims that would be infringed only as a
482
+ consequence of further modification of the contributor version. For
483
+ purposes of this definition, "control" includes the right to grant
484
+ patent sublicenses in a manner consistent with the requirements of
485
+ this License.
486
+
487
+ Each contributor grants you a non-exclusive, worldwide, royalty-free
488
+ patent license under the contributor's essential patent claims, to
489
+ make, use, sell, offer for sale, import and otherwise run, modify and
490
+ propagate the contents of its contributor version.
491
+
492
+ In the following three paragraphs, a "patent license" is any express
493
+ agreement or commitment, however denominated, not to enforce a patent
494
+ (such as an express permission to practice a patent or covenant not to
495
+ sue for patent infringement). To "grant" such a patent license to a
496
+ party means to make such an agreement or commitment not to enforce a
497
+ patent against the party.
498
+
499
+ If you convey a covered work, knowingly relying on a patent license,
500
+ and the Corresponding Source of the work is not available for anyone
501
+ to copy, free of charge and under the terms of this License, through a
502
+ publicly available network server or other readily accessible means,
503
+ then you must either (1) cause the Corresponding Source to be so
504
+ available, or (2) arrange to deprive yourself of the benefit of the
505
+ patent license for this particular work, or (3) arrange, in a manner
506
+ consistent with the requirements of this License, to extend the patent
507
+ license to downstream recipients. "Knowingly relying" means you have
508
+ actual knowledge that, but for the patent license, your conveying the
509
+ covered work in a country, or your recipient's use of the covered work
510
+ in a country, would infringe one or more identifiable patents in that
511
+ country that you have reason to believe are valid.
512
+
513
+ If, pursuant to or in connection with a single transaction or
514
+ arrangement, you convey, or propagate by procuring conveyance of, a
515
+ covered work, and grant a patent license to some of the parties
516
+ receiving the covered work authorizing them to use, propagate, modify
517
+ or convey a specific copy of the covered work, then the patent license
518
+ you grant is automatically extended to all recipients of the covered
519
+ work and works based on it.
520
+
521
+ A patent license is "discriminatory" if it does not include within
522
+ the scope of its coverage, prohibits the exercise of, or is
523
+ conditioned on the non-exercise of one or more of the rights that are
524
+ specifically granted under this License. You may not convey a covered
525
+ work if you are a party to an arrangement with a third party that is
526
+ in the business of distributing software, under which you make payment
527
+ to the third party based on the extent of your activity of conveying
528
+ the work, and under which the third party grants, to any of the
529
+ parties who would receive the covered work from you, a discriminatory
530
+ patent license (a) in connection with copies of the covered work
531
+ conveyed by you (or copies made from those copies), or (b) primarily
532
+ for and in connection with specific products or compilations that
533
+ contain the covered work, unless you entered into that arrangement,
534
+ or that patent license was granted, prior to 28 March 2007.
535
+
536
+ Nothing in this License shall be construed as excluding or limiting
537
+ any implied license or other defenses to infringement that may
538
+ otherwise be available to you under applicable patent law.
539
+
540
+ 12. No Surrender of Others' Freedom.
541
+
542
+ If conditions are imposed on you (whether by court order, agreement or
543
+ otherwise) that contradict the conditions of this License, they do not
544
+ excuse you from the conditions of this License. If you cannot convey a
545
+ covered work so as to satisfy simultaneously your obligations under this
546
+ License and any other pertinent obligations, then as a consequence you may
547
+ not convey it at all. For example, if you agree to terms that obligate you
548
+ to collect a royalty for further conveying from those to whom you convey
549
+ the Program, the only way you could satisfy both those terms and this
550
+ License would be to refrain entirely from conveying the Program.
551
+
552
+ 13. Use with the GNU Affero General Public License.
553
+
554
+ Notwithstanding any other provision of this License, you have
555
+ permission to link or combine any covered work with a work licensed
556
+ under version 3 of the GNU Affero General Public License into a single
557
+ combined work, and to convey the resulting work. The terms of this
558
+ License will continue to apply to the part which is the covered work,
559
+ but the special requirements of the GNU Affero General Public License,
560
+ section 13, concerning interaction through a network will apply to the
561
+ combination as such.
562
+
563
+ 14. Revised Versions of this License.
564
+
565
+ The Free Software Foundation may publish revised and/or new versions of
566
+ the GNU General Public License from time to time. Such new versions will
567
+ be similar in spirit to the present version, but may differ in detail to
568
+ address new problems or concerns.
569
+
570
+ Each version is given a distinguishing version number. If the
571
+ Program specifies that a certain numbered version of the GNU General
572
+ Public License "or any later version" applies to it, you have the
573
+ option of following the terms and conditions either of that numbered
574
+ version or of any later version published by the Free Software
575
+ Foundation. If the Program does not specify a version number of the
576
+ GNU General Public License, you may choose any version ever published
577
+ by the Free Software Foundation.
578
+
579
+ If the Program specifies that a proxy can decide which future
580
+ versions of the GNU General Public License can be used, that proxy's
581
+ public statement of acceptance of a version permanently authorizes you
582
+ to choose that version for the Program.
583
+
584
+ Later license versions may give you additional or different
585
+ permissions. However, no additional obligations are imposed on any
586
+ author or copyright holder as a result of your choosing to follow a
587
+ later version.
588
+
589
+ 15. Disclaimer of Warranty.
590
+
591
+ THERE IS NO WARRANTY FOR THE PROGRAM, TO THE EXTENT PERMITTED BY
592
+ APPLICABLE LAW. EXCEPT WHEN OTHERWISE STATED IN WRITING THE COPYRIGHT
593
+ HOLDERS AND/OR OTHER PARTIES PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY
594
+ OF ANY KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, BUT NOT LIMITED TO,
595
+ THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
596
+ PURPOSE. THE ENTIRE RISK AS TO THE QUALITY AND PERFORMANCE OF THE PROGRAM
597
+ IS WITH YOU. SHOULD THE PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF
598
+ ALL NECESSARY SERVICING, REPAIR OR CORRECTION.
599
+
600
+ 16. Limitation of Liability.
601
+
602
+ IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING
603
+ WILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MODIFIES AND/OR CONVEYS
604
+ THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES, INCLUDING ANY
605
+ GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING OUT OF THE
606
+ USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED TO LOSS OF
607
+ DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY YOU OR THIRD
608
+ PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER PROGRAMS),
609
+ EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE POSSIBILITY OF
610
+ SUCH DAMAGES.
611
+
612
+ 17. Interpretation of Sections 15 and 16.
613
+
614
+ If the disclaimer of warranty and limitation of liability provided
615
+ above cannot be given local legal effect according to their terms,
616
+ reviewing courts shall apply local law that most closely approximates
617
+ an absolute waiver of all civil liability in connection with the
618
+ Program, unless a warranty or assumption of liability accompanies a
619
+ copy of the Program in return for a fee.
620
+
621
+ END OF TERMS AND CONDITIONS
622
+
623
+ How to Apply These Terms to Your New Programs
624
+
625
+ If you develop a new program, and you want it to be of the greatest
626
+ possible use to the public, the best way to achieve this is to make it
627
+ free software which everyone can redistribute and change under these terms.
628
+
629
+ To do so, attach the following notices to the program. It is safest
630
+ to attach them to the start of each source file to most effectively
631
+ state the exclusion of warranty; and each file should have at least
632
+ the "copyright" line and a pointer to where the full notice is found.
633
+
634
+ <one line to give the program's name and a brief idea of what it does.>
635
+ Copyright (C) <year> <name of author>
636
+
637
+ This program is free software: you can redistribute it and/or modify
638
+ it under the terms of the GNU General Public License as published by
639
+ the Free Software Foundation, either version 3 of the License, or
640
+ (at your option) any later version.
641
+
642
+ This program is distributed in the hope that it will be useful,
643
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
644
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
645
+ GNU General Public License for more details.
646
+
647
+ You should have received a copy of the GNU General Public License
648
+ along with this program. If not, see <https://www.gnu.org/licenses/>.
649
+
650
+ Also add information on how to contact you by electronic and paper mail.
651
+
652
+ If the program does terminal interaction, make it output a short
653
+ notice like this when it starts in an interactive mode:
654
+
655
+ <program> Copyright (C) <year> <name of author>
656
+ This program comes with ABSOLUTELY NO WARRANTY; for details type `show w'.
657
+ This is free software, and you are welcome to redistribute it
658
+ under certain conditions; type `show c' for details.
659
+
660
+ The hypothetical commands `show w' and `show c' should show the appropriate
661
+ parts of the General Public License. Of course, your program's commands
662
+ might be different; for a GUI interface, you would use an "about box".
663
+
664
+ You should also get your employer (if you work as a programmer) or school,
665
+ if any, to sign a "copyright disclaimer" for the program, if necessary.
666
+ For more information on this, and how to apply and follow the GNU GPL, see
667
+ <https://www.gnu.org/licenses/>.
668
+
669
+ The GNU General Public License does not permit incorporating your program
670
+ into proprietary programs. If your program is a subroutine library, you
671
+ may consider it more useful to permit linking proprietary applications with
672
+ the library. If this is what you want to do, use the GNU Lesser General
673
+ Public License instead of this License. But first, please read
674
+ <https://www.gnu.org/licenses/why-not-lgpl.html>.
README.md ADDED
@@ -0,0 +1,92 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ # [*mEye*](https://www.pupillometry.it): A Deep Learning Tool for Pupillometry
2
+
3
+ > ⭐ MEYE is available on **MATLAB**! Check it out [here](matlab/README.md)
4
+
5
+ > Check out [pupillometry.it](https://www.pupillometry.it) for a ready-to-use web-based mEye pupillometry tool!
6
+
7
+ This branch provides the Python code to make predictions and train/finetune models.
8
+ If you are interested in the code of the pupillometry web app, check out the `gh-pages` branch.
9
+
10
+ ## Requirements
11
+ You need a Python 3 environment with the following packages installed:
12
+
13
+ - tensorflow >= 2.4
14
+ - imageio, imageio-ffmpeg
15
+ - scipy
16
+ - tqdm
17
+
18
+ If you want to train models, you also need
19
+
20
+ - adabelief_tf >= 0.2.1
21
+ - pandas
22
+ - sklearn
23
+
24
+ We provide a [Dockerfile](./Dockerfile) for building an image with docker.
25
+
26
+ ## Make Predictions with Pretrained Models
27
+
28
+ You can make predictions with pretrained models on pre-recorded videos or webcam streams.
29
+
30
+ 1. Download the [pretrained model](https://github.com/fabiocarrara/meye/releases/download/v0.1.1/meye-2022-01-24.h5). If you want to use the [old model](https://github.com/fabiocarrara/meye/releases/download/v0.1/meye-segmentation_i128_s4_c1_f16_g1_a-relu.hdf5), check out version [`v0.1` of this branch](https://github.com/fabiocarrara/meye/tree/v0.1). See available models in [Releases](https://github.com/fabiocarrara/meye/releases).
31
+ 2. Check out the `pupillometry-offline-videos.ipynb` notebook for a complete example of pupillometry data analysis.
32
+ 3. In alternative, we provide also the `predict.py` script that implements the basic loop to make predictions on video streams. E.g.:
33
+
34
+ - ```bash
35
+ # input: webcam (default)
36
+ # prediction roi: biggest central square crop (default)
37
+ # outputs: predictions.mp4, predictions.csv (default)
38
+ predict.py path/to/model
39
+ ```
40
+
41
+ - ```bash
42
+ # input: video file
43
+ # prediction roi: left=80, top=80, right=208, bottom=208
44
+ # outputs: video_with_predictions.mp4, pupil_metrics.csv
45
+ predict.py path/to/model path/to/video.mp4 -rl 80 -rt 80 -rr 208 -rb 208 -ov video_with_predictions.mp4 -oc pupil_metrics.csv
46
+ ```
47
+ - ```bash
48
+ # check all parameters with
49
+ predict.py -h
50
+ ```
51
+
52
+ ## Training Models
53
+
54
+ 1. Download our dataset ([NN_human_mouse_eyes.zip](https://doi.org/10.5281/zenodo.4488164), 246.4 MB) or prepare your dataset following our dataset's structure.
55
+ > If you need to annotate your dataset, check out [pLabeler](https://github.com/LeonardoLupori/pLabeler), a MATLAB software for labeling pupil images.
56
+
57
+ The dataset should be placed in `data/<dataset_name>`.
58
+
59
+ 2. If you are using a custom dataset, edit `train.py` to perform the train/validation/test split of your data.
60
+
61
+ 3. Train with default parameters:
62
+ ```bash
63
+ python train.py -d data/<dataset_name>
64
+ ```
65
+
66
+ - For a list of available parameters, run
67
+ ```bash
68
+ python train.py -h
69
+ ```
70
+
71
+ ## MATLAB support
72
+ Starting from MATLAB version 2021b, MEYE is also available for use on MATLAB!
73
+ A fully functional class and a tutorial for its use is available [here](matlab/README.md)!
74
+
75
+
76
+ ## References
77
+
78
+ ### Dataset
79
+ [![DOI](https://zenodo.org/badge/DOI/10.5281/zenodo.4488164.svg)](https://doi.org/10.5281/zenodo.4488164)
80
+
81
+ If you use our dataset, please cite:
82
+
83
+ @dataset{raffaele_mazziotti_2021_4488164,
84
+ author = {Raffaele Mazziotti and Fabio Carrara and Aurelia Viglione and Lupori Leonardo and Lo Verde Luca and Benedetto Alessandro and Ricci Giulia and Sagona Giulia and Amato Giuseppe and Pizzorusso Tommaso},
85
+ title = {{Human and Mouse Eyes for Pupil Semantic Segmentation}},
86
+ month = feb,
87
+ year = 2021,
88
+ publisher = {Zenodo},
89
+ version = {1.0},
90
+ doi = {10.5281/zenodo.4488164},
91
+ url = {https://doi.org/10.5281/zenodo.4488164}
92
+ }
convert_to_onnx.sh ADDED
@@ -0,0 +1,8 @@
 
 
 
 
 
 
 
 
 
1
+ #/bin/bash
2
+
3
+ # prerequisite:
4
+ # pip install tf2onnx
5
+
6
+ RUNDIR="path/to/run/dir"
7
+
8
+ python -m tf2onnx.convert --saved-model $RUNDIR/best_savedmodel --output $RUNDIR/best_model.onnx
convert_to_tfjs.py ADDED
@@ -0,0 +1,58 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import argparse
2
+ import os
3
+ os.sys.path += ['expman']
4
+ import json
5
+ from glob import glob
6
+ import subprocess
7
+ from tqdm import tqdm
8
+ import expman
9
+
10
+
11
+ def main(args):
12
+
13
+ variants = (
14
+ ('' , []),
15
+ # ('_qf16', ['--quantize_float16', '*']),
16
+ # ('_qu16', ['--quantize_uint16' , '*']),
17
+ # ('_qu8' , ['--quantize_uint8' , '*']),
18
+ )
19
+
20
+ converted_models = []
21
+
22
+ exps = expman.gather(args.run).filter(args.filter)
23
+ for exp_name, exp in tqdm(exps.items()):
24
+ # ckpt = exp.path_to('best_model.h5')
25
+ # ckpt = ckpt if os.path.exists(ckpt) else exp.path_to('last_model.h5')
26
+ ckpt = exp.path_to('best_savedmodel/')
27
+
28
+ for suffix, extra_args in variants:
29
+ name = exp_name + suffix
30
+ out = os.path.join(args.output, name) if args.output else exp.path_to(f'tfjs_graph{suffix}')
31
+
32
+ if args.force or not os.path.exists(out):
33
+ os.makedirs(out, exist_ok=True)
34
+ cmd = ['tensorflowjs_converter',
35
+ '--input_format', 'tf_saved_model',
36
+ '--output_format', 'tfjs_graph_model'] + extra_args + [ckpt, out]
37
+ subprocess.call(cmd)
38
+
39
+ converted_models.append(name)
40
+
41
+ js_output = 'models = ' + json.dumps(converted_models)
42
+ if args.output:
43
+ js_filename = os.path.join(args.output, 'models.js')
44
+ with open(js_filename, 'w') as f:
45
+ f.write(js_output)
46
+ else:
47
+ print(js_output)
48
+
49
+
50
+ if __name__ == '__main__':
51
+ parser = argparse.ArgumentParser(description='Convert runs to tfjs')
52
+ parser.add_argument('-f', '--filter', default={}, type=expman.exp_filter)
53
+ parser.add_argument('run', default='runs/')
54
+ parser.add_argument('--output', help='output dir for models, defaults to run dir')
55
+ parser.add_argument('--force', action='store_true', default=False)
56
+
57
+ args = parser.parse_args()
58
+ main(args)
data/NN_human_mouse_eyes/annotation/annotations.csv ADDED
The diff for this file is too large to render. See raw diff
 
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0001.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0002.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0003.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0004.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0005.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0006.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0007.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0008.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0009.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0010.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0011.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0012.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0013.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0014.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0015.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0016.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0017.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0018.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0019.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0020.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0021.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0022.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0023.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0024.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0025.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0026.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0027.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0028.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0029.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0040.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0041.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0042.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0043.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0044.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0045.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0046.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0047.png ADDED
data/NN_human_mouse_eyes/annotation/png/2P-fc2_save_20200806_1_0048.png ADDED
dataloader.py ADDED
@@ -0,0 +1,240 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ # dataloader.py
2
+
3
+ import os
4
+ import math
5
+ import pandas as pd
6
+ import tensorflow as tf
7
+ from functools import partial
8
+
9
+ def validate_data_files(data_df):
10
+ """Validate all files exist before starting training"""
11
+ missing_files = []
12
+
13
+ for _, row in data_df.iterrows():
14
+ if not os.path.exists(row['filename']):
15
+ missing_files.append(row['filename'])
16
+ if not os.path.exists(row['target']):
17
+ missing_files.append(row['target'])
18
+
19
+ if missing_files:
20
+ print("Missing files:")
21
+ for f in missing_files[:10]:
22
+ print(f" {f}")
23
+ if len(missing_files) > 10:
24
+ print(f" ... and {len(missing_files)-10} more")
25
+ raise FileNotFoundError(f"Found {len(missing_files)} missing files")
26
+
27
+ def _get_pupil_position(pmap, datum, x_shape):
28
+ total_mass = tf.reduce_sum(pmap)
29
+ if total_mass > 0:
30
+ shape = tf.shape(pmap)
31
+ h, w = shape[0], shape[1]
32
+ ii, jj = tf.meshgrid(tf.range(h), tf.range(w), indexing='ij')
33
+ y = tf.reduce_sum(tf.cast(ii, 'float32') * pmap) / total_mass
34
+ x = tf.reduce_sum(tf.cast(jj, 'float32') * pmap) / total_mass
35
+ return tf.stack((y, x))
36
+
37
+ if 'roi_x' in datum and 'roi_y' in datum and 'roi_w' in datum:
38
+ roi_x = tf.cast(datum['roi_x'], 'float32')
39
+ roi_y = tf.cast(datum['roi_y'], 'float32')
40
+ half = tf.cast(datum['roi_w'] / 2, 'float32')
41
+ result = tf.stack((roi_y + half, roi_x + half))
42
+ else: # fallback to center of the image
43
+ result = tf.cast(tf.stack((x_shape[0] / 2, x_shape[1] / 2)), dtype='float32')
44
+
45
+ return result
46
+
47
+ def additional_augmentations(image, mask, p=0.3):
48
+ # Keep original image with probability 1-p
49
+ if tf.random.uniform([]) > p:
50
+ return image, mask
51
+
52
+ # Random noise augmentation
53
+ if tf.random.uniform([]) < 0.3:
54
+ noise = tf.random.normal(tf.shape(image), mean=0, stddev=0.05)
55
+ image = image + noise
56
+ image = tf.clip_by_value(image, 0, 1)
57
+
58
+ # Random blur using gaussian filter
59
+ if tf.random.uniform([]) < 0.3:
60
+ kernel_size = 3
61
+ sigma = tf.random.uniform([], 0, 1.0)
62
+ x = tf.range(-kernel_size // 2 + 1, kernel_size // 2 + 1, dtype=tf.float32)
63
+ gaussian = tf.exp(-(x ** 2) / (2 * sigma ** 2))
64
+ gaussian = gaussian / tf.reduce_sum(gaussian)
65
+ gaussian = tf.reshape(gaussian, [kernel_size, 1])
66
+ gaussian_kernel = gaussian @ tf.transpose(gaussian)
67
+ gaussian_kernel = tf.reshape(gaussian_kernel, [kernel_size, kernel_size, 1, 1])
68
+ image = tf.nn.conv2d(tf.expand_dims(image, 0), gaussian_kernel,
69
+ strides=[1,1,1,1], padding='SAME')[0]
70
+
71
+ # Random sharpening
72
+ if tf.random.uniform([]) < 0.3:
73
+ kernel = tf.constant([[-1,-1,-1], [-1,9,-1], [-1,-1,-1]], dtype=tf.float32)
74
+ kernel = tf.reshape(kernel, [3, 3, 1, 1])
75
+ image = tf.nn.conv2d(tf.expand_dims(image, 0), kernel,
76
+ strides=[1,1,1,1], padding='SAME')[0]
77
+ image = tf.clip_by_value(image, 0, 1)
78
+
79
+ return image, mask
80
+
81
+ @tf.function
82
+ def load_datum(datum, x_shape=(128, 128, 1), augment=False):
83
+ try:
84
+ x = tf.io.read_file(datum['filename'])
85
+ y = tf.io.read_file(datum['target'])
86
+
87
+ # HWC [0,1] float32
88
+ x = tf.io.decode_image(x, channels=1, dtype='float32', expand_animations=False)
89
+ y = tf.io.decode_image(y, dtype='float32', expand_animations=False)
90
+
91
+ # Get image dimensions
92
+ h = tf.shape(x)[0]
93
+ w = tf.shape(x)[1]
94
+
95
+ # Extract pupil information
96
+ pupil_map = y[:, :, 0]
97
+ pupil_area = tf.reduce_sum(pupil_map)
98
+ pupil_pos_yx = _get_pupil_position(pupil_map, datum, x_shape)
99
+
100
+ # Target size we want to achieve
101
+ target_size = tf.minimum(tf.minimum(h, w), x_shape[0])
102
+
103
+ if not augment:
104
+ # Calculate center crop
105
+ h_start = (h - target_size) // 2
106
+ w_start = (w - target_size) // 2
107
+ else:
108
+ # Random crop within safe bounds considering pupil position
109
+ h_start = tf.random.uniform([], 0, h - target_size + 1, dtype=tf.int32)
110
+ w_start = tf.random.uniform([], 0, w - target_size + 1, dtype=tf.int32)
111
+
112
+ # Perform crop
113
+ x = tf.image.crop_to_bounding_box(x, h_start, w_start, target_size, target_size)
114
+ y = tf.image.crop_to_bounding_box(y, h_start, w_start, target_size, target_size)
115
+
116
+ if augment:
117
+ # Rotation with arbitrary angles
118
+ k = tf.random.uniform([], 0, 4, dtype=tf.int32)
119
+ x = tf.image.rot90(x, k=k)
120
+ y = tf.image.rot90(y, k=k)
121
+
122
+ # Flips with pupil position consideration
123
+ if tf.random.uniform([]) < 0.5:
124
+ x = tf.image.flip_left_right(x)
125
+ y = tf.image.flip_left_right(y)
126
+
127
+ if tf.random.uniform([]) < 0.5:
128
+ x = tf.image.flip_up_down(x)
129
+ y = tf.image.flip_up_down(y)
130
+
131
+ # Apply additional augmentations
132
+ x, y = additional_augmentations(x, y)
133
+
134
+ # Calculate pupil visibility after transformation
135
+ new_pupil_map = y[:, :, 0]
136
+ new_pupil_area = tf.reduce_sum(new_pupil_map)
137
+ eye = (new_pupil_area / pupil_area) if pupil_area > 0 else 0.
138
+
139
+ # Process eye and blink information
140
+ datum_eye = tf.cast(datum['eye'], 'float32')
141
+ datum_blink = tf.cast(datum['blink'], 'float32')
142
+
143
+ # Handle blink cases
144
+ if datum_eye == 0:
145
+ datum_blink = 0.
146
+ if (datum_eye == 1) & (datum_blink == 0):
147
+ datum_eye = eye
148
+
149
+ # Resize if needed
150
+ if target_size != x_shape[0]:
151
+ x = tf.image.resize(x, [x_shape[0], x_shape[0]])
152
+ y = tf.image.resize(y, [x_shape[0], x_shape[0]])
153
+
154
+ y = y[:, :, :1]
155
+ y2 = tf.stack((datum_eye, datum_blink))
156
+
157
+ return x, y, y2
158
+
159
+ except Exception as e:
160
+ print(f"Error processing datum: {str(e)}")
161
+ raise
162
+
163
+ def get_loader(dataframe, batch_size=8, shuffle=False, **kwargs):
164
+ categories = dataframe.exp.values
165
+ dataset = tf.data.Dataset.from_tensor_slices(dict(dataframe))
166
+
167
+ if shuffle:
168
+ dataset = dataset.shuffle(1000)
169
+
170
+ dataset = dataset.map(
171
+ partial(load_datum, **kwargs),
172
+ num_parallel_calls=tf.data.AUTOTUNE,
173
+ deterministic=not shuffle
174
+ )
175
+ dataset = dataset.batch(batch_size)
176
+
177
+ def _pack_targets(*ins):
178
+ inputs = ins[0]
179
+ targets = {'mask': ins[1], 'tags': ins[2]}
180
+ return [inputs, targets]
181
+
182
+ dataset = dataset.map(
183
+ _pack_targets,
184
+ num_parallel_calls=tf.data.AUTOTUNE,
185
+ deterministic=not shuffle
186
+ )
187
+ dataset = dataset.prefetch(tf.data.AUTOTUNE)
188
+ return dataset, categories
189
+
190
+ def load_datasets(dataset_dirs):
191
+ def _load_and_prepare_annotations(dataset_dir):
192
+ # Normalize path
193
+ dataset_dir = os.path.normpath(dataset_dir)
194
+
195
+ data_path = os.path.join(dataset_dir, 'annotation', 'annotations.csv')
196
+ if not os.path.exists(data_path):
197
+ raise FileNotFoundError(f"Annotations file not found: {data_path}")
198
+
199
+ data = pd.read_csv(data_path)
200
+
201
+ # Create directories if they don't exist
202
+ png_dir = os.path.join(dataset_dir, 'annotation', 'png')
203
+ full_frames_dir = os.path.join(dataset_dir, 'fullFrames')
204
+ os.makedirs(png_dir, exist_ok=True)
205
+ os.makedirs(full_frames_dir, exist_ok=True)
206
+
207
+ # Filter valid files before creating paths
208
+ valid_data = data[data.apply(lambda x: os.path.exists(os.path.join(full_frames_dir, os.path.basename(x['filename']))), axis=1)]
209
+
210
+ if len(valid_data) == 0:
211
+ raise ValueError(f"No valid image files found in {full_frames_dir}")
212
+
213
+ # Create target paths
214
+ valid_data['target'] = valid_data['filename'].apply(
215
+ lambda x: os.path.join(png_dir, os.path.splitext(os.path.basename(x))[0] + '.png')
216
+ )
217
+
218
+ valid_data['filename'] = valid_data['filename'].apply(
219
+ lambda x: os.path.join(full_frames_dir, os.path.basename(x))
220
+ )
221
+
222
+ return valid_data
223
+
224
+ all_data = []
225
+ for d in dataset_dirs:
226
+ try:
227
+ dataset = _load_and_prepare_annotations(d)
228
+ all_data.append(dataset)
229
+ except Exception as e:
230
+ print(f"Error loading dataset from {d}: {str(e)}")
231
+ continue
232
+
233
+ if not all_data:
234
+ raise ValueError("No valid datasets found in any of the provided directories")
235
+
236
+ dataset = pd.concat(all_data)
237
+ dataset['sub'] = dataset['sub'].astype(str)
238
+
239
+ print(f"Found {len(dataset)} valid image pairs")
240
+ return dataset
evaluate.py ADDED
@@ -0,0 +1,339 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ # -*- coding: utf-8 -*-
2
+ """ MEye: Semantic Segmentation """
3
+
4
+ import argparse
5
+ import os
6
+ os.sys.path += ['expman', 'models/deeplab']
7
+ import expman
8
+
9
+ import matplotlib
10
+ matplotlib.use('Agg')
11
+ import matplotlib.pyplot as plt
12
+
13
+ import numpy as np
14
+ import pandas as pd
15
+ import tensorflow as tf
16
+ from tensorflow.keras import backend as K
17
+ from tensorflow.keras.models import load_model
18
+ from sklearn.model_selection import train_test_split
19
+ from sklearn.metrics import classification_report, roc_curve, auc, precision_recall_curve, average_precision_score
20
+ from adabelief_tf import AdaBeliefOptimizer
21
+ from glob import glob
22
+ from tqdm import tqdm
23
+ from PIL import Image
24
+
25
+ from deeplabv3p.models.deeplabv3p_mobilenetv3 import hard_swish
26
+ from dataloader import get_loader, load_datasets
27
+ from utils import visualize, visualizable
28
+
29
+ from tensorflow.python.framework.convert_to_constants import convert_variables_to_constants_v2_as_graph
30
+
31
+ def iou_coef(y_true, y_pred):
32
+ y_true = tf.cast(y_true, tf.float32)
33
+ y_pred = tf.cast(y_pred, tf.float32)
34
+ intersection = K.sum(K.abs(y_true * y_pred), axis=[1, 2, 3])
35
+ union = K.sum(y_true, axis=[1, 2, 3]) + K.sum(y_pred, axis=[1, 2, 3]) - intersection
36
+ return K.mean((intersection + 1e-6) / (union + 1e-6))
37
+
38
+ def dice_coef(y_true, y_pred):
39
+ y_true = tf.cast(y_true, tf.float32)
40
+ y_pred = tf.cast(y_pred, tf.float32)
41
+ intersection = K.sum(K.abs(y_true * y_pred), axis=[1, 2, 3])
42
+ return K.mean((2. * intersection + 1e-6) / (K.sum(y_true, axis=[1, 2, 3]) + K.sum(y_pred, axis=[1, 2, 3]) + 1e-6))
43
+
44
+ def boundary_loss(y_true, y_pred):
45
+ y_true = tf.cast(y_true, tf.float32)
46
+ y_pred = tf.cast(y_pred, tf.float32)
47
+ dy_true, dx_true = tf.image.image_gradients(y_true)
48
+ dy_pred, dx_pred = tf.image.image_gradients(y_pred)
49
+ loss = tf.reduce_mean(tf.abs(dy_pred - dy_true) + tf.abs(dx_pred - dx_true))
50
+ return loss * 0.5
51
+
52
+ def enhanced_binary_crossentropy(y_true, y_pred):
53
+ y_true = tf.cast(y_true, tf.float32)
54
+ y_pred = tf.cast(y_pred, tf.float32)
55
+ bce = tf.keras.losses.binary_crossentropy(y_true, y_pred)
56
+ boundary = boundary_loss(y_true, y_pred)
57
+ return bce + boundary
58
+
59
+ def _filter_by_closeness(a, eps=10e-3):
60
+ keep = []
61
+ prev = np.array([-1, -1])
62
+ for row in a.drop('thr', axis=1).values:
63
+ if (np.abs(prev - row) > eps).any():
64
+ keep.append(True)
65
+ prev = row
66
+ else:
67
+ keep.append(False)
68
+ return a[keep]
69
+
70
+ def _weighted_roc_pr(y_true, y_scores, label, outdir, simplify=False):
71
+ # Convert to binary labels using 0.5 threshold for true values
72
+ y_true_binary = (y_true > 0.5).astype(int)
73
+
74
+ npos = y_true_binary.sum()
75
+ nneg = len(y_true_binary) - npos
76
+ pos_weight = nneg / npos if npos > 0 else 1.0
77
+ print(label, 'Tot:', len(y_true), 'P:', npos, 'N:', nneg, 'N/P:', pos_weight)
78
+ sample_weight = np.where(y_true_binary, pos_weight, 1)
79
+
80
+ fpr, tpr, thr = roc_curve(y_true_binary, y_scores, sample_weight=sample_weight)
81
+ auc_score = auc(fpr, tpr)
82
+ print(label, 'AuROC:', auc_score)
83
+
84
+ roc_metrics = pd.Series({'npos': npos, 'nneg': nneg, 'nneg_over_npos': pos_weight, 'roc_auc': auc_score})
85
+ roc_metrics_file = os.path.join(outdir, '{}_roc_metrics.csv'.format(label))
86
+ roc_metrics.to_csv(roc_metrics_file, index=False)
87
+
88
+ roc = pd.DataFrame({'fpr': fpr, 'tpr': tpr, 'thr': thr})
89
+ if simplify:
90
+ full_roc_file = os.path.join(outdir, '{}_roc_curve_full.csv.gz'.format(label))
91
+ roc.to_csv(full_roc_file, index=False)
92
+ roc = _filter_by_closeness(roc)
93
+
94
+ roc_file = os.path.join(outdir, '{}_roc_curve.csv'.format(label))
95
+ roc.to_csv(roc_file, index=False)
96
+
97
+ roc.plot(x='fpr', y='tpr', xlim=(0, 1), ylim=(0, 1))
98
+ roc_plot_file = os.path.join(outdir, '{}_roc.pdf'.format(label))
99
+ plt.savefig(roc_plot_file)
100
+ plt.close()
101
+
102
+ precision, recall, thr = precision_recall_curve(y_true_binary, y_scores, sample_weight=sample_weight)
103
+ f1_score = 2 * precision * recall / (precision + recall)
104
+ pr_auc = auc(recall, precision)
105
+
106
+ pr_metrics = pd.Series({'npos': npos, 'nneg': nneg, 'nneg_over_npos': pos_weight, 'pr_auc': pr_auc})
107
+ pr_metrics_file = os.path.join(outdir, '{}_pr_metrics.csv'.format(label))
108
+ pr_metrics.to_csv(pr_metrics_file, index=False)
109
+
110
+ thr = np.append(thr, [thr[-1]])
111
+ pr = pd.DataFrame({'precision': precision, 'recall': recall, 'f1_score': f1_score, 'thr': thr})
112
+ if simplify:
113
+ full_pr_file = os.path.join(outdir, '{}_pr_curve_full.csv.gz'.format(label))
114
+ pr.to_csv(full_pr_file, index=False)
115
+ pr = _filter_by_closeness(pr)
116
+
117
+ pr_file = os.path.join(outdir, '{}_pr_curve.csv'.format(label))
118
+ pr.to_csv(pr_file, index=False)
119
+
120
+ pr.plot(x='recall', y='precision', xlim=(0, 1), ylim=(0, 1))
121
+ pr_plot_file = os.path.join(outdir, '{}_pr.pdf'.format(label))
122
+ plt.savefig(pr_plot_file)
123
+ plt.close()
124
+
125
+ print(label, 'AuPR:', pr_auc, 'AvgP:', average_precision_score(y_true_binary, y_scores, sample_weight=sample_weight))
126
+
127
+ def get_flops(model):
128
+ concrete = tf.function(lambda inputs: model(inputs))
129
+ concrete_func = concrete.get_concrete_function(
130
+ [tf.TensorSpec([1, *inputs.shape[1:]]) for inputs in model.inputs])
131
+ frozen_func, graph_def = convert_variables_to_constants_v2_as_graph(concrete_func)
132
+ with tf.Graph().as_default() as graph:
133
+ tf.graph_util.import_graph_def(graph_def, name='')
134
+ run_meta = tf.compat.v1.RunMetadata()
135
+ opts = tf.compat.v1.profiler.ProfileOptionBuilder.float_operation()
136
+ flops = tf.compat.v1.profiler.profile(graph=graph, run_meta=run_meta, cmd="op", options=opts)
137
+
138
+ tf.compat.v1.reset_default_graph()
139
+ return flops.total_float_ops
140
+
141
+ def evaluate(exp, force=False):
142
+ ckpt_path = exp.path_to('best_model.keras')
143
+
144
+ # If the model doesn't exist and we're not forcing evaluation, just return
145
+ if not os.path.exists(ckpt_path):
146
+ print(f"Model file not found: {ckpt_path}")
147
+ return
148
+
149
+ custom_objects = {
150
+ 'AdaBeliefOptimizer': AdaBeliefOptimizer,
151
+ 'iou_coef': iou_coef,
152
+ 'dice_coef': dice_coef,
153
+ 'hard_swish': hard_swish,
154
+ 'enhanced_binary_crossentropy': enhanced_binary_crossentropy,
155
+ 'boundary_loss': boundary_loss
156
+ }
157
+
158
+ try:
159
+ model = tf.keras.models.load_model(ckpt_path, custom_objects=custom_objects)
160
+ except Exception as e:
161
+ print(f"Error loading model: {str(e)}")
162
+ return
163
+
164
+ # get flops
165
+ flop_params_path = exp.path_to('flops_nparams.csv')
166
+ if force or not os.path.exists(flop_params_path):
167
+ model.compile()
168
+ # Add .keras extension to tmp_model
169
+ tf.keras.models.save_model(model, 'tmp_model.keras', overwrite=True, include_optimizer=False)
170
+ stripped_model = tf.keras.models.load_model('tmp_model.keras')
171
+ flops = get_flops(stripped_model)
172
+ nparams = stripped_model.count_params()
173
+ del stripped_model
174
+ # Clean up temporary model file
175
+ if os.path.exists('tmp_model.keras'):
176
+ os.remove('tmp_model.keras')
177
+ print('FLOPS:', flops)
178
+ print('#PARAMS:', nparams)
179
+ pd.DataFrame({'flops': flops, 'nparams': nparams}, index=[0]).to_csv(flop_params_path)
180
+
181
+ model.compile(
182
+ loss='binary_crossentropy',
183
+ metrics={'mask': [iou_coef, dice_coef], 'tags': 'binary_accuracy'}
184
+ )
185
+
186
+ params = exp.params
187
+ np.random.seed(params.seed)
188
+ tf.random.set_seed(params.seed)
189
+
190
+ data = load_datasets(params.data)
191
+
192
+ # TRAIN/VAL/TEST SPLIT
193
+ if params.split == 'subjects': # by SUBJECTS
194
+ test_subjects = (3, 4, 19, 38, 45, 46, 51, 52)
195
+ test_data = data[data['sub'].isin(test_subjects)]
196
+ elif params.split == 'random': # 70-20-10 %
197
+ _, valtest_data = train_test_split(data, test_size=.3, shuffle=True)
198
+ _, test_data = train_test_split(valtest_data, test_size=.33)
199
+
200
+ x_shape = (params.resolution, params.resolution, 1)
201
+ test_gen, test_categories = get_loader(test_data, batch_size=1, x_shape=x_shape)
202
+
203
+ prediction_dir = exp.path_to('test_pred')
204
+ os.makedirs(prediction_dir, exist_ok=True)
205
+
206
+ loss_per_sample = None
207
+
208
+ def _get_test_predictions(test_gen, model):
209
+ x_masks = []
210
+ y_masks, y_tags = [], []
211
+ pred_masks, pred_tags = [], []
212
+ loss_per_sample = []
213
+
214
+ for x, y in tqdm(test_gen, desc='TEST'):
215
+ try:
216
+ # Generate predictions first
217
+ p_mask, p_tags = model.predict_on_batch(x)
218
+
219
+ # Append predictions and ground truth
220
+ pred_masks.append(p_mask)
221
+ pred_tags.append(p_tags)
222
+ y_masks.append(y['mask'].numpy())
223
+ y_tags.append(y['tags'].numpy())
224
+ x_masks.append(x.numpy())
225
+
226
+ # Calculate losses manually if needed
227
+ sample_loss = model.compiled_loss(
228
+ {'mask': y['mask'], 'tags': y['tags']},
229
+ {'mask': p_mask, 'tags': p_tags}
230
+ )
231
+ loss_per_sample.append(sample_loss.numpy())
232
+
233
+ except Exception as e:
234
+ print(f"Error processing batch: {str(e)}")
235
+ continue
236
+
237
+ # Check if we have any successful predictions
238
+ if not pred_masks or not pred_tags or not y_masks or not y_tags or not x_masks:
239
+ raise ValueError("No predictions were collected - all batches failed")
240
+
241
+ try:
242
+ loss_per_sample = np.array(loss_per_sample)
243
+ pred_masks = np.concatenate(pred_masks)
244
+ pred_tags = np.concatenate(pred_tags)
245
+ y_masks = np.concatenate(y_masks)
246
+ y_tags = np.concatenate(y_tags)
247
+ x_masks = np.concatenate(x_masks)
248
+ except Exception as e:
249
+ print(f"Error concatenating results: {str(e)}")
250
+ raise
251
+
252
+ return loss_per_sample, x_masks, y_masks, y_tags, pred_masks, pred_tags
253
+
254
+ mask_metrics_path = exp.path_to('test_pred/mask_metrics.csv')
255
+ if force or not os.path.exists(mask_metrics_path):
256
+ if loss_per_sample is None:
257
+ loss_per_sample, x_masks, y_masks, y_tags, pred_masks, pred_tags = _get_test_predictions(test_gen, model)
258
+
259
+ thrs = np.linspace(0, 1, 101)
260
+ ious = [iou_coef(y_masks, pred_masks).numpy() for thr in thrs]
261
+ dices = [dice_coef(y_masks, pred_masks).numpy() for thr in thrs]
262
+
263
+ best_thr = max(zip(dices, thrs))[1]
264
+
265
+ mask_metrics = pd.DataFrame({'iou': ious, 'dice': dices, 'thr': thrs})
266
+ print(mask_metrics.max(axis=0))
267
+ mask_metrics.to_csv(mask_metrics_path)
268
+ else:
269
+ mask_metrics = pd.read_csv(mask_metrics_path, index_col=0)
270
+ best_thr = mask_metrics.loc[mask_metrics.dice.idxmax(), 'thr']
271
+
272
+ if force:
273
+ if loss_per_sample is None:
274
+ loss_per_sample, x_masks, y_masks, y_tags, pred_masks, pred_tags = _get_test_predictions(test_gen, model)
275
+ _weighted_roc_pr(y_tags[:, 0], pred_tags[:, 0], 'all_eye', prediction_dir)
276
+ _weighted_roc_pr(y_tags[:, 1], pred_tags[:, 1], 'all_blink', prediction_dir)
277
+
278
+ filenames = ('top_samples.png', 'bottom_samples.png', 'random_samples.png')
279
+ if force or any(not os.path.exists(os.path.join(prediction_dir, f)) for f in filenames):
280
+ if loss_per_sample is None:
281
+ loss_per_sample, x_masks, y_masks, y_tags, pred_masks, pred_tags = _get_test_predictions(test_gen, model)
282
+
283
+ k = 5
284
+ best_selector = []
285
+ worst_selector = []
286
+ random_selector = []
287
+
288
+ idx = np.arange(len(test_data))
289
+ for cat in np.unique(test_categories):
290
+ cat_outdir = os.path.join(prediction_dir, cat)
291
+ os.makedirs(cat_outdir, exist_ok=True)
292
+
293
+ selector = test_categories == cat
294
+ _weighted_roc_pr(y_tags[selector, 0], pred_tags[selector, 0], '{}_eye'.format(cat), cat_outdir)
295
+ _weighted_roc_pr(y_tags[selector, 1], pred_tags[selector, 1], '{}_blink'.format(cat), cat_outdir)
296
+
297
+ # Fix the indexing here - loss_per_sample is 1D
298
+ cat_losses = loss_per_sample[selector]
299
+ rank = cat_losses.argsort()
300
+ topk, bottomk = rank[:k], rank[-k:]
301
+
302
+ best_selector += idx[selector][topk].tolist()
303
+ worst_selector += idx[selector][bottomk].tolist()
304
+ random_selector += np.random.choice(idx[selector], k, replace=False).tolist()
305
+
306
+ # topk-bottomk images
307
+ selectors = (best_selector, worst_selector, random_selector)
308
+ for selector, outfile in zip(selectors, filenames):
309
+ combined_m = np.concatenate((pred_masks[selector], y_masks[selector]), axis=-1)[:, :, :, ::-1]
310
+ combined_t = np.concatenate((pred_tags[selector], y_tags[selector]), axis=-1)
311
+ combined_y = (combined_m, combined_t)
312
+ out = os.path.join(prediction_dir, outfile)
313
+ visualize(x_masks[selector], combined_y, out=out, thr=best_thr, n_cols=k, width=10)
314
+
315
+ for i, (xi, yi_mask) in enumerate(zip(x_masks[selector], combined_m)):
316
+ img = visualizable(xi, yi_mask, thr=best_thr)
317
+ img = (img * 255).astype(np.uint8)
318
+ out = os.path.join(prediction_dir, outfile[:-4])
319
+ os.makedirs(out, exist_ok=True)
320
+ out = os.path.join(out, f'{i}.png')
321
+ Image.fromarray(img).save(out)
322
+
323
+ def main(args):
324
+ try:
325
+ for exp in expman.gather(args.run).filter(args.filter):
326
+ print(exp)
327
+ evaluate(exp, force=args.force)
328
+ except Exception as e:
329
+ print(f"Error in main: {str(e)}")
330
+ raise
331
+
332
+ if __name__ == '__main__':
333
+ parser = argparse.ArgumentParser(description='Evaluate Run')
334
+ parser.add_argument('run', help='Run(s) directory')
335
+ parser.add_argument('-f', '--filter', default={}, type=expman.exp_filter)
336
+ parser.add_argument('--force', default=False, action='store_true', help='Force metrics recomputation')
337
+
338
+ args = parser.parse_args()
339
+ main(args)
models/__pycache__/unet.cpython-311.pyc ADDED
Binary file (5.46 kB). View file
 
models/deeplab/inference/tflite/stb_image.h ADDED
The diff for this file is too large to render. See raw diff
 
models/deeplab/inference/tflite/stb_image_resize.h ADDED
The diff for this file is too large to render. See raw diff
 
models/deeplab/inference/tflite/stb_image_write.h ADDED
@@ -0,0 +1,1622 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ /* stb_image_write - v1.13 - public domain - http://nothings.org/stb/stb_image_write.h
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
+ If using a modern Microsoft Compiler, non-safe versions of CRT calls may cause
14
+ compilation warnings or even errors. To avoid this, also before #including,
15
+
16
+ #define STBI_MSC_SECURE_CRT
17
+
18
+ ABOUT:
19
+
20
+ This header file is a library for writing images to C stdio or a callback.
21
+
22
+ The PNG output is not optimal; it is 20-50% larger than the file
23
+ written by a decent optimizing implementation; though providing a custom
24
+ zlib compress function (see STBIW_ZLIB_COMPRESS) can mitigate that.
25
+ This library is designed for source code compactness and simplicity,
26
+ not optimal image file size or run-time performance.
27
+
28
+ BUILDING:
29
+
30
+ You can #define STBIW_ASSERT(x) before the #include to avoid using assert.h.
31
+ You can #define STBIW_MALLOC(), STBIW_REALLOC(), and STBIW_FREE() to replace
32
+ malloc,realloc,free.
33
+ You can #define STBIW_MEMMOVE() to replace memmove()
34
+ You can #define STBIW_ZLIB_COMPRESS to use a custom zlib-style compress function
35
+ for PNG compression (instead of the builtin one), it must have the following signature:
36
+ unsigned char * my_compress(unsigned char *data, int data_len, int *out_len, int quality);
37
+ The returned data will be freed with STBIW_FREE() (free() by default),
38
+ so it must be heap allocated with STBIW_MALLOC() (malloc() by default),
39
+
40
+ UNICODE:
41
+
42
+ If compiling for Windows and you wish to use Unicode filenames, compile
43
+ with
44
+ #define STBIW_WINDOWS_UTF8
45
+ and pass utf8-encoded filenames. Call stbiw_convert_wchar_to_utf8 to convert
46
+ Windows wchar_t filenames to utf8.
47
+
48
+ USAGE:
49
+
50
+ There are five functions, one for each image file format:
51
+
52
+ int stbi_write_png(char const *filename, int w, int h, int comp, const void *data, int stride_in_bytes);
53
+ int stbi_write_bmp(char const *filename, int w, int h, int comp, const void *data);
54
+ int stbi_write_tga(char const *filename, int w, int h, int comp, const void *data);
55
+ int stbi_write_jpg(char const *filename, int w, int h, int comp, const void *data, int quality);
56
+ int stbi_write_hdr(char const *filename, int w, int h, int comp, const float *data);
57
+
58
+ void stbi_flip_vertically_on_write(int flag); // flag is non-zero to flip data vertically
59
+
60
+ There are also five equivalent functions that use an arbitrary write function. You are
61
+ expected to open/close your file-equivalent before and after calling these:
62
+
63
+ 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);
64
+ int stbi_write_bmp_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const void *data);
65
+ int stbi_write_tga_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const void *data);
66
+ int stbi_write_hdr_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const float *data);
67
+ int stbi_write_jpg_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data, int quality);
68
+
69
+ where the callback is:
70
+ void stbi_write_func(void *context, void *data, int size);
71
+
72
+ You can configure it with these global variables:
73
+ int stbi_write_tga_with_rle; // defaults to true; set to 0 to disable RLE
74
+ int stbi_write_png_compression_level; // defaults to 8; set to higher for more compression
75
+ int stbi_write_force_png_filter; // defaults to -1; set to 0..5 to force a filter mode
76
+
77
+
78
+ You can define STBI_WRITE_NO_STDIO to disable the file variant of these
79
+ functions, so the library will not use stdio.h at all. However, this will
80
+ also disable HDR writing, because it requires stdio for formatted output.
81
+
82
+ Each function returns 0 on failure and non-0 on success.
83
+
84
+ The functions create an image file defined by the parameters. The image
85
+ is a rectangle of pixels stored from left-to-right, top-to-bottom.
86
+ Each pixel contains 'comp' channels of data stored interleaved with 8-bits
87
+ per channel, in the following order: 1=Y, 2=YA, 3=RGB, 4=RGBA. (Y is
88
+ monochrome color.) The rectangle is 'w' pixels wide and 'h' pixels tall.
89
+ The *data pointer points to the first byte of the top-left-most pixel.
90
+ For PNG, "stride_in_bytes" is the distance in bytes from the first byte of
91
+ a row of pixels to the first byte of the next row of pixels.
92
+
93
+ PNG creates output files with the same number of components as the input.
94
+ The BMP format expands Y to RGB in the file format and does not
95
+ output alpha.
96
+
97
+ PNG supports writing rectangles of data even when the bytes storing rows of
98
+ data are not consecutive in memory (e.g. sub-rectangles of a larger image),
99
+ by supplying the stride between the beginning of adjacent rows. The other
100
+ formats do not. (Thus you cannot write a native-format BMP through the BMP
101
+ writer, both because it is in BGR order and because it may have padding
102
+ at the end of the line.)
103
+
104
+ PNG allows you to set the deflate compression level by setting the global
105
+ variable 'stbi_write_png_compression_level' (it defaults to 8).
106
+
107
+ HDR expects linear float data. Since the format is always 32-bit rgb(e)
108
+ data, alpha (if provided) is discarded, and for monochrome data it is
109
+ replicated across all three channels.
110
+
111
+ TGA supports RLE or non-RLE compressed data. To use non-RLE-compressed
112
+ data, set the global variable 'stbi_write_tga_with_rle' to 0.
113
+
114
+ JPEG does ignore alpha channels in input data; quality is between 1 and 100.
115
+ Higher quality looks better but results in a bigger image.
116
+ JPEG baseline (no JPEG progressive).
117
+
118
+ CREDITS:
119
+
120
+
121
+ Sean Barrett - PNG/BMP/TGA
122
+ Baldur Karlsson - HDR
123
+ Jean-Sebastien Guay - TGA monochrome
124
+ Tim Kelsey - misc enhancements
125
+ Alan Hickman - TGA RLE
126
+ Emmanuel Julien - initial file IO callback implementation
127
+ Jon Olick - original jo_jpeg.cpp code
128
+ Daniel Gibson - integrate JPEG, allow external zlib
129
+ Aarni Koskela - allow choosing PNG filter
130
+
131
+ bugfixes:
132
+ github:Chribba
133
+ Guillaume Chereau
134
+ github:jry2
135
+ github:romigrou
136
+ Sergio Gonzalez
137
+ Jonas Karlsson
138
+ Filip Wasil
139
+ Thatcher Ulrich
140
+ github:poppolopoppo
141
+ Patrick Boettcher
142
+ github:xeekworx
143
+ Cap Petschulat
144
+ Simon Rodriguez
145
+ Ivan Tikhonov
146
+ github:ignotion
147
+ Adam Schackart
148
+
149
+ LICENSE
150
+
151
+ See end of file for license information.
152
+
153
+ */
154
+
155
+ #ifndef INCLUDE_STB_IMAGE_WRITE_H
156
+ #define INCLUDE_STB_IMAGE_WRITE_H
157
+
158
+ #include <stdlib.h>
159
+
160
+ // if STB_IMAGE_WRITE_STATIC causes problems, try defining STBIWDEF to 'inline' or 'static inline'
161
+ #ifndef STBIWDEF
162
+ #ifdef STB_IMAGE_WRITE_STATIC
163
+ #define STBIWDEF static
164
+ #else
165
+ #ifdef __cplusplus
166
+ #define STBIWDEF extern "C"
167
+ #else
168
+ #define STBIWDEF extern
169
+ #endif
170
+ #endif
171
+ #endif
172
+
173
+ #ifndef STB_IMAGE_WRITE_STATIC // C++ forbids static forward declarations
174
+ extern int stbi_write_tga_with_rle;
175
+ extern int stbi_write_png_compression_level;
176
+ extern int stbi_write_force_png_filter;
177
+ #endif
178
+
179
+ #ifndef STBI_WRITE_NO_STDIO
180
+ STBIWDEF int stbi_write_png(char const *filename, int w, int h, int comp, const void *data, int stride_in_bytes);
181
+ STBIWDEF int stbi_write_bmp(char const *filename, int w, int h, int comp, const void *data);
182
+ STBIWDEF int stbi_write_tga(char const *filename, int w, int h, int comp, const void *data);
183
+ STBIWDEF int stbi_write_hdr(char const *filename, int w, int h, int comp, const float *data);
184
+ STBIWDEF int stbi_write_jpg(char const *filename, int x, int y, int comp, const void *data, int quality);
185
+
186
+ #ifdef STBI_WINDOWS_UTF8
187
+ STBIWDEF int stbiw_convert_wchar_to_utf8(char *buffer, size_t bufferlen, const wchar_t* input);
188
+ #endif
189
+ #endif
190
+
191
+ typedef void stbi_write_func(void *context, void *data, int size);
192
+
193
+ 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);
194
+ STBIWDEF int stbi_write_bmp_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const void *data);
195
+ STBIWDEF int stbi_write_tga_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const void *data);
196
+ STBIWDEF int stbi_write_hdr_to_func(stbi_write_func *func, void *context, int w, int h, int comp, const float *data);
197
+ STBIWDEF int stbi_write_jpg_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data, int quality);
198
+
199
+ STBIWDEF void stbi_flip_vertically_on_write(int flip_boolean);
200
+
201
+ #endif//INCLUDE_STB_IMAGE_WRITE_H
202
+
203
+ #ifdef STB_IMAGE_WRITE_IMPLEMENTATION
204
+
205
+ #ifdef _WIN32
206
+ #ifndef _CRT_SECURE_NO_WARNINGS
207
+ #define _CRT_SECURE_NO_WARNINGS
208
+ #endif
209
+ #ifndef _CRT_NONSTDC_NO_DEPRECATE
210
+ #define _CRT_NONSTDC_NO_DEPRECATE
211
+ #endif
212
+ #endif
213
+
214
+ #ifndef STBI_WRITE_NO_STDIO
215
+ #include <stdio.h>
216
+ #endif // STBI_WRITE_NO_STDIO
217
+
218
+ #include <stdarg.h>
219
+ #include <stdlib.h>
220
+ #include <string.h>
221
+ #include <math.h>
222
+
223
+ #if defined(STBIW_MALLOC) && defined(STBIW_FREE) && (defined(STBIW_REALLOC) || defined(STBIW_REALLOC_SIZED))
224
+ // ok
225
+ #elif !defined(STBIW_MALLOC) && !defined(STBIW_FREE) && !defined(STBIW_REALLOC) && !defined(STBIW_REALLOC_SIZED)
226
+ // ok
227
+ #else
228
+ #error "Must define all or none of STBIW_MALLOC, STBIW_FREE, and STBIW_REALLOC (or STBIW_REALLOC_SIZED)."
229
+ #endif
230
+
231
+ #ifndef STBIW_MALLOC
232
+ #define STBIW_MALLOC(sz) malloc(sz)
233
+ #define STBIW_REALLOC(p,newsz) realloc(p,newsz)
234
+ #define STBIW_FREE(p) free(p)
235
+ #endif
236
+
237
+ #ifndef STBIW_REALLOC_SIZED
238
+ #define STBIW_REALLOC_SIZED(p,oldsz,newsz) STBIW_REALLOC(p,newsz)
239
+ #endif
240
+
241
+
242
+ #ifndef STBIW_MEMMOVE
243
+ #define STBIW_MEMMOVE(a,b,sz) memmove(a,b,sz)
244
+ #endif
245
+
246
+
247
+ #ifndef STBIW_ASSERT
248
+ #include <assert.h>
249
+ #define STBIW_ASSERT(x) assert(x)
250
+ #endif
251
+
252
+ #define STBIW_UCHAR(x) (unsigned char) ((x) & 0xff)
253
+
254
+ #ifdef STB_IMAGE_WRITE_STATIC
255
+ static int stbi__flip_vertically_on_write=0;
256
+ static int stbi_write_png_compression_level = 8;
257
+ static int stbi_write_tga_with_rle = 1;
258
+ static int stbi_write_force_png_filter = -1;
259
+ #else
260
+ int stbi_write_png_compression_level = 8;
261
+ int stbi__flip_vertically_on_write=0;
262
+ int stbi_write_tga_with_rle = 1;
263
+ int stbi_write_force_png_filter = -1;
264
+ #endif
265
+
266
+ STBIWDEF void stbi_flip_vertically_on_write(int flag)
267
+ {
268
+ stbi__flip_vertically_on_write = flag;
269
+ }
270
+
271
+ typedef struct
272
+ {
273
+ stbi_write_func *func;
274
+ void *context;
275
+ } stbi__write_context;
276
+
277
+ // initialize a callback-based context
278
+ static void stbi__start_write_callbacks(stbi__write_context *s, stbi_write_func *c, void *context)
279
+ {
280
+ s->func = c;
281
+ s->context = context;
282
+ }
283
+
284
+ #ifndef STBI_WRITE_NO_STDIO
285
+
286
+ static void stbi__stdio_write(void *context, void *data, int size)
287
+ {
288
+ fwrite(data,1,size,(FILE*) context);
289
+ }
290
+
291
+ #if defined(_MSC_VER) && defined(STBI_WINDOWS_UTF8)
292
+ #ifdef __cplusplus
293
+ #define STBIW_EXTERN extern "C"
294
+ #else
295
+ #define STBIW_EXTERN extern
296
+ #endif
297
+ STBIW_EXTERN __declspec(dllimport) int __stdcall MultiByteToWideChar(unsigned int cp, unsigned long flags, const char *str, int cbmb, wchar_t *widestr, int cchwide);
298
+ 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);
299
+
300
+ STBIWDEF int stbiw_convert_wchar_to_utf8(char *buffer, size_t bufferlen, const wchar_t* input)
301
+ {
302
+ return WideCharToMultiByte(65001 /* UTF8 */, 0, input, -1, buffer, (int) bufferlen, NULL, NULL);
303
+ }
304
+ #endif
305
+
306
+ static FILE *stbiw__fopen(char const *filename, char const *mode)
307
+ {
308
+ FILE *f;
309
+ #if defined(_MSC_VER) && defined(STBI_WINDOWS_UTF8)
310
+ wchar_t wMode[64];
311
+ wchar_t wFilename[1024];
312
+ if (0 == MultiByteToWideChar(65001 /* UTF8 */, 0, filename, -1, wFilename, sizeof(wFilename)))
313
+ return 0;
314
+
315
+ if (0 == MultiByteToWideChar(65001 /* UTF8 */, 0, mode, -1, wMode, sizeof(wMode)))
316
+ return 0;
317
+
318
+ #if _MSC_VER >= 1400
319
+ if (0 != _wfopen_s(&f, wFilename, wMode))
320
+ f = 0;
321
+ #else
322
+ f = _wfopen(wFilename, wMode);
323
+ #endif
324
+
325
+ #elif defined(_MSC_VER) && _MSC_VER >= 1400
326
+ if (0 != fopen_s(&f, filename, mode))
327
+ f=0;
328
+ #else
329
+ f = fopen(filename, mode);
330
+ #endif
331
+ return f;
332
+ }
333
+
334
+ static int stbi__start_write_file(stbi__write_context *s, const char *filename)
335
+ {
336
+ FILE *f = stbiw__fopen(filename, "wb");
337
+ stbi__start_write_callbacks(s, stbi__stdio_write, (void *) f);
338
+ return f != NULL;
339
+ }
340
+
341
+ static void stbi__end_write_file(stbi__write_context *s)
342
+ {
343
+ fclose((FILE *)s->context);
344
+ }
345
+
346
+ #endif // !STBI_WRITE_NO_STDIO
347
+
348
+ typedef unsigned int stbiw_uint32;
349
+ typedef int stb_image_write_test[sizeof(stbiw_uint32)==4 ? 1 : -1];
350
+
351
+ static void stbiw__writefv(stbi__write_context *s, const char *fmt, va_list v)
352
+ {
353
+ while (*fmt) {
354
+ switch (*fmt++) {
355
+ case ' ': break;
356
+ case '1': { unsigned char x = STBIW_UCHAR(va_arg(v, int));
357
+ s->func(s->context,&x,1);
358
+ break; }
359
+ case '2': { int x = va_arg(v,int);
360
+ unsigned char b[2];
361
+ b[0] = STBIW_UCHAR(x);
362
+ b[1] = STBIW_UCHAR(x>>8);
363
+ s->func(s->context,b,2);
364
+ break; }
365
+ case '4': { stbiw_uint32 x = va_arg(v,int);
366
+ unsigned char b[4];
367
+ b[0]=STBIW_UCHAR(x);
368
+ b[1]=STBIW_UCHAR(x>>8);
369
+ b[2]=STBIW_UCHAR(x>>16);
370
+ b[3]=STBIW_UCHAR(x>>24);
371
+ s->func(s->context,b,4);
372
+ break; }
373
+ default:
374
+ STBIW_ASSERT(0);
375
+ return;
376
+ }
377
+ }
378
+ }
379
+
380
+ static void stbiw__writef(stbi__write_context *s, const char *fmt, ...)
381
+ {
382
+ va_list v;
383
+ va_start(v, fmt);
384
+ stbiw__writefv(s, fmt, v);
385
+ va_end(v);
386
+ }
387
+
388
+ static void stbiw__putc(stbi__write_context *s, unsigned char c)
389
+ {
390
+ s->func(s->context, &c, 1);
391
+ }
392
+
393
+ static void stbiw__write3(stbi__write_context *s, unsigned char a, unsigned char b, unsigned char c)
394
+ {
395
+ unsigned char arr[3];
396
+ arr[0] = a; arr[1] = b; arr[2] = c;
397
+ s->func(s->context, arr, 3);
398
+ }
399
+
400
+ static void stbiw__write_pixel(stbi__write_context *s, int rgb_dir, int comp, int write_alpha, int expand_mono, unsigned char *d)
401
+ {
402
+ unsigned char bg[3] = { 255, 0, 255}, px[3];
403
+ int k;
404
+
405
+ if (write_alpha < 0)
406
+ s->func(s->context, &d[comp - 1], 1);
407
+
408
+ switch (comp) {
409
+ case 2: // 2 pixels = mono + alpha, alpha is written separately, so same as 1-channel case
410
+ case 1:
411
+ if (expand_mono)
412
+ stbiw__write3(s, d[0], d[0], d[0]); // monochrome bmp
413
+ else
414
+ s->func(s->context, d, 1); // monochrome TGA
415
+ break;
416
+ case 4:
417
+ if (!write_alpha) {
418
+ // composite against pink background
419
+ for (k = 0; k < 3; ++k)
420
+ px[k] = bg[k] + ((d[k] - bg[k]) * d[3]) / 255;
421
+ stbiw__write3(s, px[1 - rgb_dir], px[1], px[1 + rgb_dir]);
422
+ break;
423
+ }
424
+ /* FALLTHROUGH */
425
+ case 3:
426
+ stbiw__write3(s, d[1 - rgb_dir], d[1], d[1 + rgb_dir]);
427
+ break;
428
+ }
429
+ if (write_alpha > 0)
430
+ s->func(s->context, &d[comp - 1], 1);
431
+ }
432
+
433
+ 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)
434
+ {
435
+ stbiw_uint32 zero = 0;
436
+ int i,j, j_end;
437
+
438
+ if (y <= 0)
439
+ return;
440
+
441
+ if (stbi__flip_vertically_on_write)
442
+ vdir *= -1;
443
+
444
+ if (vdir < 0) {
445
+ j_end = -1; j = y-1;
446
+ } else {
447
+ j_end = y; j = 0;
448
+ }
449
+
450
+ for (; j != j_end; j += vdir) {
451
+ for (i=0; i < x; ++i) {
452
+ unsigned char *d = (unsigned char *) data + (j*x+i)*comp;
453
+ stbiw__write_pixel(s, rgb_dir, comp, write_alpha, expand_mono, d);
454
+ }
455
+ s->func(s->context, &zero, scanline_pad);
456
+ }
457
+ }
458
+
459
+ 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, ...)
460
+ {
461
+ if (y < 0 || x < 0) {
462
+ return 0;
463
+ } else {
464
+ va_list v;
465
+ va_start(v, fmt);
466
+ stbiw__writefv(s, fmt, v);
467
+ va_end(v);
468
+ stbiw__write_pixels(s,rgb_dir,vdir,x,y,comp,data,alpha,pad, expand_mono);
469
+ return 1;
470
+ }
471
+ }
472
+
473
+ static int stbi_write_bmp_core(stbi__write_context *s, int x, int y, int comp, const void *data)
474
+ {
475
+ int pad = (-x*3) & 3;
476
+ return stbiw__outfile(s,-1,-1,x,y,comp,1,(void *) data,0,pad,
477
+ "11 4 22 4" "4 44 22 444444",
478
+ 'B', 'M', 14+40+(x*3+pad)*y, 0,0, 14+40, // file header
479
+ 40, x,y, 1,24, 0,0,0,0,0,0); // bitmap header
480
+ }
481
+
482
+ STBIWDEF int stbi_write_bmp_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data)
483
+ {
484
+ stbi__write_context s;
485
+ stbi__start_write_callbacks(&s, func, context);
486
+ return stbi_write_bmp_core(&s, x, y, comp, data);
487
+ }
488
+
489
+ #ifndef STBI_WRITE_NO_STDIO
490
+ STBIWDEF int stbi_write_bmp(char const *filename, int x, int y, int comp, const void *data)
491
+ {
492
+ stbi__write_context s;
493
+ if (stbi__start_write_file(&s,filename)) {
494
+ int r = stbi_write_bmp_core(&s, x, y, comp, data);
495
+ stbi__end_write_file(&s);
496
+ return r;
497
+ } else
498
+ return 0;
499
+ }
500
+ #endif //!STBI_WRITE_NO_STDIO
501
+
502
+ static int stbi_write_tga_core(stbi__write_context *s, int x, int y, int comp, void *data)
503
+ {
504
+ int has_alpha = (comp == 2 || comp == 4);
505
+ int colorbytes = has_alpha ? comp-1 : comp;
506
+ int format = colorbytes < 2 ? 3 : 2; // 3 color channels (RGB/RGBA) = 2, 1 color channel (Y/YA) = 3
507
+
508
+ if (y < 0 || x < 0)
509
+ return 0;
510
+
511
+ if (!stbi_write_tga_with_rle) {
512
+ return stbiw__outfile(s, -1, -1, x, y, comp, 0, (void *) data, has_alpha, 0,
513
+ "111 221 2222 11", 0, 0, format, 0, 0, 0, 0, 0, x, y, (colorbytes + has_alpha) * 8, has_alpha * 8);
514
+ } else {
515
+ int i,j,k;
516
+ int jend, jdir;
517
+
518
+ 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);
519
+
520
+ if (stbi__flip_vertically_on_write) {
521
+ j = 0;
522
+ jend = y;
523
+ jdir = 1;
524
+ } else {
525
+ j = y-1;
526
+ jend = -1;
527
+ jdir = -1;
528
+ }
529
+ for (; j != jend; j += jdir) {
530
+ unsigned char *row = (unsigned char *) data + j * x * comp;
531
+ int len;
532
+
533
+ for (i = 0; i < x; i += len) {
534
+ unsigned char *begin = row + i * comp;
535
+ int diff = 1;
536
+ len = 1;
537
+
538
+ if (i < x - 1) {
539
+ ++len;
540
+ diff = memcmp(begin, row + (i + 1) * comp, comp);
541
+ if (diff) {
542
+ const unsigned char *prev = begin;
543
+ for (k = i + 2; k < x && len < 128; ++k) {
544
+ if (memcmp(prev, row + k * comp, comp)) {
545
+ prev += comp;
546
+ ++len;
547
+ } else {
548
+ --len;
549
+ break;
550
+ }
551
+ }
552
+ } else {
553
+ for (k = i + 2; k < x && len < 128; ++k) {
554
+ if (!memcmp(begin, row + k * comp, comp)) {
555
+ ++len;
556
+ } else {
557
+ break;
558
+ }
559
+ }
560
+ }
561
+ }
562
+
563
+ if (diff) {
564
+ unsigned char header = STBIW_UCHAR(len - 1);
565
+ s->func(s->context, &header, 1);
566
+ for (k = 0; k < len; ++k) {
567
+ stbiw__write_pixel(s, -1, comp, has_alpha, 0, begin + k * comp);
568
+ }
569
+ } else {
570
+ unsigned char header = STBIW_UCHAR(len - 129);
571
+ s->func(s->context, &header, 1);
572
+ stbiw__write_pixel(s, -1, comp, has_alpha, 0, begin);
573
+ }
574
+ }
575
+ }
576
+ }
577
+ return 1;
578
+ }
579
+
580
+ STBIWDEF int stbi_write_tga_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data)
581
+ {
582
+ stbi__write_context s;
583
+ stbi__start_write_callbacks(&s, func, context);
584
+ return stbi_write_tga_core(&s, x, y, comp, (void *) data);
585
+ }
586
+
587
+ #ifndef STBI_WRITE_NO_STDIO
588
+ STBIWDEF int stbi_write_tga(char const *filename, int x, int y, int comp, const void *data)
589
+ {
590
+ stbi__write_context s;
591
+ if (stbi__start_write_file(&s,filename)) {
592
+ int r = stbi_write_tga_core(&s, x, y, comp, (void *) data);
593
+ stbi__end_write_file(&s);
594
+ return r;
595
+ } else
596
+ return 0;
597
+ }
598
+ #endif
599
+
600
+ // *************************************************************************************************
601
+ // Radiance RGBE HDR writer
602
+ // by Baldur Karlsson
603
+
604
+ #define stbiw__max(a, b) ((a) > (b) ? (a) : (b))
605
+
606
+ static void stbiw__linear_to_rgbe(unsigned char *rgbe, float *linear)
607
+ {
608
+ int exponent;
609
+ float maxcomp = stbiw__max(linear[0], stbiw__max(linear[1], linear[2]));
610
+
611
+ if (maxcomp < 1e-32f) {
612
+ rgbe[0] = rgbe[1] = rgbe[2] = rgbe[3] = 0;
613
+ } else {
614
+ float normalize = (float) frexp(maxcomp, &exponent) * 256.0f/maxcomp;
615
+
616
+ rgbe[0] = (unsigned char)(linear[0] * normalize);
617
+ rgbe[1] = (unsigned char)(linear[1] * normalize);
618
+ rgbe[2] = (unsigned char)(linear[2] * normalize);
619
+ rgbe[3] = (unsigned char)(exponent + 128);
620
+ }
621
+ }
622
+
623
+ static void stbiw__write_run_data(stbi__write_context *s, int length, unsigned char databyte)
624
+ {
625
+ unsigned char lengthbyte = STBIW_UCHAR(length+128);
626
+ STBIW_ASSERT(length+128 <= 255);
627
+ s->func(s->context, &lengthbyte, 1);
628
+ s->func(s->context, &databyte, 1);
629
+ }
630
+
631
+ static void stbiw__write_dump_data(stbi__write_context *s, int length, unsigned char *data)
632
+ {
633
+ unsigned char lengthbyte = STBIW_UCHAR(length);
634
+ STBIW_ASSERT(length <= 128); // inconsistent with spec but consistent with official code
635
+ s->func(s->context, &lengthbyte, 1);
636
+ s->func(s->context, data, length);
637
+ }
638
+
639
+ static void stbiw__write_hdr_scanline(stbi__write_context *s, int width, int ncomp, unsigned char *scratch, float *scanline)
640
+ {
641
+ unsigned char scanlineheader[4] = { 2, 2, 0, 0 };
642
+ unsigned char rgbe[4];
643
+ float linear[3];
644
+ int x;
645
+
646
+ scanlineheader[2] = (width&0xff00)>>8;
647
+ scanlineheader[3] = (width&0x00ff);
648
+
649
+ /* skip RLE for images too small or large */
650
+ if (width < 8 || width >= 32768) {
651
+ for (x=0; x < width; x++) {
652
+ switch (ncomp) {
653
+ case 4: /* fallthrough */
654
+ case 3: linear[2] = scanline[x*ncomp + 2];
655
+ linear[1] = scanline[x*ncomp + 1];
656
+ linear[0] = scanline[x*ncomp + 0];
657
+ break;
658
+ default:
659
+ linear[0] = linear[1] = linear[2] = scanline[x*ncomp + 0];
660
+ break;
661
+ }
662
+ stbiw__linear_to_rgbe(rgbe, linear);
663
+ s->func(s->context, rgbe, 4);
664
+ }
665
+ } else {
666
+ int c,r;
667
+ /* encode into scratch buffer */
668
+ for (x=0; x < width; x++) {
669
+ switch(ncomp) {
670
+ case 4: /* fallthrough */
671
+ case 3: linear[2] = scanline[x*ncomp + 2];
672
+ linear[1] = scanline[x*ncomp + 1];
673
+ linear[0] = scanline[x*ncomp + 0];
674
+ break;
675
+ default:
676
+ linear[0] = linear[1] = linear[2] = scanline[x*ncomp + 0];
677
+ break;
678
+ }
679
+ stbiw__linear_to_rgbe(rgbe, linear);
680
+ scratch[x + width*0] = rgbe[0];
681
+ scratch[x + width*1] = rgbe[1];
682
+ scratch[x + width*2] = rgbe[2];
683
+ scratch[x + width*3] = rgbe[3];
684
+ }
685
+
686
+ s->func(s->context, scanlineheader, 4);
687
+
688
+ /* RLE each component separately */
689
+ for (c=0; c < 4; c++) {
690
+ unsigned char *comp = &scratch[width*c];
691
+
692
+ x = 0;
693
+ while (x < width) {
694
+ // find first run
695
+ r = x;
696
+ while (r+2 < width) {
697
+ if (comp[r] == comp[r+1] && comp[r] == comp[r+2])
698
+ break;
699
+ ++r;
700
+ }
701
+ if (r+2 >= width)
702
+ r = width;
703
+ // dump up to first run
704
+ while (x < r) {
705
+ int len = r-x;
706
+ if (len > 128) len = 128;
707
+ stbiw__write_dump_data(s, len, &comp[x]);
708
+ x += len;
709
+ }
710
+ // if there's a run, output it
711
+ if (r+2 < width) { // same test as what we break out of in search loop, so only true if we break'd
712
+ // find next byte after run
713
+ while (r < width && comp[r] == comp[x])
714
+ ++r;
715
+ // output run up to r
716
+ while (x < r) {
717
+ int len = r-x;
718
+ if (len > 127) len = 127;
719
+ stbiw__write_run_data(s, len, comp[x]);
720
+ x += len;
721
+ }
722
+ }
723
+ }
724
+ }
725
+ }
726
+ }
727
+
728
+ static int stbi_write_hdr_core(stbi__write_context *s, int x, int y, int comp, float *data)
729
+ {
730
+ if (y <= 0 || x <= 0 || data == NULL)
731
+ return 0;
732
+ else {
733
+ // Each component is stored separately. Allocate scratch space for full output scanline.
734
+ unsigned char *scratch = (unsigned char *) STBIW_MALLOC(x*4);
735
+ int i, len;
736
+ char buffer[128];
737
+ char header[] = "#?RADIANCE\n# Written by stb_image_write.h\nFORMAT=32-bit_rle_rgbe\n";
738
+ s->func(s->context, header, sizeof(header)-1);
739
+
740
+ #ifdef __STDC_WANT_SECURE_LIB__
741
+ len = sprintf_s(buffer, sizeof(buffer), "EXPOSURE= 1.0000000000000\n\n-Y %d +X %d\n", y, x);
742
+ #else
743
+ len = sprintf(buffer, "EXPOSURE= 1.0000000000000\n\n-Y %d +X %d\n", y, x);
744
+ #endif
745
+ s->func(s->context, buffer, len);
746
+
747
+ for(i=0; i < y; i++)
748
+ stbiw__write_hdr_scanline(s, x, comp, scratch, data + comp*x*(stbi__flip_vertically_on_write ? y-1-i : i));
749
+ STBIW_FREE(scratch);
750
+ return 1;
751
+ }
752
+ }
753
+
754
+ STBIWDEF int stbi_write_hdr_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const float *data)
755
+ {
756
+ stbi__write_context s;
757
+ stbi__start_write_callbacks(&s, func, context);
758
+ return stbi_write_hdr_core(&s, x, y, comp, (float *) data);
759
+ }
760
+
761
+ #ifndef STBI_WRITE_NO_STDIO
762
+ STBIWDEF int stbi_write_hdr(char const *filename, int x, int y, int comp, const float *data)
763
+ {
764
+ stbi__write_context s;
765
+ if (stbi__start_write_file(&s,filename)) {
766
+ int r = stbi_write_hdr_core(&s, x, y, comp, (float *) data);
767
+ stbi__end_write_file(&s);
768
+ return r;
769
+ } else
770
+ return 0;
771
+ }
772
+ #endif // STBI_WRITE_NO_STDIO
773
+
774
+
775
+ //////////////////////////////////////////////////////////////////////////////
776
+ //
777
+ // PNG writer
778
+ //
779
+
780
+ #ifndef STBIW_ZLIB_COMPRESS
781
+ // stretchy buffer; stbiw__sbpush() == vector<>::push_back() -- stbiw__sbcount() == vector<>::size()
782
+ #define stbiw__sbraw(a) ((int *) (a) - 2)
783
+ #define stbiw__sbm(a) stbiw__sbraw(a)[0]
784
+ #define stbiw__sbn(a) stbiw__sbraw(a)[1]
785
+
786
+ #define stbiw__sbneedgrow(a,n) ((a)==0 || stbiw__sbn(a)+n >= stbiw__sbm(a))
787
+ #define stbiw__sbmaybegrow(a,n) (stbiw__sbneedgrow(a,(n)) ? stbiw__sbgrow(a,n) : 0)
788
+ #define stbiw__sbgrow(a,n) stbiw__sbgrowf((void **) &(a), (n), sizeof(*(a)))
789
+
790
+ #define stbiw__sbpush(a, v) (stbiw__sbmaybegrow(a,1), (a)[stbiw__sbn(a)++] = (v))
791
+ #define stbiw__sbcount(a) ((a) ? stbiw__sbn(a) : 0)
792
+ #define stbiw__sbfree(a) ((a) ? STBIW_FREE(stbiw__sbraw(a)),0 : 0)
793
+
794
+ static void *stbiw__sbgrowf(void **arr, int increment, int itemsize)
795
+ {
796
+ int m = *arr ? 2*stbiw__sbm(*arr)+increment : increment+1;
797
+ void *p = STBIW_REALLOC_SIZED(*arr ? stbiw__sbraw(*arr) : 0, *arr ? (stbiw__sbm(*arr)*itemsize + sizeof(int)*2) : 0, itemsize * m + sizeof(int)*2);
798
+ STBIW_ASSERT(p);
799
+ if (p) {
800
+ if (!*arr) ((int *) p)[1] = 0;
801
+ *arr = (void *) ((int *) p + 2);
802
+ stbiw__sbm(*arr) = m;
803
+ }
804
+ return *arr;
805
+ }
806
+
807
+ static unsigned char *stbiw__zlib_flushf(unsigned char *data, unsigned int *bitbuffer, int *bitcount)
808
+ {
809
+ while (*bitcount >= 8) {
810
+ stbiw__sbpush(data, STBIW_UCHAR(*bitbuffer));
811
+ *bitbuffer >>= 8;
812
+ *bitcount -= 8;
813
+ }
814
+ return data;
815
+ }
816
+
817
+ static int stbiw__zlib_bitrev(int code, int codebits)
818
+ {
819
+ int res=0;
820
+ while (codebits--) {
821
+ res = (res << 1) | (code & 1);
822
+ code >>= 1;
823
+ }
824
+ return res;
825
+ }
826
+
827
+ static unsigned int stbiw__zlib_countm(unsigned char *a, unsigned char *b, int limit)
828
+ {
829
+ int i;
830
+ for (i=0; i < limit && i < 258; ++i)
831
+ if (a[i] != b[i]) break;
832
+ return i;
833
+ }
834
+
835
+ static unsigned int stbiw__zhash(unsigned char *data)
836
+ {
837
+ stbiw_uint32 hash = data[0] + (data[1] << 8) + (data[2] << 16);
838
+ hash ^= hash << 3;
839
+ hash += hash >> 5;
840
+ hash ^= hash << 4;
841
+ hash += hash >> 17;
842
+ hash ^= hash << 25;
843
+ hash += hash >> 6;
844
+ return hash;
845
+ }
846
+
847
+ #define stbiw__zlib_flush() (out = stbiw__zlib_flushf(out, &bitbuf, &bitcount))
848
+ #define stbiw__zlib_add(code,codebits) \
849
+ (bitbuf |= (code) << bitcount, bitcount += (codebits), stbiw__zlib_flush())
850
+ #define stbiw__zlib_huffa(b,c) stbiw__zlib_add(stbiw__zlib_bitrev(b,c),c)
851
+ // default huffman tables
852
+ #define stbiw__zlib_huff1(n) stbiw__zlib_huffa(0x30 + (n), 8)
853
+ #define stbiw__zlib_huff2(n) stbiw__zlib_huffa(0x190 + (n)-144, 9)
854
+ #define stbiw__zlib_huff3(n) stbiw__zlib_huffa(0 + (n)-256,7)
855
+ #define stbiw__zlib_huff4(n) stbiw__zlib_huffa(0xc0 + (n)-280,8)
856
+ #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))
857
+ #define stbiw__zlib_huffb(n) ((n) <= 143 ? stbiw__zlib_huff1(n) : stbiw__zlib_huff2(n))
858
+
859
+ #define stbiw__ZHASH 16384
860
+
861
+ #endif // STBIW_ZLIB_COMPRESS
862
+
863
+ STBIWDEF unsigned char * stbi_zlib_compress(unsigned char *data, int data_len, int *out_len, int quality)
864
+ {
865
+ #ifdef STBIW_ZLIB_COMPRESS
866
+ // user provided a zlib compress implementation, use that
867
+ return STBIW_ZLIB_COMPRESS(data, data_len, out_len, quality);
868
+ #else // use builtin
869
+ 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 };
870
+ 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 };
871
+ 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 };
872
+ 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 };
873
+ unsigned int bitbuf=0;
874
+ int i,j, bitcount=0;
875
+ unsigned char *out = NULL;
876
+ unsigned char ***hash_table = (unsigned char***) STBIW_MALLOC(stbiw__ZHASH * sizeof(char**));
877
+ if (hash_table == NULL)
878
+ return NULL;
879
+ if (quality < 5) quality = 5;
880
+
881
+ stbiw__sbpush(out, 0x78); // DEFLATE 32K window
882
+ stbiw__sbpush(out, 0x5e); // FLEVEL = 1
883
+ stbiw__zlib_add(1,1); // BFINAL = 1
884
+ stbiw__zlib_add(1,2); // BTYPE = 1 -- fixed huffman
885
+
886
+ for (i=0; i < stbiw__ZHASH; ++i)
887
+ hash_table[i] = NULL;
888
+
889
+ i=0;
890
+ while (i < data_len-3) {
891
+ // hash next 3 bytes of data to be compressed
892
+ int h = stbiw__zhash(data+i)&(stbiw__ZHASH-1), best=3;
893
+ unsigned char *bestloc = 0;
894
+ unsigned char **hlist = hash_table[h];
895
+ int n = stbiw__sbcount(hlist);
896
+ for (j=0; j < n; ++j) {
897
+ if (hlist[j]-data > i-32768) { // if entry lies within window
898
+ int d = stbiw__zlib_countm(hlist[j], data+i, data_len-i);
899
+ if (d >= best) { best=d; bestloc=hlist[j]; }
900
+ }
901
+ }
902
+ // when hash table entry is too long, delete half the entries
903
+ if (hash_table[h] && stbiw__sbn(hash_table[h]) == 2*quality) {
904
+ STBIW_MEMMOVE(hash_table[h], hash_table[h]+quality, sizeof(hash_table[h][0])*quality);
905
+ stbiw__sbn(hash_table[h]) = quality;
906
+ }
907
+ stbiw__sbpush(hash_table[h],data+i);
908
+
909
+ if (bestloc) {
910
+ // "lazy matching" - check match at *next* byte, and if it's better, do cur byte as literal
911
+ h = stbiw__zhash(data+i+1)&(stbiw__ZHASH-1);
912
+ hlist = hash_table[h];
913
+ n = stbiw__sbcount(hlist);
914
+ for (j=0; j < n; ++j) {
915
+ if (hlist[j]-data > i-32767) {
916
+ int e = stbiw__zlib_countm(hlist[j], data+i+1, data_len-i-1);
917
+ if (e > best) { // if next match is better, bail on current match
918
+ bestloc = NULL;
919
+ break;
920
+ }
921
+ }
922
+ }
923
+ }
924
+
925
+ if (bestloc) {
926
+ int d = (int) (data+i - bestloc); // distance back
927
+ STBIW_ASSERT(d <= 32767 && best <= 258);
928
+ for (j=0; best > lengthc[j+1]-1; ++j);
929
+ stbiw__zlib_huff(j+257);
930
+ if (lengtheb[j]) stbiw__zlib_add(best - lengthc[j], lengtheb[j]);
931
+ for (j=0; d > distc[j+1]-1; ++j);
932
+ stbiw__zlib_add(stbiw__zlib_bitrev(j,5),5);
933
+ if (disteb[j]) stbiw__zlib_add(d - distc[j], disteb[j]);
934
+ i += best;
935
+ } else {
936
+ stbiw__zlib_huffb(data[i]);
937
+ ++i;
938
+ }
939
+ }
940
+ // write out final bytes
941
+ for (;i < data_len; ++i)
942
+ stbiw__zlib_huffb(data[i]);
943
+ stbiw__zlib_huff(256); // end of block
944
+ // pad with 0 bits to byte boundary
945
+ while (bitcount)
946
+ stbiw__zlib_add(0,1);
947
+
948
+ for (i=0; i < stbiw__ZHASH; ++i)
949
+ (void) stbiw__sbfree(hash_table[i]);
950
+ STBIW_FREE(hash_table);
951
+
952
+ {
953
+ // compute adler32 on input
954
+ unsigned int s1=1, s2=0;
955
+ int blocklen = (int) (data_len % 5552);
956
+ j=0;
957
+ while (j < data_len) {
958
+ for (i=0; i < blocklen; ++i) { s1 += data[j+i]; s2 += s1; }
959
+ s1 %= 65521; s2 %= 65521;
960
+ j += blocklen;
961
+ blocklen = 5552;
962
+ }
963
+ stbiw__sbpush(out, STBIW_UCHAR(s2 >> 8));
964
+ stbiw__sbpush(out, STBIW_UCHAR(s2));
965
+ stbiw__sbpush(out, STBIW_UCHAR(s1 >> 8));
966
+ stbiw__sbpush(out, STBIW_UCHAR(s1));
967
+ }
968
+ *out_len = stbiw__sbn(out);
969
+ // make returned pointer freeable
970
+ STBIW_MEMMOVE(stbiw__sbraw(out), out, *out_len);
971
+ return (unsigned char *) stbiw__sbraw(out);
972
+ #endif // STBIW_ZLIB_COMPRESS
973
+ }
974
+
975
+ static unsigned int stbiw__crc32(unsigned char *buffer, int len)
976
+ {
977
+ #ifdef STBIW_CRC32
978
+ return STBIW_CRC32(buffer, len);
979
+ #else
980
+ static unsigned int crc_table[256] =
981
+ {
982
+ 0x00000000, 0x77073096, 0xEE0E612C, 0x990951BA, 0x076DC419, 0x706AF48F, 0xE963A535, 0x9E6495A3,
983
+ 0x0eDB8832, 0x79DCB8A4, 0xE0D5E91E, 0x97D2D988, 0x09B64C2B, 0x7EB17CBD, 0xE7B82D07, 0x90BF1D91,
984
+ 0x1DB71064, 0x6AB020F2, 0xF3B97148, 0x84BE41DE, 0x1ADAD47D, 0x6DDDE4EB, 0xF4D4B551, 0x83D385C7,
985
+ 0x136C9856, 0x646BA8C0, 0xFD62F97A, 0x8A65C9EC, 0x14015C4F, 0x63066CD9, 0xFA0F3D63, 0x8D080DF5,
986
+ 0x3B6E20C8, 0x4C69105E, 0xD56041E4, 0xA2677172, 0x3C03E4D1, 0x4B04D447, 0xD20D85FD, 0xA50AB56B,
987
+ 0x35B5A8FA, 0x42B2986C, 0xDBBBC9D6, 0xACBCF940, 0x32D86CE3, 0x45DF5C75, 0xDCD60DCF, 0xABD13D59,
988
+ 0x26D930AC, 0x51DE003A, 0xC8D75180, 0xBFD06116, 0x21B4F4B5, 0x56B3C423, 0xCFBA9599, 0xB8BDA50F,
989
+ 0x2802B89E, 0x5F058808, 0xC60CD9B2, 0xB10BE924, 0x2F6F7C87, 0x58684C11, 0xC1611DAB, 0xB6662D3D,
990
+ 0x76DC4190, 0x01DB7106, 0x98D220BC, 0xEFD5102A, 0x71B18589, 0x06B6B51F, 0x9FBFE4A5, 0xE8B8D433,
991
+ 0x7807C9A2, 0x0F00F934, 0x9609A88E, 0xE10E9818, 0x7F6A0DBB, 0x086D3D2D, 0x91646C97, 0xE6635C01,
992
+ 0x6B6B51F4, 0x1C6C6162, 0x856530D8, 0xF262004E, 0x6C0695ED, 0x1B01A57B, 0x8208F4C1, 0xF50FC457,
993
+ 0x65B0D9C6, 0x12B7E950, 0x8BBEB8EA, 0xFCB9887C, 0x62DD1DDF, 0x15DA2D49, 0x8CD37CF3, 0xFBD44C65,
994
+ 0x4DB26158, 0x3AB551CE, 0xA3BC0074, 0xD4BB30E2, 0x4ADFA541, 0x3DD895D7, 0xA4D1C46D, 0xD3D6F4FB,
995
+ 0x4369E96A, 0x346ED9FC, 0xAD678846, 0xDA60B8D0, 0x44042D73, 0x33031DE5, 0xAA0A4C5F, 0xDD0D7CC9,
996
+ 0x5005713C, 0x270241AA, 0xBE0B1010, 0xC90C2086, 0x5768B525, 0x206F85B3, 0xB966D409, 0xCE61E49F,
997
+ 0x5EDEF90E, 0x29D9C998, 0xB0D09822, 0xC7D7A8B4, 0x59B33D17, 0x2EB40D81, 0xB7BD5C3B, 0xC0BA6CAD,
998
+ 0xEDB88320, 0x9ABFB3B6, 0x03B6E20C, 0x74B1D29A, 0xEAD54739, 0x9DD277AF, 0x04DB2615, 0x73DC1683,
999
+ 0xE3630B12, 0x94643B84, 0x0D6D6A3E, 0x7A6A5AA8, 0xE40ECF0B, 0x9309FF9D, 0x0A00AE27, 0x7D079EB1,
1000
+ 0xF00F9344, 0x8708A3D2, 0x1E01F268, 0x6906C2FE, 0xF762575D, 0x806567CB, 0x196C3671, 0x6E6B06E7,
1001
+ 0xFED41B76, 0x89D32BE0, 0x10DA7A5A, 0x67DD4ACC, 0xF9B9DF6F, 0x8EBEEFF9, 0x17B7BE43, 0x60B08ED5,
1002
+ 0xD6D6A3E8, 0xA1D1937E, 0x38D8C2C4, 0x4FDFF252, 0xD1BB67F1, 0xA6BC5767, 0x3FB506DD, 0x48B2364B,
1003
+ 0xD80D2BDA, 0xAF0A1B4C, 0x36034AF6, 0x41047A60, 0xDF60EFC3, 0xA867DF55, 0x316E8EEF, 0x4669BE79,
1004
+ 0xCB61B38C, 0xBC66831A, 0x256FD2A0, 0x5268E236, 0xCC0C7795, 0xBB0B4703, 0x220216B9, 0x5505262F,
1005
+ 0xC5BA3BBE, 0xB2BD0B28, 0x2BB45A92, 0x5CB36A04, 0xC2D7FFA7, 0xB5D0CF31, 0x2CD99E8B, 0x5BDEAE1D,
1006
+ 0x9B64C2B0, 0xEC63F226, 0x756AA39C, 0x026D930A, 0x9C0906A9, 0xEB0E363F, 0x72076785, 0x05005713,
1007
+ 0x95BF4A82, 0xE2B87A14, 0x7BB12BAE, 0x0CB61B38, 0x92D28E9B, 0xE5D5BE0D, 0x7CDCEFB7, 0x0BDBDF21,
1008
+ 0x86D3D2D4, 0xF1D4E242, 0x68DDB3F8, 0x1FDA836E, 0x81BE16CD, 0xF6B9265B, 0x6FB077E1, 0x18B74777,
1009
+ 0x88085AE6, 0xFF0F6A70, 0x66063BCA, 0x11010B5C, 0x8F659EFF, 0xF862AE69, 0x616BFFD3, 0x166CCF45,
1010
+ 0xA00AE278, 0xD70DD2EE, 0x4E048354, 0x3903B3C2, 0xA7672661, 0xD06016F7, 0x4969474D, 0x3E6E77DB,
1011
+ 0xAED16A4A, 0xD9D65ADC, 0x40DF0B66, 0x37D83BF0, 0xA9BCAE53, 0xDEBB9EC5, 0x47B2CF7F, 0x30B5FFE9,
1012
+ 0xBDBDF21C, 0xCABAC28A, 0x53B39330, 0x24B4A3A6, 0xBAD03605, 0xCDD70693, 0x54DE5729, 0x23D967BF,
1013
+ 0xB3667A2E, 0xC4614AB8, 0x5D681B02, 0x2A6F2B94, 0xB40BBE37, 0xC30C8EA1, 0x5A05DF1B, 0x2D02EF8D
1014
+ };
1015
+
1016
+ unsigned int crc = ~0u;
1017
+ int i;
1018
+ for (i=0; i < len; ++i)
1019
+ crc = (crc >> 8) ^ crc_table[buffer[i] ^ (crc & 0xff)];
1020
+ return ~crc;
1021
+ #endif
1022
+ }
1023
+
1024
+ #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)
1025
+ #define stbiw__wp32(data,v) stbiw__wpng4(data, (v)>>24,(v)>>16,(v)>>8,(v));
1026
+ #define stbiw__wptag(data,s) stbiw__wpng4(data, s[0],s[1],s[2],s[3])
1027
+
1028
+ static void stbiw__wpcrc(unsigned char **data, int len)
1029
+ {
1030
+ unsigned int crc = stbiw__crc32(*data - len - 4, len+4);
1031
+ stbiw__wp32(*data, crc);
1032
+ }
1033
+
1034
+ static unsigned char stbiw__paeth(int a, int b, int c)
1035
+ {
1036
+ int p = a + b - c, pa = abs(p-a), pb = abs(p-b), pc = abs(p-c);
1037
+ if (pa <= pb && pa <= pc) return STBIW_UCHAR(a);
1038
+ if (pb <= pc) return STBIW_UCHAR(b);
1039
+ return STBIW_UCHAR(c);
1040
+ }
1041
+
1042
+ // @OPTIMIZE: provide an option that always forces left-predict or paeth predict
1043
+ 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)
1044
+ {
1045
+ static int mapping[] = { 0,1,2,3,4 };
1046
+ static int firstmap[] = { 0,1,0,5,6 };
1047
+ int *mymap = (y != 0) ? mapping : firstmap;
1048
+ int i;
1049
+ int type = mymap[filter_type];
1050
+ unsigned char *z = pixels + stride_bytes * (stbi__flip_vertically_on_write ? height-1-y : y);
1051
+ int signed_stride = stbi__flip_vertically_on_write ? -stride_bytes : stride_bytes;
1052
+
1053
+ if (type==0) {
1054
+ memcpy(line_buffer, z, width*n);
1055
+ return;
1056
+ }
1057
+
1058
+ // first loop isn't optimized since it's just one pixel
1059
+ for (i = 0; i < n; ++i) {
1060
+ switch (type) {
1061
+ case 1: line_buffer[i] = z[i]; break;
1062
+ case 2: line_buffer[i] = z[i] - z[i-signed_stride]; break;
1063
+ case 3: line_buffer[i] = z[i] - (z[i-signed_stride]>>1); break;
1064
+ case 4: line_buffer[i] = (signed char) (z[i] - stbiw__paeth(0,z[i-signed_stride],0)); break;
1065
+ case 5: line_buffer[i] = z[i]; break;
1066
+ case 6: line_buffer[i] = z[i]; break;
1067
+ }
1068
+ }
1069
+ switch (type) {
1070
+ case 1: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - z[i-n]; break;
1071
+ case 2: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - z[i-signed_stride]; break;
1072
+ case 3: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - ((z[i-n] + z[i-signed_stride])>>1); break;
1073
+ 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;
1074
+ case 5: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - (z[i-n]>>1); break;
1075
+ case 6: for (i=n; i < width*n; ++i) line_buffer[i] = z[i] - stbiw__paeth(z[i-n], 0,0); break;
1076
+ }
1077
+ }
1078
+
1079
+ STBIWDEF unsigned char *stbi_write_png_to_mem(const unsigned char *pixels, int stride_bytes, int x, int y, int n, int *out_len)
1080
+ {
1081
+ int force_filter = stbi_write_force_png_filter;
1082
+ int ctype[5] = { -1, 0, 4, 2, 6 };
1083
+ unsigned char sig[8] = { 137,80,78,71,13,10,26,10 };
1084
+ unsigned char *out,*o, *filt, *zlib;
1085
+ signed char *line_buffer;
1086
+ int j,zlen;
1087
+
1088
+ if (stride_bytes == 0)
1089
+ stride_bytes = x * n;
1090
+
1091
+ if (force_filter >= 5) {
1092
+ force_filter = -1;
1093
+ }
1094
+
1095
+ filt = (unsigned char *) STBIW_MALLOC((x*n+1) * y); if (!filt) return 0;
1096
+ line_buffer = (signed char *) STBIW_MALLOC(x * n); if (!line_buffer) { STBIW_FREE(filt); return 0; }
1097
+ for (j=0; j < y; ++j) {
1098
+ int filter_type;
1099
+ if (force_filter > -1) {
1100
+ filter_type = force_filter;
1101
+ stbiw__encode_png_line((unsigned char*)(pixels), stride_bytes, x, y, j, n, force_filter, line_buffer);
1102
+ } else { // Estimate the best filter by running through all of them:
1103
+ int best_filter = 0, best_filter_val = 0x7fffffff, est, i;
1104
+ for (filter_type = 0; filter_type < 5; filter_type++) {
1105
+ stbiw__encode_png_line((unsigned char*)(pixels), stride_bytes, x, y, j, n, filter_type, line_buffer);
1106
+
1107
+ // Estimate the entropy of the line using this filter; the less, the better.
1108
+ est = 0;
1109
+ for (i = 0; i < x*n; ++i) {
1110
+ est += abs((signed char) line_buffer[i]);
1111
+ }
1112
+ if (est < best_filter_val) {
1113
+ best_filter_val = est;
1114
+ best_filter = filter_type;
1115
+ }
1116
+ }
1117
+ if (filter_type != best_filter) { // If the last iteration already got us the best filter, don't redo it
1118
+ stbiw__encode_png_line((unsigned char*)(pixels), stride_bytes, x, y, j, n, best_filter, line_buffer);
1119
+ filter_type = best_filter;
1120
+ }
1121
+ }
1122
+ // when we get here, filter_type contains the filter type, and line_buffer contains the data
1123
+ filt[j*(x*n+1)] = (unsigned char) filter_type;
1124
+ STBIW_MEMMOVE(filt+j*(x*n+1)+1, line_buffer, x*n);
1125
+ }
1126
+ STBIW_FREE(line_buffer);
1127
+ zlib = stbi_zlib_compress(filt, y*( x*n+1), &zlen, stbi_write_png_compression_level);
1128
+ STBIW_FREE(filt);
1129
+ if (!zlib) return 0;
1130
+
1131
+ // each tag requires 12 bytes of overhead
1132
+ out = (unsigned char *) STBIW_MALLOC(8 + 12+13 + 12+zlen + 12);
1133
+ if (!out) return 0;
1134
+ *out_len = 8 + 12+13 + 12+zlen + 12;
1135
+
1136
+ o=out;
1137
+ STBIW_MEMMOVE(o,sig,8); o+= 8;
1138
+ stbiw__wp32(o, 13); // header length
1139
+ stbiw__wptag(o, "IHDR");
1140
+ stbiw__wp32(o, x);
1141
+ stbiw__wp32(o, y);
1142
+ *o++ = 8;
1143
+ *o++ = STBIW_UCHAR(ctype[n]);
1144
+ *o++ = 0;
1145
+ *o++ = 0;
1146
+ *o++ = 0;
1147
+ stbiw__wpcrc(&o,13);
1148
+
1149
+ stbiw__wp32(o, zlen);
1150
+ stbiw__wptag(o, "IDAT");
1151
+ STBIW_MEMMOVE(o, zlib, zlen);
1152
+ o += zlen;
1153
+ STBIW_FREE(zlib);
1154
+ stbiw__wpcrc(&o, zlen);
1155
+
1156
+ stbiw__wp32(o,0);
1157
+ stbiw__wptag(o, "IEND");
1158
+ stbiw__wpcrc(&o,0);
1159
+
1160
+ STBIW_ASSERT(o == out + *out_len);
1161
+
1162
+ return out;
1163
+ }
1164
+
1165
+ #ifndef STBI_WRITE_NO_STDIO
1166
+ STBIWDEF int stbi_write_png(char const *filename, int x, int y, int comp, const void *data, int stride_bytes)
1167
+ {
1168
+ FILE *f;
1169
+ int len;
1170
+ unsigned char *png = stbi_write_png_to_mem((const unsigned char *) data, stride_bytes, x, y, comp, &len);
1171
+ if (png == NULL) return 0;
1172
+
1173
+ f = stbiw__fopen(filename, "wb");
1174
+ if (!f) { STBIW_FREE(png); return 0; }
1175
+ fwrite(png, 1, len, f);
1176
+ fclose(f);
1177
+ STBIW_FREE(png);
1178
+ return 1;
1179
+ }
1180
+ #endif
1181
+
1182
+ 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)
1183
+ {
1184
+ int len;
1185
+ unsigned char *png = stbi_write_png_to_mem((const unsigned char *) data, stride_bytes, x, y, comp, &len);
1186
+ if (png == NULL) return 0;
1187
+ func(context, png, len);
1188
+ STBIW_FREE(png);
1189
+ return 1;
1190
+ }
1191
+
1192
+
1193
+ /* ***************************************************************************
1194
+ *
1195
+ * JPEG writer
1196
+ *
1197
+ * This is based on Jon Olick's jo_jpeg.cpp:
1198
+ * public domain Simple, Minimalistic JPEG writer - http://www.jonolick.com/code.html
1199
+ */
1200
+
1201
+ 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,
1202
+ 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 };
1203
+
1204
+ static void stbiw__jpg_writeBits(stbi__write_context *s, int *bitBufP, int *bitCntP, const unsigned short *bs) {
1205
+ int bitBuf = *bitBufP, bitCnt = *bitCntP;
1206
+ bitCnt += bs[1];
1207
+ bitBuf |= bs[0] << (24 - bitCnt);
1208
+ while(bitCnt >= 8) {
1209
+ unsigned char c = (bitBuf >> 16) & 255;
1210
+ stbiw__putc(s, c);
1211
+ if(c == 255) {
1212
+ stbiw__putc(s, 0);
1213
+ }
1214
+ bitBuf <<= 8;
1215
+ bitCnt -= 8;
1216
+ }
1217
+ *bitBufP = bitBuf;
1218
+ *bitCntP = bitCnt;
1219
+ }
1220
+
1221
+ static void stbiw__jpg_DCT(float *d0p, float *d1p, float *d2p, float *d3p, float *d4p, float *d5p, float *d6p, float *d7p) {
1222
+ float d0 = *d0p, d1 = *d1p, d2 = *d2p, d3 = *d3p, d4 = *d4p, d5 = *d5p, d6 = *d6p, d7 = *d7p;
1223
+ float z1, z2, z3, z4, z5, z11, z13;
1224
+
1225
+ float tmp0 = d0 + d7;
1226
+ float tmp7 = d0 - d7;
1227
+ float tmp1 = d1 + d6;
1228
+ float tmp6 = d1 - d6;
1229
+ float tmp2 = d2 + d5;
1230
+ float tmp5 = d2 - d5;
1231
+ float tmp3 = d3 + d4;
1232
+ float tmp4 = d3 - d4;
1233
+
1234
+ // Even part
1235
+ float tmp10 = tmp0 + tmp3; // phase 2
1236
+ float tmp13 = tmp0 - tmp3;
1237
+ float tmp11 = tmp1 + tmp2;
1238
+ float tmp12 = tmp1 - tmp2;
1239
+
1240
+ d0 = tmp10 + tmp11; // phase 3
1241
+ d4 = tmp10 - tmp11;
1242
+
1243
+ z1 = (tmp12 + tmp13) * 0.707106781f; // c4
1244
+ d2 = tmp13 + z1; // phase 5
1245
+ d6 = tmp13 - z1;
1246
+
1247
+ // Odd part
1248
+ tmp10 = tmp4 + tmp5; // phase 2
1249
+ tmp11 = tmp5 + tmp6;
1250
+ tmp12 = tmp6 + tmp7;
1251
+
1252
+ // The rotator is modified from fig 4-8 to avoid extra negations.
1253
+ z5 = (tmp10 - tmp12) * 0.382683433f; // c6
1254
+ z2 = tmp10 * 0.541196100f + z5; // c2-c6
1255
+ z4 = tmp12 * 1.306562965f + z5; // c2+c6
1256
+ z3 = tmp11 * 0.707106781f; // c4
1257
+
1258
+ z11 = tmp7 + z3; // phase 5
1259
+ z13 = tmp7 - z3;
1260
+
1261
+ *d5p = z13 + z2; // phase 6
1262
+ *d3p = z13 - z2;
1263
+ *d1p = z11 + z4;
1264
+ *d7p = z11 - z4;
1265
+
1266
+ *d0p = d0; *d2p = d2; *d4p = d4; *d6p = d6;
1267
+ }
1268
+
1269
+ static void stbiw__jpg_calcBits(int val, unsigned short bits[2]) {
1270
+ int tmp1 = val < 0 ? -val : val;
1271
+ val = val < 0 ? val-1 : val;
1272
+ bits[1] = 1;
1273
+ while(tmp1 >>= 1) {
1274
+ ++bits[1];
1275
+ }
1276
+ bits[0] = val & ((1<<bits[1])-1);
1277
+ }
1278
+
1279
+ static int stbiw__jpg_processDU(stbi__write_context *s, int *bitBuf, int *bitCnt, float *CDU, float *fdtbl, int DC, const unsigned short HTDC[256][2], const unsigned short HTAC[256][2]) {
1280
+ const unsigned short EOB[2] = { HTAC[0x00][0], HTAC[0x00][1] };
1281
+ const unsigned short M16zeroes[2] = { HTAC[0xF0][0], HTAC[0xF0][1] };
1282
+ int dataOff, i, diff, end0pos;
1283
+ int DU[64];
1284
+
1285
+ // DCT rows
1286
+ for(dataOff=0; dataOff<64; dataOff+=8) {
1287
+ stbiw__jpg_DCT(&CDU[dataOff], &CDU[dataOff+1], &CDU[dataOff+2], &CDU[dataOff+3], &CDU[dataOff+4], &CDU[dataOff+5], &CDU[dataOff+6], &CDU[dataOff+7]);
1288
+ }
1289
+ // DCT columns
1290
+ for(dataOff=0; dataOff<8; ++dataOff) {
1291
+ stbiw__jpg_DCT(&CDU[dataOff], &CDU[dataOff+8], &CDU[dataOff+16], &CDU[dataOff+24], &CDU[dataOff+32], &CDU[dataOff+40], &CDU[dataOff+48], &CDU[dataOff+56]);
1292
+ }
1293
+ // Quantize/descale/zigzag the coefficients
1294
+ for(i=0; i<64; ++i) {
1295
+ float v = CDU[i]*fdtbl[i];
1296
+ // DU[stbiw__jpg_ZigZag[i]] = (int)(v < 0 ? ceilf(v - 0.5f) : floorf(v + 0.5f));
1297
+ // ceilf() and floorf() are C99, not C89, but I /think/ they're not needed here anyway?
1298
+ DU[stbiw__jpg_ZigZag[i]] = (int)(v < 0 ? v - 0.5f : v + 0.5f);
1299
+ }
1300
+
1301
+ // Encode DC
1302
+ diff = DU[0] - DC;
1303
+ if (diff == 0) {
1304
+ stbiw__jpg_writeBits(s, bitBuf, bitCnt, HTDC[0]);
1305
+ } else {
1306
+ unsigned short bits[2];
1307
+ stbiw__jpg_calcBits(diff, bits);
1308
+ stbiw__jpg_writeBits(s, bitBuf, bitCnt, HTDC[bits[1]]);
1309
+ stbiw__jpg_writeBits(s, bitBuf, bitCnt, bits);
1310
+ }
1311
+ // Encode ACs
1312
+ end0pos = 63;
1313
+ for(; (end0pos>0)&&(DU[end0pos]==0); --end0pos) {
1314
+ }
1315
+ // end0pos = first element in reverse order !=0
1316
+ if(end0pos == 0) {
1317
+ stbiw__jpg_writeBits(s, bitBuf, bitCnt, EOB);
1318
+ return DU[0];
1319
+ }
1320
+ for(i = 1; i <= end0pos; ++i) {
1321
+ int startpos = i;
1322
+ int nrzeroes;
1323
+ unsigned short bits[2];
1324
+ for (; DU[i]==0 && i<=end0pos; ++i) {
1325
+ }
1326
+ nrzeroes = i-startpos;
1327
+ if ( nrzeroes >= 16 ) {
1328
+ int lng = nrzeroes>>4;
1329
+ int nrmarker;
1330
+ for (nrmarker=1; nrmarker <= lng; ++nrmarker)
1331
+ stbiw__jpg_writeBits(s, bitBuf, bitCnt, M16zeroes);
1332
+ nrzeroes &= 15;
1333
+ }
1334
+ stbiw__jpg_calcBits(DU[i], bits);
1335
+ stbiw__jpg_writeBits(s, bitBuf, bitCnt, HTAC[(nrzeroes<<4)+bits[1]]);
1336
+ stbiw__jpg_writeBits(s, bitBuf, bitCnt, bits);
1337
+ }
1338
+ if(end0pos != 63) {
1339
+ stbiw__jpg_writeBits(s, bitBuf, bitCnt, EOB);
1340
+ }
1341
+ return DU[0];
1342
+ }
1343
+
1344
+ static int stbi_write_jpg_core(stbi__write_context *s, int width, int height, int comp, const void* data, int quality) {
1345
+ // Constants that don't pollute global namespace
1346
+ 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};
1347
+ static const unsigned char std_dc_luminance_values[] = {0,1,2,3,4,5,6,7,8,9,10,11};
1348
+ 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};
1349
+ static const unsigned char std_ac_luminance_values[] = {
1350
+ 0x01,0x02,0x03,0x00,0x04,0x11,0x05,0x12,0x21,0x31,0x41,0x06,0x13,0x51,0x61,0x07,0x22,0x71,0x14,0x32,0x81,0x91,0xa1,0x08,
1351
+ 0x23,0x42,0xb1,0xc1,0x15,0x52,0xd1,0xf0,0x24,0x33,0x62,0x72,0x82,0x09,0x0a,0x16,0x17,0x18,0x19,0x1a,0x25,0x26,0x27,0x28,
1352
+ 0x29,0x2a,0x34,0x35,0x36,0x37,0x38,0x39,0x3a,0x43,0x44,0x45,0x46,0x47,0x48,0x49,0x4a,0x53,0x54,0x55,0x56,0x57,0x58,0x59,
1353
+ 0x5a,0x63,0x64,0x65,0x66,0x67,0x68,0x69,0x6a,0x73,0x74,0x75,0x76,0x77,0x78,0x79,0x7a,0x83,0x84,0x85,0x86,0x87,0x88,0x89,
1354
+ 0x8a,0x92,0x93,0x94,0x95,0x96,0x97,0x98,0x99,0x9a,0xa2,0xa3,0xa4,0xa5,0xa6,0xa7,0xa8,0xa9,0xaa,0xb2,0xb3,0xb4,0xb5,0xb6,
1355
+ 0xb7,0xb8,0xb9,0xba,0xc2,0xc3,0xc4,0xc5,0xc6,0xc7,0xc8,0xc9,0xca,0xd2,0xd3,0xd4,0xd5,0xd6,0xd7,0xd8,0xd9,0xda,0xe1,0xe2,
1356
+ 0xe3,0xe4,0xe5,0xe6,0xe7,0xe8,0xe9,0xea,0xf1,0xf2,0xf3,0xf4,0xf5,0xf6,0xf7,0xf8,0xf9,0xfa
1357
+ };
1358
+ 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};
1359
+ static const unsigned char std_dc_chrominance_values[] = {0,1,2,3,4,5,6,7,8,9,10,11};
1360
+ 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};
1361
+ static const unsigned char std_ac_chrominance_values[] = {
1362
+ 0x00,0x01,0x02,0x03,0x11,0x04,0x05,0x21,0x31,0x06,0x12,0x41,0x51,0x07,0x61,0x71,0x13,0x22,0x32,0x81,0x08,0x14,0x42,0x91,
1363
+ 0xa1,0xb1,0xc1,0x09,0x23,0x33,0x52,0xf0,0x15,0x62,0x72,0xd1,0x0a,0x16,0x24,0x34,0xe1,0x25,0xf1,0x17,0x18,0x19,0x1a,0x26,
1364
+ 0x27,0x28,0x29,0x2a,0x35,0x36,0x37,0x38,0x39,0x3a,0x43,0x44,0x45,0x46,0x47,0x48,0x49,0x4a,0x53,0x54,0x55,0x56,0x57,0x58,
1365
+ 0x59,0x5a,0x63,0x64,0x65,0x66,0x67,0x68,0x69,0x6a,0x73,0x74,0x75,0x76,0x77,0x78,0x79,0x7a,0x82,0x83,0x84,0x85,0x86,0x87,
1366
+ 0x88,0x89,0x8a,0x92,0x93,0x94,0x95,0x96,0x97,0x98,0x99,0x9a,0xa2,0xa3,0xa4,0xa5,0xa6,0xa7,0xa8,0xa9,0xaa,0xb2,0xb3,0xb4,
1367
+ 0xb5,0xb6,0xb7,0xb8,0xb9,0xba,0xc2,0xc3,0xc4,0xc5,0xc6,0xc7,0xc8,0xc9,0xca,0xd2,0xd3,0xd4,0xd5,0xd6,0xd7,0xd8,0xd9,0xda,
1368
+ 0xe2,0xe3,0xe4,0xe5,0xe6,0xe7,0xe8,0xe9,0xea,0xf2,0xf3,0xf4,0xf5,0xf6,0xf7,0xf8,0xf9,0xfa
1369
+ };
1370
+ // Huffman tables
1371
+ 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}};
1372
+ 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}};
1373
+ static const unsigned short YAC_HT[256][2] = {
1374
+ {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},
1375
+ {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},
1376
+ {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},
1377
+ {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},
1378
+ {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},
1379
+ {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},
1380
+ {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},
1381
+ {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},
1382
+ {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},
1383
+ {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},
1384
+ {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},
1385
+ {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},
1386
+ {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},
1387
+ {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},
1388
+ {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},
1389
+ {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}
1390
+ };
1391
+ static const unsigned short UVAC_HT[256][2] = {
1392
+ {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},
1393
+ {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},
1394
+ {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},
1395
+ {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},
1396
+ {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},
1397
+ {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},
1398
+ {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},
1399
+ {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},
1400
+ {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},
1401
+ {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},
1402
+ {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},
1403
+ {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},
1404
+ {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},
1405
+ {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},
1406
+ {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},
1407
+ {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}
1408
+ };
1409
+ 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,
1410
+ 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};
1411
+ 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,
1412
+ 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};
1413
+ static const float aasf[] = { 1.0f * 2.828427125f, 1.387039845f * 2.828427125f, 1.306562965f * 2.828427125f, 1.175875602f * 2.828427125f,
1414
+ 1.0f * 2.828427125f, 0.785694958f * 2.828427125f, 0.541196100f * 2.828427125f, 0.275899379f * 2.828427125f };
1415
+
1416
+ int row, col, i, k;
1417
+ float fdtbl_Y[64], fdtbl_UV[64];
1418
+ unsigned char YTable[64], UVTable[64];
1419
+
1420
+ if(!data || !width || !height || comp > 4 || comp < 1) {
1421
+ return 0;
1422
+ }
1423
+
1424
+ quality = quality ? quality : 90;
1425
+ quality = quality < 1 ? 1 : quality > 100 ? 100 : quality;
1426
+ quality = quality < 50 ? 5000 / quality : 200 - quality * 2;
1427
+
1428
+ for(i = 0; i < 64; ++i) {
1429
+ int uvti, yti = (YQT[i]*quality+50)/100;
1430
+ YTable[stbiw__jpg_ZigZag[i]] = (unsigned char) (yti < 1 ? 1 : yti > 255 ? 255 : yti);
1431
+ uvti = (UVQT[i]*quality+50)/100;
1432
+ UVTable[stbiw__jpg_ZigZag[i]] = (unsigned char) (uvti < 1 ? 1 : uvti > 255 ? 255 : uvti);
1433
+ }
1434
+
1435
+ for(row = 0, k = 0; row < 8; ++row) {
1436
+ for(col = 0; col < 8; ++col, ++k) {
1437
+ fdtbl_Y[k] = 1 / (YTable [stbiw__jpg_ZigZag[k]] * aasf[row] * aasf[col]);
1438
+ fdtbl_UV[k] = 1 / (UVTable[stbiw__jpg_ZigZag[k]] * aasf[row] * aasf[col]);
1439
+ }
1440
+ }
1441
+
1442
+ // Write Headers
1443
+ {
1444
+ 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 };
1445
+ static const unsigned char head2[] = { 0xFF,0xDA,0,0xC,3,1,0,2,0x11,3,0x11,0,0x3F,0 };
1446
+ const unsigned char head1[] = { 0xFF,0xC0,0,0x11,8,(unsigned char)(height>>8),STBIW_UCHAR(height),(unsigned char)(width>>8),STBIW_UCHAR(width),
1447
+ 3,1,0x11,0,2,0x11,1,3,0x11,1,0xFF,0xC4,0x01,0xA2,0 };
1448
+ s->func(s->context, (void*)head0, sizeof(head0));
1449
+ s->func(s->context, (void*)YTable, sizeof(YTable));
1450
+ stbiw__putc(s, 1);
1451
+ s->func(s->context, UVTable, sizeof(UVTable));
1452
+ s->func(s->context, (void*)head1, sizeof(head1));
1453
+ s->func(s->context, (void*)(std_dc_luminance_nrcodes+1), sizeof(std_dc_luminance_nrcodes)-1);
1454
+ s->func(s->context, (void*)std_dc_luminance_values, sizeof(std_dc_luminance_values));
1455
+ stbiw__putc(s, 0x10); // HTYACinfo
1456
+ s->func(s->context, (void*)(std_ac_luminance_nrcodes+1), sizeof(std_ac_luminance_nrcodes)-1);
1457
+ s->func(s->context, (void*)std_ac_luminance_values, sizeof(std_ac_luminance_values));
1458
+ stbiw__putc(s, 1); // HTUDCinfo
1459
+ s->func(s->context, (void*)(std_dc_chrominance_nrcodes+1), sizeof(std_dc_chrominance_nrcodes)-1);
1460
+ s->func(s->context, (void*)std_dc_chrominance_values, sizeof(std_dc_chrominance_values));
1461
+ stbiw__putc(s, 0x11); // HTUACinfo
1462
+ s->func(s->context, (void*)(std_ac_chrominance_nrcodes+1), sizeof(std_ac_chrominance_nrcodes)-1);
1463
+ s->func(s->context, (void*)std_ac_chrominance_values, sizeof(std_ac_chrominance_values));
1464
+ s->func(s->context, (void*)head2, sizeof(head2));
1465
+ }
1466
+
1467
+ // Encode 8x8 macroblocks
1468
+ {
1469
+ static const unsigned short fillBits[] = {0x7F, 7};
1470
+ const unsigned char *imageData = (const unsigned char *)data;
1471
+ int DCY=0, DCU=0, DCV=0;
1472
+ int bitBuf=0, bitCnt=0;
1473
+ // comp == 2 is grey+alpha (alpha is ignored)
1474
+ int ofsG = comp > 2 ? 1 : 0, ofsB = comp > 2 ? 2 : 0;
1475
+ int x, y, pos;
1476
+ for(y = 0; y < height; y += 8) {
1477
+ for(x = 0; x < width; x += 8) {
1478
+ float YDU[64], UDU[64], VDU[64];
1479
+ for(row = y, pos = 0; row < y+8; ++row) {
1480
+ // row >= height => use last input row
1481
+ int clamped_row = (row < height) ? row : height - 1;
1482
+ int base_p = (stbi__flip_vertically_on_write ? (height-1-clamped_row) : clamped_row)*width*comp;
1483
+ for(col = x; col < x+8; ++col, ++pos) {
1484
+ float r, g, b;
1485
+ // if col >= width => use pixel from last input column
1486
+ int p = base_p + ((col < width) ? col : (width-1))*comp;
1487
+
1488
+ r = imageData[p+0];
1489
+ g = imageData[p+ofsG];
1490
+ b = imageData[p+ofsB];
1491
+ YDU[pos]=+0.29900f*r+0.58700f*g+0.11400f*b-128;
1492
+ UDU[pos]=-0.16874f*r-0.33126f*g+0.50000f*b;
1493
+ VDU[pos]=+0.50000f*r-0.41869f*g-0.08131f*b;
1494
+ }
1495
+ }
1496
+
1497
+ DCY = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, YDU, fdtbl_Y, DCY, YDC_HT, YAC_HT);
1498
+ DCU = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, UDU, fdtbl_UV, DCU, UVDC_HT, UVAC_HT);
1499
+ DCV = stbiw__jpg_processDU(s, &bitBuf, &bitCnt, VDU, fdtbl_UV, DCV, UVDC_HT, UVAC_HT);
1500
+ }
1501
+ }
1502
+
1503
+ // Do the bit alignment of the EOI marker
1504
+ stbiw__jpg_writeBits(s, &bitBuf, &bitCnt, fillBits);
1505
+ }
1506
+
1507
+ // EOI
1508
+ stbiw__putc(s, 0xFF);
1509
+ stbiw__putc(s, 0xD9);
1510
+
1511
+ return 1;
1512
+ }
1513
+
1514
+ STBIWDEF int stbi_write_jpg_to_func(stbi_write_func *func, void *context, int x, int y, int comp, const void *data, int quality)
1515
+ {
1516
+ stbi__write_context s;
1517
+ stbi__start_write_callbacks(&s, func, context);
1518
+ return stbi_write_jpg_core(&s, x, y, comp, (void *) data, quality);
1519
+ }
1520
+
1521
+
1522
+ #ifndef STBI_WRITE_NO_STDIO
1523
+ STBIWDEF int stbi_write_jpg(char const *filename, int x, int y, int comp, const void *data, int quality)
1524
+ {
1525
+ stbi__write_context s;
1526
+ if (stbi__start_write_file(&s,filename)) {
1527
+ int r = stbi_write_jpg_core(&s, x, y, comp, data, quality);
1528
+ stbi__end_write_file(&s);
1529
+ return r;
1530
+ } else
1531
+ return 0;
1532
+ }
1533
+ #endif
1534
+
1535
+ #endif // STB_IMAGE_WRITE_IMPLEMENTATION
1536
+
1537
+ /* Revision history
1538
+ 1.10 (2019-02-07)
1539
+ support utf8 filenames in Windows; fix warnings and platform ifdefs
1540
+ 1.09 (2018-02-11)
1541
+ fix typo in zlib quality API, improve STB_I_W_STATIC in C++
1542
+ 1.08 (2018-01-29)
1543
+ add stbi__flip_vertically_on_write, external zlib, zlib quality, choose PNG filter
1544
+ 1.07 (2017-07-24)
1545
+ doc fix
1546
+ 1.06 (2017-07-23)
1547
+ writing JPEG (using Jon Olick's code)
1548
+ 1.05 ???
1549
+ 1.04 (2017-03-03)
1550
+ monochrome BMP expansion
1551
+ 1.03 ???
1552
+ 1.02 (2016-04-02)
1553
+ avoid allocating large structures on the stack
1554
+ 1.01 (2016-01-16)
1555
+ STBIW_REALLOC_SIZED: support allocators with no realloc support
1556
+ avoid race-condition in crc initialization
1557
+ minor compile issues
1558
+ 1.00 (2015-09-14)
1559
+ installable file IO function
1560
+ 0.99 (2015-09-13)
1561
+ warning fixes; TGA rle support
1562
+ 0.98 (2015-04-08)
1563
+ added STBIW_MALLOC, STBIW_ASSERT etc
1564
+ 0.97 (2015-01-18)
1565
+ fixed HDR asserts, rewrote HDR rle logic
1566
+ 0.96 (2015-01-17)
1567
+ add HDR output
1568
+ fix monochrome BMP
1569
+ 0.95 (2014-08-17)
1570
+ add monochrome TGA output
1571
+ 0.94 (2014-05-31)
1572
+ rename private functions to avoid conflicts with stb_image.h
1573
+ 0.93 (2014-05-27)
1574
+ warning fixes
1575
+ 0.92 (2010-08-01)
1576
+ casts to unsigned char to fix warnings
1577
+ 0.91 (2010-07-17)
1578
+ first public release
1579
+ 0.90 first internal release
1580
+ */
1581
+
1582
+ /*
1583
+ ------------------------------------------------------------------------------
1584
+ This software is available under 2 licenses -- choose whichever you prefer.
1585
+ ------------------------------------------------------------------------------
1586
+ ALTERNATIVE A - MIT License
1587
+ Copyright (c) 2017 Sean Barrett
1588
+ Permission is hereby granted, free of charge, to any person obtaining a copy of
1589
+ this software and associated documentation files (the "Software"), to deal in
1590
+ the Software without restriction, including without limitation the rights to
1591
+ use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies
1592
+ of the Software, and to permit persons to whom the Software is furnished to do
1593
+ so, subject to the following conditions:
1594
+ The above copyright notice and this permission notice shall be included in all
1595
+ copies or substantial portions of the Software.
1596
+ THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
1597
+ IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
1598
+ FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
1599
+ AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
1600
+ LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
1601
+ OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
1602
+ SOFTWARE.
1603
+ ------------------------------------------------------------------------------
1604
+ ALTERNATIVE B - Public Domain (www.unlicense.org)
1605
+ This is free and unencumbered software released into the public domain.
1606
+ Anyone is free to copy, modify, publish, use, compile, sell, or distribute this
1607
+ software, either in source code form or as a compiled binary, for any purpose,
1608
+ commercial or non-commercial, and by any means.
1609
+ In jurisdictions that recognize copyright laws, the author or authors of this
1610
+ software dedicate any and all copyright interest in the software to the public
1611
+ domain. We make this dedication for the benefit of the public at large and to
1612
+ the detriment of our heirs and successors. We intend this dedication to be an
1613
+ overt act of relinquishment in perpetuity of all present and future rights to
1614
+ this software under copyright law.
1615
+ THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
1616
+ IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
1617
+ FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
1618
+ AUTHORS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN
1619
+ ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
1620
+ WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
1621
+ ------------------------------------------------------------------------------
1622
+ */