[49] | 1 | #include "supmorph.hpp" |
---|
| 2 | #include "specs.hpp" |
---|
| 3 | #include "timage.hpp" |
---|
| 4 | #include "timing.hpp" |
---|
| 5 | #include "filter.hpp" |
---|
| 6 | #include "video.hpp" |
---|
| 7 | #include "jrand.hpp" |
---|
| 8 | |
---|
| 9 | #define p_swap(x,y) { x^=y; y^=x; x^=y; } |
---|
| 10 | #define p_dist(x1,y1,x2,y2) (((int)(x1)-(int)x2)*((int)(x1)-(int)x2)+ \ |
---|
| 11 | ((int)(y1)-(int)y2)*((int)(y1)-(int)y2)) |
---|
| 12 | |
---|
| 13 | super_morph::super_morph(trans_image *hint1, trans_image *hint2, |
---|
| 14 | int aneal_steps, void (*stat_fun)(int)) |
---|
| 15 | { |
---|
| 16 | int x,y,w1=hint1->width(), |
---|
| 17 | h1=hint1->height(), |
---|
| 18 | w2=hint2->width(), |
---|
| 19 | h2=hint2->height(); |
---|
| 20 | if (w1>w2) w=w1; else w=w2; |
---|
| 21 | if (h1>h2) h=h1; else h=h2; |
---|
| 22 | unsigned char *dp; |
---|
| 23 | |
---|
| 24 | /************************ count up the number of hints ****************************/ |
---|
| 25 | unsigned short hints1[256],hints2[256]; |
---|
| 26 | memset(hints1,0,256*2); |
---|
| 27 | memset(hints2,0,256*2); |
---|
| 28 | |
---|
| 29 | dp=hint1->t_data(); |
---|
| 30 | for (y=0;y<h1;y++) |
---|
| 31 | { |
---|
| 32 | x=0; |
---|
| 33 | while (x<w1) |
---|
| 34 | { |
---|
| 35 | x+=*dp; // skip over space |
---|
| 36 | dp++; |
---|
| 37 | if (x<w1) |
---|
| 38 | { |
---|
| 39 | int rl=*(dp++); |
---|
| 40 | while (rl--) { hints1[*(dp++)]++; x++; } |
---|
| 41 | } |
---|
| 42 | } |
---|
| 43 | } |
---|
| 44 | |
---|
| 45 | // hint2 image2 |
---|
| 46 | dp=hint2->t_data(); |
---|
| 47 | for (y=0;y<h2;y++) |
---|
| 48 | { |
---|
| 49 | x=0; |
---|
| 50 | while (x<w2) |
---|
| 51 | { |
---|
| 52 | x+=*dp; // skip over space |
---|
| 53 | dp++; |
---|
| 54 | if (x<w2) |
---|
| 55 | { |
---|
| 56 | int rl=*(dp++); |
---|
| 57 | while (rl--) { hints2[*(dp++)]++; x++; } |
---|
| 58 | } |
---|
| 59 | } |
---|
| 60 | } |
---|
| 61 | |
---|
| 62 | |
---|
| 63 | /****************** now sum up hints and alloc memory for movers *************/ |
---|
| 64 | unsigned short start1[256],start2[256], |
---|
| 65 | totals[256], |
---|
| 66 | start=0,*h1p=hints1,*h2p=hints2; |
---|
| 67 | unsigned char hint_color[256],total_hints=0; |
---|
| 68 | |
---|
| 69 | for (y=0;y<256;y++,h1p++,h2p++) |
---|
| 70 | if (*h1p) |
---|
| 71 | { |
---|
| 72 | if (*h2p==0) |
---|
| 73 | { |
---|
| 74 | t=0; |
---|
| 75 | return ; |
---|
| 76 | } |
---|
| 77 | start1[y]=start2[y]=start; // specify start of hint range |
---|
| 78 | if (*h1p>*h2p) |
---|
| 79 | totals[y]=*h1p; |
---|
| 80 | else totals[y]=*h2p; |
---|
| 81 | start+=totals[y]; |
---|
| 82 | hint_color[total_hints++]=y; |
---|
| 83 | } |
---|
| 84 | |
---|
| 85 | t=start; |
---|
| 86 | movers=(unsigned char *)jmalloc(t*4,"morph movers"); |
---|
| 87 | |
---|
| 88 | |
---|
| 89 | /**************** Now scan the images again setup hints *********************/ |
---|
| 90 | dp=hint1->t_data(); |
---|
| 91 | for (y=0;y<h1;y++) |
---|
| 92 | { |
---|
| 93 | x=0; |
---|
| 94 | while (x<w1) |
---|
| 95 | { |
---|
| 96 | x+=*dp; // skip over space |
---|
| 97 | dp++; |
---|
| 98 | if (x<w1) |
---|
| 99 | { |
---|
| 100 | int rl=*(dp++); |
---|
| 101 | while (rl--) |
---|
| 102 | { |
---|
| 103 | int maddr=(start1[*(dp++)]++)*4; |
---|
| 104 | movers[(maddr++)]=x; |
---|
| 105 | movers[maddr]=y; |
---|
| 106 | x++; |
---|
| 107 | } |
---|
| 108 | } |
---|
| 109 | } |
---|
| 110 | } |
---|
| 111 | |
---|
| 112 | dp=hint2->t_data(); |
---|
| 113 | for (y=0;y<h2;y++) |
---|
| 114 | { |
---|
| 115 | x=0; |
---|
| 116 | while (x<w2) |
---|
| 117 | { |
---|
| 118 | x+=*dp; // skip over space |
---|
| 119 | dp++; |
---|
| 120 | if (x<w2) |
---|
| 121 | { |
---|
| 122 | int rl=*(dp++); |
---|
| 123 | while (rl--) |
---|
| 124 | { |
---|
| 125 | int maddr=(start2[*(dp++)]++)*4+2; |
---|
| 126 | movers[(maddr++)]=x; |
---|
| 127 | movers[maddr]=y; |
---|
| 128 | x++; |
---|
| 129 | } |
---|
| 130 | } |
---|
| 131 | } |
---|
| 132 | } |
---|
| 133 | |
---|
| 134 | /********* if hint sizes don't match duplicate the smaller until sizes are equal **********/ |
---|
| 135 | for (start=0,x=0;x<total_hints;x++) |
---|
| 136 | { |
---|
| 137 | y=hint_color[x]; |
---|
| 138 | int dups; |
---|
| 139 | for (dp=movers+start1[y]*4,dups=totals[y]-hints1[y];dups;dups--) |
---|
| 140 | { |
---|
| 141 | *dp=*(dp-4); dp++; // copy previous x,y position |
---|
| 142 | *dp=*(dp-4); dp++; |
---|
| 143 | dp+=2; |
---|
| 144 | } |
---|
| 145 | start1[y]-=2*totals[y]-hints1[y]; // set the start back to the begining of hint range |
---|
| 146 | } |
---|
| 147 | |
---|
| 148 | for (start=0,x=0;x<total_hints;x++) |
---|
| 149 | { |
---|
| 150 | y=hint_color[x]; |
---|
| 151 | int dups; |
---|
| 152 | for (dp=movers+start2[y]*4+2,dups=totals[y]-hints2[y];dups;dups--) |
---|
| 153 | { |
---|
| 154 | *dp=*(dp-4); dp++; // copy previous x,y position |
---|
| 155 | *dp=*(dp-4); dp++; |
---|
| 156 | dp+=2; |
---|
| 157 | } |
---|
| 158 | start2[y]-=hints2[y]; // set the start back to the begining of hint range |
---|
| 159 | } |
---|
| 160 | |
---|
| 161 | |
---|
| 162 | /******* Now apply simulated annealing to solve for a smaller total distance ********/ |
---|
| 163 | int rand_on=0; |
---|
| 164 | for (y=0;y<aneal_steps;y++) |
---|
| 165 | { |
---|
| 166 | if (stat_fun) |
---|
| 167 | stat_fun(y); |
---|
| 168 | dp=movers; |
---|
| 169 | for (x=0;x<total_hints;x++) |
---|
| 170 | { |
---|
| 171 | int hc=hint_color[x]; |
---|
| 172 | int a,z=totals[hc]; |
---|
| 173 | unsigned char *range_start=dp; |
---|
| 174 | for (a=0;a<z;a++,dp+=4) |
---|
| 175 | { |
---|
| 176 | unsigned char *swap=range_start+(rtable[((rand_on++)&(RAND_TABLE_SIZE-1))]%z)*4; |
---|
| 177 | int d_old=p_dist(dp[0],dp[1],dp[2],dp[3])+p_dist(swap[0],swap[1],swap[2],swap[3]); |
---|
| 178 | int d_new=p_dist(dp[0],dp[1],swap[2],swap[3])+p_dist(swap[0],swap[1],dp[2],dp[3]); |
---|
| 179 | if (d_new<d_old) |
---|
| 180 | { |
---|
| 181 | unsigned char s; |
---|
| 182 | s=swap[2]; swap[2]=dp[2]; dp[2]=s; |
---|
| 183 | s=swap[3]; swap[3]=dp[3]; dp[3]=s; |
---|
| 184 | } |
---|
| 185 | } |
---|
| 186 | } |
---|
| 187 | } |
---|
| 188 | } |
---|
| 189 | |
---|
| 190 | smorph_player::smorph_player(super_morph *m, palette *pal, image *i1, image *i2, int frames, int dir) |
---|
| 191 | { |
---|
| 192 | int i,x1,y1,x2,y2; |
---|
| 193 | unsigned char *d=m->movers,*paddr=(unsigned char *)pal->addr(),*pa; |
---|
| 194 | stepper *p; |
---|
| 195 | p=steps=(stepper *)jmalloc(sizeof(stepper)*m->t,"smorph steps"); |
---|
| 196 | f_left=frames; |
---|
| 197 | frames--; |
---|
| 198 | t=m->t; |
---|
| 199 | w=m->w; h=m->h; |
---|
| 200 | |
---|
| 201 | for (i=0;i<t;i++,p++) |
---|
| 202 | { |
---|
| 203 | x1=*(d++); |
---|
| 204 | y1=*(d++); |
---|
| 205 | x2=*(d++); |
---|
| 206 | y2=*(d++); |
---|
| 207 | |
---|
| 208 | unsigned char r1,g1,b1,r2,g2,b2; |
---|
| 209 | pa=paddr+(int)(*(i1->scan_line(y1)+x1))*3; |
---|
| 210 | r1=*(pa++); |
---|
| 211 | g1=*(pa++); |
---|
| 212 | b1=*(pa++); |
---|
| 213 | |
---|
| 214 | pa=paddr+(int)(*(i2->scan_line(y2)+x2))*3; |
---|
| 215 | r2=*(pa++); |
---|
| 216 | g2=*(pa++); |
---|
| 217 | b2=*(pa++); |
---|
| 218 | |
---|
| 219 | p->r=r1<<16; |
---|
| 220 | p->g=g1<<16; |
---|
| 221 | p->b=b1<<16; |
---|
| 222 | |
---|
| 223 | p->dr=(long)(((int)r2-(int)r1)<<16)/frames; |
---|
| 224 | p->dg=(long)(((int)g2-(int)g1)<<16)/frames; |
---|
| 225 | p->db=(long)(((int)b2-(int)b1)<<16)/frames; |
---|
| 226 | |
---|
| 227 | if (dir<0) |
---|
| 228 | { |
---|
| 229 | x1=w-x1-1; |
---|
| 230 | x2=w-x2-1; |
---|
| 231 | } |
---|
| 232 | p->dx=((x2-x1)<<16)/frames; |
---|
| 233 | p->dy=((y2-y1)<<16)/frames; |
---|
| 234 | p->x=x1<<16; |
---|
| 235 | p->y=y1<<16; |
---|
| 236 | } |
---|
| 237 | |
---|
| 238 | hole=(unsigned char *)jmalloc(w*h,"hole image"); |
---|
| 239 | } |
---|
| 240 | |
---|
| 241 | |
---|
| 242 | |
---|
| 243 | |
---|
| 244 | |
---|
| 245 | int smorph_player::show(image *screen, int x, int y, color_filter *fil, palette *pal, |
---|
| 246 | int blur_threshold) |
---|
| 247 | { |
---|
| 248 | if (f_left) |
---|
| 249 | { |
---|
| 250 | int i,px,py,ix,iy; |
---|
| 251 | short x1,y1,x2,y2; |
---|
| 252 | screen->get_clip(x1,y1,x2,y2); |
---|
| 253 | screen->add_dirty(x,y,x+w-1,y+h-1); |
---|
| 254 | stepper *ss; |
---|
| 255 | memset(hole,0,w*h); |
---|
| 256 | unsigned char *paddr=(unsigned char *)pal->addr(); |
---|
| 257 | for (ss=steps,i=0;i<t;i++,ss++) |
---|
| 258 | { |
---|
| 259 | ix=(ss->x>>(16)); |
---|
| 260 | iy=(ss->y>>(16)); |
---|
| 261 | px=ix+x; |
---|
| 262 | py=iy+y; |
---|
| 263 | if (px>=x1 && px<=x2 && py>=y1 && py<=y2) |
---|
| 264 | { |
---|
| 265 | hole[ix+iy*w]=*(screen->scan_line(py)+px)=fil->lookup_color(ss->r>>(19), |
---|
| 266 | ss->g>>(19), |
---|
| 267 | ss->b>>(19)); |
---|
| 268 | } |
---|
| 269 | ss->x+=ss->dx; |
---|
| 270 | ss->y+=ss->dy; |
---|
| 271 | ss->r+=ss->dr; |
---|
| 272 | ss->g+=ss->dg; |
---|
| 273 | ss->b+=ss->db; |
---|
| 274 | } |
---|
| 275 | f_left--; |
---|
| 276 | if (!f_left) // skip hole fills and smoothing on last frame |
---|
| 277 | return 1; |
---|
| 278 | |
---|
| 279 | unsigned char *ll=hole+1,*tl=hole+w+1,*nl=hole+w*2+1; |
---|
| 280 | for (iy=1;iy<h-1;iy++) // now scan the for holes to fill |
---|
| 281 | { |
---|
| 282 | for (ix=1;ix<w-1;ix++,ll++,tl++,nl++) |
---|
| 283 | { |
---|
| 284 | if (x+ix>=x1 && x+ix<=x2 && y+iy>=y1 && y+iy<=y2) |
---|
| 285 | { |
---|
| 286 | int t=0; |
---|
| 287 | unsigned char *pa; |
---|
| 288 | int r=0,g=0,b=0; |
---|
| 289 | /* if (*(tl-1)) t++; |
---|
| 290 | if (*(tl+1)) t++; |
---|
| 291 | if (*ll) t++; |
---|
| 292 | if (*nl) t++;*/ |
---|
| 293 | |
---|
| 294 | if (*(tl-1)) { t++; pa=paddr+(*(tl-1))*3; r+=*(pa++); g+=*(pa++); b+=*(pa++); } |
---|
| 295 | if (*(tl+1)) { t++; pa=paddr+(*(tl+1))*3; r+=*(pa++); g+=*(pa++); b+=*(pa++); } |
---|
| 296 | if (*(ll)) { t++; pa=paddr+(*ll)*3; r+=*(pa++); g+=*(pa++); b+=*(pa++); } |
---|
| 297 | if (*(nl)) { t++; pa=paddr+(*nl)*3; r+=*(pa++); g+=*(pa++); b+=*(pa++); } |
---|
| 298 | |
---|
| 299 | if (*tl) |
---|
| 300 | { |
---|
| 301 | if (t) |
---|
| 302 | { |
---|
| 303 | pa=paddr+(*tl)*3; |
---|
| 304 | r/=t; g/=t; b/=t; |
---|
| 305 | int dist=((int)(*pa)-r)*((int)(*pa)-r); pa++; |
---|
| 306 | dist+=((int)(*pa)-g)*((int)(*pa)-g); pa++; |
---|
| 307 | dist+=((int)(*pa)-b)*((int)(*pa)-b); |
---|
| 308 | if (dist>blur_threshold) |
---|
| 309 | *(tl)=*(screen->scan_line(y+iy)+x+ix)=fil->lookup_color(r>>3,g>>3,b>>3); |
---|
| 310 | } else *(tl)=*(screen->scan_line(y+iy)+x+ix)=0; // kill single pixels |
---|
| 311 | } |
---|
| 312 | else if (t>=3) |
---|
| 313 | *(tl)=*(screen->scan_line(y+iy)+x+ix)=fil->lookup_color((r/t)>>3,(g/t)>>3,(b/t)>>3); |
---|
| 314 | } |
---|
| 315 | } |
---|
| 316 | ll+=2; |
---|
| 317 | tl+=2; |
---|
| 318 | nl+=2; |
---|
| 319 | } |
---|
| 320 | return 1; |
---|
| 321 | } else return 0; |
---|
| 322 | } |
---|
| 323 | |
---|
| 324 | |
---|
| 325 | |
---|
| 326 | |
---|
| 327 | /*void free_up_memory() { printf("you're screwed\n"); } |
---|
| 328 | |
---|
| 329 | main(int argc, char **argv) |
---|
| 330 | { |
---|
| 331 | image_init(); |
---|
| 332 | jrand_init(); |
---|
| 333 | FILE *fp=fopen("art/mrphmask.spe","rb"); |
---|
| 334 | spec_directory sd(fp); |
---|
| 335 | image *h1=new image(sd.find("20 h"),fp), |
---|
| 336 | *h2=new image(sd.find("1h"),fp), |
---|
| 337 | *i1=new image(sd.find("20"),fp), |
---|
| 338 | *i2=new image(sd.find("1"),fp); |
---|
| 339 | palette *pal=new palette(sd.find(SPEC_PALETTE),fp); |
---|
| 340 | color_filter *fil=new color_filter(sd.find(SPEC_COLOR_TABLE),fp); |
---|
| 341 | |
---|
| 342 | int steps=atoi(argv[1]); |
---|
| 343 | if (steps<2) steps=50; |
---|
| 344 | trans_image *hh1=new trans_image(h1,"hint1"),*hh2=new trans_image(h2,"hint2"); |
---|
| 345 | |
---|
| 346 | time_marker time1; |
---|
| 347 | super_morph sm(hh1,hh2,steps); |
---|
| 348 | int frames=atoi(argv[2]); |
---|
| 349 | if (frames<2) frames=16; |
---|
| 350 | smorph_player sp(&sm,pal,i1,i2,frames,-1); |
---|
| 351 | |
---|
| 352 | time_marker time2; |
---|
| 353 | printf("time = %lf\n",time2.diff_time(&time1)); |
---|
| 354 | |
---|
| 355 | set_mode(19,argc,argv); |
---|
| 356 | pal->load(); |
---|
| 357 | i1->put_image(screen,30,30); |
---|
| 358 | update_dirty(screen); |
---|
| 359 | sleep(2); |
---|
| 360 | while (sp.show(screen,30,30,fil,pal)) |
---|
| 361 | { update_dirty(screen); |
---|
| 362 | screen->bar(30,30,30+sp.w,30+sp.h,0); |
---|
| 363 | } |
---|
| 364 | sleep(2); |
---|
| 365 | close_graphics(); |
---|
| 366 | }*/ |
---|
| 367 | |
---|
| 368 | |
---|