0N/A/*
2362N/A * Copyright (c) 1997, 2003, Oracle and/or its affiliates. All rights reserved.
0N/A * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER.
0N/A *
0N/A * This code is free software; you can redistribute it and/or modify it
0N/A * under the terms of the GNU General Public License version 2 only, as
2362N/A * published by the Free Software Foundation. Oracle designates this
0N/A * particular file as subject to the "Classpath" exception as provided
2362N/A * by Oracle in the LICENSE file that accompanied this code.
0N/A *
0N/A * This code is distributed in the hope that it will be useful, but WITHOUT
0N/A * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
0N/A * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
0N/A * version 2 for more details (a copy is included in the LICENSE file that
0N/A * accompanied this code).
0N/A *
0N/A * You should have received a copy of the GNU General Public License version
0N/A * 2 along with this work; if not, write to the Free Software Foundation,
0N/A * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.
0N/A *
2362N/A * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA
2362N/A * or visit www.oracle.com if you need additional information or have any
2362N/A * questions.
0N/A */
0N/A
0N/A
0N/A/*
0N/A * FUNCTION
0N/A * mlib_ImageAffine_u8_1ch_nn
0N/A * mlib_ImageAffine_u8_2ch_nn
0N/A * mlib_ImageAffine_u8_3ch_nn
0N/A * mlib_ImageAffine_u8_4ch_nn
0N/A * mlib_ImageAffine_s16_1ch_nn
0N/A * mlib_ImageAffine_s16_2ch_nn
0N/A * mlib_ImageAffine_s16_3ch_nn
0N/A * mlib_ImageAffine_s16_4ch_nn
0N/A * - image affine transformation with Nearest Neighbor filtering
0N/A * SYNOPSIS
0N/A * mlib_status mlib_ImageAffine_[u8|s16]_?ch_nn(mlib_s32 *leftEdges,
0N/A * mlib_s32 *rightEdges,
0N/A * mlib_s32 *xStarts,
0N/A * mlib_s32 *yStarts,
0N/A * mlib_s32 *sides,
0N/A * mlib_u8 *dstData,
0N/A * mlib_u8 **lineAddr,
0N/A * mlib_s32 dstYStride,
0N/A * mlib_s32 is_affine)
0N/A *
0N/A * ARGUMENTS
0N/A * leftEdges array[dstHeight] of xLeft coordinates
0N/A * RightEdges array[dstHeight] of xRight coordinates
0N/A * xStarts array[dstHeight] of xStart * 65536 coordinates
0N/A * yStarts array[dstHeight] of yStart * 65536 coordinates
0N/A * sides output array[4]. sides[0] is yStart, sides[1] is yFinish,
0N/A * sides[2] is dx * 65536, sides[3] is dy * 65536
0N/A * dstData pointer to the first pixel on (yStart - 1) line
0N/A * lineAddr array[srcHeight] of pointers to the first pixel on
0N/A * the corresponding lines
0N/A * dstYStride stride of destination image
0N/A * is_affine indicator (Affine - GridWarp)
0N/A *
0N/A * DESCRIPTION
0N/A * The functions step along the lines from xLeft to xRight and get the
0N/A * nearest pixel values as being with the following coordinates
0N/A * ((xStart - (i - xLeft) * dx) >> 16, (yStart - (i - xLeft) * dy) >> 16)
0N/A *
0N/A */
0N/A
0N/A#include "mlib_ImageAffine.h"
0N/A
0N/A/***************************************************************/
0N/A#undef DTYPE
0N/A#define DTYPE mlib_u8
0N/A
0N/Amlib_status mlib_ImageAffine_u8_1ch_nn(mlib_affine_param *param)
0N/A{
0N/A DECLAREVAR_NN();
0N/A DTYPE *dstLineEnd;
0N/A
0N/A for (j = yStart; j <= yFinish; j++) {
0N/A DTYPE pix0;
0N/A
0N/A CLIP(1);
0N/A dstLineEnd = (DTYPE *) dstData + xRight;
0N/A
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/A for (; dstPixelPtr <= dstLineEnd; dstPixelPtr++) {
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc);
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A pix0 = srcPixelPtr[xSrc];
0N/A dstPixelPtr[0] = pix0;
0N/A }
0N/A }
0N/A
0N/A return MLIB_SUCCESS;
0N/A}
0N/A
0N/A/***************************************************************/
0N/Amlib_status mlib_ImageAffine_u8_2ch_nn(mlib_affine_param *param)
0N/A{
0N/A DECLAREVAR_NN();
0N/A DTYPE *dstLineEnd;
0N/A
0N/A for (j = yStart; j <= yFinish; j++) {
0N/A DTYPE pix0, pix1;
0N/A
0N/A CLIP(2);
0N/A dstLineEnd = (DTYPE *) dstData + 2 * xRight;
0N/A
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 2 * xSrc;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/A for (; dstPixelPtr < dstLineEnd; dstPixelPtr += 2) {
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 2 * xSrc;
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A }
0N/A
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A }
0N/A
0N/A return MLIB_SUCCESS;
0N/A}
0N/A
0N/A/***************************************************************/
0N/Amlib_status mlib_ImageAffine_u8_3ch_nn(mlib_affine_param *param)
0N/A{
0N/A DECLAREVAR_NN();
0N/A DTYPE *dstLineEnd;
0N/A
0N/A for (j = yStart; j <= yFinish; j++) {
0N/A DTYPE pix0, pix1, pix2;
0N/A
0N/A CLIP(3);
0N/A dstLineEnd = (DTYPE *) dstData + 3 * xRight;
0N/A
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 3 * xSrc;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A pix2 = srcPixelPtr[2];
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/A for (; dstPixelPtr < dstLineEnd; dstPixelPtr += 3) {
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 3 * xSrc;
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A dstPixelPtr[2] = pix2;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A pix2 = srcPixelPtr[2];
0N/A }
0N/A
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A dstPixelPtr[2] = pix2;
0N/A }
0N/A
0N/A return MLIB_SUCCESS;
0N/A}
0N/A
0N/A/***************************************************************/
0N/Amlib_status mlib_ImageAffine_u8_4ch_nn(mlib_affine_param *param)
0N/A{
0N/A DECLAREVAR_NN();
0N/A DTYPE *dstLineEnd;
0N/A
0N/A for (j = yStart; j <= yFinish; j++) {
0N/A DTYPE pix0, pix1, pix2, pix3;
0N/A CLIP(4);
0N/A dstLineEnd = (DTYPE *) dstData + 4 * xRight;
0N/A
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 4 * xSrc;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A pix2 = srcPixelPtr[2];
0N/A pix3 = srcPixelPtr[3];
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A for (; dstPixelPtr < dstLineEnd; dstPixelPtr += 4) {
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 4 * xSrc;
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A dstPixelPtr[2] = pix2;
0N/A dstPixelPtr[3] = pix3;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A pix2 = srcPixelPtr[2];
0N/A pix3 = srcPixelPtr[3];
0N/A }
0N/A
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A dstPixelPtr[2] = pix2;
0N/A dstPixelPtr[3] = pix3;
0N/A }
0N/A
0N/A return MLIB_SUCCESS;
0N/A}
0N/A
0N/A/***************************************************************/
0N/A#undef DTYPE
0N/A#define DTYPE mlib_u16
0N/A
0N/Amlib_status mlib_ImageAffine_s16_1ch_nn(mlib_affine_param *param)
0N/A{
0N/A DECLAREVAR_NN();
0N/A DTYPE *dstLineEnd;
0N/A
0N/A for (j = yStart; j <= yFinish; j++) {
0N/A mlib_s32 pix0;
0N/A
0N/A CLIP(1);
0N/A dstLineEnd = (DTYPE *) dstData + xRight;
0N/A
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc);
0N/A pix0 = srcPixelPtr[xSrc];
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A for (; dstPixelPtr < dstLineEnd; dstPixelPtr++) {
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc);
0N/A dstPixelPtr[0] = pix0;
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A pix0 = srcPixelPtr[xSrc];
0N/A }
0N/A
0N/A dstPixelPtr[0] = pix0;
0N/A }
0N/A
0N/A return MLIB_SUCCESS;
0N/A}
0N/A
0N/A/***************************************************************/
0N/Amlib_status mlib_ImageAffine_s16_2ch_nn(mlib_affine_param *param)
0N/A{
0N/A DECLAREVAR_NN();
0N/A DTYPE *dstLineEnd;
0N/A
0N/A for (j = yStart; j <= yFinish; j++) {
0N/A mlib_s32 pix0, pix1;
0N/A
0N/A CLIP(2);
0N/A dstLineEnd = (DTYPE *) dstData + 2 * xRight;
0N/A
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 2 * xSrc;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/A for (; dstPixelPtr < dstLineEnd; dstPixelPtr += 2) {
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 2 * xSrc;
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A }
0N/A
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A }
0N/A
0N/A return MLIB_SUCCESS;
0N/A}
0N/A
0N/A/***************************************************************/
0N/Amlib_status mlib_ImageAffine_s16_3ch_nn(mlib_affine_param *param)
0N/A{
0N/A DECLAREVAR_NN();
0N/A DTYPE *dstLineEnd;
0N/A
0N/A for (j = yStart; j <= yFinish; j++) {
0N/A mlib_s32 pix0, pix1, pix2;
0N/A
0N/A CLIP(3);
0N/A dstLineEnd = (DTYPE *) dstData + 3 * xRight;
0N/A
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 3 * xSrc;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A pix2 = srcPixelPtr[2];
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/A for (; dstPixelPtr < dstLineEnd; dstPixelPtr += 3) {
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 3 * xSrc;
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A dstPixelPtr[2] = pix2;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A pix2 = srcPixelPtr[2];
0N/A }
0N/A
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A dstPixelPtr[2] = pix2;
0N/A }
0N/A
0N/A return MLIB_SUCCESS;
0N/A}
0N/A
0N/A/***************************************************************/
0N/Amlib_status mlib_ImageAffine_s16_4ch_nn(mlib_affine_param *param)
0N/A{
0N/A DECLAREVAR_NN();
0N/A DTYPE *dstLineEnd;
0N/A
0N/A for (j = yStart; j <= yFinish; j++) {
0N/A mlib_s32 pix0, pix1, pix2, pix3;
0N/A CLIP(4);
0N/A dstLineEnd = (DTYPE *) dstData + 4 * xRight;
0N/A
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 4 * xSrc;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A pix2 = srcPixelPtr[2];
0N/A pix3 = srcPixelPtr[3];
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A for (; dstPixelPtr < dstLineEnd; dstPixelPtr += 4) {
0N/A srcPixelPtr = MLIB_POINTER_GET(lineAddr, ySrc) + 4 * xSrc;
0N/A ySrc = MLIB_POINTER_SHIFT(Y);
0N/A Y += dY;
0N/A xSrc = X >> MLIB_SHIFT;
0N/A X += dX;
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A dstPixelPtr[2] = pix2;
0N/A dstPixelPtr[3] = pix3;
0N/A pix0 = srcPixelPtr[0];
0N/A pix1 = srcPixelPtr[1];
0N/A pix2 = srcPixelPtr[2];
0N/A pix3 = srcPixelPtr[3];
0N/A }
0N/A
0N/A dstPixelPtr[0] = pix0;
0N/A dstPixelPtr[1] = pix1;
0N/A dstPixelPtr[2] = pix2;
0N/A dstPixelPtr[3] = pix3;
0N/A }
0N/A
0N/A return MLIB_SUCCESS;
0N/A}
0N/A
0N/A/***************************************************************/