Вырезка части изображения по двум точкам - C++

bool TiffFile::Read(const char *fileName, int x1, int y1, int x2, int y2)
{
    tif=TIFFOpen(fileName,"r");
    if(tif)
     {
        int compression=0;
        int planar=0;
        TIFFGetField(tif,TIFFTAG_COMPRESSION,&compression);
        TIFFGetField(tif,TIFFTAG_PLANARCONFIG,&planar);
        if((compression==1)&&(planar==PLANARCONFIG_CONTIG))
        {
            if(mat!=0)
            {
                delete[] mat;
                height=0;
                width=0;
                depth=0;
                channels=0;
            }
            TIFFGetField(tif,TIFFTAG_IMAGELENGTH,&height);
            TIFFGetField(tif,TIFFTAG_IMAGEWIDTH,&width);
            TIFFGetField(tif,TIFFTAG_BITSPERSAMPLE,&depth);
            TIFFGetField(tif,TIFFTAG_SAMPLESPERPIXEL,&channels);
            if((x1<width&&x1>=0)&&(x2<width&&x2>=0)&&(y1<height&&y1>=0)&&(y2<height&&y2>=0)&&(x2>x1)&&(y2>y1))
            {
                mat=new char[(x2-x1)*(y2-y1)*(depth/8)*channels];
                //cout<<"size of "<<(x2-x1)*(y2-y1)*(depth/8)*channels<<endl;
                char *line=new char[width*(depth/8)*channels];
                for(int i=y1;i<y2;i++)
                {
                    TIFFReadScanline(tif,line,i);
                    for(int j=x1;j<x2;j++)
                    {
                        mat[(i-y1)*(x2-x1)+j-x1]=line[j];
                    }
                }
                width=x2-x1;
                height=y2-y1;
                TIFFClose(tif);
                return true;
            }
            else
            {
                cout<<"Point not in the image area."<<endl;
                return false;
            }
        }
        else
        {
            cout<<"Sorry.Given type of image is not supported.Refer to Lukyanau Maxim"<<endl;
            return false;
        }
     }
     else
    {
         cout<<"Can not open image"<<endl;
         return false;
    }
}