Parent Directory | Revision Log
Revision 677 -
(view)
(download)
Original Path: trunk/xvidcore/src/quant/adapt_quant.c
1 : | edgomez | 443 | /***************************************************************************** |
2 : | * | ||
3 : | * XVID MPEG-4 VIDEO CODEC | ||
4 : | * - Adaptive quantization functions - | ||
5 : | * | ||
6 : | suxen_drol | 499 | * Copyright(C) 2002 Peter Ross <pross@xvid.org> |
7 : | edgomez | 605 | * 2002 Christoph Lampert <gruel@web.de> |
8 : | edgomez | 443 | * |
9 : | edgomez | 653 | * This file is part of XviD, a free MPEG-4 video encoder/decoder |
10 : | edgomez | 443 | * |
11 : | edgomez | 653 | * XviD is free software; you can redistribute it and/or modify it |
12 : | * under the terms of the GNU General Public License as published by | ||
13 : | * the Free Software Foundation; either version 2 of the License, or | ||
14 : | edgomez | 443 | * (at your option) any later version. |
15 : | * | ||
16 : | * This program is distributed in the hope that it will be useful, | ||
17 : | edgomez | 653 | * but WITHOUT ANY WARRANTY; without even the implied warranty of |
18 : | edgomez | 443 | * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
19 : | * GNU General Public License for more details. | ||
20 : | * | ||
21 : | * You should have received a copy of the GNU General Public License | ||
22 : | edgomez | 653 | * along with this program; if not, write to the Free Software |
23 : | edgomez | 443 | * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA |
24 : | * | ||
25 : | edgomez | 653 | * Under section 8 of the GNU General Public License, the copyright |
26 : | * holders of XVID explicitly forbid distribution in the following | ||
27 : | * countries: | ||
28 : | edgomez | 443 | * |
29 : | edgomez | 653 | * - Japan |
30 : | * - United States of America | ||
31 : | * | ||
32 : | * Linking XviD statically or dynamically with other modules is making a | ||
33 : | * combined work based on XviD. Thus, the terms and conditions of the | ||
34 : | * GNU General Public License cover the whole combination. | ||
35 : | * | ||
36 : | * As a special exception, the copyright holders of XviD give you | ||
37 : | * permission to link XviD with independent modules that communicate with | ||
38 : | * XviD solely through the VFW1.1 and DShow interfaces, regardless of the | ||
39 : | * license terms of these independent modules, and to copy and distribute | ||
40 : | * the resulting combined work under terms of your choice, provided that | ||
41 : | * every copy of the combined work is accompanied by a complete copy of | ||
42 : | * the source code of XviD (the version of XviD used to produce the | ||
43 : | * combined work), being distributed under the terms of the GNU General | ||
44 : | * Public License plus this exception. An independent module is a module | ||
45 : | * which is not derived from or based on XviD. | ||
46 : | * | ||
47 : | * Note that people who make modified versions of XviD are not obligated | ||
48 : | * to grant this special exception for their modified versions; it is | ||
49 : | * their choice whether to do so. The GNU General Public License gives | ||
50 : | * permission to release a modified version without this exception; this | ||
51 : | * exception also makes it possible to release a modified version which | ||
52 : | * carries forward this exception. | ||
53 : | * | ||
54 : | edgomez | 677 | * $Id: adapt_quant.c,v 1.13 2002-11-26 23:44:11 edgomez Exp $ |
55 : | edgomez | 653 | * |
56 : | edgomez | 443 | ****************************************************************************/ |
57 : | |||
58 : | Isibaar | 3 | #include "../portab.h" |
59 : | #include "adapt_quant.h" | ||
60 : | |||
61 : | edgomez | 195 | #include <stdlib.h> /* free, malloc */ |
62 : | edgomez | 18 | |
63 : | Isibaar | 3 | #define MAX(a,b) (((a) > (b)) ? (a) : (b)) |
64 : | #define RDIFF(a,b) ((int)(a+0.5)-(int)(b+0.5)) | ||
65 : | |||
66 : | edgomez | 443 | /***************************************************************************** |
67 : | * Functions | ||
68 : | ****************************************************************************/ | ||
69 : | |||
70 : | edgomez | 195 | int |
71 : | normalize_quantizer_field(float *in, | ||
72 : | int *out, | ||
73 : | int num, | ||
74 : | int min_quant, | ||
75 : | int max_quant) | ||
76 : | Isibaar | 3 | { |
77 : | int i; | ||
78 : | int finished; | ||
79 : | edgomez | 195 | |
80 : | do { | ||
81 : | Isibaar | 3 | finished = 1; |
82 : | edgomez | 195 | for (i = 1; i < num; i++) { |
83 : | if (RDIFF(in[i], in[i - 1]) > 2) { | ||
84 : | Isibaar | 3 | in[i] -= (float) 0.5; |
85 : | finished = 0; | ||
86 : | edgomez | 195 | } else if (RDIFF(in[i], in[i - 1]) < -2) { |
87 : | in[i - 1] -= (float) 0.5; | ||
88 : | Isibaar | 3 | finished = 0; |
89 : | } | ||
90 : | edgomez | 195 | |
91 : | if (in[i] > max_quant) { | ||
92 : | edgomez | 18 | in[i] = (float) max_quant; |
93 : | finished = 0; | ||
94 : | } | ||
95 : | edgomez | 195 | if (in[i] < min_quant) { |
96 : | edgomez | 18 | in[i] = (float) min_quant; |
97 : | finished = 0; | ||
98 : | } | ||
99 : | edgomez | 195 | if (in[i - 1] > max_quant) { |
100 : | in[i - 1] = (float) max_quant; | ||
101 : | edgomez | 18 | finished = 0; |
102 : | } | ||
103 : | edgomez | 195 | if (in[i - 1] < min_quant) { |
104 : | in[i - 1] = (float) min_quant; | ||
105 : | edgomez | 18 | finished = 0; |
106 : | } | ||
107 : | Isibaar | 3 | } |
108 : | edgomez | 195 | } while (!finished); |
109 : | |||
110 : | Isibaar | 3 | out[0] = 0; |
111 : | for (i = 1; i < num; i++) | ||
112 : | edgomez | 195 | out[i] = RDIFF(in[i], in[i - 1]); |
113 : | |||
114 : | Isibaar | 3 | return (int) (in[0] + 0.5); |
115 : | } | ||
116 : | |||
117 : | edgomez | 195 | int |
118 : | adaptive_quantization(unsigned char *buf, | ||
119 : | int stride, | ||
120 : | int *intquant, | ||
121 : | int framequant, | ||
122 : | int min_quant, | ||
123 : | int max_quant, | ||
124 : | int mb_width, | ||
125 : | edgomez | 677 | int mb_height) /* no qstride because normalization */ |
126 : | Isibaar | 3 | { |
127 : | edgomez | 195 | int i, j, k, l; |
128 : | |||
129 : | edgomez | 391 | float *quant; |
130 : | Isibaar | 3 | unsigned char *ptr; |
131 : | float *val; | ||
132 : | Isibaar | 357 | float global = 0.; |
133 : | uint32_t mid_range = 0; | ||
134 : | Isibaar | 3 | |
135 : | Isibaar | 357 | const float DarkAmpl = 14 / 2; |
136 : | const float BrightAmpl = 10 / 2; | ||
137 : | const float DarkThres = 70; | ||
138 : | const float BrightThres = 200; | ||
139 : | Isibaar | 3 | |
140 : | Isibaar | 357 | const float GlobalDarkThres = 60; |
141 : | const float GlobalBrightThres = 170; | ||
142 : | Isibaar | 3 | |
143 : | Isibaar | 357 | const float MidRangeThres = 20; |
144 : | const float UpperLimit = 200; | ||
145 : | const float LowerLimit = 25; | ||
146 : | Isibaar | 3 | |
147 : | Isibaar | 357 | |
148 : | edgomez | 391 | if (!(quant = (float *) malloc(mb_width * mb_height * sizeof(float)))) |
149 : | return(-1); | ||
150 : | Isibaar | 3 | |
151 : | edgomez | 496 | if(!(val = (float *) malloc(mb_width * mb_height * sizeof(float)))) { |
152 : | free(quant); | ||
153 : | edgomez | 391 | return(-1); |
154 : | edgomez | 496 | } |
155 : | Isibaar | 3 | |
156 : | edgomez | 195 | for (k = 0; k < mb_height; k++) { |
157 : | edgomez | 677 | for (l = 0; l < mb_width; l++) /* do this for all macroblocks individually */ |
158 : | Isibaar | 3 | { |
159 : | edgomez | 195 | quant[k * mb_width + l] = (float) framequant; |
160 : | |||
161 : | edgomez | 677 | /* calculate luminance-masking */ |
162 : | ptr = &buf[16 * k * stride + 16 * l]; /* address of MB */ | ||
163 : | Isibaar | 3 | |
164 : | edgomez | 195 | val[k * mb_width + l] = 0.; |
165 : | |||
166 : | for (i = 0; i < 16; i++) | ||
167 : | for (j = 0; j < 16; j++) | ||
168 : | val[k * mb_width + l] += ptr[i * stride + j]; | ||
169 : | val[k * mb_width + l] /= 256.; | ||
170 : | Isibaar | 357 | global +=val[k * mb_width + l]; |
171 : | |||
172 : | if ((val[k * mb_width + l] > LowerLimit) && | ||
173 : | (val[k * mb_width + l] < UpperLimit)) | ||
174 : | mid_range++; | ||
175 : | Isibaar | 3 | } |
176 : | } | ||
177 : | |||
178 : | Isibaar | 357 | global /=mb_width * mb_height; |
179 : | Isibaar | 3 | |
180 : | Isibaar | 357 | if (((global <GlobalBrightThres) &&(global >GlobalDarkThres)) |
181 : | || (mid_range < MidRangeThres)) { | ||
182 : | for (k = 0; k < mb_height; k++) { | ||
183 : | edgomez | 677 | for (l = 0; l < mb_width; l++) /* do this for all macroblocks individually */ |
184 : | Isibaar | 357 | { |
185 : | if (val[k * mb_width + l] < DarkThres) | ||
186 : | quant[k * mb_width + l] += | ||
187 : | DarkAmpl * (DarkThres - | ||
188 : | val[k * mb_width + l]) / DarkThres; | ||
189 : | else if (val[k * mb_width + l] > BrightThres) | ||
190 : | quant[k * mb_width + l] += | ||
191 : | BrightAmpl * (val[k * mb_width + l] - | ||
192 : | BrightThres) / (255 - BrightThres); | ||
193 : | } | ||
194 : | Isibaar | 3 | } |
195 : | } | ||
196 : | edgomez | 495 | |
197 : | i = normalize_quantizer_field(quant, intquant, | ||
198 : | mb_width * mb_height, | ||
199 : | min_quant, max_quant); | ||
200 : | |||
201 : | Isibaar | 3 | free(val); |
202 : | edgomez | 391 | free(quant); |
203 : | edgomez | 495 | |
204 : | return(i); | ||
205 : | |||
206 : | edgomez | 18 | } |
No admin address has been configured | ViewVC Help |
Powered by ViewVC 1.0.4 |