00001
00002
00003
00004
00005
00006
00007
00008 #include <string.h>
00009 #include <stdlib.h>
00010 #include <cmath>
00011
00012 #include "susan.h"
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
00069
00070
00071
00072
00073
00074
00075
00076
00077
00078
00079
00080
00081
00082
00083
00084
00085
00086
00087
00088
00089
00090
00091
00092
00093
00094
00095
00096
00097
00098
00099
00100
00101
00102
00103
00104
00105
00106
00107
00108
00109
00110
00111
00112
00113
00114
00115
00116
00117
00118
00119
00120
00121
00122
00123
00124
00125
00126
00127
00128
00129
00130
00131
00132
00133
00134
00135
00136
00137
00138
00139
00140
00141
00142
00143
00144
00145
00146
00147
00148
00149
00150
00151
00152
00153
00154
00155
00156
00157
00158
00159
00160
00161
00162
00163
00164
00165
00166
00167
00168
00169
00170
00171
00172
00173
00174
00175
00176
00177
00178
00179
00180
00181
00182
00183
00184
00185
00186
00187
00188
00189
00190
00191
00192
00193
00194
00195
00196
00197
00198
00199
00200
00201
00202
00203
00204
00205
00206
00207
00208
00209
00210
00211
00212
00213
00214
00215
00216
00217
00218
00219
00220
00221
00222
00223
00224
00225
00226
00227
00228
00229
00230
00231
00232
00233
00234
00235
00236
00237
00238
00239
00240
00241
00242
00243
00244
00245
00246
00247
00248
00249
00250
00251
00252
00253
00254
00255
00256
00257
00258
00259
00260
00261
00262
00263
00264
00265
00266
00267
00268
00269
00270
00271
00272
00273
00274
00275
00276
00277
00278
00279
00280
00281
00282
00283
00284
00285
00286
00287
00288
00289
00290
00291
00292
00293
00294
00295
00296 namespace DualCoding {
00297
00298
00299
00300
00301
00302
00303
00304 void setup_brightness_lut(uchar **bp, int thresh, int form)
00305 {
00306 int k;
00307 float temp;
00308
00309 *bp=(unsigned char *)malloc(516);
00310 *bp=*bp+258;
00311
00312 for(k=-256;k<257;k++)
00313 {
00314 temp=((float)k)/((float)thresh);
00315 temp=temp*temp;
00316 if (form==6)
00317 temp=temp*temp*temp;
00318 temp=100*std::exp(-temp);
00319 *(*bp+k)= (uchar)temp;
00320 }
00321 }
00322
00323
00324
00325
00326
00327
00328
00329
00330
00331
00332 void susan_principle(uchar *in, int *r, uchar *bp, int max_no,
00333 int x_size, int y_size)
00334 {
00335 int i, j, n;
00336 uchar *p,*cp;
00337
00338 memset (r,0,x_size * y_size * sizeof(int));
00339
00340
00341 for (i=3;i<y_size-3;i++)
00342 for (j=3;j<x_size-3;j++)
00343 {
00344 n=100;
00345 p=in + (i-3)*x_size + j - 1;
00346 cp=bp + in[i*x_size+j];
00347
00348 n+=*(cp-*p++);
00349 n+=*(cp-*p++);
00350 n+=*(cp-*p);
00351 p+=x_size-3;
00352
00353 n+=*(cp-*p++);
00354 n+=*(cp-*p++);
00355 n+=*(cp-*p++);
00356 n+=*(cp-*p++);
00357 n+=*(cp-*p);
00358 p+=x_size-5;
00359
00360 n+=*(cp-*p++);
00361 n+=*(cp-*p++);
00362 n+=*(cp-*p++);
00363 n+=*(cp-*p++);
00364 n+=*(cp-*p++);
00365 n+=*(cp-*p++);
00366 n+=*(cp-*p);
00367 p+=x_size-6;
00368
00369 n+=*(cp-*p++);
00370 n+=*(cp-*p++);
00371 n+=*(cp-*p);
00372 p+=2;
00373 n+=*(cp-*p++);
00374 n+=*(cp-*p++);
00375 n+=*(cp-*p);
00376 p+=x_size-6;
00377
00378 n+=*(cp-*p++);
00379 n+=*(cp-*p++);
00380 n+=*(cp-*p++);
00381 n+=*(cp-*p++);
00382 n+=*(cp-*p++);
00383 n+=*(cp-*p++);
00384 n+=*(cp-*p);
00385 p+=x_size-5;
00386
00387 n+=*(cp-*p++);
00388 n+=*(cp-*p++);
00389 n+=*(cp-*p++);
00390 n+=*(cp-*p++);
00391 n+=*(cp-*p);
00392 p+=x_size-3;
00393
00394 n+=*(cp-*p++);
00395 n+=*(cp-*p++);
00396 n+=*(cp-*p);
00397
00398 if (n<=max_no)
00399 r[i*x_size+j] = max_no - n;
00400 }
00401 }
00402
00403
00404
00405
00406
00407
00408
00409
00410
00411
00412 void susan_thin(int *r, uchar *mid, int x_size, int y_size)
00413 {
00414 int l[9], centre,
00415 b01, b12, b21, b10,
00416 p1, p2, p3, p4,
00417 b00, b02, b20, b22,
00418 m, n, a=0, b=0, x, y, i, j;
00419 uchar *mp;
00420
00421 for (i=4;i<y_size-4;i++)
00422 for (j=4;j<x_size-4;j++)
00423 if (mid[i*x_size+j]<8)
00424 {
00425 centre = r[i*x_size+j];
00426
00427
00428 mp=mid + (i-1)*x_size + j-1;
00429
00430 n = (*mp<8) +
00431 (*(mp+1)<8) +
00432 (*(mp+2)<8) +
00433 (*(mp+x_size)<8) +
00434 (*(mp+x_size+2)<8) +
00435 (*(mp+x_size+x_size)<8) +
00436 (*(mp+x_size+x_size+1)<8) +
00437 (*(mp+x_size+x_size+2)<8);
00438
00439
00440
00441
00442 if (n==0)
00443 mid[i*x_size+j]=100;
00444
00445
00446
00447
00448
00449
00450 if ( (n==1) && (mid[i*x_size+j]<6) )
00451 {
00452
00453
00454
00455
00456
00457
00458 l[0]=r[(i-1)*x_size+j-1]; l[1]=r[(i-1)*x_size+j]; l[2]=r[(i-1)*x_size+j+1];
00459 l[3]=r[(i )*x_size+j-1]; l[4]=0; l[5]=r[(i )*x_size+j+1];
00460 l[6]=r[(i+1)*x_size+j-1]; l[7]=r[(i+1)*x_size+j]; l[8]=r[(i+1)*x_size+j+1];
00461
00462 if (mid[(i-1)*x_size+j-1]<8) { l[0]=0; l[1]=0; l[3]=0; l[2]*=2;
00463 l[6]*=2; l[5]*=3; l[7]*=3; l[8]*=4; }
00464 else { if (mid[(i-1)*x_size+j]<8) { l[1]=0; l[0]=0; l[2]=0; l[3]*=2;
00465 l[5]*=2; l[6]*=3; l[8]*=3; l[7]*=4; }
00466 else { if (mid[(i-1)*x_size+j+1]<8) { l[2]=0; l[1]=0; l[5]=0; l[0]*=2;
00467 l[8]*=2; l[3]*=3; l[7]*=3; l[6]*=4; }
00468 else { if (mid[(i)*x_size+j-1]<8) { l[3]=0; l[0]=0; l[6]=0; l[1]*=2;
00469 l[7]*=2; l[2]*=3; l[8]*=3; l[5]*=4; }
00470 else { if (mid[(i)*x_size+j+1]<8) { l[5]=0; l[2]=0; l[8]=0; l[1]*=2;
00471 l[7]*=2; l[0]*=3; l[6]*=3; l[3]*=4; }
00472 else { if (mid[(i+1)*x_size+j-1]<8) { l[6]=0; l[3]=0; l[7]=0; l[0]*=2;
00473 l[8]*=2; l[1]*=3; l[5]*=3; l[2]*=4; }
00474 else { if (mid[(i+1)*x_size+j]<8) { l[7]=0; l[6]=0; l[8]=0; l[3]*=2;
00475 l[5]*=2; l[0]*=3; l[2]*=3; l[1]*=4; }
00476 else { if (mid[(i+1)*x_size+j+1]<8) { l[8]=0; l[5]=0; l[7]=0; l[6]*=2;
00477 l[2]*=2; l[1]*=3; l[3]*=3; l[0]*=4; } }}}}}}}
00478
00479 m=0;
00480 for(y=0; y<3; y++)
00481 for(x=0; x<3; x++)
00482 if (l[y+y+y+x]>m) { m=l[y+y+y+x]; a=y; b=x; }
00483
00484 if (m>0)
00485 {
00486 if (mid[i*x_size+j]<4)
00487 mid[(i+a-1)*x_size+j+b-1] = 4;
00488 else
00489 mid[(i+a-1)*x_size+j+b-1] = mid[i*x_size+j]+1;
00490 if ( (a+a+b) < 3 )
00491 {
00492 i+=a-1;
00493 j+=b-2;
00494 if (i<4) i=4;
00495 if (j<4) j=4;
00496 }
00497 }
00498 }
00499
00500
00501
00502
00503 if (n==2)
00504 {
00505
00506 b00 = mid[(i-1)*x_size+j-1]<8;
00507 b02 = mid[(i-1)*x_size+j+1]<8;
00508 b20 = mid[(i+1)*x_size+j-1]<8;
00509 b22 = mid[(i+1)*x_size+j+1]<8;
00510 if ( ((b00+b02+b20+b22)==2) && ((b00|b22)&(b02|b20)))
00511 {
00512
00513
00514
00515 if (b00)
00516 {
00517 if (b02) { x=0; y=-1; }
00518 else { x=-1; y=0; }
00519 }
00520 else
00521 {
00522 if (b02) { x=1; y=0; }
00523 else { x=0; y=1; }
00524 }
00525 if (((float)r[(i+y)*x_size+j+x]/(float)centre) > 0.7)
00526 {
00527 if ( ( (x==0) && (mid[(i+(2*y))*x_size+j]>7) && (mid[(i+(2*y))*x_size+j-1]>7) && (mid[(i+(2*y))*x_size+j+1]>7) ) ||
00528 ( (y==0) && (mid[(i)*x_size+j+(2*x)]>7) && (mid[(i+1)*x_size+j+(2*x)]>7) && (mid[(i-1)*x_size+j+(2*x)]>7) ) )
00529 {
00530 mid[(i)*x_size+j]=100;
00531 mid[(i+y)*x_size+j+x]=3;
00532 }
00533 }
00534 }
00535 else
00536 {
00537 b01 = mid[(i-1)*x_size+j ]<8;
00538 b12 = mid[(i )*x_size+j+1]<8;
00539 b21 = mid[(i+1)*x_size+j ]<8;
00540 b10 = mid[(i )*x_size+j-1]<8;
00541
00542
00543 #ifdef IGNORETHIS
00544 if ( (b00&b01)|(b00&b10)|(b02&b01)|(b02&b12)|(b20&b10)|(b20&b21)|(b22&b21)|(b22&b12) )
00545 {
00546
00547
00548
00549 if ( ((b01)&(mid[(i-2)*x_size+j-1]>7)&(mid[(i-2)*x_size+j]>7)&(mid[(i-2)*x_size+j+1]>7)&
00550 ((b00&((2*r[(i-1)*x_size+j+1])>centre))|(b02&((2*r[(i-1)*x_size+j-1])>centre)))) |
00551 ((b10)&(mid[(i-1)*x_size+j-2]>7)&(mid[(i)*x_size+j-2]>7)&(mid[(i+1)*x_size+j-2]>7)&
00552 ((b00&((2*r[(i+1)*x_size+j-1])>centre))|(b20&((2*r[(i-1)*x_size+j-1])>centre)))) |
00553 ((b12)&(mid[(i-1)*x_size+j+2]>7)&(mid[(i)*x_size+j+2]>7)&(mid[(i+1)*x_size+j+2]>7)&
00554 ((b02&((2*r[(i+1)*x_size+j+1])>centre))|(b22&((2*r[(i-1)*x_size+j+1])>centre)))) |
00555 ((b21)&(mid[(i+2)*x_size+j-1]>7)&(mid[(i+2)*x_size+j]>7)&(mid[(i+2)*x_size+j+1]>7)&
00556 ((b20&((2*r[(i+1)*x_size+j+1])>centre))|(b22&((2*r[(i+1)*x_size+j-1])>centre)))) )
00557 {
00558 mid[(i)*x_size+j]=100;
00559 if (b10&b20) j-=2;
00560 if (b00|b01|b02) { i--; j-=2; }
00561 }
00562 }
00563 #endif
00564
00565
00566 if ( ((b01+b12+b21+b10)==2) && ((b10|b12)&(b01|b21)) &&
00567 ((b01&((mid[(i-2)*x_size+j-1]<8)|(mid[(i-2)*x_size+j+1]<8)))|(b10&((mid[(i-1)*x_size+j-2]<8)|(mid[(i+1)*x_size+j-2]<8)))|
00568 (b12&((mid[(i-1)*x_size+j+2]<8)|(mid[(i+1)*x_size+j+2]<8)))|(b21&((mid[(i+2)*x_size+j-1]<8)|(mid[(i+2)*x_size+j+1]<8)))) )
00569 {
00570
00571
00572
00573 mid[(i)*x_size+j]=100;
00574 i--;
00575 j-=2;
00576 if (i<4) i=4;
00577 if (j<4) j=4;
00578 }
00579 }
00580 }
00581
00582
00583
00584
00585 if (n>2)
00586 {
00587 b01 = mid[(i-1)*x_size+j ]<8;
00588 b12 = mid[(i )*x_size+j+1]<8;
00589 b21 = mid[(i+1)*x_size+j ]<8;
00590 b10 = mid[(i )*x_size+j-1]<8;
00591 if((b01+b12+b21+b10)>1)
00592 {
00593 b00 = mid[(i-1)*x_size+j-1]<8;
00594 b02 = mid[(i-1)*x_size+j+1]<8;
00595 b20 = mid[(i+1)*x_size+j-1]<8;
00596 b22 = mid[(i+1)*x_size+j+1]<8;
00597 p1 = b00 | b01;
00598 p2 = b02 | b12;
00599 p3 = b22 | b21;
00600 p4 = b20 | b10;
00601
00602 if( ((p1 + p2 + p3 + p4) - ((b01 & p2)+(b12 & p3)+(b21 & p4)+(b10 & p1))) < 2)
00603 {
00604 mid[(i)*x_size+j]=100;
00605 i--;
00606 j-=2;
00607 if (i<4) i=4;
00608 if (j<4) j=4;
00609 }
00610 }
00611 }
00612
00613 }
00614 }
00615
00616
00617
00618
00619
00620
00621
00622
00623
00624 void susan_edges_internal(uchar *in, int *r, uchar *mid, uchar *bp,
00625 int max_no, int x_size, int y_size)
00626 {
00627 float z;
00628 int do_symmetry, i, j, m, n, a, b, x, y, w;
00629 uchar c,*p,*cp;
00630
00631 memset (r,0,x_size * y_size * sizeof(int));
00632
00633 for (i=3;i<y_size-3;i++)
00634 for (j=3;j<x_size-3;j++)
00635 {
00636 n=100;
00637 p=in + (i-3)*x_size + j - 1;
00638 cp=bp + in[i*x_size+j];
00639
00640 n+=*(cp-*p++);
00641 n+=*(cp-*p++);
00642 n+=*(cp-*p);
00643 p+=x_size-3;
00644
00645 n+=*(cp-*p++);
00646 n+=*(cp-*p++);
00647 n+=*(cp-*p++);
00648 n+=*(cp-*p++);
00649 n+=*(cp-*p);
00650 p+=x_size-5;
00651
00652 n+=*(cp-*p++);
00653 n+=*(cp-*p++);
00654 n+=*(cp-*p++);
00655 n+=*(cp-*p++);
00656 n+=*(cp-*p++);
00657 n+=*(cp-*p++);
00658 n+=*(cp-*p);
00659 p+=x_size-6;
00660
00661 n+=*(cp-*p++);
00662 n+=*(cp-*p++);
00663 n+=*(cp-*p);
00664 p+=2;
00665 n+=*(cp-*p++);
00666 n+=*(cp-*p++);
00667 n+=*(cp-*p);
00668 p+=x_size-6;
00669
00670 n+=*(cp-*p++);
00671 n+=*(cp-*p++);
00672 n+=*(cp-*p++);
00673 n+=*(cp-*p++);
00674 n+=*(cp-*p++);
00675 n+=*(cp-*p++);
00676 n+=*(cp-*p);
00677 p+=x_size-5;
00678
00679 n+=*(cp-*p++);
00680 n+=*(cp-*p++);
00681 n+=*(cp-*p++);
00682 n+=*(cp-*p++);
00683 n+=*(cp-*p);
00684 p+=x_size-3;
00685
00686 n+=*(cp-*p++);
00687 n+=*(cp-*p++);
00688 n+=*(cp-*p);
00689
00690 if (n<=max_no)
00691 r[i*x_size+j] = max_no - n;
00692 }
00693
00694 for (i=4;i<y_size-4;i++)
00695 for (j=4;j<x_size-4;j++)
00696 {
00697 if (r[i*x_size+j]>0)
00698 {
00699 m=r[i*x_size+j];
00700 n=max_no - m;
00701 cp=bp + in[i*x_size+j];
00702
00703 if (n>600)
00704 {
00705 p=in + (i-3)*x_size + j - 1;
00706 x=0;y=0;
00707
00708 c=*(cp-*p++);x-=c;y-=3*c;
00709 c=*(cp-*p++);y-=3*c;
00710 c=*(cp-*p);x+=c;y-=3*c;
00711 p+=x_size-3;
00712
00713 c=*(cp-*p++);x-=2*c;y-=2*c;
00714 c=*(cp-*p++);x-=c;y-=2*c;
00715 c=*(cp-*p++);y-=2*c;
00716 c=*(cp-*p++);x+=c;y-=2*c;
00717 c=*(cp-*p);x+=2*c;y-=2*c;
00718 p+=x_size-5;
00719
00720 c=*(cp-*p++);x-=3*c;y-=c;
00721 c=*(cp-*p++);x-=2*c;y-=c;
00722 c=*(cp-*p++);x-=c;y-=c;
00723 c=*(cp-*p++);y-=c;
00724 c=*(cp-*p++);x+=c;y-=c;
00725 c=*(cp-*p++);x+=2*c;y-=c;
00726 c=*(cp-*p);x+=3*c;y-=c;
00727 p+=x_size-6;
00728
00729 c=*(cp-*p++);x-=3*c;
00730 c=*(cp-*p++);x-=2*c;
00731 c=*(cp-*p);x-=c;
00732 p+=2;
00733 c=*(cp-*p++);x+=c;
00734 c=*(cp-*p++);x+=2*c;
00735 c=*(cp-*p);x+=3*c;
00736 p+=x_size-6;
00737
00738 c=*(cp-*p++);x-=3*c;y+=c;
00739 c=*(cp-*p++);x-=2*c;y+=c;
00740 c=*(cp-*p++);x-=c;y+=c;
00741 c=*(cp-*p++);y+=c;
00742 c=*(cp-*p++);x+=c;y+=c;
00743 c=*(cp-*p++);x+=2*c;y+=c;
00744 c=*(cp-*p);x+=3*c;y+=c;
00745 p+=x_size-5;
00746
00747 c=*(cp-*p++);x-=2*c;y+=2*c;
00748 c=*(cp-*p++);x-=c;y+=2*c;
00749 c=*(cp-*p++);y+=2*c;
00750 c=*(cp-*p++);x+=c;y+=2*c;
00751 c=*(cp-*p);x+=2*c;y+=2*c;
00752 p+=x_size-3;
00753
00754 c=*(cp-*p++);x-=c;y+=3*c;
00755 c=*(cp-*p++);y+=3*c;
00756 c=*(cp-*p);x+=c;y+=3*c;
00757
00758 z = std::sqrt((float)((x*x) + (y*y)));
00759 if (z > (0.9*(float)n))
00760 {
00761 do_symmetry=0;
00762 if (x==0)
00763 z=1000000.0f;
00764 else
00765 z=((float)y) / ((float)x);
00766 if (z < 0) { z=-z; w=-1; }
00767 else w=1;
00768 if (z < 0.5) { a=0; b=1; }
00769 else { if (z > 2.0) { a=1; b=0; }
00770 else { if (w>0) { a=1; b=1; }
00771 else { a=-1; b=1; }}}
00772 if ( (m > r[(i+a)*x_size+j+b]) && (m >= r[(i-a)*x_size+j-b]) &&
00773 (m > r[(i+(2*a))*x_size+j+(2*b)]) && (m >= r[(i-(2*a))*x_size+j-(2*b)]) )
00774 mid[i*x_size+j] = 1;
00775 }
00776 else
00777 do_symmetry=1;
00778 }
00779 else
00780 do_symmetry=1;
00781
00782 if (do_symmetry==1)
00783 {
00784 p=in + (i-3)*x_size + j - 1;
00785 x=0; y=0; w=0;
00786
00787
00788
00789
00790
00791 c=*(cp-*p++);x+=c;y+=9*c;w+=3*c;
00792 c=*(cp-*p++);y+=9*c;
00793 c=*(cp-*p);x+=c;y+=9*c;w-=3*c;
00794 p+=x_size-3;
00795
00796 c=*(cp-*p++);x+=4*c;y+=4*c;w+=4*c;
00797 c=*(cp-*p++);x+=c;y+=4*c;w+=2*c;
00798 c=*(cp-*p++);y+=4*c;
00799 c=*(cp-*p++);x+=c;y+=4*c;w-=2*c;
00800 c=*(cp-*p);x+=4*c;y+=4*c;w-=4*c;
00801 p+=x_size-5;
00802
00803 c=*(cp-*p++);x+=9*c;y+=c;w+=3*c;
00804 c=*(cp-*p++);x+=4*c;y+=c;w+=2*c;
00805 c=*(cp-*p++);x+=c;y+=c;w+=c;
00806 c=*(cp-*p++);y+=c;
00807 c=*(cp-*p++);x+=c;y+=c;w-=c;
00808 c=*(cp-*p++);x+=4*c;y+=c;w-=2*c;
00809 c=*(cp-*p);x+=9*c;y+=c;w-=3*c;
00810 p+=x_size-6;
00811
00812 c=*(cp-*p++);x+=9*c;
00813 c=*(cp-*p++);x+=4*c;
00814 c=*(cp-*p);x+=c;
00815 p+=2;
00816 c=*(cp-*p++);x+=c;
00817 c=*(cp-*p++);x+=4*c;
00818 c=*(cp-*p);x+=9*c;
00819 p+=x_size-6;
00820
00821 c=*(cp-*p++);x+=9*c;y+=c;w-=3*c;
00822 c=*(cp-*p++);x+=4*c;y+=c;w-=2*c;
00823 c=*(cp-*p++);x+=c;y+=c;w-=c;
00824 c=*(cp-*p++);y+=c;
00825 c=*(cp-*p++);x+=c;y+=c;w+=c;
00826 c=*(cp-*p++);x+=4*c;y+=c;w+=2*c;
00827 c=*(cp-*p);x+=9*c;y+=c;w+=3*c;
00828 p+=x_size-5;
00829
00830 c=*(cp-*p++);x+=4*c;y+=4*c;w-=4*c;
00831 c=*(cp-*p++);x+=c;y+=4*c;w-=2*c;
00832 c=*(cp-*p++);y+=4*c;
00833 c=*(cp-*p++);x+=c;y+=4*c;w+=2*c;
00834 c=*(cp-*p);x+=4*c;y+=4*c;w+=4*c;
00835 p+=x_size-3;
00836
00837 c=*(cp-*p++);x+=c;y+=9*c;w-=3*c;
00838 c=*(cp-*p++);y+=9*c;
00839 c=*(cp-*p);x+=c;y+=9*c;w+=3*c;
00840
00841 if (y==0)
00842 z = 1000000.0f;
00843 else
00844 z = ((float)x) / ((float)y);
00845 if (z < 0.5) { a=0; b=1; }
00846 else { if (z > 2.0) { a=1; b=0; }
00847 else { if (w>0) { a=-1; b=1; }
00848 else { a=1; b=1; }}}
00849 if ( (m > r[(i+a)*x_size+j+b]) && (m >= r[(i-a)*x_size+j-b]) &&
00850 (m > r[(i+(2*a))*x_size+j+(2*b)]) && (m >= r[(i-(2*a))*x_size+j-(2*b)]) )
00851 mid[i*x_size+j] = 2;
00852 }
00853 }
00854 }
00855 }
00856
00857
00858
00859
00860
00861
00862
00863 void edge_draw(uchar *in, uchar *mid, int x_size, int y_size, int drawing_mode)
00864 {
00865 int i;
00866 uchar *inp, *midp;
00867
00868 if (drawing_mode==0)
00869 {
00870
00871 midp=mid;
00872 for (i=0; i<x_size*y_size; i++)
00873 {
00874 if (*midp<8)
00875 {
00876 inp = in + (midp - mid) - x_size - 1;
00877 *inp++=255; *inp++=255; *inp=255; inp+=x_size-2;
00878 *inp++=255; ++inp; *inp=255; inp+=x_size-2;
00879 *inp++=255; *inp++=255; *inp=255;
00880 }
00881 midp++;
00882 }
00883 }
00884
00885
00886 midp=mid;
00887 for (i=0; i<x_size*y_size; i++)
00888 {
00889 if (*midp<8)
00890
00891 *(in + (midp - mid)) = 1;
00892 midp++;
00893 }
00894 }
00895
00896 }