loadsave.cpp 30.2 KB
Newer Older
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
/*M///////////////////////////////////////////////////////////////////////////////////////
//
//  IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING.
//
//  By downloading, copying, installing or using the software you agree to this license.
//  If you do not agree to this license, do not download, install,
//  copy or use the software.
//
//
//                        Intel License Agreement
//                For Open Source Computer Vision Library
//
// Copyright (C) 2000, Intel Corporation, all rights reserved.
// Third party copyrights are property of their respective owners.
//
// Redistribution and use in source and binary forms, with or without modification,
// are permitted provided that the following conditions are met:
//
//   * Redistribution's of source code must retain the above copyright notice,
//     this list of conditions and the following disclaimer.
//
//   * Redistribution's in binary form must reproduce the above copyright notice,
//     this list of conditions and the following disclaimer in the documentation
//     and/or other materials provided with the distribution.
//
//   * The name of Intel Corporation may not be used to endorse or promote products
//     derived from this software without specific prior written permission.
//
// This software is provided by the copyright holders and contributors "as is" and
// any express or implied warranties, including, but not limited to, the implied
// warranties of merchantability and fitness for a particular purpose are disclaimed.
// In no event shall the Intel Corporation or contributors be liable for any direct,
// indirect, incidental, special, exemplary, or consequential damages
// (including, but not limited to, procurement of substitute goods or services;
// loss of use, data, or profits; or business interruption) however caused
// and on any theory of liability, whether in contract, strict liability,
// or tort (including negligence or otherwise) arising in any way out of
// the use of this software, even if advised of the possibility of such damage.
//
//M*/

//
//  Loading and saving IPL images.
//

#include "precomp.hpp"
#include "grfmts.hpp"
48 49
#include "utils.hpp"
#include "exif.hpp"
50 51
#undef min
#undef max
52
#include <iostream>
53
#include <fstream>
54 55 56 57

/****************************************************************************************\
*                                      Image Codecs                                      *
\****************************************************************************************/
58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78

