| | |
| | | #include "image.h" |
| | | #include "utils.h" |
| | | #include <stdio.h> |
| | | #include <math.h> |
| | | |
| | | int windows = 0; |
| | | |
| | |
| | | return copy; |
| | | } |
| | | |
| | | |
| | | void show_image(image p, char *name) |
| | | { |
| | | int x,y,k; |
| | |
| | | return out; |
| | | } |
| | | |
| | | void rotate_image(image m) |
| | | image rotate_image(image im, float rad) |
| | | { |
| | | int i,j; |
| | | for(j = 0; j < m.c; ++j){ |
| | | for(i = 0; i < m.h*m.w/2; ++i){ |
| | | float swap = m.data[j*m.h*m.w + i]; |
| | | m.data[j*m.h*m.w + i] = m.data[j*m.h*m.w + (m.h*m.w-1 - i)]; |
| | | m.data[j*m.h*m.w + (m.h*m.w-1 - i)] = swap; |
| | | int x, y, c; |
| | | float cx = im.w/2.; |
| | | float cy = im.h/2.; |
| | | image rot = make_image(im.w, im.h, im.c); |
| | | for(c = 0; c < im.c; ++c){ |
| | | for(y = 0; y < im.h; ++y){ |
| | | for(x = 0; x < im.w; ++x){ |
| | | float rx = cos(rad)*(x-cx) - sin(rad)*(y-cy) + cx; |
| | | float ry = sin(rad)*(x-cx) + cos(rad)*(y-cy) + cy; |
| | | float val = billinear_interpolate(im, rx, ry, c); |
| | | set_pixel(rot, x, y, c, val); |
| | | } |
| | | } |
| | | } |
| | | return rot; |
| | | } |
| | | |
| | | void translate_image(image m, float s) |
| | |
| | | for(i = 0; i < w; ++i){ |
| | | int r = j + dy; |
| | | int c = i + dx; |
| | | float val = 128; |
| | | float val = 0; |
| | | if (r >= 0 && r < im.h && c >= 0 && c < im.w) { |
| | | val = get_pixel(im, c, r, k); |
| | | } |
| | |
| | | return cropped; |
| | | } |
| | | |
| | | image grayscale_image(image im) |
| | | { |
| | | assert(im.c == 3); |
| | | int i, j, k; |
| | | image gray = make_image(im.w, im.h, im.c); |
| | | float scale[] = {0.114, 0.587, 0.299}; |
| | | for(k = 0; k < im.c; ++k){ |
| | | for(j = 0; j < im.h; ++j){ |
| | | for(i = 0; i < im.w; ++i){ |
| | | gray.data[i+im.w*j] += scale[k]*get_pixel(im, i, j, k); |
| | | } |
| | | } |
| | | } |
| | | memcpy(gray.data + im.w*im.h*1, gray.data, sizeof(float)*im.w*im.h); |
| | | memcpy(gray.data + im.w*im.h*2, gray.data, sizeof(float)*im.w*im.h); |
| | | return gray; |
| | | } |
| | | |
| | | image blend_image(image fore, image back, float alpha) |
| | | { |
| | | assert(fore.w == back.w && fore.h == back.h && fore.c == back.c); |
| | | image blend = make_image(fore.w, fore.h, fore.c); |
| | | int i, j, k; |
| | | for(k = 0; k < fore.c; ++k){ |
| | | for(j = 0; j < fore.h; ++j){ |
| | | for(i = 0; i < fore.w; ++i){ |
| | | float val = alpha * get_pixel(fore, i, j, k) + |
| | | (1 - alpha)* get_pixel(back, i, j, k); |
| | | set_pixel(blend, i, j, k, val); |
| | | } |
| | | } |
| | | } |
| | | return blend; |
| | | } |
| | | |
| | | image saturate_image(image im, float sat) |
| | | { |
| | | image gray = grayscale_image(im); |
| | | image blend = blend_image(im, gray, sat); |
| | | free_image(gray); |
| | | return blend; |
| | | } |
| | | |
| | | image brightness_image(image im, float b) |
| | | { |
| | | image bright = make_image(im.w, im.h, im.c); |
| | | } |
| | | |
| | | float billinear_interpolate(image im, float x, float y, int c) |
| | | { |
| | | int ix = (int) x; |
| | | int iy = (int) y; |
| | | int ix = (int) floorf(x); |
| | | int iy = (int) floorf(y); |
| | | |
| | | float dx = x - ix; |
| | | float dy = y - iy; |
| | | |
| | | float val = (1-dy) * (1-dx) * get_pixel_extend(im, ix, iy, c) + |
| | | dy * (1-dx) * get_pixel_extend(im, ix, iy+1, c) + |
| | | (1-dy) * dx * get_pixel_extend(im, ix+1, iy, c) + |
| | | dy * dx * get_pixel_extend(im, ix+1, iy+1, c); |
| | | dy * (1-dx) * get_pixel_extend(im, ix, iy+1, c) + |
| | | (1-dy) * dx * get_pixel_extend(im, ix+1, iy, c) + |
| | | dy * dx * get_pixel_extend(im, ix+1, iy+1, c); |
| | | return val; |
| | | } |
| | | |
| | |
| | | void test_resize(char *filename) |
| | | { |
| | | image im = load_image(filename, 0,0); |
| | | translate_image(im, -128); |
| | | image small = resize_image(im, 65, 63); |
| | | image big = resize_image(im, 513, 512); |
| | | image crop = crop_image(im, 50, 10, 100, 100); |
| | | image crop2 = crop_image(im, -30, -50, 291, 400); |
| | | image rot = rotate_image(big, .02); |
| | | image rot2 = rotate_image(big, 3.14159265/2.); |
| | | image test = rotate_image(im, .6); |
| | | image gray = grayscale_image(im); |
| | | image sat = saturate_image(im, 2); |
| | | image sat2 = saturate_image(im, .5); |
| | | show_image(im, "original"); |
| | | show_image(small, "smaller"); |
| | | show_image(big, "bigger"); |
| | | show_image(crop, "crop"); |
| | | show_image(crop2, "crop2"); |
| | | show_image(gray, "gray"); |
| | | show_image(sat, "sat"); |
| | | show_image(sat2, "sat2"); |
| | | /* |
| | | show_image(small, "smaller"); |
| | | show_image(big, "bigger"); |
| | | show_image(crop, "crop"); |
| | | show_image(crop2, "crop2"); |
| | | show_image(rot, "rot"); |
| | | show_image(rot2, "rot2"); |
| | | show_image(test, "test"); |
| | | */ |
| | | cvWaitKey(0); |
| | | } |
| | | |