矩阵论还真有用：图片缩放旋转

#include <stdio.h>
#include <stdlib.h>
#include <unistd.h>
#include <math.h>


#define SCREEN_W 800
#define SCREEN_H 480


//#pragma pack(1)


typedef struct{
    short  bfType;     //bfType(2字节)，这里恒定等于&H4D42,ASCII字符'BM'
    unsigned int bfSize;     //文件大小，以4字节为单位
    short  bfReserve1; //备用
    short  bfReserve2; //备用
    unsigned int bfoffBits;  //数据区在文件中的位置偏移量
}__attribute__((packed))BITMAPFILEHEADER;   //文件头结构体，14字节


typedef struct{
    long bitSize;    //位图信息头大小 本结构体大小 40btye
    long biWidth;    //图象宽度，像素单位
    long biHeight;   //图象高度，像素单位
    short  biPlanes;   //位平面树=1
    short  biBitCount; //单位像素的位数，表示bmp图片的颜色位数，即24位图、32位图
    long biCompression;  //图片的压缩属性，bmp图片是不压缩的，等于0
    long biSizeImage;    //表示bmp图片数据区的大小，当上一个属性biCompression等于0时，这里的值可以省略不填
    long biXPlosPerMeter;//水平分辨率，可省略
    long biYPlosPerMeter;//垂直分辨率，可省略
    long biClrUsed;      //表示使用了多少个颜色索引表，一般biBitCount属性小于16才会用到，
                                  //等于0时表示有2^biBitCount个颜色索引表
    long biClrImportant; //表示有多少个重要的颜色，等于0时表示所有颜色都很重要
}BITMAPINFOHEADER;       //位图信息头，40字节


typedef struct tagRGBQUAD {
    unsigned char rgbBlue;     // 蓝色的亮度(值范围为0-255)
    unsigned char rgbGreen;    // 绿色的亮度(值范围为0-255)
    unsigned char rgbRed;      // 红色的亮度(值范围为0-255)
    unsigned char rgbReserved; // 保留，必须为0
} RGBQUAD;


typedef struct tagBITMAPINFO {
    BITMAPINFOHEADER bmiHeader; // 位图信息头
    RGBQUAD bmiColors[1];       // 颜色表
} BITMAPINFO;


//#pragma pack()


struct pic_obj {
    int nxp,nyp;
    unsigned bmap[256];
    unsigned char *data;        
};

struct pic_obj pics[2] = {0};
//将图片读入内存
int init_pics(const char *file, int num)
{
    int i, j;
    int nread;
    FILE *fp;
    char buf[1024*4];
    BITMAPFILEHEADER header;
    BITMAPINFOHEADER info;
    
    fp = fopen(file, "rb");
    if(fp == NULL){
        printf("open pic failed!\n");
    }

    nread = fread(&header, 1, sizeof(BITMAPFILEHEADER), fp);
    if(nread != sizeof(BITMAPFILEHEADER)){
        printf("read header fail\n");
    }
    
    nread = fread(&info, 1, sizeof(BITMAPINFOHEADER), fp);
    if(nread != sizeof(BITMAPINFOHEADER)){
        printf("read info fail\n");
    }

    nread = fread(buf, 1, sizeof(RGBQUAD)*256, fp);    
    if(nread != sizeof(RGBQUAD)*256){
        printf("read pic map!\n");
    }

    pics[num].nxp = info.biWidth;
    pics[num].nyp = info.biHeight;

    unsigned int value;
    RGBQUAD *rgbmap = buf;

    pics[num].data = (unsigned char *)malloc(header.bfSize - header.bfoffBits);

    nread = fread(pics[num].data, 1, header.bfSize - header.bfoffBits, fp); 
    if(nread != header.bfSize - header.bfoffBits){
        printf("read pic data failed!\n");
    }

    fclose(fp);

    return 0;
}

//x = x*c - y*s
//y = x*s + y*c

unsigned char gbuf[10*1024*1024] = {0,};
float cbuf[2000]={0,};
float sbuf[2000]={0,};
float scbuf[2000]={0,};

int main(int argc,char *argv[])
{
    int i,j;    

    int ang = atoi(argv[1]);
    float scal = (float)atoi(argv[2])/100;
    init_pics("./in.bmp",0);
    
    float c = cos(ang*3.1415926/180);
    float s = sin(ang*3.1415926/180);
    int tmp = sizeof(sbuf)/sizeof(sbuf[0])/2;


    for(i=0;i<sizeof(cbuf)/sizeof(cbuf[0]);i++){
        cbuf[i] = (i-tmp)*c;
        sbuf[i] = (i-tmp)*s;
        scbuf[i] = (i-tmp)*scal;
    }    

    printf("%f,%f,%f\n",c,s,scal);
    printf("%f,%f,%f\n",cbuf[50],sbuf[50],scal);
    printf("%f,%f,%f\n",cbuf[100],sbuf[100],scal);
    
    int w = pics[0].nxp;
    int h = pics[0].nyp;
    
    char *dat = malloc(w*h);
    memset(dat,0,w*h);
    
    int w2 = w/2;
    int h2 = h/2;    
    
    int jj = -1*w;
    unsigned char *ptr[1000]={0,}; 
    for(i=0; i < h; i++){
        ptr[i] = &pics[0].data[i*w];
    }

    //c=0;
    //s=1;

    for(i=0; i < h; i++){
        jj += w;
        for(j=0; j < w; j++){
            int xx,yy,xc,yc;
            
            //xx = j - w/2;
            //yy = h/2 - i;
            //xc = xx*c - yy*s;
            //yc = xx*s + yy*c;
            //xc = xc*scal + w/2;
            //yc = h/2 + yc*scal;
            
            xx = j - w2;
            yy = h2 - i;
            xc = cbuf[xx+tmp] - sbuf[yy+tmp];
            yc = sbuf[xx+tmp] + cbuf[yy+tmp];
            xc = scbuf[xc+tmp] + w2;
            yc = h2 + scbuf[yc+tmp];


            if(xc >= 0 && xc < w && yc >= 0 && yc < h){
                //dat[(h-yc)*w+xc] = pics[0].data[i*w+j];
                //dat[jj+j] = pics[0].data[(h-yc)*w+xc];
                dat[jj+j] = ptr[h-yc-1][xc];
            }
            
        }
    }
    
    #if 0
    for(i=1; i < h-1; i++){
        for(j=1; j < w-1; j++){
            if(dat[i*w+j] == 0){
                //dat[i*800+j] = (dat[i*800+j+1]+dat[i*800+j-1]+
                //dat[(i+1)*800+j]+dat[(i-1)*800+j])/4;
                dat[i*w+j] = dat[(i+1)*w+j];
            }            
        }
    }
    #endif
    
    FILE *outfp = fopen("out.bmp","w+");
    FILE *infp = fopen("in.bmp","r+");
    
    fread(gbuf,sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+sizeof(RGBQUAD)*256,1,infp);
    fwrite(gbuf,sizeof(BITMAPFILEHEADER)+sizeof(BITMAPINFOHEADER)+sizeof(RGBQUAD)*256,1,outfp);
    
    fwrite(dat,w*h,1,outfp);

    free(dat);
    fclose(infp);
    fclose(outfp);

}