pikl 0.1.0

data/ext/pikl/pikl_io.h

@@ -0,0 +1,27 @@
+#ifndef _LIB_PIKL_IO_
+#define _LIB_PIKL_IO_
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+
+#include "pikl.h"
+#include "pikl_private.h"
+#include "pikl_jpeg.h"
+#include "pikl_bitmap.h"
+#include "pikl_png.h"
+
+#define PKL_HEADER_SIZE 8
+struct PKL_HEADER {
+	PKL_FORMAT format;
+	int size;
+	unsigned char mark[PKL_HEADER_SIZE];
+};
+
+static struct PKL_HEADER pklheader[] = {
+	{ PKL_FORMAT_JPEG,   2, {0xFF, 0xD8} },
+	{ PKL_FORMAT_PNG,    8, {0x89, 0x50, 0x4E, 0x47, 0x0D, 0x0A, 0x1A, 0x0A} },
+	{ PKL_FORMAT_BITMAP, 2, {0x42, 0x4D} },
+};
+
+#endif

data/ext/pikl/pikl_jpeg.c

@@ -0,0 +1,161 @@
+#include "pikl_jpeg.h"
+
+static void pkljpeg_error_exit(j_common_ptr cinfo);
+static int color_p2j(int color);
+static int color_j2p(int color);
+static int compress_p2j(int level);
+
+//=============================================================================
+// pkljpeg_error_exit
+//=============================================================================
+static void pkljpeg_error_exit(j_common_ptr cinfo)
+{
+	struct err_mgr *jp_err = (struct err_mgr *)cinfo->err;
+	longjmp(jp_err->setjmp_buffer, 1);
+}
+
+//=============================================================================
+// load_jpeg
+//=============================================================================
+int load_jpeg(PKLImage pkl, FILE *image)
+{
+	struct jpeg_decompress_struct cinfo;
+	struct err_mgr jerr;
+	unsigned char *p;
+	
+	cinfo.err = jpeg_std_error(&jerr.pub);
+	jerr.pub.error_exit = pkljpeg_error_exit;
+	
+	if( setjmp(jerr.setjmp_buffer) ){
+		jpeg_destroy_decompress(&cinfo);
+		return(1);
+	}
+
+	jpeg_create_decompress(&cinfo);
+	jpeg_stdio_src(&cinfo, image);
+	jpeg_read_header(&cinfo, TRUE);
+	jpeg_start_decompress(&cinfo);
+	
+	pkl->width   = cinfo.output_width;
+	pkl->height  = cinfo.output_height;
+	pkl->color   = color_j2p(cinfo.out_color_space);
+	pkl->channel = cinfo.output_components;
+	if(pkl->color == PKL_UNKNOWN) return(1);
+	
+	pkl->image = malloc(cinfo.output_height * cinfo.output_width * cinfo.output_components);
+	if(!pkl->image) return(1);
+	
+	p = pkl->image;
+	while(cinfo.output_scanline < cinfo.output_height){
+		jpeg_read_scanlines(&cinfo, &p, 1);
+		p += (cinfo.output_width * cinfo.output_components);
+	}
+	jpeg_finish_decompress(&cinfo);
+	jpeg_destroy_decompress(&cinfo);
+
+	return(0);
+}
+
+//=============================================================================
+// save_jpg
+//=============================================================================
+int save_jpeg(PKLImage pkl, FILE *image)
+{
+	struct jpeg_compress_struct oinfo;
+	struct err_mgr jerr;
+	JSAMPROW   row_pointer[1];
+	int i;
+
+	if(pkl->color == PKL_UNKNOWN) return(1);
+	
+	oinfo.err = jpeg_std_error(&jerr.pub);
+	jerr.pub.error_exit = pkljpeg_error_exit;
+	
+	if( setjmp(jerr.setjmp_buffer) ){
+		jpeg_destroy_compress(&oinfo);
+		return(1);
+	}
+
+	jpeg_create_compress(&oinfo);
+	jpeg_stdio_dest(&oinfo, image);
+
+	oinfo.image_width      = pkl->width;
+	oinfo.image_height     = pkl->height;
+	oinfo.in_color_space   = color_p2j(pkl->color);
+	oinfo.input_components = pkl->channel;
+	jpeg_set_defaults(&oinfo);
+	jpeg_set_quality(&oinfo, compress_p2j(pkl->compress), TRUE);
+	oinfo.write_JFIF_header= TRUE;
+	oinfo.density_unit     = 1;			//(1)dots/inch, (2)dots/cm
+	
+	jpeg_start_compress(&oinfo, TRUE);
+	
+	for(i=0; i<pkl->height ; i++){
+		row_pointer[0] = &pkl->image[i*pkl->width*pkl->channel];
+		jpeg_write_scanlines(&oinfo, row_pointer, 1);
+	}
+	
+	jpeg_finish_compress(&oinfo);
+	jpeg_destroy_compress(&oinfo);
+	
+
+	return(0);
+}
+
+//=============================================================================
+// color_p2j
+//=============================================================================
+static int color_p2j(int color)
+{
+	switch(color){
+		case PKL_GRAYSCALE:
+			return JCS_GRAYSCALE;
+		case PKL_RGB:
+			return JCS_RGB;
+		//case PKL_YCbCr:
+		//	return JCS_YCbCr;
+		case PKL_CMYK:
+			return JCS_CMYK;
+		//case PKL_YCCK:
+		//	return JCS_YCCK;
+		//case PKL_RGBA:
+		default:
+			return JCS_UNKNOWN;
+	}
+	return JCS_UNKNOWN;
+}
+
+//=============================================================================
+// color_j2p
+//=============================================================================
+static int color_j2p(int color)
+{
+	switch(color){
+		case JCS_GRAYSCALE:
+			return PKL_GRAYSCALE;
+		case JCS_RGB:
+			return PKL_RGB;
+		//case JCS_YCbCr:
+		//	return PKL_CMYK;
+		case JCS_CMYK:
+			return PKL_CMYK;
+		//case JCS_YCCK:
+		//	return PKL_YCCK;
+		default:
+			return PKL_UNKNOWN;
+	}
+	return PKL_UNKNOWN;
+}
+
+//=============================================================================
+// compress_p2j
+//=============================================================================
+static int compress_p2j(int level)
+{
+	int comp;
+	
+	if(level < 0) return(90);
+	comp = 100 - (level*10);
+	if(comp <= 0) return(1);
+	return(comp);
+}