namespace cv {

// TODO Add runtime configuration
#define CV_IO_MAX_IMAGE_PARAMS (50)
#define CV_IO_MAX_IMAGE_WIDTH (1<<20)
#define CV_IO_MAX_IMAGE_HEIGHT (1<<20)
#define CV_IO_MAX_IMAGE_PIXELS (1<<30) // 1 Gigapixel

static Size validateInputImageSize(const Size& size)
{
    CV_Assert(size.width > 0);
    CV_Assert(size.width <= CV_IO_MAX_IMAGE_WIDTH);
    CV_Assert(size.height > 0);
    CV_Assert(size.height <= CV_IO_MAX_IMAGE_HEIGHT);
    uint64 pixels = (uint64)size.width * (uint64)size.height;
    CV_Assert(pixels <= CV_IO_MAX_IMAGE_PIXELS);
    return size;
}


79 80 81 82 83 84 85 86 87 88 89 90 91
namespace {

class ByteStreamBuffer: public std::streambuf
{
public:
    ByteStreamBuffer(char* base, size_t length)
    {
        setg(base, base, base + length);
    }

protected:
    virtual pos_type seekoff( off_type offset,
                              std::ios_base::seekdir dir,
92
                              std::ios_base::openmode ) CV_OVERRIDE
93
    {
B
Brian Armstrong 已提交
94
        char* whence = eback();
95 96
        if (dir == std::ios_base::cur)
        {
B
Brian Armstrong 已提交
97
            whence = gptr();
98 99 100
        }
        else if (dir == std::ios_base::end)
        {
B
Brian Armstrong 已提交
101
            whence = egptr();
102
        }
B
Brian Armstrong 已提交
103
        char* to = whence + offset;
104 105

        // check limits
B
Brian Armstrong 已提交
106
        if (to >= eback() && to <= egptr())
107
        {
B
Brian Armstrong 已提交
108
            setg(eback(), to, egptr());
109 110 111 112 113 114 115 116 117
            return gptr() - eback();
        }

        return -1;
    }
};

}

118 119 120 121 122
/**
 * @struct ImageCodecInitializer
 *
 * Container which stores the registered codecs to be used by OpenCV
*/
123 124
struct ImageCodecInitializer
{
125 126 127
    /**
     * Default Constructor for the ImageCodeInitializer
    */
128 129
    ImageCodecInitializer()
    {
130
        /// BMP Support
R
Roman Donchenko 已提交
131 132
        decoders.push_back( makePtr<BmpDecoder>() );
        encoders.push_back( makePtr<BmpEncoder>() );
133

134
    #ifdef HAVE_IMGCODEC_HDR
135 136
        decoders.push_back( makePtr<HdrDecoder>() );
        encoders.push_back( makePtr<HdrEncoder>() );
137
    #endif
138
    #ifdef HAVE_JPEG
R
Roman Donchenko 已提交
139 140
        decoders.push_back( makePtr<JpegDecoder>() );
        encoders.push_back( makePtr<JpegEncoder>() );
141 142
    #endif
    #ifdef HAVE_WEBP
R
Roman Donchenko 已提交
143 144
        decoders.push_back( makePtr<WebPDecoder>() );
        encoders.push_back( makePtr<WebPEncoder>() );
145
    #endif
146
    #ifdef HAVE_IMGCODEC_SUNRASTER
R
Roman Donchenko 已提交
147 148
        decoders.push_back( makePtr<SunRasterDecoder>() );
        encoders.push_back( makePtr<SunRasterEncoder>() );
149 150
    #endif
    #ifdef HAVE_IMGCODEC_PXM
R
Roman Donchenko 已提交
151
        decoders.push_back( makePtr<PxMDecoder>() );
152 153 154 155
        encoders.push_back( makePtr<PxMEncoder>(PXM_TYPE_AUTO) );
        encoders.push_back( makePtr<PxMEncoder>(PXM_TYPE_PBM) );
        encoders.push_back( makePtr<PxMEncoder>(PXM_TYPE_PGM) );
        encoders.push_back( makePtr<PxMEncoder>(PXM_TYPE_PPM) );
156 157 158
        decoders.push_back( makePtr<PAMDecoder>() );
        encoders.push_back( makePtr<PAMEncoder>() );
    #endif
159
    #ifdef HAVE_TIFF
R
Roman Donchenko 已提交
160 161
        decoders.push_back( makePtr<TiffDecoder>() );
        encoders.push_back( makePtr<TiffEncoder>() );
162
    #endif
163
    #ifdef HAVE_PNG
R
Roman Donchenko 已提交
164 165
        decoders.push_back( makePtr<PngDecoder>() );
        encoders.push_back( makePtr<PngEncoder>() );
166
    #endif
167 168 169
    #ifdef HAVE_GDCM
        decoders.push_back( makePtr<DICOMDecoder>() );
    #endif
170
    #ifdef HAVE_JASPER
R
Roman Donchenko 已提交
171 172
        decoders.push_back( makePtr<Jpeg2KDecoder>() );
        encoders.push_back( makePtr<Jpeg2KEncoder>() );
173 174
    #endif
    #ifdef HAVE_OPENEXR
R
Roman Donchenko 已提交
175 176
        decoders.push_back( makePtr<ExrDecoder>() );
        encoders.push_back( makePtr<ExrEncoder>() );
177
    #endif
178 179 180 181 182

    #ifdef HAVE_GDAL
        /// Attach the GDAL Decoder
        decoders.push_back( makePtr<GdalDecoder>() );
    #endif/*HAVE_GDAL*/
183 184
    }

185 186
    std::vector<ImageDecoder> decoders;
    std::vector<ImageEncoder> encoders;
187 188 189
};

static ImageCodecInitializer codecs;
190

191 192 193 194 195 196 197 198 199
/**
 * Find the decoders
 *
 * @param[in] filename File to search
 *
 * @return Image decoder to parse image file.
*/
static ImageDecoder findDecoder( const String& filename ) {

200
    size_t i, maxlen = 0;
201 202

    /// iterate through list of registered codecs
203
    for( i = 0; i < codecs.decoders.size(); i++ )
204
    {
205
        size_t len = codecs.decoders[i]->signatureLength();
206 207 208
        maxlen = std::max(maxlen, len);
    }

209
    /// Open the file
210
    FILE* f= fopen( filename.c_str(), "rb" );
211 212

    /// in the event of a failure, return an empty image decoder
213 214
    if( !f )
        return ImageDecoder();
215 216

    // read the file signature
217
    String signature(maxlen, ' ');
218
    maxlen = fread( (void*)signature.c_str(), 1, maxlen, f );
219 220 221
    fclose(f);
    signature = signature.substr(0, maxlen);

222
    /// compare signature against all decoders
223
    for( i = 0; i < codecs.decoders.size(); i++ )
224
    {
225 226
        if( codecs.decoders[i]->checkSignature(signature) )
            return codecs.decoders[i]->newDecoder();
227 228
    }

229
    /// If no decoder was found, return base type
230 231 232
    return ImageDecoder();
}

233
static ImageDecoder findDecoder( const Mat& buf )
234 235 236 237 238 239
{
    size_t i, maxlen = 0;

    if( buf.rows*buf.cols < 1 || !buf.isContinuous() )
        return ImageDecoder();

240
    for( i = 0; i < codecs.decoders.size(); i++ )
241
    {
242
        size_t len = codecs.decoders[i]->signatureLength();
243 244 245
        maxlen = std::max(maxlen, len);
    }

246
    String signature(maxlen, ' ');
247 248
    size_t bufSize = buf.rows*buf.cols*buf.elemSize();
    maxlen = std::min(maxlen, bufSize);
249
    memcpy( (void*)signature.c_str(), buf.data, maxlen );
250

251
    for( i = 0; i < codecs.decoders.size(); i++ )
252
    {
253 254
        if( codecs.decoders[i]->checkSignature(signature) )
            return codecs.decoders[i]->newDecoder();
255 256 257 258 259
    }

    return ImageDecoder();
}

260
static ImageEncoder findEncoder( const String& _ext )
261 262 263 264 265 266 267 268
{
    if( _ext.size() <= 1 )
        return ImageEncoder();

    const char* ext = strrchr( _ext.c_str(), '.' );
    if( !ext )
        return ImageEncoder();
    int len = 0;
J
Julien Nabet 已提交
269
    for( ext++; len < 128 && isalnum(ext[len]); len++ )
270 271
        ;

272
    for( size_t i = 0; i < codecs.encoders.size(); i++ )
273
    {
274
        String description = codecs.encoders[i]->getDescription();
275 276 277 278 279 280 281 282
        const char* descr = strchr( description.c_str(), '(' );

        while( descr )
        {
            descr = strchr( descr + 1, '.' );
            if( !descr )
                break;
            int j = 0;
J
Julien Nabet 已提交
283
            for( descr++; j < len && isalnum(descr[j]) ; j++ )
284 285 286 287 288 289 290
            {
                int c1 = tolower(ext[j]);
                int c2 = tolower(descr[j]);
                if( c1 != c2 )
                    break;
            }
            if( j == len && !isalnum(descr[j]))
291
                return codecs.encoders[i]->newEncoder();
292 293 294 295 296 297 298
            descr += j;
        }
    }

    return ImageEncoder();
}

299

300 301
enum { LOAD_CVMAT=0, LOAD_IMAGE=1, LOAD_MAT=2 };

302
static void ExifTransform(int orientation, Mat& img)
A
Arkadiusz Raj 已提交
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
{
    switch( orientation )
    {
        case    IMAGE_ORIENTATION_TL: //0th row == visual top, 0th column == visual left-hand side
            //do nothing, the image already has proper orientation
            break;
        case    IMAGE_ORIENTATION_TR: //0th row == visual top, 0th column == visual right-hand side
            flip(img, img, 1); //flip horizontally
            break;
        case    IMAGE_ORIENTATION_BR: //0th row == visual bottom, 0th column == visual right-hand side
            flip(img, img, -1);//flip both horizontally and vertically
            break;
        case    IMAGE_ORIENTATION_BL: //0th row == visual bottom, 0th column == visual left-hand side
            flip(img, img, 0); //flip vertically
            break;
        case    IMAGE_ORIENTATION_LT: //0th row == visual left-hand side, 0th column == visual top
            transpose(img, img);
            break;
        case    IMAGE_ORIENTATION_RT: //0th row == visual right-hand side, 0th column == visual top
            transpose(img, img);
            flip(img, img, 1); //flip horizontally
            break;
        case    IMAGE_ORIENTATION_RB: //0th row == visual right-hand side, 0th column == visual bottom
            transpose(img, img);
            flip(img, img, -1); //flip both horizontally and vertically
            break;
        case    IMAGE_ORIENTATION_LB: //0th row == visual left-hand side, 0th column == visual bottom
            transpose(img, img);
            flip(img, img, 0); //flip vertically
            break;
        default:
            //by default the image read has normal (JPEG_ORIENTATION_TL) orientation
            break;
    }
}
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
static void ApplyExifOrientation(const String& filename, Mat& img)
{
    int orientation = IMAGE_ORIENTATION_TL;

    if (filename.size() > 0)
    {
        std::ifstream stream( filename.c_str(), std::ios_base::in | std::ios_base::binary );
        ExifReader reader( stream );
        if( reader.parse() )
        {
            ExifEntry_t entry = reader.getTag( ORIENTATION );
            if (entry.tag != INVALID_TAG)
            {
                orientation = entry.field_u16; //orientation is unsigned short, so check field_u16
            }
        }
        stream.close();
    }

    ExifTransform(orientation, img);
}

static void ApplyExifOrientation(const Mat& buf, Mat& img)
{
    int orientation = IMAGE_ORIENTATION_TL;

    if( buf.isContinuous() )
    {
        ByteStreamBuffer bsb( reinterpret_cast<char*>(buf.data), buf.total() * buf.elemSize() );
        std::istream stream( &bsb );
        ExifReader reader( stream );
        if( reader.parse() )
        {
            ExifEntry_t entry = reader.getTag( ORIENTATION );
            if (entry.tag != INVALID_TAG)
            {
                orientation = entry.field_u16; //orientation is unsigned short, so check field_u16
            }
        }
    }

    ExifTransform(orientation, img);
}

383 384 385 386 387 388 389 390 391 392
/**
 * Read an image into memory and return the information
 *
 * @param[in] filename File to load
 * @param[in] flags Flags
 * @param[in] hdrtype { LOAD_CVMAT=0,
 *                      LOAD_IMAGE=1,
 *                      LOAD_MAT=2
 *                    }
 * @param[in] mat Reference to C++ Mat object (If LOAD_MAT)
S
Suleyman TURKMEN 已提交
393
 * @param[in] scale_denom Scale value
394 395
 *
*/
396
static void*
397
imread_( const String& filename, int flags, int hdrtype, Mat* mat=0 )
398
{
399 400
    CV_Assert(mat || hdrtype != LOAD_MAT); // mat is required in LOAD_MAT case

401 402 403 404
    IplImage* image = 0;
    CvMat *matrix = 0;
    Mat temp, *data = &temp;

405 406 407 408
    /// Search for the relevant decoder to handle the imagery
    ImageDecoder decoder;

#ifdef HAVE_GDAL
M
Maksim Shabunin 已提交
409
    if(flags != IMREAD_UNCHANGED && (flags & IMREAD_LOAD_GDAL) == IMREAD_LOAD_GDAL ){
410 411 412
        decoder = GdalDecoder().newDecoder();
    }else{
#endif
413
        decoder = findDecoder( filename );
414 415 416 417 418 419
#ifdef HAVE_GDAL
    }
#endif

    /// if no decoder was found, return nothing.
    if( !decoder ){
420
        return 0;
421 422
    }

423 424 425
    int scale_denom = 1;
    if( flags > IMREAD_LOAD_GDAL )
    {
S
Suleyman TURKMEN 已提交
426
    if( flags & IMREAD_REDUCED_GRAYSCALE_2 )
427
        scale_denom = 2;
S
Suleyman TURKMEN 已提交
428
    else if( flags & IMREAD_REDUCED_GRAYSCALE_4 )
429
        scale_denom = 4;
S
Suleyman TURKMEN 已提交
430
    else if( flags & IMREAD_REDUCED_GRAYSCALE_8 )
431 432 433
        scale_denom = 8;
    }

S
Suleyman TURKMEN 已提交
434 435 436
    /// set the scale_denom in the driver
    decoder->setScale( scale_denom );

437
    /// set the filename in the driver
438
    decoder->setSource( filename );
439

M
Maksim Shabunin 已提交
440
    CV_TRY
441 442 443 444 445
    {
        // read the header to make sure it succeeds
        if( !decoder->readHeader() )
            return 0;
    }
M
Maksim Shabunin 已提交
446
    CV_CATCH (cv::Exception, e)
447 448
    {
        std::cerr << "imread_('" << filename << "'): can't read header: " << e.what() << std::endl << std::flush;
449
        return 0;
450
    }
M
Maksim Shabunin 已提交
451
    CV_CATCH_ALL
452 453 454 455 456
    {
        std::cerr << "imread_('" << filename << "'): can't read header: unknown exception" << std::endl << std::flush;
        return 0;
    }

457 458

    // established the required input image size
459
    Size size = validateInputImageSize(Size(decoder->width(), decoder->height()));
460

461
    // grab the decoded type
462
    int type = decoder->type();
463
    if( (flags & IMREAD_LOAD_GDAL) != IMREAD_LOAD_GDAL && flags != IMREAD_UNCHANGED )
464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479
    {
        if( (flags & CV_LOAD_IMAGE_ANYDEPTH) == 0 )
            type = CV_MAKETYPE(CV_8U, CV_MAT_CN(type));

        if( (flags & CV_LOAD_IMAGE_COLOR) != 0 ||
           ((flags & CV_LOAD_IMAGE_ANYCOLOR) != 0 && CV_MAT_CN(type) > 1) )
            type = CV_MAKETYPE(CV_MAT_DEPTH(type), 3);
        else
            type = CV_MAKETYPE(CV_MAT_DEPTH(type), 1);
    }

    if( hdrtype == LOAD_CVMAT || hdrtype == LOAD_MAT )
    {
        if( hdrtype == LOAD_CVMAT )
        {
            matrix = cvCreateMat( size.height, size.width, type );
480
            temp = cvarrToMat( matrix );
481 482 483 484 485 486 487 488 489 490
        }
        else
        {
            mat->create( size.height, size.width, type );
            data = mat;
        }
    }
    else
    {
        image = cvCreateImage( size, cvIplDepth(type), CV_MAT_CN(type) );
491
        temp = cvarrToMat( image );
492 493
    }

494
    // read the image data
495
    bool success = false;
M
Maksim Shabunin 已提交
496
    CV_TRY
497 498 499 500
    {
        if (decoder->readData(*data))
            success = true;
    }
M
Maksim Shabunin 已提交
501
    CV_CATCH (cv::Exception, e)
502 503 504
    {
        std::cerr << "imread_('" << filename << "'): can't read data: " << e.what() << std::endl << std::flush;
    }
M
Maksim Shabunin 已提交
505
    CV_CATCH_ALL
506 507 508 509
    {
        std::cerr << "imread_('" << filename << "'): can't read data: unknown exception" << std::endl << std::flush;
    }
    if (!success)
510 511 512 513 514 515 516 517
    {
        cvReleaseImage( &image );
        cvReleaseMat( &matrix );
        if( mat )
            mat->release();
        return 0;
    }

518
    if( decoder->setScale( scale_denom ) > 1 ) // if decoder is JpegDecoder then decoder->setScale always returns 1
S
Suleyman TURKMEN 已提交
519
    {
520
        resize( *mat, *mat, Size( size.width / scale_denom, size.height / scale_denom ), 0, 0, INTER_LINEAR_EXACT);
S
Suleyman TURKMEN 已提交
521 522
    }

523 524 525 526
    return hdrtype == LOAD_CVMAT ? (void*)matrix :
        hdrtype == LOAD_IMAGE ? (void*)image : (void*)mat;
}

527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542

/**
* Read an image into memory and return the information
*
* @param[in] filename File to load
* @param[in] flags Flags
* @param[in] mats Reference to C++ vector<Mat> object to hold the images
*
*/
static bool
imreadmulti_(const String& filename, int flags, std::vector<Mat>& mats)
{
    /// Search for the relevant decoder to handle the imagery
    ImageDecoder decoder;

#ifdef HAVE_GDAL
M
Maksim Shabunin 已提交
543
    if (flags != IMREAD_UNCHANGED && (flags & IMREAD_LOAD_GDAL) == IMREAD_LOAD_GDAL){
544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561
        decoder = GdalDecoder().newDecoder();
    }
    else{
#endif
        decoder = findDecoder(filename);
#ifdef HAVE_GDAL
    }
#endif

    /// if no decoder was found, return nothing.
    if (!decoder){
        return 0;
    }

    /// set the filename in the driver
    decoder->setSource(filename);

    // read the header to make sure it succeeds
M
Maksim Shabunin 已提交
562
    CV_TRY
563 564 565 566 567
    {
        // read the header to make sure it succeeds
        if( !decoder->readHeader() )
            return 0;
    }
M
Maksim Shabunin 已提交
568
    CV_CATCH (cv::Exception, e)
569 570
    {
        std::cerr << "imreadmulti_('" << filename << "'): can't read header: " << e.what() << std::endl << std::flush;
571
        return 0;
572
    }
M
Maksim Shabunin 已提交
573
    CV_CATCH_ALL
574 575 576 577
    {
        std::cerr << "imreadmulti_('" << filename << "'): can't read header: unknown exception" << std::endl << std::flush;
        return 0;
    }
578 579 580 581 582

    for (;;)
    {
        // grab the decoded type
        int type = decoder->type();
583
        if( (flags & IMREAD_LOAD_GDAL) != IMREAD_LOAD_GDAL && flags != IMREAD_UNCHANGED )
584 585 586 587 588 589 590 591 592 593 594
        {
            if ((flags & CV_LOAD_IMAGE_ANYDEPTH) == 0)
                type = CV_MAKETYPE(CV_8U, CV_MAT_CN(type));

            if ((flags & CV_LOAD_IMAGE_COLOR) != 0 ||
                ((flags & CV_LOAD_IMAGE_ANYCOLOR) != 0 && CV_MAT_CN(type) > 1))
                type = CV_MAKETYPE(CV_MAT_DEPTH(type), 3);
            else
                type = CV_MAKETYPE(CV_MAT_DEPTH(type), 1);
        }

595 596 597
        // established the required input image size
        Size size = validateInputImageSize(Size(decoder->width(), decoder->height()));

598
        // read the image data
599 600
        Mat mat(size.height, size.width, type);
        bool success = false;
M
Maksim Shabunin 已提交
601
        CV_TRY
602
        {
603 604 605
            if (decoder->readData(mat))
                success = true;
        }
M
Maksim Shabunin 已提交
606
        CV_CATCH (cv::Exception, e)
607 608 609
        {
            std::cerr << "imreadmulti_('" << filename << "'): can't read data: " << e.what() << std::endl << std::flush;
        }
M
Maksim Shabunin 已提交
610
        CV_CATCH_ALL
611 612 613 614
        {
            std::cerr << "imreadmulti_('" << filename << "'): can't read data: unknown exception" << std::endl << std::flush;
        }
        if (!success)
615
            break;
616 617 618 619 620

        // optionally rotate the data if EXIF' orientation flag says so
        if( (flags & IMREAD_IGNORE_ORIENTATION) == 0 && flags != IMREAD_UNCHANGED )
        {
            ApplyExifOrientation(filename, mat);
621 622 623 624 625 626 627 628 629 630 631 632
        }

        mats.push_back(mat);
        if (!decoder->nextPage())
        {
            break;
        }
    }

    return !mats.empty();
}

633 634 635 636 637 638 639 640
/**
 * Read an image
 *
 *  This function merely calls the actual implementation above and returns itself.
 *
 * @param[in] filename File to load
 * @param[in] flags Flags you wish to set.
*/
641
Mat imread( const String& filename, int flags )
642
{
643 644
    CV_TRACE_FUNCTION();

645
    /// create the basic container
646
    Mat img;
647 648

    /// load the data
649
    imread_( filename, flags, LOAD_MAT, &img );
650

651
    /// optionally rotate the data if EXIF' orientation flag says so
652
    if( !img.empty() && (flags & IMREAD_IGNORE_ORIENTATION) == 0 && flags != IMREAD_UNCHANGED )
653
    {
A
Arkadiusz Raj 已提交
654
        ApplyExifOrientation(filename, img);
655 656
    }

657
    /// return a reference to the data
658 659 660
    return img;
}

661 662 663 664 665 666 667 668 669 670 671 672
/**
* Read a multi-page image
*
*  This function merely calls the actual implementation above and returns itself.
*
* @param[in] filename File to load
* @param[in] mats Reference to C++ vector<Mat> object to hold the images
* @param[in] flags Flags you wish to set.
*
*/
bool imreadmulti(const String& filename, std::vector<Mat>& mats, int flags)
{
673 674
    CV_TRACE_FUNCTION();

675 676 677
    return imreadmulti_(filename, flags, mats);
}

678
static bool imwrite_( const String& filename, const std::vector<Mat>& img_vec,
679
                      const std::vector<int>& params, bool flipv )
680
{
681 682
    bool isMultiImg = img_vec.size() > 1;
    std::vector<Mat> write_vec;
683 684

    ImageEncoder encoder = findEncoder( filename );
R
Roman Donchenko 已提交
685
    if( !encoder )
686 687
        CV_Error( CV_StsError, "could not find a writer for the specified extension" );

688
    for (size_t page = 0; page < img_vec.size(); page++)
689
    {
690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707
        Mat image = img_vec[page];
        CV_Assert( image.channels() == 1 || image.channels() == 3 || image.channels() == 4 );

        Mat temp;
        if( !encoder->isFormatSupported(image.depth()) )
        {
            CV_Assert( encoder->isFormatSupported(CV_8U) );
            image.convertTo( temp, CV_8U );
            image = temp;
        }

        if( flipv )
        {
            flip(image, temp, 0);
            image = temp;
        }

        write_vec.push_back(image);
708 709 710
    }

    encoder->setDestination( filename );
711
    CV_Assert(params.size() <= CV_IO_MAX_IMAGE_PARAMS*2);
712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727
    bool code = false;
    try
    {
        if (!isMultiImg)
            code = encoder->write( write_vec[0], params );
        else
            code = encoder->writemulti( write_vec, params ); //to be implemented
    }
    catch (const cv::Exception& e)
    {
        std::cerr << "imwrite_('" << filename << "'): can't write data: " << e.what() << std::endl << std::flush;
    }
    catch (...)
    {
        std::cerr << "imwrite_('" << filename << "'): can't write data: unknown exception" << std::endl << std::flush;
    }
728 729 730 731 732

    //    CV_Assert( code );
    return code;
}

733
bool imwrite( const String& filename, InputArray _img,
734
              const std::vector<int>& params )
735
{
736
    CV_TRACE_FUNCTION();
737
    std::vector<Mat> img_vec;
738
    if (_img.isMatVector() || _img.isUMatVector())
739
        _img.getMatVector(img_vec);
740
    else
741
        img_vec.push_back(_img.getMat());
742

743
    CV_Assert(!img_vec.empty());
744
    return imwrite_(filename, img_vec, params, false);
745 746 747 748 749
}

static void*
imdecode_( const Mat& buf, int flags, int hdrtype, Mat* mat=0 )
{
750
    CV_Assert(!buf.empty() && buf.isContinuous());
751 752 753
    IplImage* image = 0;
    CvMat *matrix = 0;
    Mat temp, *data = &temp;
754
    String filename;
755 756

    ImageDecoder decoder = findDecoder(buf);
R
Roman Donchenko 已提交
757
    if( !decoder )
758 759 760 761
        return 0;

    if( !decoder->setSource(buf) )
    {
R
Roy Reapor 已提交
762
        filename = tempfile();
763
        FILE* f = fopen( filename.c_str(), "wb" );
764 765 766
        if( !f )
            return 0;
        size_t bufSize = buf.cols*buf.rows*buf.elemSize();
767 768 769 770 771 772 773 774 775
        if( fwrite( buf.ptr(), 1, bufSize, f ) != bufSize )
        {
            fclose( f );
            CV_Error( CV_StsError, "failed to write image data to temporary file" );
        }
        if( fclose(f) != 0 )
        {
            CV_Error( CV_StsError, "failed to write image data to temporary file" );
        }
776 777 778
        decoder->setSource(filename);
    }

779
    bool success = false;
M
Maksim Shabunin 已提交
780
    CV_TRY
781 782 783 784
    {
        if (decoder->readHeader())
            success = true;
    }
M
Maksim Shabunin 已提交
785
    CV_CATCH (cv::Exception, e)
786 787 788
    {
        std::cerr << "imdecode_('" << filename << "'): can't read header: " << e.what() << std::endl << std::flush;
    }
M
Maksim Shabunin 已提交
789
    CV_CATCH_ALL
790 791 792 793
    {
        std::cerr << "imdecode_('" << filename << "'): can't read header: unknown exception" << std::endl << std::flush;
    }
    if (!success)
794
    {
795
        decoder.release();
796
        if (!filename.empty())
797
        {
798
            if (0 != remove(filename.c_str()))
799
            {
800
                std::cerr << "unable to remove temporary file:" << filename << std::endl << std::flush;
801 802
            }
        }
803 804 805
        return 0;
    }

806 807
    // established the required input image size
    Size size = validateInputImageSize(Size(decoder->width(), decoder->height()));
808 809

    int type = decoder->type();
810
    if( (flags & IMREAD_LOAD_GDAL) != IMREAD_LOAD_GDAL && flags != IMREAD_UNCHANGED )
811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840
    {
        if( (flags & CV_LOAD_IMAGE_ANYDEPTH) == 0 )
            type = CV_MAKETYPE(CV_8U, CV_MAT_CN(type));

        if( (flags & CV_LOAD_IMAGE_COLOR) != 0 ||
           ((flags & CV_LOAD_IMAGE_ANYCOLOR) != 0 && CV_MAT_CN(type) > 1) )
            type = CV_MAKETYPE(CV_MAT_DEPTH(type), 3);
        else
            type = CV_MAKETYPE(CV_MAT_DEPTH(type), 1);
    }

    if( hdrtype == LOAD_CVMAT || hdrtype == LOAD_MAT )
    {
        if( hdrtype == LOAD_CVMAT )
        {
            matrix = cvCreateMat( size.height, size.width, type );
            temp = cvarrToMat(matrix);
        }
        else
        {
            mat->create( size.height, size.width, type );
            data = mat;
        }
    }
    else
    {
        image = cvCreateImage( size, cvIplDepth(type), CV_MAT_CN(type) );
        temp = cvarrToMat(image);
    }

841
    success = false;
M
Maksim Shabunin 已提交
842
    CV_TRY
843 844 845 846
    {
        if (decoder->readData(*data))
            success = true;
    }
M
Maksim Shabunin 已提交
847
    CV_CATCH (cv::Exception, e)
848 849 850
    {
        std::cerr << "imdecode_('" << filename << "'): can't read data: " << e.what() << std::endl << std::flush;
    }
M
Maksim Shabunin 已提交
851
    CV_CATCH_ALL
852 853 854
    {
        std::cerr << "imdecode_('" << filename << "'): can't read data: unknown exception" << std::endl << std::flush;
    }
855
    decoder.release();
856
    if (!filename.empty())
857
    {
858
        if (0 != remove(filename.c_str()))
859
        {
860
            std::cerr << "unable to remove temporary file:" << filename << std::endl << std::flush;
861 862
        }
    }
863

864
    if (!success)
865 866 867 868 869 870 871 872 873 874 875 876 877
    {
        cvReleaseImage( &image );
        cvReleaseMat( &matrix );
        if( mat )
            mat->release();
        return 0;
    }

    return hdrtype == LOAD_CVMAT ? (void*)matrix :
        hdrtype == LOAD_IMAGE ? (void*)image : (void*)mat;
}


878
Mat imdecode( InputArray _buf, int flags )
879
{
880 881
    CV_TRACE_FUNCTION();

882
    Mat buf = _buf.getMat(), img;
883
    imdecode_( buf, flags, LOAD_MAT, &img );
884 885 886 887 888 889 890

    /// optionally rotate the data if EXIF' orientation flag says so
    if( !img.empty() && (flags & IMREAD_IGNORE_ORIENTATION) == 0 && flags != IMREAD_UNCHANGED )
    {
        ApplyExifOrientation(buf, img);
    }

891 892
    return img;
}
893

894 895
Mat imdecode( InputArray _buf, int flags, Mat* dst )
{
896 897
    CV_TRACE_FUNCTION();

898 899 900
    Mat buf = _buf.getMat(), img;
    dst = dst ? dst : &img;
    imdecode_( buf, flags, LOAD_MAT, dst );
901 902 903 904 905 906 907

    /// optionally rotate the data if EXIF' orientation flag says so
    if( !dst->empty() && (flags & IMREAD_IGNORE_ORIENTATION) == 0 && flags != IMREAD_UNCHANGED )
    {
        ApplyExifOrientation(buf, *dst);
    }

908 909
    return *dst;
}
910

911
bool imencode( const String& ext, InputArray _image,
912
               std::vector<uchar>& buf, const std::vector<int>& params )
913
{
914 915
    CV_TRACE_FUNCTION();

916
    Mat image = _image.getMat();
917 918 919 920 921

    int channels = image.channels();
    CV_Assert( channels == 1 || channels == 3 || channels == 4 );

    ImageEncoder encoder = findEncoder( ext );
R
Roman Donchenko 已提交
922
    if( !encoder )
923 924 925 926 927
        CV_Error( CV_StsError, "could not find encoder for the specified extension" );

    if( !encoder->isFormatSupported(image.depth()) )
    {
        CV_Assert( encoder->isFormatSupported(CV_8U) );
928
        Mat temp;
929
        image.convertTo(temp, CV_8U);
930
        image = temp;
931 932 933 934 935 936
    }

    bool code;
    if( encoder->setDestination(buf) )
    {
        code = encoder->write(image, params);
937
        encoder->throwOnEror();
938 939 940 941
        CV_Assert( code );
    }
    else
    {
942
        String filename = tempfile();
943 944
        code = encoder->setDestination(filename);
        CV_Assert( code );
945

946
        code = encoder->write(image, params);
947
        encoder->throwOnEror();
948
        CV_Assert( code );
949

950
        FILE* f = fopen( filename.c_str(), "rb" );
951 952 953 954 955 956 957
        CV_Assert(f != 0);
        fseek( f, 0, SEEK_END );
        long pos = ftell(f);
        buf.resize((size_t)pos);
        fseek( f, 0, SEEK_SET );
        buf.resize(fread( &buf[0], 1, buf.size(), f ));
        fclose(f);
958
        remove(filename.c_str());
959 960 961 962 963 964 965
    }
    return code;
}

}

