#include <stdio.h>
#include <stdlib.h>
#include <malloc.h>
#include <time.h>
#include <math.h>
#define SNOW_MATERIAL "WHITE MATTE"
#define WATER_MATERIAL "BLUE GLASS"
#define GRASS_MATERIAL "GREEN MATTE"
#define ROCK_MATERIAL "BEIGE MATTE"
#define RAND_FNC rand()
#define SIZE (unsigned long)((1 << ITER))
#define VARI (float)pow((double)VARI_CHANGE,(double)ITER)
char FTYPE;
short ITER=4;
short MATERIAL=1;
float VARI_CHANGE=2.00;
float SEA=0.40;
float SNOW=1.00;
short XSIZE=200;
short YSIZE=200;
short ZSIZE=200;
float x_pos;
float y_pos;
unsigned long count1;
unsigned long count2;
unsigned long buf_size;
unsigned long faces;
FILE *handle;
float huge *map;
float snow;
float water;
float altitude;
float min;
float max;
char filename[81];
void iterate(long,long,long,long,float,float huge *);
void write_material(FILE *,
unsigned long,
unsigned long, unsigned long,
float huge *,
float, float, float);
int write_faces(FILE *hand, char FTYPE)
{
if (FTYPE==0)
{
fprintf(hand,"Face list:\n");
for (count1=0; count1<SIZE-1; count1++)
for (count2=0; count2<SIZE-1; count2++)
{
fprintf(hand,"Face %lu: A:%lu B:%lu C:%lu AB:1 BC:1 CA:1\n",
count1*((SIZE-1)<<1)+count2*2,
count1*SIZE+count2,
count1*SIZE+count2+1,
(count1+1)*SIZE+count2);
if (MATERIAL)
write_material(hand,
count1*SIZE+count2,
count1*SIZE+count2+1,
(count1+1)*SIZE+count2,
map,water,snow,altitude);
fprintf(hand,"Face %lu: A:%lu B:%lu C:%lu AB:1 BC:1 CA:1\n",
count1*((SIZE-1)<<1)+count2*2+1,
count1*SIZE+count2+1,
(count1+1)*SIZE+count2+1,
(count1+1)*SIZE+count2);
if (MATERIAL)
write_material(hand,
count1*SIZE+count2+1,
(count1+1)*SIZE+count2+1,
(count1+1)*SIZE+count2,
map,water,snow,altitude);
}
}
return 0;
}
int write_verts(FILE *hand, char FTYPE)
{
if (FTYPE==0)
{
fprintf(hand,"Vertex list:\n");
for (count1=0; count1<SIZE; count1++)
for (count2=0; count2<SIZE; count2++)
fprintf(hand,"Vertex %lu: X: %f Y: %f Z: %f\n",
count1*SIZE+count2,
(float)(count2*XSIZE)/SIZE,(float)(count1*YSIZE)/SIZE,
(float)(map[count1*SIZE+count2]*ZSIZE)/altitude);
}
if (FTYPE==1)
{
fprintf(hand,"\tLL_DATA\t");
for (count1=0; count1<SIZE*8; count1++)
{
fprintf(hand,"db ");
for (count2=0; count2<SIZE/8-1; count2++)
fprintf(hand,"%2u,",(unsigned char)((float)((map[count1*(SIZE/8)+count2]-water)*ZSIZE)/altitude));
fprintf(hand,"%2u\n\t\t",(unsigned char)((float)((map[count1*(SIZE/8)+count2]-water)*ZSIZE)/altitude));
}
}
return 0;
}
FILE *init_file(char *fname,char FTYPE)
{
FILE *hand;
if ((hand=fopen(fname,"wt"))==0)
return 0;
if (FTYPE==0)
{
fprintf(hand,"Ambient light color: Red=0.3 Green=0.3 Blue=0.3\n\n");
fprintf(hand,"Named object: \"Fractal Mtn\"\n");
fprintf(hand,"Tri-mesh, Vertices: %lu Faces: %lu\n",buf_size,faces);
}
if (FTYPE==1)
{
fprintf(hand,";;\n");
fprintf(hand,";; Lord Logics Landscape Mesh File\n");
fprintf(hand,";;\n");
fprintf(hand,".fardata\tLL_MESH\n");
}
return hand;
}
main(int argc, char *argv[])
{
strcpy(filename,"");
for (count1=1; count1<argc; count1++)
{
if ((argv[count1][0]!='\\') && (argv[count1][0]!='/') && (argv[count1][0]!='-'))
{
printf("3D Fractal Landscape Creator v1.5 [12-08-92]\n");
printf("Copyright (C) 1992-93 James P. Ketrenos\n\n");
printf("Valid Command Line Options:\n");
printf(" -nFileame.ext [default: FRACTAL.??? ]\n");
printf(" -iIterations [default: 4 ]\n");
printf(" -mMaterials [default: 1 (ON) ]\n");
printf(" -wWaterLevel [default: 0.40 ]\n");
printf(" -sSnowLevel [default: 1.00 ]\n");
printf(" -vVariable [default: 2.00 ]\n");
printf(" -xSize [default: 200 ]\n");
printf(" -ySize [default: 200 ]\n");
printf(" -zSize [default: 200 ]\n");
printf(" -tFileType [default: 3D Studio .ASC]\n\n");
printf("Filename.ext consits of any valid DOS filename. Up to 80 characters\n");
printf(" may be used to include path etc. The file written is a\n");
printf(" 3D Studio ASCII Text file and can get quite large.\n\n");
printf("Iterations is an integer that can range from 1 to 15. This is\n");
printf(" used to determine the number of splits performed.\n\n");
printf("Materials is used to determine whether or not materials will\n");
printf(" assigned to faces or not. 1=ON 0=OFF.\n\n");
printf("<< Press any key to continue >>");
if (getch()==27)
{
printf("\r \r");
exit(0);
}
printf("\r \r");
printf("WaterLevel is the percentage of the landform to make below the SEA\n");
printf(" LEVEL. All land below this will be leveled to SEA. It\n");
printf(" is entered as a decimal. (1.00=All .. 0.25=1/4 etc.)\n\n");
printf("SnowLevel is the percentage of the landform to make below the SNOW\n");
printf(" LEVEL. All land above this will be WHITE. This makes no\n");
printf(" plateau adjustments. (1.00=All .. 0.25=1/4 etc.)\n\n");
printf("Variable is the variable degredation of variance. Basically, what\n");
printf(" this does is determine how much the variance will decrease\n");
printf(" with every land split. Variable>1 will model \"life\".\n");
printf(" 0<Variable<=1 will resemble \"chaos\". Decent values can\n");
printf(" be found between 2.00 and 5.00.\n\n");
printf("Size is the length of the box containing the landscape along\n");
printf(" perspective axis (X, Y, or Z). It is better to keep the\n");
printf(" values semi-large (in the hundreds) to keep from having\n");
printf(" distortions due to round-off errors.\n\n");
printf("<< Press any key to continue >>");
if (getch()==27)
{
printf("\r \r");
exit(0);
}
printf("\r \r");
printf("Type is the filetype to create. Types are:\n");
printf(" 0 - 3D Studio .ASC [default]\n");
printf(" 1 - Lord Logics Landscape Data Mesh File .INC\n");
printf(" Note: Iterations = 8, ZSize = 256\n");
printf(" XSize = 256, and YSize = 256\n");
printf(" Materials are OFF,\n");
printf(" Snow is OFF Water is at 0.35\n\n");
printf("Any questions, comments, suggestions, donations, or information is welcome.\n");
printf("Mail to:\n\n");
printf(" James P. Ketrenos\n");
printf(" 885 S.W. 84th Ct.\n");
printf(" Portland, OR 97225\n\n");
printf("Or e-mail at: ketrenoj@ucs.orst.edu up till 06-01-93.\n\n");
printf("If you would like to see the data output to a different file format\n");
printf("besides the 3D Studio ASCII file, then please send me a note saying\n");
printf("so along with a description of the file format and how it is made up.\n");
printf("Primarily I need to know how to define vertices and triangle faces.\n");
exit(0);
}
switch(argv[count1][1]<'a' ? (argv[count1][1]-'A'+'a') : argv[count1][1])
{
case 'n':strcpy(filename,&argv[count1][2]); break;
case 'i':ITER=atoi(&argv[count1][2]); break;
case 'm':MATERIAL=atoi(&argv[count1][2]); break;
case 'w':SEA=atof(&argv[count1][2]); break;
case 's':SNOW=atof(&argv[count1][2]); break;
case 'v':VARI_CHANGE=atof(&argv[count1][2]); break;
case 'x':XSIZE=atoi(&argv[count1][2]); break;
case 'y':YSIZE=atoi(&argv[count1][2]); break;
case 'z':ZSIZE=atoi(&argv[count1][2]); break;
case 't':FTYPE=atoi(&argv[count1][2]); break;
default:argv[count1][0]='?';
count1--;
}
}
if (filename[0]==0)
switch(FTYPE)
{
case 0:strcpy(filename,"FRACTAL.ASC"); break;
case 1:strcpy(filename,"FRACTAL.INC"); break;
}
if (FTYPE==1)
{
MATERIAL=0;
SEA=0.35;
SNOW=1.00;
ITER=8;
ZSIZE=256;
XSIZE=YSIZE=256;
}
printf("-=[ Fractal Landscape Creator ]=-\n");
printf(" -=[ Version 1.5 ]=-\n\n");
printf("[ - Fractal Parameters - ]\n");
printf("[ Matrix Splits : %-5u ]",ITER);
if ((ITER>0) && (ITER<16))
printf("\n");
else
{
printf(" Invalid range!\n");
exit(-1);
}
printf("[ Matrix Base : %-4lu by %-4lu ]\n",SIZE,SIZE);
printf("[ Matrix Sea Level : %-4.2f ]",SEA);
if ((SEA>=0.0) && (SEA<=1.0))
printf("\n");
else
{
printf(" Invalid range!\n");
exit(-1);
}
printf("[ Matrix Snow Level : %-4.2f ]",SNOW);
if ((SNOW>=0.0) && (SNOW<=1.0))
printf("\n");
else
{
printf(" Invalid range!\n");
exit(-1);
}
printf("[ Variance Degredation : %-5.2f ]",VARI_CHANGE);
if (VARI_CHANGE>0.0)
printf("\n");
else
{
printf(" Invalid range!\n");
exit(-1);
}
printf("[ Maximum Variance : %-7.2f ]\n",VARI);
printf("[ Surface Vertices : %-8lu ]\n",buf_size=SIZE*SIZE);
if (FTYPE!=1)
printf("[ Surface Faces : %-8lu ]\n",faces=((SIZE-1)*(SIZE-1)*2));
printf("[ Cube Containment Size :(%-4ux%-4ux%-4u) ]",XSIZE,YSIZE,ZSIZE);
if ((XSIZE<=0.0) || (YSIZE<=0.0) || (ZSIZE<=0.0))
{
printf(" Invalid range!\n");
exit(-1);
}
printf("\n\n");
printf("[ - Fractal Creation Progress - ]\n");
printf("[ Initializing Matrix : WAIT. ]");
if ((map=(float huge *)halloc(buf_size,sizeof(map[0])))==0)
{
printf(" Error!\n");
printf(" - Insufficient Memory -\n");
exit(-1);
}
for (count1=0; count1<buf_size; count1++)
map[count1]=(float)987;
/* Initialize the four corners */
srand((unsigned)time(0)); // Initialize the RNG
map[0*SIZE+0]=(float)((float)RAND_FNC/RAND_MAX) * VARI;
map[0*SIZE+(SIZE-1)]=(float)((float)RAND_FNC/RAND_MAX) * VARI;
map[(SIZE-1)*SIZE+(SIZE-1)]=(float)((float)RAND_FNC/RAND_MAX) * VARI;
map[(SIZE-1)*SIZE+0]=(float)((float)RAND_FNC/RAND_MAX) * VARI;
printf("\r[ Initializing Matrix : OK. ]\n");
printf("[ Creating Landscape : WAIT. ]");
iterate(0,0,SIZE-1,SIZE-1,(float)VARI,map);
min=map[0];
max=map[0];
for (count1=1; count1<buf_size; count1++)
{
if (map[count1]<min) min=map[count1];
if (map[count1]>max) max=map[count1];
}
altitude=max-min;
printf("\r[ Creating Landscape : OK. ]\n");
printf("[ Adding Elements : WAIT. ]");
water=(float)altitude*SEA+min;
snow=(float)altitude*SNOW+min;
if (altitude==0.0)
{
printf(" Error!\n");
printf("No Landmass present.\n");
free(map);
exit(-2);
}
for (count1=0; count1<buf_size; count1++)
if (map[count1]<=water)
map[count1]=water;
printf("\r[ Adding Elements : OK. ]\n\n");
printf("[ Initializing Mesh File: WAIT. ]");
if ((handle=init_file(filename,FTYPE))==0)
{
printf("\r[ Initializing Mesh File: ERROR!! ]\n");
free(map);
exit(-1);
}
if (strlen(filename)>12)
printf("\r[ Initializing Mesh File: ...%-12s ]\n",&filename[strlen(filename)-12]);
else
printf("\r[ Initializing Mesh File: %-12s ]\n",filename);
printf("[ Writing Vertices : WAIT. ]");
if (write_verts(handle, FTYPE)==0)
printf("\r[ Writing Vertices : OK. ]\n");
else
{
printf("\r[ Writing Vertices : ERROR! ]\n");
close(handle);
free(map);
exit(-1);
}
printf("[ Writing Faces : WAIT. ]");
if (write_faces(handle, FTYPE)==0)
printf("\r[ Writing Faces : OK. ]\n");
else
{
printf("\r[ Writing Faces : ERROR! ]\n");
close(handle);
free(map);
exit(-1);
}
printf("[ Closing Mesh File : WAIT. ]");
fclose(handle);
printf("\r[ Closing Mesh File : OK. ]\n");
free(map);
}
void iterate(long c1, long c2, long c3, long c4, float variance, float huge *map)
{
long longi;
long latit;
float height;
float change;
if ((c3 - (c1 + 1)) & 1)
change=0.5;
else
change=0.0;
longi=((c1 + c3) / 2);
latit=((c2 + c4) / 2);
if (map[c2*SIZE+longi]==987)
{
height=((map[c2*SIZE+c1] + map[c2*SIZE+c3]) / 2);
height+=(float)((float)RAND_FNC/RAND_MAX) * variance - (variance * 0.5);
map[c2*SIZE+longi]=height+((map[c2*SIZE+c1]-height)*change)/((longi+change)-c1);
}
if (map[c4*SIZE+longi]==987)
{
height=((map[c4*SIZE+c1] + map[c4*SIZE+c3]) / 2);
height+=(float)((float)RAND_FNC/RAND_MAX) * variance - (variance * 0.5);
map[c4*SIZE+longi]=height+((map[c4*SIZE+c1]-height)*change)/((longi+change)-c1);
}
if (map[latit*SIZE+c1]==987)
{
height=((map[c2*SIZE+c1] + map[c4*SIZE+c1]) / 2);
height+=(float)((float)RAND_FNC/RAND_MAX) * variance - (variance * 0.5);
map[latit*SIZE+c1]=height+((map[c2*SIZE+c1]-height)*change)/((latit+change)-c2);
}
if (map[latit*SIZE+c3]==987)
{
height=((map[c2*SIZE+c3] + map[c4*SIZE+c3]) / 2);
height+=(float)((float)RAND_FNC/RAND_MAX) * variance - (variance * 0.5);
map[latit*SIZE+c3]=height+((map[c2*SIZE+c3]-height)*change)/((latit+change)-c2);
}
height=((map[c2*SIZE+c1] + map[c2*SIZE+c3] + map[c4*SIZE+c1] + map[c4*SIZE+c3] +
map[c2*SIZE+longi] + map[c4*SIZE+longi] + map[latit*SIZE+c1] + map[latit*SIZE+c3]) / 8);
height+=(float)((float)RAND_FNC/RAND_MAX) * variance - (variance * 0.5);
if (change==0.5)
map[latit*SIZE+longi]=height+((map[c2*SIZE+c1]-height)*0.7071)/(1.4142*((longi+change)-c1));
else
map[latit*SIZE+longi]=height;
if ((longi==c1) || (longi==c3) || (latit==c2) || (latit==c4)) return;
iterate(c1, c2, longi, latit, (float)((longi-c1)*variance)/(((c3-c1)*VARI_CHANGE)/2), map);
iterate(longi, c2, c3, latit, (float)((c3-longi)*variance)/(((c3-c1)*VARI_CHANGE)/2), map);
iterate(c1, latit, longi, c4, (float)((latit-c2)*variance)/(((c4-c2)*VARI_CHANGE)/2), map);
iterate(longi, latit, c3, c4, (float)((c4-latit)*variance)/(((c4-c2)*VARI_CHANGE)/2), map);
return;
}
void write_material(FILE *handle,unsigned long c1, unsigned long c2, unsigned long c3,float huge *map,
float water, float snow, float altitude)
{
float average;
average=(map[c1] + map[c2] + map[c3]) / 3;
if (average >= snow)
fprintf(handle,"Material:\"%s\"\n",SNOW_MATERIAL);
else
if ((map[c1]==water) && (map[c2]==water) && (map[c3]==water))
fprintf(handle,"Material:\"%s\"\n",WATER_MATERIAL);
else
{
if ( ((fabs(map[c1]-average)*(ZSIZE/altitude)*2*SIZE/YSIZE)>4.5) ||
((fabs(map[c2]-average)*(ZSIZE/altitude)*2*SIZE/YSIZE)>4.5) ||
((fabs(map[c3]-average)*(ZSIZE/altitude)*2*SIZE/XSIZE)>4.5))
fprintf(handle,"Material:\"%s\"\n",ROCK_MATERIAL);
else
fprintf(handle,"Material:\"%s\"\n",GRASS_MATERIAL);
}
}
