| /* |
| * Copyright (c) 2003, Oracle and/or its affiliates. All rights reserved. |
| * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER. |
| * |
| * This code is free software; you can redistribute it and/or modify it |
| * under the terms of the GNU General Public License version 2 only, as |
| * published by the Free Software Foundation. Oracle designates this |
| * particular file as subject to the "Classpath" exception as provided |
| * by Oracle in the LICENSE file that accompanied this code. |
| * |
| * This code is distributed in the hope that it will be useful, but WITHOUT |
| * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
| * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
| * version 2 for more details (a copy is included in the LICENSE file that |
| * accompanied this code). |
| * |
| * You should have received a copy of the GNU General Public License version |
| * 2 along with this work; if not, write to the Free Software Foundation, |
| * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. |
| * |
| * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA |
| * or visit www.oracle.com if you need additional information or have any |
| * questions. |
| */ |
| |
| |
| /* |
| * FUNCTION |
| * mlib_ImageClipping |
| * mlib_ImageClippingMxN |
| * Clipping for image processing in case of pixel-to-pixel |
| * square kernel filtering. Source and destination images can have |
| * different sizes, center of the source is mapped to the center of |
| * the destination image. |
| * Examples of this type of image processing are Convolve, Gradient, |
| * Dilate/Erode functions, etc. |
| * |
| * SYNOPSIS |
| * mlib_status mlib_ImageClipping(mlib_image *dst_i, |
| * mlib_image *src_i, |
| * mlib_image *dst_e, |
| * mlib_image *src_e, |
| * mlib_s32 *edg_sizes, |
| * const mlib_image *dst, |
| * const mlib_image *src, |
| * mlib_s32 ker_size) |
| * |
| * mlib_status mlib_ImageClippingMxN(mlib_image *dst_i, |
| * mlib_image *src_i, |
| * mlib_image *dst_e, |
| * mlib_image *src_e, |
| * mlib_s32 *edg_sizes, |
| * const mlib_image *dst, |
| * const mlib_image *src, |
| * mlib_s32 kw, |
| * mlib_s32 kh, |
| * mlib_s32 kw1, |
| * mlib_s32 kh1) |
| * |
| * OUTPUT ARGUMENTS |
| * dst_i Pointer to destination image of internal pixels |
| * src_i Pointer to source image of internal pixels |
| * dst_e Pointer to destination image for edge processing |
| * src_e Pointer to source image for edge processing |
| * edg_sizes Array of edge sizes |
| * |
| * INPUT ARGUMENTS |
| * dst Pointer to destination image. |
| * src Pointer to source image. |
| * ksize Size of kernel |
| * |
| * RESTRICTION |
| * The src and the dst must be images of the same type. |
| * The src and dst must have same number of channels. |
| * |
| */ |
| |
| #include "mlib_image.h" |
| #include "mlib_ImageCheck.h" |
| #include "mlib_ImageClipping.h" |
| #include "mlib_ImageCreate.h" |
| |
| /***************************************************************/ |
| mlib_status mlib_ImageClippingMxN(mlib_image *dst_i, |
| mlib_image *src_i, |
| mlib_image *dst_e, |
| mlib_image *src_e, |
| mlib_s32 *edg_sizes, |
| const mlib_image *dst, |
| const mlib_image *src, |
| mlib_s32 kw, |
| mlib_s32 kh, |
| mlib_s32 kw1, |
| mlib_s32 kh1) |
| { |
| mlib_s32 kw2 = kw - 1 - kw1; |
| mlib_s32 kh2 = kh - 1 - kh1; |
| mlib_s32 src_wid, src_hgt, dst_wid, dst_hgt; |
| mlib_s32 dx, dy, dxd, dxs, dyd, dys, wid_e, hgt_e; |
| mlib_s32 dx_l, dx_r, dy_t, dy_b, wid_i, hgt_i; |
| |
| MLIB_IMAGE_CHECK(dst); |
| MLIB_IMAGE_CHECK(src); |
| MLIB_IMAGE_TYPE_EQUAL(dst, src); |
| MLIB_IMAGE_CHAN_EQUAL(dst, src); |
| |
| dst_wid = mlib_ImageGetWidth(dst); |
| dst_hgt = mlib_ImageGetHeight(dst); |
| src_wid = mlib_ImageGetWidth(src); |
| src_hgt = mlib_ImageGetHeight(src); |
| |
| /* X clipping */ |
| dx = src_wid - dst_wid; |
| |
| if (dx > 0) { |
| dxs = (dx + 1) >> 1; |
| dxd = 0; |
| } else { |
| dxs = 0; |
| dxd = (-dx) >> 1; |
| } |
| |
| dx_l = kw1 - dxs; |
| dx_r = kw2 + dxs - dx; |
| |
| if (dx_l < 0) dx_l = 0; |
| if (dx_r < 0) dx_r = 0; |
| if (dx_r > kw2) dx_r = kw2; |
| |
| /* Y clipping */ |
| dy = src_hgt - dst_hgt; |
| |
| if (dy > 0) { |
| dys = (dy + 1) >> 1; |
| dyd = 0; |
| } else { |
| dys = 0; |
| dyd = (-dy) >> 1; |
| } |
| |
| dy_t = kh1 - dys; |
| dy_b = kh2 + dys - dy; |
| |
| if (dy_t < 0) dy_t = 0; |
| if (dy_b < 0) dy_b = 0; |
| if (dy_b > kh2) dy_b = kh2; |
| |
| /* image sizes */ |
| wid_e = (src_wid < dst_wid) ? src_wid : dst_wid; |
| hgt_e = (src_hgt < dst_hgt) ? src_hgt : dst_hgt; |
| wid_i = wid_e + (kw1 - dx_l) + (kw2 - dx_r); |
| hgt_i = hgt_e + (kh1 - dy_t) + (kh2 - dy_b); |
| |
| mlib_ImageSetSubimage(dst_i, dst, dxd - (kw1 - dx_l), dyd - (kh1 - dy_t), wid_i, hgt_i); |
| mlib_ImageSetSubimage(src_i, src, dxs - (kw1 - dx_l), dys - (kh1 - dy_t), wid_i, hgt_i); |
| |
| if (dst_e != NULL && src_e != NULL) { /* images for edge processing */ |
| mlib_ImageSetSubimage(dst_e, dst, dxd, dyd, wid_e, hgt_e); |
| mlib_ImageSetSubimage(src_e, src, dxs, dys, wid_e, hgt_e); |
| } |
| |
| if (edg_sizes != NULL) { /* save edges */ |
| edg_sizes[0] = dx_l; |
| edg_sizes[1] = dx_r; |
| edg_sizes[2] = dy_t; |
| edg_sizes[3] = dy_b; |
| } |
| |
| return MLIB_SUCCESS; |
| } |
| |
| /***************************************************************/ |
| mlib_status mlib_ImageClipping(mlib_image *dst_i, |
| mlib_image *src_i, |
| mlib_image *dst_e, |
| mlib_image *src_e, |
| mlib_s32 *edg_sizes, |
| const mlib_image *dst, |
| const mlib_image *src, |
| mlib_s32 ker_size) |
| { |
| mlib_s32 kw1 = (ker_size - 1)/2; |
| return mlib_ImageClippingMxN(dst_i, src_i, dst_e, src_e, edg_sizes, |
| dst, src, ker_size, ker_size, kw1, kw1); |
| } |
| |
| /***************************************************************/ |
