我想学写bmp

[雕公]

因为看到网上有人用代码生成图片，所以在安卓机上安装了c4droid,然而生成的图片是ppm格式的，手机上根本无法显示，找了个文件格式转换的apk，结果才转了1张就要出钱才可以继续，自己不知道也懒得去破解，找了个在线转换的网站，却转个3*2的图片都变成1024*1024大小的，我一开始还以为转换失败，后来放大才发现转了，于是觉得不满意。去网上搜索了ppm2bmp的代码，是cpp的，发现其实把ppm2bmp.cpp在c4droid里运行就可以了，其他的都可以不用。于是我就有2个cpp了，一个用来生成ppm图片，主要是用什么算法生成rgb数值，另一个就是把ppm转换成bmp，这样就可以直接在手机上看了。然后干脆把2个cpp合成一个，让它生成ppm后再转换出对应的bmp。最后，我想研究下直接生成bmp图片，rgb值还是用生成ppm里的，只是把生成ppm图片的部分换成bmp格式。对于有一点水平的人来说估计也就1小时的事，查资料文件读写和bmp格式半小时，编辑调试源文件半小时。对于我来说，要1个月吗？从把生成ppm的cpp输入手机运行开始，到现在已经1星期了，成果就是：否定了文件转换apk和在线转换网站；留下一个同时生成ppm和bmp的cpp源文件（其实是先生成ppm,再转成bmp);知道了P3和P6格式ppm是怎么回事，但是无法自己独立生成，也就是不能自己用c实现。想法：研究bmp格式，学会编程生成bmp，然后把那个2合一的cpp弄成3合一即生成bmp,生成ppm,转ppm成bmp，这样我可以拆成3个cpp,一个功能一个。为了照顾自己的能力，ppm仅限于P6的，bmp仅限于24位不压缩的。还要多久才能把bmp弄出来呢？是为序。

[雕公]

看了自己之前发的帖子，有些东西现在自己根本不懂也没心情懂，像MTK里的c,PocketC之类，我当初竟然去学，现在根本没眼看那些代码，现在安卓手机确实功能不错，但是要成为桌上电脑那样的操作系统恐怕不可能，比如电脑上我要查看ppm图片，一大堆工具，而安卓上呢，只好自己先用c4droid写cpp应付，可见安卓开发大有作为啊！

[雕公]

#include <stdio.h>
#include <math.h>
//#include <iostream.h>
#include <stdlib.h>
int main(){
char head[54]={0x42,0x4d,0x36,0,0x30,0,0,0,0,0,0x36,0,0,0,0x28,0,0,0,0,0x4,0,0,0,0x4,0,0,1,0,0x18,0,0,0,0,0,0,0,0x30,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};
int corlor[3];
FILE *f=fopen("aaaa.bmp","wb");
// /*
for(int i=0;i<54;i++){
char c;c=(char)head[i];
fwrite(&c,1,1,f);
}
// */
fwrite(head,1,54,f);
for(int i=0;i<1024;i++){
for(int j=0;j<1024;j++){

int color[3];
color[0]=i+j;
color[1]=i*j;
color[2]=log(47)*i+j;
fwrite(color,1,3,f);
}
}
char c[54];
for(int i=0;i<54;i++){
;
fread(&c,1,1,f);
printf("%x ",c[i]);}
fclose(f);
}
这个勉强写出了bmp图片，但最后想读出文件头总觉得不太对头的样子。

[雕公]

所以才想学写bmp嘛，那个代码估计是Linux系统的人弄的，所以他生成的是ppm格式，而我在安卓手机上运行，才发现安卓无法显示ppm图像，所以才想学下写bmp.

[雕公]