data/ext/pikl/pikl_jpeg.h

@@ -0,0 +1,22 @@
+#ifndef _LIB_PIKL_JPEG_
+#define _LIB_PIKL_JPEG_
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include <setjmp.h>
+#include <jpeglib.h>
+
+#include "pikl.h"
+#include "pikl_private.h"
+
+/* libjpeg error handler */
+struct err_mgr {
+	struct jpeg_error_mgr pub;
+	jmp_buf setjmp_buffer;
+};
+
+int load_jpeg(PKLImage pkl, FILE *image);
+int save_jpeg(PKLImage pkl, FILE *image);
+
+#endif

data/ext/pikl/pikl_png.c

@@ -0,0 +1,188 @@
+#include "pikl_png.h"
+
+static int color_pn2pk(int color);
+static int color_pk2pn(int color);
+static int pkl_png_channels(int color);
+static int pkl_png_depth(int color);
+static int compress_pk2pn(int level);
+
+//=============================================================================
+// load_png
+//=============================================================================
+int load_png(PKLImage pkl, FILE *image)
+{
+	png_structp png_ptr;
+	png_infop   info_ptr;
+	int i;
+	
+	png_ptr = png_create_read_struct(PNG_LIBPNG_VER_STRING, NULL, NULL, NULL);
+	if(!png_ptr) return(1);
+	
+	info_ptr = png_create_info_struct(png_ptr);
+	if(!info_ptr){
+		png_destroy_read_struct(&png_ptr, NULL, NULL);
+		return(1);
+	}
+
+	if( setjmp(png_jmpbuf(png_ptr)) ){
+		png_destroy_info_struct(png_ptr, &info_ptr);
+		png_destroy_read_struct(&png_ptr, &info_ptr, (png_infopp)NULL);
+		return(1);
+	}
+
+	png_init_io(png_ptr, image);
+	png_read_info(png_ptr, info_ptr);
+	pkl->width   = info_ptr->width;
+	pkl->height  = info_ptr->height;
+	pkl->color   = color_pn2pk(info_ptr->color_type);
+	pkl->channel = pkl_png_channels(pkl->color);
+	if(pkl->color == PKL_UNKNOWN) return(1);
+
+	pkl->image = malloc(pkl->width * pkl->height * pkl->channel);
+	if(!pkl->image) return(1);
+
+	if(info_ptr->color_type==PNG_COLOR_TYPE_RGB_ALPHA  ||
+	   info_ptr->color_type==PNG_COLOR_TYPE_GRAY_ALPHA || info_ptr->num_trans) png_set_strip_alpha(png_ptr);
+	if(info_ptr->bit_depth<8 &&
+	   info_ptr->color_type==PNG_COLOR_TYPE_GRAY) png_set_gray_1_2_4_to_8(png_ptr);
+	if(info_ptr->bit_depth==16) png_set_strip_16(png_ptr);
+	if(info_ptr->color_type==PNG_COLOR_TYPE_PALETTE) png_set_palette_to_rgb(png_ptr);
+
+	for(i=0; i<pkl->height; i++)
+		png_read_row(png_ptr, &pkl->image[i*pkl->width*pkl->channel], NULL);
+	
+	png_destroy_info_struct(png_ptr, &info_ptr);
+	png_destroy_read_struct(&png_ptr, &info_ptr, (png_infopp)NULL);
+	
+	return(0);
+}
+
+//===================================================================
+// save_png
+//===================================================================
+int save_png(PKLImage pkl, FILE *image)
+{
+	png_structp png_ptr;
+	png_infop   info_ptr;
+	int  i;
+	
+	png_ptr = png_create_write_struct(PNG_LIBPNG_VER_STRING, NULL, NULL, NULL);
+	if(!png_ptr) return(1);
+	
+	info_ptr = png_create_info_struct(png_ptr);
+	if(!info_ptr){
+		png_destroy_write_struct(&png_ptr,  (png_infopp)NULL);
+		return(1);
+	}
+	
+	if( setjmp(png_jmpbuf(png_ptr)) ){
+		png_destroy_info_struct(png_ptr, &info_ptr);
+		png_destroy_write_struct(&png_ptr, &info_ptr);
+		return(1);
+	}
+	
+	png_init_io(png_ptr, image);
+	info_ptr->width  = pkl->width;
+	info_ptr->height = pkl->height;
+	info_ptr->bit_depth  = pkl_png_depth(pkl->color);
+	info_ptr->color_type = color_pk2pn(pkl->color);
+	png_set_compression_level(png_ptr, compress_pk2pn(pkl->compress));
+	png_write_info(png_ptr, info_ptr);
+	
+	for(i=0; i<pkl->height; i++)
+		png_write_row(png_ptr, &pkl->image[i*pkl->width*pkl->channel]);
+	png_write_end(png_ptr, info_ptr);
+	png_destroy_info_struct(png_ptr, &info_ptr);
+	png_destroy_write_struct(&png_ptr, &info_ptr);
+
+	return(0);
+}
+
+//=============================================================================
+// color_pn2pk
+//=============================================================================
+static int color_pn2pk(int color)
+{
+	switch(color){
+		case PNG_COLOR_TYPE_GRAY:
+		case PNG_COLOR_TYPE_GRAY_ALPHA:
+			return(PKL_GRAYSCALE);
+		case PNG_COLOR_TYPE_PALETTE:
+		case PNG_COLOR_TYPE_RGB:
+		case PNG_COLOR_TYPE_RGB_ALPHA:
+			return(PKL_RGB);
+		default:
+			return(PKL_UNKNOWN);
+	}
+	return(PKL_UNKNOWN);
+}
+
+//=============================================================================
+// pkl_png_channels
+//=============================================================================
+static int pkl_png_channels(int color)
+{
+	switch(color){
+		//case PKL_BLACKWHITE:
+		//	return ;
+		case PKL_GRAYSCALE:
+			return 1;
+		case PKL_RGB:
+			return 3;
+		//case PKL_RGBA:
+		//case PKL_CMYK:
+		//	return 4;
+		default:
+			return 0;
+	}
+	return 0;
+}
+
+//=============================================================================
+// pkl_png_depth
+//=============================================================================
+static int pkl_png_depth(int color)
+{
+	switch(color){
+		//case PKL_BLACKWHITE:
+		//	return 1;
+		case PKL_GRAYSCALE:
+		case PKL_RGB:
+		//case PKL_CMYK:
+		//case PKL_RGBA:
+			return 8;
+		default:
+			return 0;
+	}
+	return 0;
+}
+
+//=============================================================================
+// color_pk2pn
+//=============================================================================
+static int color_pk2pn(int color)
+{
+	switch(color){
+		//case PKL_BLACKWHITE:
+		//	return PNG_COLOR_TYPE_GRAY;
+		case PKL_GRAYSCALE:
+			return PNG_COLOR_TYPE_GRAY;
+		case PKL_RGB:
+			return PNG_COLOR_TYPE_RGB;
+		//case PKL_CMYK:
+		//case PKL_RGBA:
+		default:
+			return -1;
+	}
+	return -1;
+}
+
+//=============================================================================
+// compress_pk2pn
+//=============================================================================
+static int compress_pk2pn(int level)
+{
+	if(level < 0) return(Z_BEST_COMPRESSION);
+	if(level > 9) return(9);
+	return(level);
+}