/****************************************************************************************\
966
*                         Imgcodecs loading & saving function implementation            *
967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000
\****************************************************************************************/

CV_IMPL int
cvHaveImageReader( const char* filename )
{
    cv::ImageDecoder decoder = cv::findDecoder(filename);
    return !decoder.empty();
}

CV_IMPL int cvHaveImageWriter( const char* filename )
{
    cv::ImageEncoder encoder = cv::findEncoder(filename);
    return !encoder.empty();
}

CV_IMPL IplImage*
cvLoadImage( const char* filename, int iscolor )
{
    return (IplImage*)cv::imread_(filename, iscolor, cv::LOAD_IMAGE );
}

CV_IMPL CvMat*
cvLoadImageM( const char* filename, int iscolor )
{
    return (CvMat*)cv::imread_( filename, iscolor, cv::LOAD_CVMAT );
}

CV_IMPL int
cvSaveImage( const char* filename, const CvArr* arr, const int* _params )
{
    int i = 0;
    if( _params )
    {
        for( ; _params[i] > 0; i += 2 )
1001
            CV_Assert(i < CV_IO_MAX_IMAGE_PARAMS*2); // Limit number of params for security reasons
1002 1003
    }
    return cv::imwrite_(filename, cv::cvarrToMat(arr),
1004
        i > 0 ? std::vector<int>(_params, _params+i) : std::vector<int>(),
1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031
        CV_IS_IMAGE(arr) && ((const IplImage*)arr)->origin == IPL_ORIGIN_BL );
}