#include <stdio.h>
#include <stdlib.h>
int writePPMHeader(FILE *f,char magic,int w,int h,int color){
if(f==NULL) exit(0);
if(magic=='A') fprintf(f,"P3\n");
else if(magic=='B') fprintf(f,"P6\n");
else exit(0);
fprintf(f,"%d %d\n%d\n",w,h,color);
return 0; }
int writePPMdatap3(FILE *f,unsigned char* img,int w,int h){
int i,j;
for(i=0;i<h;i++){
for(j=0;j<w;j++){
fprintf(f,"%d ",img[i*w*3+j*3]);
fprintf(f,"%d ",img[i*w*3+j*3+1]);
fprintf(f,"%d ",img[i*w*3+j*3+2]);}
fprintf(f,"\n"); } }
int writePPMdatap6(FILE *f,unsigned char* img,int w,int h){
int i,j;
for(i=0;i<h;i++)
for(j=0;j<w;j++)
fwrite(img,w*h,3,f); }
#define WIDTH 3
#define HIGHT 2
unsigned char img[WIDTH*HIGHT*3]={255,0,0,0,255,0,0,0,255,255,255,0,255,255,255,0,0,0};
int main(){
char *f1="testp3.ppm";
char *f2="testp6.ppm";
FILE *f3=fopen(f1,"w");
FILE *f6=fopen(f2,"w");
if(f3==NULL) exit(0);
if(f6==NULL) exit(0);
writePPMHeader(f3,'A',WIDTH,HIGHT,255);
writePPMdatap3(f3,img,3,2);
writePPMHeader(f6,'B',WIDTH,HIGHT,255);
writePPMdatap6(f6,img,3,2);
fclose(f3);fclose(f6);return 0;}

这个是生成ppm格式的，p3和p6都可以，实例是3*2的图片，太小了。

[赵4老师]

[日立奔腾浪潮微软松下联想]

这个...不知所云，你要生成图片，直接生成bmp好了，干嘛先ppm再转换。

[雕公]

我的bmp是偷懒固定了1024*1024大小的，那个文件头的char[54]是我读一个生成的1024*1024bmp文件得到的。我的那个直接生成bmp的源代码最后想读出生成的文件头，但是结果跟我自己写上去的文件头也就是char[54]似乎根本不相同，所以我说好像有什么不对，是fread用法不当没达到读出文件的前54字节的效果么？
下面附上网上找来的ppm2bmp.cpp，我的那个char[54]就是从它那里读出来的。
#include <iostream>
#include <string>
#include <cstdlib>
#include <fstream>
#include <iomanip>
using namespace std;

int main(int argc, char * argv[])
{
if (argc != 3)
{
cout << "please use this command like this: YuessMpeg2dec.exe src.ppm des.bmp" << endl;
}

string pFilename="MathPic.ppm"; //pFilename=argv[1];//ppm file name

string pOutFilename="MathPic.bmp";
ifstream ifile;//input ppm file
ofstream ofile;//output bmp file
ifile.open(pFilename,ios::binary);
if (!ifile) {
cout<<"open error!"<<endl;
}
ofile.open(pOutFilename,ios::binary);


string FileType;//file type-p6
int Width,Height,Count;//count is the number of pixels
int DataSize,HeadSize,FileSize;
ifile>>FileType>>Width>>Height>>Count;//get the para of the ppm file

cout<<FileType<<"!"<<Width<<"!"<<Height<<endl;

DataSize=Width*Height*3;//every pixel need 3 byte to store
HeadSize=0x36;//54 Bype
FileSize=HeadSize+DataSize;
Count=Width*Height;

cout<<Count<<endl;


unsigned char BmpHead[54];

for(int i=0;i<53;i++){
BmpHead[i]=0;
}

BmpHead[0]=0x42;
BmpHead[1]=0x4D;//type

cout<<FileSize<<endl;

int SizeNum=2;//the size of bmp file
while (FileSize!=0) {
BmpHead[SizeNum++]=FileSize % 256;
cout<<FileSize%256<<endl;
FileSize = FileSize / 256;
cout<<FileSize<<endl;
}

cout<<BmpHead[2]<<BmpHead[3]<<BmpHead[4]<<BmpHead[5]<<endl;

BmpHead[0x0A]=0x36;//data begin here
BmpHead[0x0E]=0x28;//size of bitmap information head

int FileWidthCount = 0x12;//width
while (Width != 0)
{

BmpHead[FileWidthCount++] = Width % 256;
Width = Width / 256;

}
int FileHeightCount = 0x16;//height
while (Height != 0)
{

BmpHead[FileHeightCount++] = Height % 256;
Height = Height / 256;
}

BmpHead[0x1A]=0x1;//device
BmpHead[0x1C]=0x18;//every pixel need 3 byte


int FileSizeCount=0x22;
while (DataSize!=0) {
BmpHead[FileSizeCount++]=DataSize%256;
DataSize=DataSize/256;
}
cout<<"下面是54字节的文件头信息："<<endl;
for (int i=0; i<54; i++)
{ofile << BmpHead[i];
int t=BmpHead[i]&255; cout<<i<<":"<<t<<endl;
}
cout<<"上面是54字节的文件头信息"<<endl;
char blue,green,red;

ifile.ignore();
cout<<"ok"<<endl;
for(int j = Count;j >=1 ;j--){
ifile.get(blue);
ifile.get(green);
ifile.get(red);

ofile<<hex;
ofile.put(red);
ofile.put(green);
ofile.put(blue);
}

ifile.close();
ofile.close();

return 0;

}
我改了下因为他是要在windows下变成xx.exe在命令行下执行的，我就直接把参数固定直接写进源文件，用c4droid在安卓手机上运行。
最后我想说一下一件很坑爹的事情，就是我的那个直接生成bmp的文件我把它的rgb函数的内容换成生成曼德波罗集的，结果呢，如果不在写像素rgb值的函数里加上命令行输出一些东西的语句，它就自己运行了一下就停止，生成的bmp根本达不到3M,明显是半途而费；而加了输出的语句后，因为命令行那里一直在输出，所以能一直运行到结束，然而因为我不以为然的觉得最后的fclose(f);这句没什么卵用，就没加上去，结果运行了最少3分钟才生成的bmp虽然大小是3M了，命令行也显示运行完了，我在main里最后加上输出done!的提示信息，次次都打不开生成的bmp,后来加上了fclose(f)后，运行3分钟多生成的图像终于可以看见了，确实是曼德波罗那个分形图像，不过跟网上给出的好像不太一样，颜色好像我的很单调，形状估计是因为bmp他的图形像素是从下到上从左到右输入的原因，感觉我的图像好像是要怎么翻转下，当然这没什么问题的。
下一部是想弄结构体填充头，这样就可以生成不固定1024宽高的图像了，估计有点水准的人也就不到30分钟——百度下相关信息，然后输入到元文件里去，我现在基础都不牢固，要先学习下基本的东西，理解后才能做到。

