Rev 1085 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed
Rev | Author | Line No. | Line |
---|---|---|---|
1085 | pj | 1 | /* store.c, picture output routines */ |
2 | |||
3 | /* Copyright (C) 1996, MPEG Software Simulation Group. All Rights Reserved. */ |
||
4 | |||
5 | /* |
||
6 | * Disclaimer of Warranty |
||
7 | * |
||
8 | * These software programs are available to the user without any license fee or |
||
9 | * royalty on an "as is" basis. The MPEG Software Simulation Group disclaims |
||
10 | * any and all warranties, whether express, implied, or statuary, including any |
||
11 | * implied warranties or merchantability or of fitness for a particular |
||
12 | * purpose. In no event shall the copyright-holder be liable for any |
||
13 | * incidental, punitive, or consequential damages of any kind whatsoever |
||
14 | * arising from the use of these programs. |
||
15 | * |
||
16 | * This disclaimer of warranty extends to the user of these programs and user's |
||
17 | * customers, employees, agents, transferees, successors, and assigns. |
||
18 | * |
||
19 | * The MPEG Software Simulation Group does not represent or warrant that the |
||
20 | * programs furnished hereunder are free of infringement of any third-party |
||
21 | * patents. |
||
22 | * |
||
23 | * Commercial implementations of MPEG-1 and MPEG-2 video, including shareware, |
||
24 | * are subject to royalty fees to patent holders. Many of these patents are |
||
25 | * general enough such that they are unavoidable regardless of implementation |
||
26 | * design. |
||
27 | * |
||
28 | */ |
||
29 | |||
30 | #include <stdio.h> |
||
31 | #include <stdlib.h> |
||
32 | #include <fcntl.h> |
||
33 | |||
34 | #include "config.h" |
||
35 | #include "global.h" |
||
36 | |||
37 | /* private prototypes */ |
||
38 | static void store_one _ANSI_ARGS_((char *outname, unsigned char *src[], |
||
39 | int offset, int incr, int height)); |
||
40 | static void store_yuv _ANSI_ARGS_((char *outname, unsigned char *src[], |
||
41 | int offset, int incr, int height)); |
||
42 | static void store_sif _ANSI_ARGS_((char *outname, unsigned char *src[], |
||
43 | int offset, int incr, int height)); |
||
44 | static void store_ppm_tga _ANSI_ARGS_((char *outname, unsigned char *src[], |
||
45 | int offset, int incr, int height, int tgaflag)); |
||
46 | static void store_yuv1 _ANSI_ARGS_((char *name, unsigned char *src, |
||
47 | int offset, int incr, int width, int height)); |
||
48 | static void putbyte _ANSI_ARGS_((int c)); |
||
49 | static void putword _ANSI_ARGS_((int w)); |
||
50 | static void conv422to444 _ANSI_ARGS_((unsigned char *src, unsigned char *dst)); |
||
51 | static void conv420to422 _ANSI_ARGS_((unsigned char *src, unsigned char *dst)); |
||
52 | |||
53 | #define OBFRSIZE 4096 |
||
54 | static unsigned char obfr[OBFRSIZE]; |
||
55 | static unsigned char *optr; |
||
56 | static int outfile; |
||
57 | |||
58 | /* |
||
59 | * store a picture as either one frame or two fields |
||
60 | */ |
||
61 | void Write_Frame(src,frame) |
||
62 | unsigned char *src[]; |
||
63 | int frame; |
||
64 | { |
||
65 | char outname[FILENAME_LENGTH]; |
||
66 | |||
67 | if (progressive_sequence || progressive_frame || Frame_Store_Flag) |
||
68 | { |
||
69 | /* progressive */ |
||
70 | sprintf(outname,Output_Picture_Filename,frame,'f'); |
||
71 | store_one(outname,src,0,Coded_Picture_Width,vertical_size); |
||
72 | } |
||
73 | else |
||
74 | { |
||
75 | /* interlaced */ |
||
76 | sprintf(outname,Output_Picture_Filename,frame,'a'); |
||
77 | store_one(outname,src,0,Coded_Picture_Width<<1,vertical_size>>1); |
||
78 | |||
79 | sprintf(outname,Output_Picture_Filename,frame,'b'); |
||
80 | store_one(outname,src, |
||
81 | Coded_Picture_Width,Coded_Picture_Width<<1,vertical_size>>1); |
||
82 | } |
||
83 | } |
||
84 | |||
85 | /* |
||
86 | * store one frame or one field |
||
87 | */ |
||
88 | static void store_one(outname,src,offset,incr,height) |
||
89 | char *outname; |
||
90 | unsigned char *src[]; |
||
91 | int offset, incr, height; |
||
92 | { |
||
93 | switch (Output_Type) |
||
94 | { |
||
95 | case T_YUV: |
||
96 | store_yuv(outname,src,offset,incr,height); |
||
97 | break; |
||
98 | case T_SIF: |
||
99 | store_sif(outname,src,offset,incr,height); |
||
100 | break; |
||
101 | case T_TGA: |
||
102 | store_ppm_tga(outname,src,offset,incr,height,1); |
||
103 | break; |
||
104 | case T_PPM: |
||
105 | store_ppm_tga(outname,src,offset,incr,height,0); |
||
106 | break; |
||
107 | #ifdef DISPLAY |
||
108 | case T_X11: |
||
109 | dither(src); |
||
110 | break; |
||
111 | #endif |
||
112 | default: |
||
113 | break; |
||
114 | } |
||
115 | } |
||
116 | |||
117 | /* separate headerless files for y, u and v */ |
||
118 | static void store_yuv(outname,src,offset,incr,height) |
||
119 | char *outname; |
||
120 | unsigned char *src[]; |
||
121 | int offset,incr,height; |
||
122 | { |
||
123 | int hsize; |
||
124 | char tmpname[FILENAME_LENGTH]; |
||
125 | |||
126 | hsize = horizontal_size; |
||
127 | |||
128 | sprintf(tmpname,"%s.Y",outname); |
||
129 | store_yuv1(tmpname,src[0],offset,incr,hsize,height); |
||
130 | |||
131 | if (chroma_format!=CHROMA444) |
||
132 | { |
||
133 | offset>>=1; incr>>=1; hsize>>=1; |
||
134 | } |
||
135 | |||
136 | if (chroma_format==CHROMA420) |
||
137 | { |
||
138 | height>>=1; |
||
139 | } |
||
140 | |||
141 | sprintf(tmpname,"%s.U",outname); |
||
142 | store_yuv1(tmpname,src[1],offset,incr,hsize,height); |
||
143 | |||
144 | sprintf(tmpname,"%s.V",outname); |
||
145 | store_yuv1(tmpname,src[2],offset,incr,hsize,height); |
||
146 | } |
||
147 | |||
148 | /* auxiliary routine */ |
||
149 | static void store_yuv1(name,src,offset,incr,width,height) |
||
150 | char *name; |
||
151 | unsigned char *src; |
||
152 | int offset,incr,width,height; |
||
153 | { |
||
154 | int i, j; |
||
155 | unsigned char *p; |
||
156 | |||
157 | if (!Quiet_Flag) |
||
158 | fprintf(stderr,"saving %s\n",name); |
||
159 | |||
160 | if ((outfile = open(name,O_CREAT|O_TRUNC|O_WRONLY|O_BINARY,0666))==-1) |
||
161 | { |
||
162 | sprintf(Error_Text,"Couldn't create %s\n",name); |
||
163 | Error(Error_Text); |
||
164 | } |
||
165 | |||
166 | optr=obfr; |
||
167 | |||
168 | for (i=0; i<height; i++) |
||
169 | { |
||
170 | p = src + offset + incr*i; |
||
171 | for (j=0; j<width; j++) |
||
172 | putbyte(*p++); |
||
173 | } |
||
174 | |||
175 | if (optr!=obfr) |
||
176 | write(outfile,obfr,optr-obfr); |
||
177 | |||
178 | close(outfile); |
||
179 | } |
||
180 | |||
181 | /* |
||
182 | * store as headerless file in U,Y,V,Y format |
||
183 | */ |
||
184 | static void store_sif (outname,src,offset,incr,height) |
||
185 | char *outname; |
||
186 | unsigned char *src[]; |
||
187 | int offset, incr, height; |
||
188 | { |
||
189 | int i,j; |
||
190 | unsigned char *py, *pu, *pv; |
||
191 | static unsigned char *u422, *v422; |
||
192 | |||
193 | if (chroma_format==CHROMA444) |
||
194 | Error("4:4:4 not supported for SIF format"); |
||
195 | |||
196 | if (chroma_format==CHROMA422) |
||
197 | { |
||
198 | u422 = src[1]; |
||
199 | v422 = src[2]; |
||
200 | } |
||
201 | else |
||
202 | { |
||
203 | if (!u422) |
||
204 | { |
||
205 | if (!(u422 = (unsigned char *)malloc((Coded_Picture_Width>>1) |
||
206 | *Coded_Picture_Height))) |
||
207 | Error("malloc failed"); |
||
208 | if (!(v422 = (unsigned char *)malloc((Coded_Picture_Width>>1) |
||
209 | *Coded_Picture_Height))) |
||
210 | Error("malloc failed"); |
||
211 | } |
||
212 | |||
213 | conv420to422(src[1],u422); |
||
214 | conv420to422(src[2],v422); |
||
215 | } |
||
216 | |||
217 | strcat(outname,".SIF"); |
||
218 | |||
219 | if (!Quiet_Flag) |
||
220 | fprintf(stderr,"saving %s\n",outname); |
||
221 | |||
222 | if ((outfile = open(outname,O_CREAT|O_TRUNC|O_WRONLY|O_BINARY,0666))==-1) |
||
223 | { |
||
224 | sprintf(Error_Text,"Couldn't create %s\n",outname); |
||
225 | Error(Error_Text); |
||
226 | } |
||
227 | |||
228 | optr = obfr; |
||
229 | |||
230 | for (i=0; i<height; i++) |
||
231 | { |
||
232 | py = src[0] + offset + incr*i; |
||
233 | pu = u422 + (offset>>1) + (incr>>1)*i; |
||
234 | pv = v422 + (offset>>1) + (incr>>1)*i; |
||
235 | |||
236 | for (j=0; j<horizontal_size; j+=2) |
||
237 | { |
||
238 | putbyte(*pu++); |
||
239 | putbyte(*py++); |
||
240 | putbyte(*pv++); |
||
241 | putbyte(*py++); |
||
242 | } |
||
243 | } |
||
244 | |||
245 | if (optr!=obfr) |
||
246 | write(outfile,obfr,optr-obfr); |
||
247 | |||
248 | close(outfile); |
||
249 | } |
||
250 | |||
251 | /* |
||
252 | * store as PPM (PBMPLUS) or uncompressed Truevision TGA ('Targa') file |
||
253 | */ |
||
254 | static void store_ppm_tga(outname,src,offset,incr,height,tgaflag) |
||
255 | char *outname; |
||
256 | unsigned char *src[]; |
||
257 | int offset, incr, height; |
||
258 | int tgaflag; |
||
259 | { |
||
260 | int i, j; |
||
261 | int y, u, v, r, g, b; |
||
262 | int crv, cbu, cgu, cgv; |
||
263 | unsigned char *py, *pu, *pv; |
||
264 | static unsigned char tga24[14] = {0,0,2,0,0,0,0, 0,0,0,0,0,24,32}; |
||
265 | char header[FILENAME_LENGTH]; |
||
266 | static unsigned char *u422, *v422, *u444, *v444; |
||
267 | |||
268 | if (chroma_format==CHROMA444) |
||
269 | { |
||
270 | u444 = src[1]; |
||
271 | v444 = src[2]; |
||
272 | } |
||
273 | else |
||
274 | { |
||
275 | if (!u444) |
||
276 | { |
||
277 | if (chroma_format==CHROMA420) |
||
278 | { |
||
279 | if (!(u422 = (unsigned char *)malloc((Coded_Picture_Width>>1) |
||
280 | *Coded_Picture_Height))) |
||
281 | Error("malloc failed"); |
||
282 | if (!(v422 = (unsigned char *)malloc((Coded_Picture_Width>>1) |
||
283 | *Coded_Picture_Height))) |
||
284 | Error("malloc failed"); |
||
285 | } |
||
286 | |||
287 | if (!(u444 = (unsigned char *)malloc(Coded_Picture_Width |
||
288 | *Coded_Picture_Height))) |
||
289 | Error("malloc failed"); |
||
290 | |||
291 | if (!(v444 = (unsigned char *)malloc(Coded_Picture_Width |
||
292 | *Coded_Picture_Height))) |
||
293 | Error("malloc failed"); |
||
294 | } |
||
295 | |||
296 | if (chroma_format==CHROMA420) |
||
297 | { |
||
298 | conv420to422(src[1],u422); |
||
299 | conv420to422(src[2],v422); |
||
300 | conv422to444(u422,u444); |
||
301 | conv422to444(v422,v444); |
||
302 | } |
||
303 | else |
||
304 | { |
||
305 | conv422to444(src[1],u444); |
||
306 | conv422to444(src[2],v444); |
||
307 | } |
||
308 | } |
||
309 | |||
310 | strcat(outname,tgaflag ? ".tga" : ".ppm"); |
||
311 | |||
312 | if (!Quiet_Flag) |
||
313 | fprintf(stderr,"saving %s\n",outname); |
||
314 | |||
315 | if ((outfile = open(outname,O_CREAT|O_TRUNC|O_WRONLY|O_BINARY,0666))==-1) |
||
316 | { |
||
317 | sprintf(Error_Text,"Couldn't create %s\n",outname); |
||
318 | Error(Error_Text); |
||
319 | } |
||
320 | |||
321 | optr = obfr; |
||
322 | |||
323 | if (tgaflag) |
||
324 | { |
||
325 | /* TGA header */ |
||
326 | for (i=0; i<12; i++) |
||
327 | putbyte(tga24[i]); |
||
328 | |||
329 | putword(horizontal_size); putword(height); |
||
330 | putbyte(tga24[12]); putbyte(tga24[13]); |
||
331 | } |
||
332 | else |
||
333 | { |
||
334 | /* PPM header */ |
||
335 | sprintf(header,"P6\n%d %d\n255\n",horizontal_size,height); |
||
336 | |||
337 | for (i=0; header[i]!=0; i++) |
||
338 | putbyte(header[i]); |
||
339 | } |
||
340 | |||
341 | /* matrix coefficients */ |
||
342 | crv = Inverse_Table_6_9[matrix_coefficients][0]; |
||
343 | cbu = Inverse_Table_6_9[matrix_coefficients][1]; |
||
344 | cgu = Inverse_Table_6_9[matrix_coefficients][2]; |
||
345 | cgv = Inverse_Table_6_9[matrix_coefficients][3]; |
||
346 | |||
347 | for (i=0; i<height; i++) |
||
348 | { |
||
349 | py = src[0] + offset + incr*i; |
||
350 | pu = u444 + offset + incr*i; |
||
351 | pv = v444 + offset + incr*i; |
||
352 | |||
353 | for (j=0; j<horizontal_size; j++) |
||
354 | { |
||
355 | u = *pu++ - 128; |
||
356 | v = *pv++ - 128; |
||
357 | y = 76309 * (*py++ - 16); /* (255/219)*65536 */ |
||
358 | r = Clip[(y + crv*v + 32768)>>16]; |
||
359 | g = Clip[(y - cgu*u - cgv*v + 32768)>>16]; |
||
360 | b = Clip[(y + cbu*u + 32786)>>16]; |
||
361 | |||
362 | if (tgaflag) |
||
363 | { |
||
364 | putbyte(b); putbyte(g); putbyte(r); |
||
365 | } |
||
366 | else |
||
367 | { |
||
368 | putbyte(r); putbyte(g); putbyte(b); |
||
369 | } |
||
370 | } |
||
371 | } |
||
372 | |||
373 | if (optr!=obfr) |
||
374 | write(outfile,obfr,optr-obfr); |
||
375 | |||
376 | close(outfile); |
||
377 | } |
||
378 | |||
379 | static void putbyte(c) |
||
380 | int c; |
||
381 | { |
||
382 | *optr++ = c; |
||
383 | |||
384 | if (optr == obfr+OBFRSIZE) |
||
385 | { |
||
386 | write(outfile,obfr,OBFRSIZE); |
||
387 | optr = obfr; |
||
388 | } |
||
389 | } |
||
390 | |||
391 | static void putword(w) |
||
392 | int w; |
||
393 | { |
||
394 | putbyte(w); putbyte(w>>8); |
||
395 | } |
||
396 | |||
397 | /* horizontal 1:2 interpolation filter */ |
||
398 | static void conv422to444(src,dst) |
||
399 | unsigned char *src,*dst; |
||
400 | { |
||
401 | int i, i2, w, j, im3, im2, im1, ip1, ip2, ip3; |
||
402 | |||
403 | w = Coded_Picture_Width>>1; |
||
404 | |||
405 | if (base.MPEG2_Flag) |
||
406 | { |
||
407 | for (j=0; j<Coded_Picture_Height; j++) |
||
408 | { |
||
409 | for (i=0; i<w; i++) |
||
410 | { |
||
411 | i2 = i<<1; |
||
412 | im2 = (i<2) ? 0 : i-2; |
||
413 | im1 = (i<1) ? 0 : i-1; |
||
414 | ip1 = (i<w-1) ? i+1 : w-1; |
||
415 | ip2 = (i<w-2) ? i+2 : w-1; |
||
416 | ip3 = (i<w-3) ? i+3 : w-1; |
||
417 | |||
418 | /* FIR filter coefficients (*256): 21 0 -52 0 159 256 159 0 -52 0 21 */ |
||
419 | /* even samples (0 0 256 0 0) */ |
||
420 | dst[i2] = src[i]; |
||
421 | |||
422 | /* odd samples (21 -52 159 159 -52 21) */ |
||
423 | dst[i2+1] = Clip[(int)(21*(src[im2]+src[ip3]) |
||
424 | -52*(src[im1]+src[ip2]) |
||
425 | +159*(src[i]+src[ip1])+128)>>8]; |
||
426 | } |
||
427 | src+= w; |
||
428 | dst+= Coded_Picture_Width; |
||
429 | } |
||
430 | } |
||
431 | else |
||
432 | { |
||
433 | for (j=0; j<Coded_Picture_Height; j++) |
||
434 | { |
||
435 | for (i=0; i<w; i++) |
||
436 | { |
||
437 | |||
438 | i2 = i<<1; |
||
439 | im3 = (i<3) ? 0 : i-3; |
||
440 | im2 = (i<2) ? 0 : i-2; |
||
441 | im1 = (i<1) ? 0 : i-1; |
||
442 | ip1 = (i<w-1) ? i+1 : w-1; |
||
443 | ip2 = (i<w-2) ? i+2 : w-1; |
||
444 | ip3 = (i<w-3) ? i+3 : w-1; |
||
445 | |||
446 | /* FIR filter coefficients (*256): 5 -21 70 228 -37 11 */ |
||
447 | dst[i2] = Clip[(int)( 5*src[im3] |
||
448 | -21*src[im2] |
||
449 | +70*src[im1] |
||
450 | +228*src[i] |
||
451 | -37*src[ip1] |
||
452 | +11*src[ip2]+128)>>8]; |
||
453 | |||
454 | dst[i2+1] = Clip[(int)( 5*src[ip3] |
||
455 | -21*src[ip2] |
||
456 | +70*src[ip1] |
||
457 | +228*src[i] |
||
458 | -37*src[im1] |
||
459 | +11*src[im2]+128)>>8]; |
||
460 | } |
||
461 | src+= w; |
||
462 | dst+= Coded_Picture_Width; |
||
463 | } |
||
464 | } |
||
465 | } |
||
466 | |||
467 | /* vertical 1:2 interpolation filter */ |
||
468 | static void conv420to422(src,dst) |
||
469 | unsigned char *src,*dst; |
||
470 | { |
||
471 | int w, h, i, j, j2; |
||
472 | int jm6, jm5, jm4, jm3, jm2, jm1, jp1, jp2, jp3, jp4, jp5, jp6, jp7; |
||
473 | |||
474 | w = Coded_Picture_Width>>1; |
||
475 | h = Coded_Picture_Height>>1; |
||
476 | |||
477 | if (progressive_frame) |
||
478 | { |
||
479 | /* intra frame */ |
||
480 | for (i=0; i<w; i++) |
||
481 | { |
||
482 | for (j=0; j<h; j++) |
||
483 | { |
||
484 | j2 = j<<1; |
||
485 | jm3 = (j<3) ? 0 : j-3; |
||
486 | jm2 = (j<2) ? 0 : j-2; |
||
487 | jm1 = (j<1) ? 0 : j-1; |
||
488 | jp1 = (j<h-1) ? j+1 : h-1; |
||
489 | jp2 = (j<h-2) ? j+2 : h-1; |
||
490 | jp3 = (j<h-3) ? j+3 : h-1; |
||
491 | |||
492 | /* FIR filter coefficients (*256): 5 -21 70 228 -37 11 */ |
||
493 | /* New FIR filter coefficients (*256): 3 -16 67 227 -32 7 */ |
||
494 | dst[w*j2] = Clip[(int)( 3*src[w*jm3] |
||
495 | -16*src[w*jm2] |
||
496 | +67*src[w*jm1] |
||
497 | +227*src[w*j] |
||
498 | -32*src[w*jp1] |
||
499 | +7*src[w*jp2]+128)>>8]; |
||
500 | |||
501 | dst[w*(j2+1)] = Clip[(int)( 3*src[w*jp3] |
||
502 | -16*src[w*jp2] |
||
503 | +67*src[w*jp1] |
||
504 | +227*src[w*j] |
||
505 | -32*src[w*jm1] |
||
506 | +7*src[w*jm2]+128)>>8]; |
||
507 | } |
||
508 | src++; |
||
509 | dst++; |
||
510 | } |
||
511 | } |
||
512 | else |
||
513 | { |
||
514 | /* intra field */ |
||
515 | for (i=0; i<w; i++) |
||
516 | { |
||
517 | for (j=0; j<h; j+=2) |
||
518 | { |
||
519 | j2 = j<<1; |
||
520 | |||
521 | /* top field */ |
||
522 | jm6 = (j<6) ? 0 : j-6; |
||
523 | jm4 = (j<4) ? 0 : j-4; |
||
524 | jm2 = (j<2) ? 0 : j-2; |
||
525 | jp2 = (j<h-2) ? j+2 : h-2; |
||
526 | jp4 = (j<h-4) ? j+4 : h-2; |
||
527 | jp6 = (j<h-6) ? j+6 : h-2; |
||
528 | |||
529 | /* Polyphase FIR filter coefficients (*256): 2 -10 35 242 -18 5 */ |
||
530 | /* New polyphase FIR filter coefficients (*256): 1 -7 30 248 -21 5 */ |
||
531 | dst[w*j2] = Clip[(int)( 1*src[w*jm6] |
||
532 | -7*src[w*jm4] |
||
533 | +30*src[w*jm2] |
||
534 | +248*src[w*j] |
||
535 | -21*src[w*jp2] |
||
536 | +5*src[w*jp4]+128)>>8]; |
||
537 | |||
538 | /* Polyphase FIR filter coefficients (*256): 11 -38 192 113 -30 8 */ |
||
539 | /* New polyphase FIR filter coefficients (*256):7 -35 194 110 -24 4 */ |
||
540 | dst[w*(j2+2)] = Clip[(int)( 7*src[w*jm4] |
||
541 | -35*src[w*jm2] |
||
542 | +194*src[w*j] |
||
543 | +110*src[w*jp2] |
||
544 | -24*src[w*jp4] |
||
545 | +4*src[w*jp6]+128)>>8]; |
||
546 | |||
547 | /* bottom field */ |
||
548 | jm5 = (j<5) ? 1 : j-5; |
||
549 | jm3 = (j<3) ? 1 : j-3; |
||
550 | jm1 = (j<1) ? 1 : j-1; |
||
551 | jp1 = (j<h-1) ? j+1 : h-1; |
||
552 | jp3 = (j<h-3) ? j+3 : h-1; |
||
553 | jp5 = (j<h-5) ? j+5 : h-1; |
||
554 | jp7 = (j<h-7) ? j+7 : h-1; |
||
555 | |||
556 | /* Polyphase FIR filter coefficients (*256): 11 -38 192 113 -30 8 */ |
||
557 | /* New polyphase FIR filter coefficients (*256):7 -35 194 110 -24 4 */ |
||
558 | dst[w*(j2+1)] = Clip[(int)( 7*src[w*jp5] |
||
559 | -35*src[w*jp3] |
||
560 | +194*src[w*jp1] |
||
561 | +110*src[w*jm1] |
||
562 | -24*src[w*jm3] |
||
563 | +4*src[w*jm5]+128)>>8]; |
||
564 | |||
565 | dst[w*(j2+3)] = Clip[(int)( 1*src[w*jp7] |
||
566 | -7*src[w*jp5] |
||
567 | +30*src[w*jp3] |
||
568 | +248*src[w*jp1] |
||
569 | -21*src[w*jm1] |
||
570 | +5*src[w*jm3]+128)>>8]; |
||
571 | } |
||
572 | src++; |
||
573 | dst++; |
||
574 | } |
||
575 | } |
||
576 | } |