1 |
/************************************************************************** |
/***************************************************************************** |
2 |
* |
* |
3 |
* XVID MPEG-4 VIDEO CODEC |
* XVID MPEG-4 VIDEO CODEC |
4 |
* image stuff |
* - image module - |
5 |
|
* |
6 |
|
* Copyright(C) 2002 Peter Ross <pross@xvid.org> |
7 |
* |
* |
8 |
* This program is an implementation of a part of one or more MPEG-4 |
* This program is an implementation of a part of one or more MPEG-4 |
9 |
* Video tools as specified in ISO/IEC 14496-2 standard. Those intending |
* Video tools as specified in ISO/IEC 14496-2 standard. Those intending |
26 |
* |
* |
27 |
* You should have received a copy of the GNU General Public License |
* You should have received a copy of the GNU General Public License |
28 |
* along with this program; if not, write to the Free Software |
* along with this program; if not, write to the Free Software |
29 |
* Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. |
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA |
|
* |
|
|
*************************************************************************/ |
|
|
|
|
|
/************************************************************************** |
|
|
* |
|
|
* History: |
|
30 |
* |
* |
31 |
* 05.10.2002 support for interpolated images in qpel mode - Isibaar |
****************************************************************************/ |
|
* 01.05.2002 BFRAME image-based u,v interpolation |
|
|
* 22.04.2002 added some B-frame support |
|
|
* 14.04.2002 added image_dump_yuvpgm(), added image_mad() |
|
|
* XVID_CSP_USER input support |
|
|
* 09.04.2002 PSNR calculations - Isibaar |
|
|
* 06.04.2002 removed interlaced edging from U,V blocks (as per spec) |
|
|
* 26.03.2002 interlacing support (field-based edging in set_edges) |
|
|
* 26.01.2002 rgb555, rgb565 |
|
|
* 07.01.2001 commented u,v interpolation (not required for uv-block-based) |
|
|
* 23.12.2001 removed #ifdefs, added function pointers + init_common() |
|
|
* 22.12.2001 cpu #ifdefs |
|
|
* 19.12.2001 image_dump(); useful for debugging |
|
|
* 6.12.2001 inital version; (c)2001 peter ross <pross@cs.rmit.edu.au> |
|
|
* |
|
|
*************************************************************************/ |
|
32 |
|
|
33 |
#include <stdlib.h> |
#include <stdlib.h> |
34 |
#include <string.h> // memcpy, memset |
#include <string.h> // memcpy, memset |
142 |
uint32_t edged_width, |
uint32_t edged_width, |
143 |
uint32_t edged_height, |
uint32_t edged_height, |
144 |
uint32_t width, |
uint32_t width, |
145 |
uint32_t height) |
uint32_t height, |
146 |
|
uint32_t interlacing) |
147 |
{ |
{ |
148 |
const uint32_t edged_width2 = edged_width / 2; |
const uint32_t edged_width2 = edged_width / 2; |
149 |
const uint32_t width2 = width / 2; |
const uint32_t width2 = width / 2; |
156 |
src = image->y; |
src = image->y; |
157 |
|
|
158 |
for (i = 0; i < EDGE_SIZE; i++) { |
for (i = 0; i < EDGE_SIZE; i++) { |
159 |
|
/* // if interlacing, edges contain top-most data from each field |
160 |
|
if (interlacing && (i & 1)) { |
161 |
|
memset(dst, *(src + edged_width), EDGE_SIZE); |
162 |
|
memcpy(dst + EDGE_SIZE, src + edged_width, width); |
163 |
|
memset(dst + edged_width - EDGE_SIZE, |
164 |
|
*(src + edged_width + width - 1), EDGE_SIZE); |
165 |
|
} else {*/ |
166 |
memset(dst, *src, EDGE_SIZE); |
memset(dst, *src, EDGE_SIZE); |
167 |
memcpy(dst + EDGE_SIZE, src, width); |
memcpy(dst + EDGE_SIZE, src, width); |
168 |
memset(dst + edged_width - EDGE_SIZE, *(src + width - 1), |
memset(dst + edged_width - EDGE_SIZE, *(src + width - 1), |
169 |
EDGE_SIZE); |
EDGE_SIZE); |
170 |
|
/*}*/ |
171 |
dst += edged_width; |
dst += edged_width; |
172 |
} |
} |
173 |
|
|
180 |
|
|
181 |
src -= edged_width; |
src -= edged_width; |
182 |
for (i = 0; i < EDGE_SIZE; i++) { |
for (i = 0; i < EDGE_SIZE; i++) { |
183 |
|
/* // if interlacing, edges contain bottom-most data from each field |
184 |
|
if (interlacing && !(i & 1)) { |
185 |
|
memset(dst, *(src - edged_width), EDGE_SIZE); |
186 |
|
memcpy(dst + EDGE_SIZE, src - edged_width, width); |
187 |
|
memset(dst + edged_width - EDGE_SIZE, |
188 |
|
*(src - edged_width + width - 1), EDGE_SIZE); |
189 |
|
} else {*/ |
190 |
memset(dst, *src, EDGE_SIZE); |
memset(dst, *src, EDGE_SIZE); |
191 |
memcpy(dst + EDGE_SIZE, src, width); |
memcpy(dst + EDGE_SIZE, src, width); |
192 |
memset(dst + edged_width - EDGE_SIZE, *(src + width - 1), |
memset(dst + edged_width - EDGE_SIZE, *(src + width - 1), |
193 |
EDGE_SIZE); |
EDGE_SIZE); |
194 |
|
/*}*/ |
195 |
dst += edged_width; |
dst += edged_width; |
196 |
} |
} |
197 |
|
|
260 |
IMAGE * refhv, |
IMAGE * refhv, |
261 |
uint32_t edged_width, |
uint32_t edged_width, |
262 |
uint32_t edged_height, |
uint32_t edged_height, |
|
uint32_t quarterpel, |
|
263 |
uint32_t rounding) |
uint32_t rounding) |
264 |
{ |
{ |
265 |
const uint32_t offset = EDGE_SIZE2 * (edged_width + 1); // we only interpolate half of the edge area |
const uint32_t offset = EDGE_SIZE * (edged_width + 1); |
266 |
const uint32_t stride_add = 7 * edged_width; |
const uint32_t stride_add = 7 * edged_width; |
267 |
/* |
|
268 |
#ifdef BFRAMES |
#ifdef BFRAMES |
269 |
const uint32_t edged_width2 = edged_width / 2; |
const uint32_t edged_width2 = edged_width / 2; |
270 |
const uint32_t edged_height2 = edged_height / 2; |
const uint32_t edged_height2 = edged_height / 2; |
271 |
const uint32_t offset2 = EDGE_SIZE2 * (edged_width2 + 1); |
const uint32_t offset2 = EDGE_SIZE2 * (edged_width2 + 1); |
272 |
const uint32_t stride_add2 = 7 * edged_width2; |
const uint32_t stride_add2 = 7 * edged_width2; |
273 |
#endif |
#endif |
274 |
*/ |
|
275 |
uint8_t *n_ptr, *h_ptr, *v_ptr, *hv_ptr; |
uint8_t *n_ptr, *h_ptr, *v_ptr, *hv_ptr; |
276 |
uint32_t x, y; |
uint32_t x, y; |
277 |
|
|
286 |
v_ptr -= offset; |
v_ptr -= offset; |
287 |
hv_ptr -= offset; |
hv_ptr -= offset; |
288 |
|
|
289 |
if(quarterpel) { |
for (y = 0; y < edged_height; y = y + 8) { |
290 |
|
for (x = 0; x < edged_width; x = x + 8) { |
|
for (y = 0; y < (edged_height - EDGE_SIZE); y += 8) { |
|
|
for (x = 0; x < (edged_width - EDGE_SIZE); x += 8) { |
|
|
interpolate8x8_6tap_lowpass_h(h_ptr, n_ptr, edged_width, rounding); |
|
|
interpolate8x8_6tap_lowpass_v(v_ptr, n_ptr, edged_width, rounding); |
|
|
|
|
|
n_ptr += 8; |
|
|
h_ptr += 8; |
|
|
v_ptr += 8; |
|
|
} |
|
|
|
|
|
n_ptr += EDGE_SIZE; |
|
|
h_ptr += EDGE_SIZE; |
|
|
v_ptr += EDGE_SIZE; |
|
|
|
|
|
h_ptr += stride_add; |
|
|
v_ptr += stride_add; |
|
|
n_ptr += stride_add; |
|
|
} |
|
|
|
|
|
h_ptr = refh->y; |
|
|
h_ptr -= offset; |
|
|
|
|
|
for (y = 0; y < (edged_height - EDGE_SIZE); y = y + 8) { |
|
|
for (x = 0; x < (edged_width - EDGE_SIZE); x = x + 8) { |
|
|
interpolate8x8_6tap_lowpass_v(hv_ptr, h_ptr, edged_width, rounding); |
|
|
hv_ptr += 8; |
|
|
h_ptr += 8; |
|
|
} |
|
|
|
|
|
hv_ptr += EDGE_SIZE2; |
|
|
h_ptr += EDGE_SIZE2; |
|
|
|
|
|
hv_ptr += stride_add; |
|
|
h_ptr += stride_add; |
|
|
} |
|
|
} |
|
|
else { |
|
|
|
|
|
for (y = 0; y < (edged_height - EDGE_SIZE); y += 8) { |
|
|
for (x = 0; x < (edged_width - EDGE_SIZE); x += 8) { |
|
291 |
interpolate8x8_halfpel_h(h_ptr, n_ptr, edged_width, rounding); |
interpolate8x8_halfpel_h(h_ptr, n_ptr, edged_width, rounding); |
292 |
interpolate8x8_halfpel_v(v_ptr, n_ptr, edged_width, rounding); |
interpolate8x8_halfpel_v(v_ptr, n_ptr, edged_width, rounding); |
293 |
interpolate8x8_halfpel_hv(hv_ptr, n_ptr, edged_width, rounding); |
interpolate8x8_halfpel_hv(hv_ptr, n_ptr, edged_width, rounding); |
297 |
v_ptr += 8; |
v_ptr += 8; |
298 |
hv_ptr += 8; |
hv_ptr += 8; |
299 |
} |
} |
|
|
|
|
h_ptr += EDGE_SIZE; |
|
|
v_ptr += EDGE_SIZE; |
|
|
hv_ptr += EDGE_SIZE; |
|
|
n_ptr += EDGE_SIZE; |
|
|
|
|
300 |
h_ptr += stride_add; |
h_ptr += stride_add; |
301 |
v_ptr += stride_add; |
v_ptr += stride_add; |
302 |
hv_ptr += stride_add; |
hv_ptr += stride_add; |
303 |
n_ptr += stride_add; |
n_ptr += stride_add; |
304 |
} |
} |
305 |
} |
|
|
/* |
|
306 |
#ifdef BFRAMES |
#ifdef BFRAMES |
307 |
n_ptr = refn->u; |
n_ptr = refn->u; |
308 |
h_ptr = refh->u; |
h_ptr = refh->u; |
314 |
v_ptr -= offset2; |
v_ptr -= offset2; |
315 |
hv_ptr -= offset2; |
hv_ptr -= offset2; |
316 |
|
|
317 |
for (y = 0; y < edged_height2; y += 8) { |
for (y = 0; y < edged_height2; y = y + 8) { |
318 |
for (x = 0; x < edged_width2; x += 8) { |
for (x = 0; x < edged_width2; x = x + 8) { |
319 |
interpolate8x8_halfpel_h(h_ptr, n_ptr, edged_width2, rounding); |
interpolate8x8_halfpel_h(h_ptr, n_ptr, edged_width2, rounding); |
320 |
interpolate8x8_halfpel_v(v_ptr, n_ptr, edged_width2, rounding); |
interpolate8x8_halfpel_v(v_ptr, n_ptr, edged_width2, rounding); |
321 |
interpolate8x8_halfpel_hv(hv_ptr, n_ptr, edged_width2, rounding); |
interpolate8x8_halfpel_hv(hv_ptr, n_ptr, edged_width2, rounding); |
358 |
n_ptr += stride_add2; |
n_ptr += stride_add2; |
359 |
} |
} |
360 |
#endif |
#endif |
361 |
*/ |
|
362 |
/* |
/* |
363 |
interpolate_halfpel_h( |
interpolate_halfpel_h( |
364 |
refh->y - offset, |
refh->y - offset, |
512 |
height = -height; |
height = -height; |
513 |
} |
} |
514 |
|
|
|
// --- xvid 2.1 compatiblity patch --- |
|
|
// --- remove when xvid_dec_frame->stride equals real stride |
|
|
if ((csp & ~XVID_CSP_VFLIP) == XVID_CSP_RGB555 || |
|
|
(csp & ~XVID_CSP_VFLIP) == XVID_CSP_RGB565 || |
|
|
(csp & ~XVID_CSP_VFLIP) == XVID_CSP_YUY2 || |
|
|
(csp & ~XVID_CSP_VFLIP) == XVID_CSP_YVYU || |
|
|
(csp & ~XVID_CSP_VFLIP) == XVID_CSP_UYVY) |
|
|
{ |
|
|
dst_stride *= 2; |
|
|
} |
|
|
else if ((csp & ~XVID_CSP_VFLIP) == XVID_CSP_RGB24) |
|
|
{ |
|
|
dst_stride *= 3; |
|
|
} |
|
|
else if ((csp & ~XVID_CSP_VFLIP) == XVID_CSP_RGB32 || |
|
|
(csp & ~XVID_CSP_VFLIP) == XVID_CSP_ABGR || |
|
|
(csp & ~XVID_CSP_VFLIP) == XVID_CSP_RGBA) |
|
|
{ |
|
|
dst_stride *= 4; |
|
|
} |
|
|
// ^--- xvid 2.1 compatiblity fix ---^ |
|
|
|
|
|
|
|
515 |
switch (csp & ~XVID_CSP_VFLIP) { |
switch (csp & ~XVID_CSP_VFLIP) { |
516 |
case XVID_CSP_RGB555: |
case XVID_CSP_RGB555: |
517 |
yv12_to_rgb555(dst, dst_stride, image->y, image->u, image->v, |
yv12_to_rgb555(dst, dst_stride, image->y, image->u, image->v, |
533 |
edged_width, edged_width / 2, width, height); |
edged_width, edged_width / 2, width, height); |
534 |
return 0; |
return 0; |
535 |
|
|
|
case XVID_CSP_ABGR: |
|
|
yv12_to_abgr(dst, dst_stride, image->y, image->u, image->v, |
|
|
edged_width, edged_width / 2, width, height); |
|
|
return 0; |
|
|
|
|
|
case XVID_CSP_RGBA: |
|
|
yv12_to_rgba(dst, dst_stride, image->y, image->u, image->v, |
|
|
edged_width, edged_width / 2, width, height); |
|
|
return 0; |
|
|
|
|
536 |
case XVID_CSP_I420: |
case XVID_CSP_I420: |
537 |
yv12_to_yuv(dst, dst_stride, image->y, image->u, image->v, edged_width, |
yv12_to_yuv(dst, dst_stride, image->y, image->u, image->v, edged_width, |
538 |
edged_width / 2, width, height); |
edged_width / 2, width, height); |