[雕公]

// NOTE: compile with g++ filename.cpp -std=c++11

#include <iostream>
#include <cmath>
#include <cstdlib>
#define DIM 1024
#define DM1 (DIM-1)
#define _sq(x) ((x)*(x)) // square
#define _cb(x) abs((x)*(x)*(x)) // absolute value of cube
#define _cr(x) (unsigned char)(pow((x),1.0/3.0)) // cube root

unsigned char GR(int,int);
unsigned char BL(int,int);
unsigned char RD(int,int);

unsigned char RD(int i,int j){
// YOUR CODE HERE
// return (char)(_sq(cos(atan2(j-512,i-512)/2))*255);
float x=0,y=0;int k;for(k=0;k++<256;){float a=x*x-y*y+(i-768.0)/512;y=2*x*y+(j-512.0)/512;x=a;if(x*x+y*y>4)break;}return log(k)*47;
}
unsigned char GR(int i,int j){
// YOUR CODE HERE
// return (char)(_sq(cos(atan2(j-512,i-512)/2))*255);
float x=0,y=0;int k;for(k=0;k++<256;){float a=x*x-y*y+(i-768.0)/512;y=2*x*y+(j-512.0)/512;x=a;if(x*x+y*y>4)break;}return log(k)*47;
}
unsigned char BL(int i,int j){
// YOUR CODE HERE
// return (char)(_sq(cos(atan2(j-512,i-512)/2))*255);
float x=0,y=0;int k;for(k=0;k++<256;){float a=x*x-y*y+(i-768.0)/512;y=2*x*y+(j-512.0)/512;x=a;if(x*x+y*y>4)break;}return 128-log(k)*23;
}

void pixel_write(int,int);
FILE *fp;
int main(){
fp = fopen("MathPic.ppm","wb");
fprintf(fp, "P6\n%d %d\n255\n", DIM, DIM);
for(int j=0;j<DIM;j++)
for(int i=0;i<DIM;i++)
pixel_write(i,j);
fclose(fp);
return 0;
}
void pixel_write(int i, int j){
static unsigned char color[3];
color[0] = RD(i,j)&255;
color[1] = GR(i,j)&255;
color[2] = BL(i,j)&255;
fwrite(color, 1, 3, fp);
}

