about summary refs log tree commit diff
path: root/converter/ppm/ppmtompeg/jpeg.c
blob: 6680404b0c44f3739abb05b7ce25ea83b477cfe6 (plain) (blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
/*===========================================================================*
 * jpeg.c
 *
 *  procedures to deal with JPEG files
 *
 * EXPORTED PROCEDURES:
 *  JMovie2JPEG
 *      ReadJPEG
 *
 *===========================================================================*/

/* COPYRIGHT INFORMATION IS AT THE END OF THIS FILE */

/*==============*
 * HEADER FILES *
 *==============*/
#define _XOPEN_SOURCE    /* Make sure stdio.h contains fileno() */
#include <unistd.h>
#include <stdio.h>
/* With the lossless jpeg patch applied to the Jpeg library
    (ftp://ftp.wizards.dupont.com/pub/ImageMagick/delegates/ljpeg-6b.tar.gz),
    the name of min_DCT_scaled_size changes to min_codec_data_unit,
    for some reason.  With this macro, we change it back.
*/
#define min_codec_data_unit min_DCT_scaled_size
#include <jpeglib.h>
#undef min_codec_data_unit

#include "netpbm/pm_config.h"
#include "netpbm/pm_c_util.h"
#include "netpbm/mallocvar.h"
#include "netpbm/pm.h"

#include "fsize.h"
#include "frame.h"

#include "jpeg.h"


/* make it happier.... */
#undef DCTSIZE2

/* jcopy_sample_rows() is an internal routine in the JPEG library, not
   meant for use by us.  We should figure out what the official interface
   for this is and use it.  The following is copied out of jpegint.h, which
   is part of the JPEG library source code.
   */
extern void jcopy_sample_rows JPP((JSAMPARRAY input_array, int source_row,
                    JSAMPARRAY output_array, int dest_row,
                    int num_rows, JDIMENSION num_cols));


#define HEADER_SIZE 607   /*JFIF header size used on output images*/


static int
minDctVScaledSize(struct jpeg_decompress_struct const cinfo) {

#if JPEG_LIB_VERSION >= 70
    return cinfo.min_DCT_v_scaled_size;
#else
    return cinfo.min_DCT_scaled_size;
#endif
}



/*=======================================================================*
 *                                                                       *
 * JMovie2JPEG                                                           *
 *                                                                       *
 *      Splits up a Parallax J_Movie into a set of JFIF image files      *
 *                                                                       *
 * RETURNS:     nothing                                                  *
 *                                                                       *
 * SIDE EFFECTS:    none                                                 *
 *                                                                       *
 *   Contributed By James Boucher(jboucher@flash.bu.edu)                 *
 *               Boston University Multimedia Communications Lab         *
 * This code was adapted from the Berkeley Playone.c and Brian Smith's   *
 * JGetFrame code after being modified on 10-7-93 by Dinesh Venkatesh    *
 * of BU.                                                                *
 *       This code converts a version 2 Parallax J_Movie into a          *
 * set of JFIF compatible JPEG images. It works for all image            *
 * sizes and qualities.                                                  *
 ************************************************************************/
void
JMovie2JPEG(const char * const infilename,
                /* input filename string */
            const char * const obase,
                /* output filename base string=>obase##.jpg */
            int          const start,
                /* first frame to be extracted */
            int          const end
                /* last frame to be extracted */
    ) {

    FILE *inFile;         /* Jmovie file pointer */
    FILE *outFile;        /* JPEG file pointer for output file */
    int fd, i;            /* input file descriptor and a counting variable*/
    char ofname[256];     /* output filename string */
    int Temp = 0, temp = 0;   /* dummy variables */
    int image_offset = 0;     /* counting variable */
    /* J_Movie header information */
    int ver_no;           /* version number - expected to be 2 */
    int fps;              /* frame rate - frames per second */
    int no_frames;        /* total number of frames in jmovie */
    int bandwidth;        /* bandwidth required for normal playback*/
    int qfactor;          /* quality factor used to scale Q matrix */
    int mapsize;          /* number of color map entries - 2^24 */
    int audio_tracks;     /* number of audio tracks ==1    */
    int audiosize;        /*number of bytes in audio tracks */
    int *inoffsets;       /* input frame offsets from start of jmovie*/
    int width;            /* image width */
    int height;           /* image height */
    int size;             /* total image size in bytes */
    char op_code;         /* jmovie op_code */
    char jpeg_size[4];        /* jpeg data size */
    static char junk[1000];   /* data sink for audio data */
    int last;             /* Number of last frame we will process */

    /* The next array represents the default JFIF header for
       quality = 100 and size = 320x240. The values are
       adjusted as the J_Movie header is read.  The default
       size of this array is set large so as to make room
       for the appending of the jpeg bitstream. It can be
       made smaller if you have a better idea of its expected size
    */
    static unsigned char inbuffer[300000] = {
        0xFF, 0xD8, 0xFF, 0xE0, 0x00, 0x10, 0x4A, 0x46,
        0x49, 0x46, 0x00, 0x01, 0x01, 0x00, 0x00, 0x01,
        0x00, 0x01, 0x00, 0x00, 0xFF, 0xC0, 0x00, 0x11,
        0x08, 0x00, 0xF0, 0x01, 0x40, 0x03, 0x01, 0x21,
        0x00, 0x02, 0x11, 0x01, 0x03, 0x11, 0x01, 0xFF,
        0xDB, 0x00, 0x84, 0x00, 0x10, 0x0B, 0x0C, 0x0E,
        0x0C, 0x0A, 0x10, 0x0E, 0x0D, 0x0E, 0x12,
        0x11, 0x10, 0x13, 0x18, 0x28, 0x1A, 0x18, 0x16,
        0x16, 0x18, 0x31, 0x23, 0x25, 0x1D, 0x28, 0x3A,
        0x33, 0x3D, 0x3C, 0x39, 0x33, 0x38, 0x37, 0x40,
        0x48, 0x5C, 0x4E, 0x40, 0x44, 0x57, 0x45, 0x37,
        0x38, 0x50, 0x6D, 0x51, 0x57, 0x5F, 0x62, 0x67,
        0x68, 0x67, 0x3E, 0x4D, 0x71, 0x79, 0x70, 0x64,
        0x78, 0x5C, 0x65, 0x67, 0x63, 0x01, 0x11, 0x12,
        0x12, 0x18, 0x15, 0x18, 0x2F, 0x1A, 0x1A, 0x2F,
        0x63, 0x42, 0x38, 0x42, 0x63, 0x63, 0x63, 0x63,
        0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0x63,
        0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0x63,
        0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0x63,
        0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0x63,
        0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0x63,
        0x63, 0x63, 0x63, 0x63, 0x63, 0x63, 0xFF, 0xC4,
        0x01, 0xA2, 0x00, 0x00, 0x01, 0x05, 0x01, 0x01,
        0x01, 0x01, 0x01, 0x01, 0x00, 0x00, 0x00, 0x00,
        0x00, 0x00, 0x00, 0x00, 0x01, 0x02, 0x03, 0x04,
        0x05, 0x06, 0x07, 0x08, 0x09, 0x0A, 0x0B,
        0x10, 0x00, 0x02, 0x01, 0x03, 0x03, 0x02, 0x04,
        0x03, 0x05, 0x05, 0x04, 0x04, 0x00, 0x00, 0x01,
        0x7D, 0x01, 0x02, 0x03, 0x00, 0x04, 0x11, 0x05,
        0x12, 0x21, 0x31, 0x41, 0x06, 0x13, 0x51, 0x61,
        0x07, 0x22, 0x71, 0x14, 0x32, 0x81, 0x91, 0xA1,
        0x08, 0x23, 0x42, 0xB1, 0xC1, 0x15, 0x52, 0xD1,
        0xF0, 0x24, 0x33, 0x62, 0x72, 0x82, 0x09,
        0x0A, 0x16, 0x17, 0x18, 0x19, 0x1A, 0x25, 0x26,
        0x27, 0x28, 0x29, 0x2A, 0x34, 0x35, 0x36, 0x37,
        0x38, 0x39, 0x3A, 0x43, 0x44, 0x45, 0x46, 0x47,
        0x48, 0x49, 0x4A, 0x53, 0x54, 0x55, 0x56, 0x57,
        0x58, 0x59, 0x5A, 0x63, 0x64, 0x65, 0x66, 0x67,
        0x68, 0x69, 0x6A, 0x73, 0x74, 0x75, 0x76, 0x77,
        0x78, 0x79, 0x7A, 0x83, 0x84, 0x85, 0x86, 0x87,
        0x88, 0x89, 0x8A, 0x92, 0x93, 0x94, 0x95, 0x96,
        0x97, 0x98, 0x99, 0x9A, 0xA2, 0xA3, 0xA4, 0xA5,
        0xA6, 0xA7, 0xA8, 0xA9, 0xAA, 0xB2, 0xB3, 0xB4,
        0xB5, 0xB6, 0xB7, 0xB8, 0xB9, 0xBA, 0xC2, 0xC3,
        0xC4, 0xC5, 0xC6, 0xC7, 0xC8, 0xC9, 0xCA, 0xD2,
        0xD3, 0xD4, 0xD5, 0xD6, 0xD7, 0xD8, 0xD9, 0xDA,
        0xE1, 0xE2, 0xE3, 0xE4, 0xE5, 0xE6, 0xE7, 0xE8,
        0xE9, 0xEA, 0xF1, 0xF2, 0xF3, 0xF4, 0xF5, 0xF6,
        0xF7, 0xF8, 0xF9, 0xFA, 0x01, 0x00, 0x03, 0x01,
        0x01, 0x01, 0x01, 0x01, 0x01, 0x01, 0x01, 0x01,
        0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x01, 0x02,
        0x03, 0x04, 0x05, 0x06, 0x07, 0x08, 0x09,
        0x0A, 0x0B, 0x11, 0x00, 0x02, 0x01, 0x02, 0x04,
        0x04, 0x03, 0x04, 0x07, 0x05, 0x04, 0x04, 0x00,
        0x01, 0x02, 0x77, 0x00, 0x01, 0x02, 0x03, 0x11,
        0x04, 0x05, 0x21, 0x31, 0x06, 0x12, 0x41, 0x51,
        0x07, 0x61, 0x71, 0x13, 0x22, 0x32, 0x81, 0x08,
        0x14, 0x42, 0x91, 0xA1, 0xB1, 0xC1, 0x09, 0x23,
        0x33, 0x52, 0xF0, 0x15, 0x62, 0x72, 0xD1,
        0x0A, 0x16, 0x24, 0x34, 0xE1, 0x25, 0xF1, 0x17,
        0x18, 0x19, 0x1A, 0x26, 0x27, 0x28, 0x29, 0x2A,
        0x35, 0x36, 0x37, 0x38, 0x39, 0x3A, 0x43, 0x44,
        0x45, 0x46, 0x47, 0x48, 0x49, 0x4A, 0x53, 0x54,
        0x55, 0x56, 0x57, 0x58, 0x59, 0x5A, 0x63, 0x64,
        0x65, 0x66, 0x67, 0x68, 0x69, 0x6A, 0x73, 0x74,
        0x75, 0x76, 0x77, 0x78, 0x79, 0x7A, 0x82, 0x83,
        0x84, 0x85, 0x86, 0x87, 0x88, 0x89, 0x8A, 0x92,
        0x93, 0x94, 0x95, 0x96, 0x97, 0x98, 0x99, 0x9A,
        0xA2, 0xA3, 0xA4, 0xA5, 0xA6, 0xA7, 0xA8, 0xA9,
        0xAA, 0xB2, 0xB3, 0xB4, 0xB5, 0xB6, 0xB7, 0xB8,
        0xB9, 0xBA, 0xC2, 0xC3, 0xC4, 0xC5, 0xC6, 0xC7,
        0xC8, 0xC9, 0xCA, 0xD2, 0xD3, 0xD4, 0xD5, 0xD6,
        0xD7, 0xD8, 0xD9, 0xDA, 0xE2, 0xE3, 0xE4, 0xE5,
        0xE6, 0xE7, 0xE8, 0xE9, 0xEA, 0xF2, 0xF3, 0xF4,
        0xF5, 0xF6, 0xF7, 0xF8, 0xF9, 0xFA, 0xFF, 0xDA,
        0x00, 0x0C, 0x03, 0x01, 0x00, 0x02, 0x11, 0x03,
        0x11, 0x00, 0x3F, 0x00

    };

    if (start > end) {
        fprintf(stderr,"bad frame numbers\n");
        exit(1);
    }

    /* open J_Movie */
    inFile = fopen(infilename, "rb");
    if (inFile == NULL) {
        perror (infilename);
        exit(1);
    }

    /* get file descriptor */
    fd = fileno(inFile);

    /* The following lines parse the jpeg_movie header and recover the */
    /* relevant information */

    fseek(inFile, 8 * sizeof(char), 0);

    if (fread(&ver_no, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading version");
        exit(1);
    }
    if (ver_no != 2) {
        perror("Unrecognized version - Quantization tables may be wrong\n");
    }
    if (fread(&fps, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading fps");
        exit(1);
    }
    if (fread (&no_frames, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading no_frames");
        exit(1);
    }

    MALLOCARRAY(inoffsets, no_frames);

    if (fread(&width, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading width");
        exit(1);
    }
    /* set image width in JFIF header */
    inbuffer[27] = (char)(0xFF & (width >> 8));
    inbuffer[28] = (char)(0xFF & width);

    if (fread(&height, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading height");
        exit(1);
    }
    /* set image height in JFIF header */
    inbuffer[25] = (char)(0xFF & (height >> 8));
    inbuffer[26] = (char)(0xFF & height);

    if (fread(&bandwidth, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading bandwidth");
        exit(1);
    }

    if (fread(&qfactor, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading qfactor");
        exit(1);
    }
    /* The default quality factor = 100, therefore, if
       our quality factor does not equal 100 we must
       scale the quantization matrices in the JFIF header
    */
    /* Note values are clipped to a max of 255 */
    if (qfactor != 100) {
        for (Temp = 44; Temp < 108; ++Temp) {
            temp= (inbuffer[Temp]*qfactor)/100;
            inbuffer[Temp] = (char)((temp<255) ? temp : 255);
        }
        for (Temp = 109; Temp < 173; ++Temp) {
            temp = (inbuffer[Temp]*qfactor)/100;
            inbuffer[Temp] = (char)((temp<255) ? temp : 255);
        }
    }

    if (fread(&mapsize, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading mapsize");
        exit(1);
    }
    if (fread (&image_offset, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading image offset");
        exit(1);
    }
    if (fread (&audio_tracks, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading audio tracks");
        exit(1);
    }

    fread(junk,sizeof(int),1,inFile);

    if (fread (&audiosize, sizeof(int), 1, inFile) != 1) {
        perror("Error in reading audiosize");
        exit(1);
    }

    fseek (inFile, image_offset, 0);

    last = MIN(end, no_frames-1);

    for (i = 0; i < no_frames; ++i)  {
        fread(&(inoffsets[i]), sizeof(int), 1, inFile);
    } /* Reads in the frame sizes into the array */

    rewind(inFile);

    /* Extract JFIF files from J_Movie */
    for (i = start; i <= last; ++i) {
        size = inoffsets[i] - inoffsets[i-1]- 5;
        lseek(fd, inoffsets[i-1], 0);
        read(fd, &(op_code), 1);
        while (op_code !=  0xffffffec) {
            read(fd,junk,audiosize);
            read(fd, &(op_code), 1);
            size = size - audiosize ;
        } /* To skip the audio bytes in each frame */
        read(fd, jpeg_size, 4);
        read(fd, &(inbuffer[607]), size);
        sprintf(ofname, "%s%d.jpg", obase, i);
        outFile = fopen(ofname, "wb");
        fwrite(inbuffer, (size+607), sizeof(char), outFile);
        fclose(outFile);
    }
    free(inoffsets);
    fclose(inFile);
}



/*===========================================================================*
 *
 * ReadJPEG  contributed by James Arthur Boucher of Boston University's
 *                                Multimedia Communications Lab
 *
 *      read a JPEG file and copy data into frame original data arrays
 *
 * RETURNS:     mf modified
 *
 * SIDE EFFECTS:    none
 *
 *===========================================================================*/
/*************************JPEG LIBRARY INTERFACE*********************/
/*
 * THE BIG PICTURE:
 *
 * The rough outline this JPEG decompression operation is:
 *
 *      allocate and initialize JPEG decompression object
 *      specify data source (eg, a file)
 *      jpeg_read_header();      obtain image dimensions and other parameters
 *      set parameters for decompression
 *      jpeg_start_decompress();
 *      while (scan lines remain to be read)
 *              jpeg_read_scanlines(...);
 *      jpeg_finish_decompress();
 *      release JPEG decompression object
 *
 */
void
ReadJPEG(MpegFrame * const mf,
         FILE *      const fp) {

    static struct jpeg_decompress_struct cinfo;
        /* The JPEG decompression parameters and pointers to
           working data (which is allocated as needed by the JPEG library).
        */
    struct jpeg_error_mgr jerr;
    JSAMPARRAY scanarray[3];
    int ci,cd,cp;
    JDIMENSION ncols[3];
    JDIMENSION nrows[3];
    jpeg_component_info *compptr;
    int buffer_height;
    int current_row[3];
    uint8_t **orig[3];
    int h_samp[3],v_samp[3];
    int max_h_samp,max_v_samp;
    int temp_h, temp_v;
    int temp;

    /* Allocate and initialize JPEG decompression object */
    cinfo.err = jpeg_std_error(&jerr);

    /*
    ** If we're reading from stdin we want to create the cinfo struct
    ** ONCE (during the first read).  This is because when reading jpeg
    ** from stdin we will not release the cinfo struct, because doing
    ** so would totally screw up the read buffer and make it impossible
    ** to read jpegs from stdin.
    ** Dave Scott (dhs), UofO, 7/19/95
    */
    {
        static int first_stdin;
        first_stdin = TRUE;  /* initial value */
        if ((fp != stdin) || first_stdin) {
            first_stdin = FALSE;
            /* Now we can initialize the JPEG decompression object. */
            jpeg_create_decompress(&cinfo);
            /* specify data source (eg, a file) */
            jpeg_stdio_src(&cinfo, fp);
        }
    }

    /* specify data source (eg, a file) */

    jpeg_stdio_src(&cinfo, fp);

    /* read file parameters with jpeg_read_header() */


    (void) jpeg_read_header(&cinfo, TRUE);
    /* We can ignore the return value from jpeg_read_header since
     *   (a) suspension is not possible with the stdio data source, and
     *   (b) we passed TRUE to reject a tables-only JPEG file as an error.
     */

    /* set parameters for decompression */
    cinfo.raw_data_out = TRUE;
    cinfo.out_color_space = JCS_YCbCr;

    /* calculate image output dimensions */
    jpeg_calc_output_dimensions(&cinfo);
    /* the above calculation will set these soon */
    /* for now we'll set them ourselves */

    /* tell mpeg_encode the size of the JPEG Image*/
    Fsize_Note(mf->id,(int)(cinfo.image_width),(int)(cinfo.image_height));

    /* Allocate memory for the raw YCbCr data to occupy*/
    Frame_AllocYCC(mf);      /*allocate space for mpeg frame*/

    /* copy pointers to array structure- this make the following
       code more compact
    */
    orig[0] = mf->orig_y;
    orig[1] = mf->orig_cb;
    orig[2] = mf->orig_cr;

    /* Note that we can use the info obtained from jpeg_read_header.
     */

    /* Start decompressor */

    jpeg_start_decompress(&cinfo);


    /* JSAMPLEs per row in output buffer  */
    /* collect component subsample values*/
    for (cp=0, compptr = cinfo.comp_info;
         cp < cinfo.num_components;
         cp++, compptr++) {
        h_samp[cp] = compptr->h_samp_factor;
        v_samp[cp] = compptr->v_samp_factor;
    }
    /* calculate max subsample values*/
    temp_h = (h_samp[0]<h_samp[1]) ? h_samp[1] : h_samp[0];
    max_h_samp = (temp_h<h_samp[2]) ? h_samp[2]:temp_h;
    temp_v = (v_samp[0]<v_samp[1]) ? v_samp[1] : v_samp[0];
    max_v_samp = (temp_v<v_samp[2]) ? v_samp[2]:temp_v;

    /* Make an 8-row-high sample array that will go away when done
       with image
    */
    buffer_height = cinfo.max_v_samp_factor * minDctVScaledSize(cinfo);

    for(cp=0,compptr = cinfo.comp_info;cp<cinfo.num_components;
        cp++,compptr++) {
        ncols[cp] = (JDIMENSION)((cinfo.image_width*compptr->h_samp_factor)/
                                 max_h_samp);

        nrows[cp] = (JDIMENSION)((buffer_height*compptr->v_samp_factor)/
                                 max_v_samp);

        scanarray[cp] = (*cinfo.mem->alloc_sarray)
            ((j_common_ptr) &cinfo, JPOOL_IMAGE, ncols[cp], nrows[cp]);
    }

    /*  while (scan lines remain to be read)
           jpeg_read_scanlines(...);
    */

    /* Here we use the library's state variable cinfo.output_scanline as the
     * loop counter, so that we don't have to keep track ourselves.
     */

    while (cinfo.output_scanline < cinfo.output_height) {

        (void) jpeg_read_raw_data(&cinfo, scanarray, buffer_height);

        /* alter subsample ratio's if necessary */
        if ((h_samp[0]==2) && (h_samp[1]==1) && (h_samp[2]==1) &&
            (v_samp[0]==2) && (v_samp[1]==1) && (v_samp[2]==1)) {
            /* we are 4:1:1 as expected by the encoder*/
        } else if((h_samp[0] == 2) && (h_samp[1] == 1) && (h_samp[2]==1) &&
                  (v_samp[0] == 1) && (v_samp[1] == 1) && (v_samp[2]==1)) {
            /* must subsample 2:1 vertically and adjust params*/
            for (ci = 1; ci < 3; ++ci) {
                for (cp = 0; cp < (buffer_height/2); cp = cp+1) {
                    for (cd = 0; cd < ncols[ci]; ++cd) {
                        temp = ((scanarray[ci][cp*2][cd] +
                                 scanarray[ci][(cp*2)+1][cd]) / 2);
                        scanarray[ci][cp][cd] = (JSAMPLE)(temp);
                    }
                }
            }
            /* only reset values the first time through*/
            if (cinfo.output_scanline == buffer_height) {
                nrows[1] = nrows[1]/2;
                nrows[2] = nrows[2]/2;
                max_v_samp = 2;
                v_samp[0] = 2;
            }
        } else
            pm_error("Invalid subsampling ratio");

        /* transfer data from jpeg buffer to MPEG frame */
        /* calculate the row we wish to output into */
        for (ci = 0, compptr = cinfo.comp_info;
             ci < cinfo.num_components;
             ++ci, ++compptr) {
            current_row[ci] =((cinfo.output_scanline - buffer_height)*
                              (v_samp[ci])/max_v_samp);

            jcopy_sample_rows(scanarray[ci],0,(JSAMPARRAY)(orig[ci]),
                              current_row[ci],nrows[ci],ncols[ci]);
        }
    }

    /* Step 7: Finish decompression */

    (void) jpeg_finish_decompress(&cinfo);
    /* We can ignore the return value since suspension is not possible
     * with the stdio data source.
     */

    /* Step 8: Release JPEG decompression object */

    /*
    ** DO NOT release the cinfo struct if we are reading from stdin, this
    ** is because the cinfo struct contains the read buffer, and the read
    ** buffer may (and almost always does) contain the end of one image and
    ** the beginning of another.  If we throw away the read buffer then
    ** we loose the beginning of the next image, and we're screwed.
    ** Dave Scott (dhs), UofO, 7/19/95
    */
    if (fp == stdin) {
        static int no_from_stdin = 0;
        ++no_from_stdin;
    } else {
        /* This is an important step since it will release a good deal
           of memory.
        */
        jpeg_destroy_decompress(&cinfo);
    }

    /* After finish_decompress, we can close the input file.  Here we
       postpone it until after no more JPEG errors are possible, so as
       to simplify the setjmp error logic above.  (Actually, I don't
       think that jpeg_destroy can do an error exit, but why assume
       anything...)
    */

    /* At this point you may want to check to see whether any corrupt-data
     * warnings occurred (test whether jerr.pub.num_warnings is nonzero).
     * If you prefer to treat corrupt data as a fatal error, override the
     * error handler's emit_message method to call error_exit on a warning.
     */
}



/*
 * SOME FINE POINTS:
 *
 * In the above loop, we ignored the return value of jpeg_read_scanlines,
 * which is the number of scanlines actually read.  We could get away with
 * this for the same reasons discussed in the compression example.  Actually
 * there is one perfectly normal situation in which jpeg_read_scanlines may
 * return fewer lines than you asked for: at the bottom of the image.  But the
 * loop above can't ask for more lines than there are in the image since it
 * reads only one line at a time.
 *
 * In some high-speed operating modes, some data copying can be saved by
 * making the buffer passed to jpeg_read_scanlines be cinfo.rec_outbuf_height
 * lines high (or a multiple thereof).  This will usually be 1, 2, or 4 lines.
 *
 * To decompress multiple images, you can repeat the whole sequence, or you
 * can keep the JPEG object around and just repeat steps 2-7.  This will
 * save a little bit of startup/shutdown time.
 *
 * As with compression, some operating modes may require temporary files.
 * On some systems you may need to set up a signal handler to ensure that
 * temporary files are deleted if the program is interrupted.
 *
 * Scanlines are returned in the same order as they appear in the JPEG file,
 * which is standardly top-to-bottom.  If you must have data supplied
 * bottom-to-top, you can use one of the virtual arrays provided by the
 * JPEG memory manager to invert the data.  See wrrle.c for an example.
 */


/*
 * Copyright (c) 1995 The Regents of the University of California.
 * All rights reserved.
 *
 * Permission to use, copy, modify, and distribute this software and its
 * documentation for any purpose, without fee, and without written agreement is
 * hereby granted, provided that the above copyright notice and the following
 * two paragraphs appear in all copies of this software.
 *
 * IN NO EVENT SHALL THE UNIVERSITY OF CALIFORNIA BE LIABLE TO ANY PARTY FOR
 * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, OR CONSEQUENTIAL DAMAGES ARISING OUT
 * OF THE USE OF THIS SOFTWARE AND ITS DOCUMENTATION, EVEN IF THE UNIVERSITY OF
 * CALIFORNIA HAS BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 *
 * THE UNIVERSITY OF CALIFORNIA SPECIFICALLY DISCLAIMS ANY WARRANTIES,
 * INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY
 * AND FITNESS FOR A PARTICULAR PURPOSE.  THE SOFTWARE PROVIDED HEREUNDER IS
 * ON AN "AS IS" BASIS, AND THE UNIVERSITY OF CALIFORNIA HAS NO OBLIGATION TO
 * PROVIDE MAINTENANCE, SUPPORT, UPDATES, ENHANCEMENTS, OR MODIFICATIONS.
 */