/* decode image stored in the buffer */
CV_IMPL IplImage*
cvDecodeImage( const CvMat* _buf, int iscolor )
{
    CV_Assert( _buf && CV_IS_MAT_CONT(_buf->type) );
    cv::Mat buf(1, _buf->rows*_buf->cols*CV_ELEM_SIZE(_buf->type), CV_8U, _buf->data.ptr);
    return (IplImage*)cv::imdecode_(buf, iscolor, cv::LOAD_IMAGE );
}

CV_IMPL CvMat*
cvDecodeImageM( const CvMat* _buf, int iscolor )
{
    CV_Assert( _buf && CV_IS_MAT_CONT(_buf->type) );
    cv::Mat buf(1, _buf->rows*_buf->cols*CV_ELEM_SIZE(_buf->type), CV_8U, _buf->data.ptr);
    return (CvMat*)cv::imdecode_(buf, iscolor, cv::LOAD_CVMAT );
}

CV_IMPL CvMat*
cvEncodeImage( const char* ext, const CvArr* arr, const int* _params )
{
    int i = 0;
    if( _params )
    {
        for( ; _params[i] > 0; i += 2 )
1032
            CV_Assert(i < CV_IO_MAX_IMAGE_PARAMS*2); // Limit number of params for security reasons
1033 1034 1035 1036 1037 1038 1039 1040
    }
    cv::Mat img = cv::cvarrToMat(arr);
    if( CV_IS_IMAGE(arr) && ((const IplImage*)arr)->origin == IPL_ORIGIN_BL )
    {
        cv::Mat temp;
        cv::flip(img, temp, 0);
        img = temp;
    }
1041
    std::vector<uchar> buf;
1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053

    bool code = cv::imencode(ext, img, buf,
        i > 0 ? std::vector<int>(_params, _params+i) : std::vector<int>() );
    if( !code )
        return 0;
    CvMat* _buf = cvCreateMat(1, (int)buf.size(), CV_8U);
    memcpy( _buf->data.ptr, &buf[0], buf.size() );

    return _buf;
}

/* End of file. */