一切的源头就是这文件，我的安卓机装上c4droid,gcc,sdl等之后，运行这文件，才发现生成的ppm竟然打不开，然后就是在网上找这种格式的信息，最后痛苦的发现，可以看这种图片的是windows下运行的程序，不是安卓；找了很久找到的一个apk可以转成bmp，却一次后就要冲钱，我只能放弃；找到一个上传后在线转换可以下载和网上预览的网站，用其间找ppm信息时发现的可以生成P3和P6格式ppm的源代码把P6格式的那个3*2大小的ppm上传转换，预览好像不堆，下载来一看3M多——似乎太大了？图片一浏览，根本没那3*2的方块，又重新弄了次再下载，还是3M还是找不到那颜色分明的3*2方块，最后放大才在中间看到了，太小不清楚。这才决定把ppm转成bmp。于是就有了ppm2bmp.cpp这个文件。
先运行生成ppm的cpp,再运行ppm2bmp.cpp，然后浏览生成的位图，觉得麻烦了点，于是就想办法把这2个cpp合2为一，一次就生成ppm和bmp。因为实际生成是先生成ppm然后再根据ppm生成bmp的，又想弄成直接把rgb函数生成的像素颜值直接灌进bmp中，于是就有了这个帖子。。
找bmp的信息，到现在稍为有点成果了。其实我看不懂ppm2bmp.cpp里那些填充bmp文件头的语句，才把填好的打印出来记住，先直接打印，有些字符无法显示，就不知道是什么。然后&255变成数字记下来，，最后才确定char[54],把数字又换成0x开头，又找了fwrite，才发现确定可以1次把头输入，不用for54次，就有了第一次可以用的直接生成bmp的c,不是cpp,cpp我很难看懂。最后是实验了生成曼德波罗，解决了一些小问题。这小问题难了我好久。
下一阶段用结构体来弄位图头，又要找书找资料学了。

[雕公]

今天最后一次回贴把曼德波罗的送上，最后一定要fclose，就它坑惨我了。
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
unsigned char r( int , int ) ;
unsigned char g( int , int ) ;
unsigned char b( int , int ) ;
int main() {
char head[54] = {0x42 ,0x4d ,0x36, 0 ,0x30,0 ,0 ,0, 0 ,0 ,0x36 ,0 ,0 ,0 ,0x28, 0 ,0, 0 ,0 ,0x4, 0 ,0 ,0 ,0x4 ,0 ,0 ,1 ,0 ,0x18,0 ,0 ,0 ,0 ,0 ,0 ,0, 0x30 ,0, 0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0 ,0, 0 ,0 ,0 ,0} ;
FILE *f ;
f = fopen( "xxxx.bmp" , "wb" ) ;
fwrite(head,1,54,f);
int d , i , j , k ;

for(i=0;i<1024;i++){
for(j=0;j<1024;j++)
{int color[3];
color[0]=r(i,j);
color[1]=g(i,j);
color[2]=b(i,j);
fwrite(color,1,3,f);
printf("%d%d",i%2,i%2?0:1);
}
}

printf("done!");
fclose(f);
}
unsigned char r( int i , int j ) {
float x=0,y=0;int k;for(k=0;k++<256;){float a=x*x-y*y+(i-768.0)/512;y=2*x*y+(j-512.0)/512;x=a;if(x*x+y*y>4)break;}return log(k)*47; }
unsigned char g( int i , int j ) {
float x=0,y=0;int k;for(k=0;k++<256;){float a=x*x-y*y+(i-768.0)/512;y=2*x*y+(j-512.0)/512;x=a;if(x*x+y*y>4)break;}return log(k)*47; }
unsigned char b (int i , int j ) {
float x=0,y=0;int k;for(k=0;k++<256;){float a=x*x-y*y+(i-768.0)/512;y=2*x*y+(j-512.0)/512;x=a;if(x*x+y*y>4)break;}return 128-log(k)*23; }


fwrite写rgb时一定要在后面printf点东西，不然手机不会运行到结束，因为计算这个分形的数据挺恐怖的，手机会罢工。我记得这个函数好像是（z=cz+c;z是复数a+bi而且a,b都不是0。c是实数且好像不是1？拼命跌代后的一堆复数z在复平面上画出来就变成分形图像。）显然这里不是用这种方法弄的，反正看不懂。
运行这个文件，你会得到分形图片，不过不是科研上可以无限放大的那种。

[日立奔腾浪潮微软松下联想]

文件头设置得不对，应该是 char head[54]={0x42,0x4d, 0x66,0,0,0, 0,0, 0,0, 0x36,0,0,0, 0x28,0,0,0, 0x4,0,0,0, 0x4,0,0,0, 1,0, 0x18,0, 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0, 0,0,0,0}; 见BITMAPFILEHEADER、BITMAPINFOHEADER两个结构的定义。