data/ext/pikl/pikl_png.h

@@ -0,0 +1,17 @@
+#ifndef _LIB_PIKL_PNG_
+#define _LIB_PIKL_PNG_
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include <setjmp.h>
+#include <png.h>
+
+#include "pikl.h"
+#include "pikl_private.h"
+
+int load_png(PKLImage pkl, FILE *image);
+int save_png(PKLImage pkl, FILE *image);
+
+#endif
+

data/ext/pikl/pikl_private.h

@@ -0,0 +1,22 @@
+#ifndef _LIB_PIKL_PRIVATE_
+#define _LIB_PIKL_PRIVATE_
+
+#define PKL_MINMAX(a,b,i) (((a>=i)?a:((b<=i)?b:i)))
+#define PKL_RGB(v) v<0?0:v>255?255:v
+#define PKL_CHANNEL 4
+#define PKL_COLOR 256
+
+//=============================================================================
+// PKLImage
+//=============================================================================
+struct _PKLImage {
+	int width;
+	int height;
+	PKL_COLOR_SPACE color;
+	int channel;
+	unsigned char *image;
+	PKL_FORMAT format;
+	int compress;
+};
+
+#endif

data/ext/pikl/pikl_resize.c

@@ -0,0 +1,338 @@
+#include "pikl_resize.h"
+
+static int resize_pa(PKLImage pkl, int width, int height);
+static int resize_nn(PKLImage pkl, int width, int height);
+static int resize_lz(PKLImage pkl, int width, int height);
+static int resize_bc(PKLImage pkl, int width, int height);
+static int resize_bl(PKLImage pkl, int width, int height);
+static int resize_lz_cr(PKLImage pkl, double scale, int ax);
+static double lz_weight(double x);
+static int get_pixel(PKLImage pkl, int x, int y, unsigned char *pix);
+
+//=================================================================================
+// pkl_resize
+//=================================================================================
+PKLExport int pkl_resize(PKLImage pkl, int width, int height, PKL_SAMPLE sample)
+{
+	// color mode check
+	
+	if(width < 1 || height < 1) return(1);
+	if(width==pkl->width && height==pkl->height) return(0);
+	
+	switch(sample){
+		case PKL_SAMPLE_NN:
+			return resize_nn(pkl, width, height);
+		case PKL_SAMPLE_BL:
+			return resize_bl(pkl, width, height);
+		case PKL_SAMPLE_BC:
+			return resize_bc(pkl, width, height);
+		case PKL_SAMPLE_PA:
+			if(width<pkl->width && height<pkl->height) return resize_pa(pkl, width, height);
+		case PKL_SAMPLE_LZ:
+			return resize_lz(pkl, width, height);
+		default:
+			return(1);
+	}
+	return(1);
+}
+
+//=============================================================================
+// resize_nn
+//=============================================================================
+static int resize_nn(PKLImage pkl, int width, int height)
+{
+	unsigned char *wrk;
+	double zx, zy;
+	int i, j, px, py;
+	
+	wrk = malloc(width*height*pkl->channel);
+	if(!wrk) return(1);
+	memset(wrk, 0, width*height*pkl->channel);
+	
+	zx = (double)width /pkl->width;
+	zy = (double)height/pkl->height;
+
+	for(i=0; i<height; i++){
+		for(j=0; j<width; j++){
+			px = (int)((double)j)/zx;
+			py = (int)((double)i)/zy;
+			if(px>=0 && px<pkl->width && py>=0 && py<pkl->height){
+				memcpy(&wrk[(i*width+j)*pkl->channel], &pkl->image[(py*pkl->width+px)*pkl->channel], pkl->channel);
+			}
+		}
+	}
+	free(pkl->image);
+	pkl->image = wrk;
+	pkl->width = width;
+	pkl->height = height;
+	return(0);
+}
+
+//=============================================================================
+// resize_bl
+//=============================================================================
+static int resize_bl(PKLImage pkl, int width, int height)
+{
+	unsigned char *wrk;
+	double zx, zy, dx, dy, px, py, f[4];
+	int i, j, k, wx, wy, sx, sy;
+	
+	wrk = malloc(width*height*pkl->channel);
+	if(!wrk) return(1);
+	memset(wrk, 0, width*height*pkl->channel);
+	
+	zx = (double)pkl->width / (double)width;
+	zy = (double)pkl->height/ (double)height;
+
+	for(i=0; i<height; i++){
+		for(j=0; j<width; j++){
+			dx = modf(zx*j, &px);
+			dy = modf(zy*i, &py);
+			f[0] = (1.0-dx)*(1.0-dy);
+			f[1] = dx*(1.0-dy);
+			f[2] = (1.0-dx)*dy;
+			f[3] = dx*dy;
+
+			for(wy=0; wy<2; wy++){
+				for(wx=0; wx<2; wx++){
+					sx = ((int)px+wx<0) ? pkl->width +((int)px+wx)%pkl->width -1 : ((int)px+wx)%pkl->width;
+					sy = ((int)py+wy<0) ? pkl->height+((int)py+wy)%pkl->height-1 : ((int)py+wy)%pkl->height;
+					for(k=0; k<pkl->channel; k++){
+						wrk[(i*width+j)*pkl->channel+k] = 
+							PKL_RGB(wrk[(i*width+j)*pkl->channel+k] + pkl->image[(sy*pkl->width+sx)*pkl->channel+k] * f[wy*2+wx]);
+					}
+				}
+			}
+		}
+	}
+	
+	free(pkl->image);
+	pkl->image = wrk;
+	pkl->width = width;
+	pkl->height = height;
+	return(0);
+}
+
+//=============================================================================
+// resize_bc
+//=============================================================================
+static int resize_bc(PKLImage pkl, int width, int height)
+{
+	int i, j, k, px, py, x0, y0;
+	double wx, wy, dx, dy, zx, zy, sx, sy;
+	double data[PKL_CHANNEL];
+	unsigned char *wrk;
+	
+	wrk = malloc(width*height*pkl->channel);
+	if(!wrk) return(1);
+	memset(wrk, 0, width*height*pkl->channel);
+
+	zx = (double)pkl->width / (double)width;
+	zy = (double)pkl->height/ (double)height;
+
+	for(i=0; i<height; i++){
+		for(j=0; j<width; j++){
+			memset(&data, 0, sizeof(data));
+			sx = zx * j;
+			sy = zy * i;
+			for(py=(int)sy-1; py<=(int)sy+2; py++){
+				for(px=(int)sx-1; px<=(int)sx+2; px++){
+					dx = fabs(sx-px);
+					dy = fabs(sy-py);
+					wx = (dx<1) ? (dx-1)*(dx*dx-dx-1) : -(dx-1)*(dx-2)*(dx-2);
+					wy = (dy<1) ? (dy-1)*(dy*dy-dy-1) : -(dy-1)*(dy-2)*(dy-2);
+					x0 = (px<0 || px>pkl->width -1) ? sx : px;
+					y0 = (py<0 || py>pkl->height-1) ? sy : py;
+					for(k=0; k<pkl->channel; k++) data[k]+=pkl->image[(y0*pkl->width+x0)*pkl->channel+k]*wx*wy;
+				}
+			}
+			for(k=0; k<pkl->channel; k++) wrk[(i*width+j)*pkl->channel+k]=PKL_RGB(data[k]);
+		}
+	}
+	free(pkl->image);
+	pkl->image = wrk;
+	pkl->width  = width;
+	pkl->height = height;
+	return(0);
+}
+
+//=============================================================================
+// resize_pa
+//=============================================================================
+static int resize_pa(PKLImage pkl, int width, int height)
+{
+	unsigned char *wrk;
+	int i, j;
+	double zx, zy;
+	
+	wrk = malloc(width * height * pkl->channel);
+	if(!wrk) return(1);
+	memset(wrk, 0, width*height*pkl->channel);
+	
+	zx = (double)pkl->width  / (double)width;
+	zy = (double)pkl->height / (double)height;
+
+	for(i=0; i<height; i++){
+		for(j=0; j<width; j++){
+			double w1, w2, h1, h2;
+			double w1b, w2b, h1b, h2b;
+			double bi, ixg[PKL_CHANNEL] = {0.0, 0.0, 0.0, 0.0};
+			int fsx ,fsy, s, t, k;
+			
+			w1b = modf(zx *  j,    &w1);
+			w2b = modf(zx * (j+1), &w2);
+			h1b = modf(zy *  i,    &h1);
+			h2b = modf(zy * (i+1), &h2);
+
+			for(s=(int)w1; s<=(int)w2; s++){
+				for(t=(int)h1; t<=(int)h2; t++){
+					bi = 1.0;
+					if(s == (int)w1) bi*=(1-w1b);
+					if(s == (int)w2) bi*=(w2b);
+					if(t == (int)h1) bi*=(1-h1b);
+					if(t == (int)h2) bi*=(h2b);
+					fsx = (s < 0) ? 0 : (s > pkl->width -1) ? pkl->width -1 : s;
+					fsy = (t < 0) ? 0 : (t > pkl->height-1) ? pkl->height-1 : t;
+					for(k=0; k<pkl->channel; k++)
+						ixg[k] += pkl->image[(fsy*pkl->width+fsx)*pkl->channel+k] * bi;
+				}
+			}
+			for(k=0; k<pkl->channel; k++)
+				wrk[(i*width+j)*pkl->channel+k] = PKL_RGB(ixg[k]/(zx*zy));
+		}
+	}
+	free(pkl->image);
+	pkl->image = wrk;
+	pkl->width  = width;
+	pkl->height = height;
+	return(0);
+}
+
+//=============================================================================
+// resize_lz
+//=============================================================================
+static int resize_lz(PKLImage pkl, int width, int height)
+{
+	double scale;
+	
+	scale = (double)width / (double)pkl->width;
+	if( resize_lz_cr(pkl, scale, 1) ) return(1);
+	
+	scale = (double)height / (double)pkl->height;
+	if( resize_lz_cr(pkl, scale, 0) ) return(1);
+	
+	return(0);
+}
+
+//=============================================================================
+// resize_lz_cr
+//=============================================================================
+static int resize_lz_cr(PKLImage pkl, double scale, int ax)
+{
+	int hsize, vsize, i, j, k, l, lMax, zx, zy, factor;
+	double *w0, *w1, T, t, F, data[PKL_CHANNEL];
+	unsigned char *wrk, v[PKL_CHANNEL], v1[PKL_CHANNEL], v2[PKL_CHANNEL];
+	int mx, Mx;
+
+	if(ax){
+		zy = hsize = (int)(0.5 + pkl->width * scale);
+		zx = vsize = pkl->height;
+	}else{
+		zx = hsize = pkl->width;
+		zy = vsize = (int)(0.5 + pkl->height * scale);
+	}
+
+	wrk = malloc(hsize*vsize*pkl->channel);
+	if(!wrk) return(1);
+	memset(wrk, 0, hsize*vsize*pkl->channel);
+
+	factor = (scale>=1) ? 1 : (int)(1.4/scale); 
+	lMax = factor << 1;
+	w0 = malloc(lMax * sizeof(double));
+	w1 = malloc(lMax * sizeof(double));
+
+	for(j=0; j<zy; j++){
+		F = modf((double)j/scale, &T);
+		for(l=0; l<lMax; l++){
+			w0[lMax-l-1] = lz_weight( ((lMax-l-1) + F)/factor );
+			w1[l]        = lz_weight( ((l+1)      - F)/factor );
+		}
+		t = 0.0;
+		for(l=0; l<lMax; l++){
+			t+=w0[l];
+			t+=w1[l];
+		}
+		t /= (double)factor;
+		for(l=0; l<lMax; l++){
+			w0[l] /= t;
+			w1[l] /= t;
+		}
+
+		if(ax){
+			for(i=0; i<zx; i++){
+				memset(data, 0, sizeof(data));
+				for(l=0; l<lMax; l++){
+					mx = (T-lMax+l+1 < 0) ? 0 : T-lMax+l+1;
+					Mx = (T+l+1 >= pkl->width) ? pkl->width-1 : T+l+1;
+					get_pixel(pkl, Mx, i, v1);
+					get_pixel(pkl, mx, i, v2);
+					for(k=0; k<pkl->channel; k++) data[k]+=(w1[l]*v1[k]+w0[lMax-l-1]*v2[k]);
+				}
+				for(k=0; k<pkl->channel; k++) v[k]=PKL_MINMAX(0, 255, 0.5+(data[k]/factor));
+				if(i>=0 && i<vsize && j>=0 && j<hsize) memcpy(&wrk[(i*hsize+j)*pkl->channel], &v, pkl->channel);
+			}
+		}else{
+			for(i=0; i<zx; i++){
+				memset(data, 0, sizeof(data));
+				for(l=0; l<lMax; l++){
+					mx = (T-lMax+l+1 < 0) ? 0 : T-lMax+l+1;
+					Mx = (T+l+1 >= pkl->height) ? pkl->height-1 : T+l+1;
+					get_pixel(pkl, i, Mx, v1);
+					get_pixel(pkl, i, mx, v2);
+					for(k=0; k<pkl->channel; k++) data[k]+=(w1[l]*v1[k]+w0[lMax-l-1]*v2[k]);
+				}
+				for(k=0; k<pkl->channel; k++) v[k]=PKL_MINMAX(0, 255, 0.5+(data[k]/factor));
+				if(j>=0 && j<vsize && i>=0 && i<hsize) memcpy(&wrk[(j*hsize+i)*pkl->channel], &v, pkl->channel);
+			}
+		}
+	}
+	free(w0);
+	free(w1);
+	free(pkl->image);
+	pkl->image  = wrk;
+	pkl->width  = hsize;
+	pkl->height = vsize;
+	return(0);
+}
+
+//=============================================================================
+// get_pixel
+//=============================================================================
+static int get_pixel(PKLImage pkl, int x, int y, unsigned char *pix)
+{
+	if(x<0 || x>=pkl->width || y<0 || y>=pkl->height) return(1);
+	memset(pix, 0, pkl->channel);
+	memcpy(pix, &pkl->image[(y*pkl->width+x)*pkl->channel], pkl->channel);
+	return(0);
+}
+
+//=============================================================================
+// lz_weight
+//=============================================================================
+static double lz_weight(double x)
+{
+	double PIx, PIx2;
+	
+	PIx  = M_PI * x;
+	PIx2 = PIx / 2.0;
+
+	if(x>=2.0 || x<=-2.0){
+		return (0.0);
+	}else
+	if(x == 0.0){
+		return (1.0);
+	}else{
+		return(sin(PIx) / PIx * sin(PIx2) / PIx2);
+	}
+}
+