/*
// Copyright 2015 2018 Intel Corporation All Rights Reserved.
//
// The source code, information and material ("Material") contained herein is
// owned by Intel Corporation or its suppliers or licensors, and title
// to such Material remains with Intel Corporation or its suppliers or
// licensors. The Material contains proprietary information of Intel
// or its suppliers and licensors. The Material is protected by worldwide
// copyright laws and treaty provisions. No part of the Material may be used,
// copied, reproduced, modified, published, uploaded, posted, transmitted,
// distributed or disclosed in any way without Intel's prior express written
// permission. No license under any patent, copyright or other intellectual
// property rights in the Material is granted to or conferred upon you,
// either expressly, by implication, inducement, estoppel or otherwise.
// Any license under such intellectual property rights must be express and
// approved by Intel in writing.
//
// Unless otherwise agreed by Intel in writing,
// you may not remove or alter this notice or any other notice embedded in
// Materials by Intel or Intel's suppliers or licensors in any way.
*/
// A simple example of perspective warping of an image
// using Intel(R) Integrated Primitives (Intel(R) IPP) functions:
// ippiGetPerspectiveTransform
// ippiWarpPerspectiveSize
// ippiWarpPerspectiveNearestInit
// ippiWarpGetBufferSize
// ippiWarpPerspectiveNearest_32f_C1R
#include <stdio.h>
#include "ipp.h"
#define WIDTH 8 /* source image width */
#define HEIGHT 8 /* source image height */
/* Next two defines are created to simplify code reading and understanding */
#define EXIT_MAIN exitLine: /* Label for Exit */
#define check_sts(st) if((st) != ippStsNoErr) goto exitLine; /* Go to Exit if Intel(R) IPP function returned status different from ippStsNoErr */
static Ipp32f pSrc[WIDTH*HEIGHT] =/* Pointers to source images */
{
1.f, 1.f, 1.f, 1.f, 1.f, 1.f, 1.f, 1.f,
1.f, 1.f, 1.f, 5.f, 5.f, 1.f, 1.f, 1.f,
1.f, 1.f, 5.f, 5.f, 5.f, 5.f, 1.f, 1.f,
1.f, 5.f, 5.f, 5.f, 5.f, 5.f, 5.f, 1.f,
1.f, 5.f, 5.f, 5.f, 5.f, 5.f, 5.f, 1.f,
1.f, 1.f, 5.f, 5.f, 5.f, 5.f, 1.f, 1.f,
1.f, 1.f, 1.f, 5.f, 5.f, 1.f, 1.f, 1.f,
1.f, 1.f, 1.f, 1.f, 1.f, 1.f, 1.f, 1.f };
static const double quad[4][2] = /* Vertex coordinates of the quadrangle */
{
{ 0.0, 0.0 },
{ 8.0, 0.0 },
{ 6.0, 8.0 },
{ 0.0, 8.0 } };
/* Results of ippMalloc() are not validated because Intel(R) IPP functions perform bad arguments check and will return an appropriate status */
int main(void)
{
IppStatus status = ippStsNoErr;
IppiWarpSpec* pSpec = NULL; /* Pointer to the specification structure for the forward transform */
IppiWarpSpec* pSpecB = NULL; /* Pointer to the specification structure for the backward transform */
IppiSize srcSize = { WIDTH, HEIGHT }; /* Size of source ROI in pixels */
IppiRect srcRoi = { 0, 0, WIDTH, HEIGHT }; /* Region of interest in the source image */
IppiRect dstRoi = { 0, 0, WIDTH, HEIGHT }; /* Region of interest in the destination image */
IppiPoint dstOffset = { 0, 0 }; /* Offset of the destination image ROI with respect to the destination image origin */
IppiPoint dstbOffset = { 0, 0 }; /* Offset of the back destination image ROI with respect to the back destination image origin */
int srcStep = WIDTH * sizeof(Ipp32f), dstStep = srcStep, dstbStep = srcStep; /* Steps, in bytes, through the source/destination/back destination images */
Ipp32f pDst[WIDTH*HEIGHT], pDstB[WIDTH*HEIGHT]; /* Pointers to destination/back destination images */
IppiSize dstSize = {WIDTH, HEIGHT}; /* Size of the destination image */
IppiSize dstbSize = {WIDTH, HEIGHT}; /* Size of the back destination image */
double coeffs[3][3] = {0};
int specSize = 0, initSize = 0, bufSize = 0; /* Buffer sizes for the forward transform */
int specSizeB = 0, initSizeB = 0, bufSizeB = 0; /* Buffer sizes for the backward transform */
Ipp8u* pBuffer = NULL, *pBufferB = NULL; /* Work buffers */
IppiBorderType borderType = ippBorderTransp;
IppiWarpDirection direction = ippWarpForward; /* Forward transformation direction */
IppiWarpDirection directionB = ippWarpBackward; /* Backward transformation direction */
check_sts( status = ippiSet_32f_C1R(2, pDst, dstStep, dstSize) )
check_sts( status = ippiSet_32f_C1R(3, pDstB, dstbStep, dstbSize) )
/* Computes perspective transform matrix */
check_sts( status = ippiGetPerspectiveTransform(srcRoi, quad, (double(*)[3])coeffs) )
/* Spec and init buffer sizes for the forward transform*/
check_sts( status = ippiWarpPerspectiveGetSize(srcSize, srcRoi, dstSize, ipp32f, coeffs, ippNearest, direction, borderType, &specSize, &initSize) )
/* Spec and init buffer sizes for the backward transform*/
check_sts( status = ippiWarpPerspectiveGetSize(dstSize, dstRoi, dstbSize, ipp32f, coeffs, ippNearest, directionB, borderType, &specSizeB, &initSizeB) )
pSpec = (IppiWarpSpec*)ippsMalloc_8u(specSize);
pSpecB = (IppiWarpSpec*)ippsMalloc_8u(specSizeB);
/* Filter initialization for the forward transform */
check_sts( status = ippiWarpPerspectiveNearestInit(srcSize, srcRoi, dstSize, ipp32f, coeffs, direction, 1, borderType, NULL, 0, pSpec) )
/* Filter initialization for the backward transform */
check_sts( status = ippiWarpPerspectiveNearestInit(dstSize, dstRoi, dstbSize, ipp32f, coeffs, directionB, 1, borderType, NULL, 0, pSpecB) )
/* Get work buffer size for the forward transform */
check_sts( status = ippiWarpGetBufferSize(pSpec, dstSize, &bufSize) )
/* Get work buffer size for the backward transform */
check_sts( status = ippiWarpGetBufferSize(pSpecB, dstSize, &bufSizeB) )
pBuffer = ippsMalloc_8u(bufSize);
pBufferB = ippsMalloc_8u(bufSizeB);
/* Warp perspective processing for the forward transform */
check_sts( status = ippiWarpPerspectiveNearest_32f_C1R(pSrc, srcStep, pDst, dstStep, dstOffset, dstSize, pSpec, pBuffer) )
/* Warp perspective processing for the backward transform */
check_sts( status = ippiWarpPerspectiveNearest_32f_C1R(pDst, dstStep, pDstB, dstbStep, dstbOffset, dstbSize, pSpecB, pBufferB) )
EXIT_MAIN
ippsFree(pSpec);
ippsFree(pSpecB);
ippsFree(pBuffer);
ippsFree(pBufferB);
printf("Exit status %d (%s)\n", (int)status, ippGetStatusString(status));
return (int)status;
}
/*
pDst
1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
1.0 1.0 1.0 5.0 5.0 1.0 1.0 1.0
1.0 1.0 1.0 5.0 5.0 1.0 1.0 1.0
1.0 1.0 5.0 5.0 5.0 5.0 1.0 1.0
1.0 5.0 5.0 5.0 5.0 5.0 5.0 1.0
1.0 5.0 5.0 5.0 5.0 5.0 5.0 2.0
1.0 1.0 5.0 5.0 5.0 5.0 1.0 2.0
1.0 1.0 1.0 5.0 5.0 1.0 1.0 2.0
pDstB
1.0 1.0 1.0 1.0 1.0 1.0 1.0 3.0
1.0 1.0 1.0 5.0 5.0 1.0 1.0 3.0
1.0 1.0 5.0 5.0 5.0 5.0 1.0 3.0
1.0 5.0 5.0 5.0 5.0 5.0 5.0 1.0
1.0 5.0 5.0 5.0 5.0 5.0 5.0 2.0
1.0 1.0 5.0 5.0 5.0 5.0 1.0 1.0
3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0
3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0
*/