ccdist.c

Computing 2D Skeleton

  */
   
 
 /* NORMALS: BASE LINE FOR SCHIZOPHRENICS */
 /*  
   strcpy(dir_name,"::CorpusCallosum:ccSchzNorms:");

   strcpy(str[0],"::CorpusCallosum:ccSchzNorms:1019_2_128.buchar");
   strcpy(str[1],"::CorpusCallosum:ccSchzNorms:1039_4_128.buchar");
   strcpy(str[2],"::CorpusCallosum:ccSchzNorms:1042_2_128.buchar");
   strcpy(str[3],"::CorpusCallosum:ccSchzNorms:1044_2_128.buchar");
   strcpy(str[4],"::CorpusCallosum:ccSchzNorms:1045_2_128.buchar");
   strcpy(str[5],"::CorpusCallosum:ccSchzNorms:1046_2_128.buchar");
   strcpy(str[6],"::CorpusCallosum:ccSchzNorms:1050_3_128.buchar");
   strcpy(str[7],"::CorpusCallosum:ccSchzNorms:1069_2_128.buchar");
   strcpy(str[8],"::CorpusCallosum:ccSchzNorms:1077_2_128.buchar");
   strcpy(str[9],"::CorpusCallosum:ccSchzNorms:1110_2_128.buchar");
   strcpy(str[10],"::CorpusCallosum:ccSchzNorms:1112_2_128.buchar");
   strcpy(str[11],"::CorpusCallosum:ccSchzNorms:1120_2_128.buchar");
   strcpy(str[12],"::CorpusCallosum:ccSchzNorms:1134_2_128.buchar");
   strcpy(str[13],"::CorpusCallosum:ccSchzNorms:1150_2_128.buchar");
   strcpy(str[14],"::CorpusCallosum:ccSchzNorms:1151_2_128.buchar");
   strcpy(str[15],"::CorpusCallosum:ccSchzNorms:1152_2_128.buchar");
   strcpy(str[16],"::CorpusCallosum:ccSchzNorms:1153_2_128.buchar");
   strcpy(str[17],"::CorpusCallosum:ccSchzNorms:1154_2_128.buchar");
   strcpy(str[18],"::CorpusCallosum:ccSchzNorms:1155_2_128.buchar");
   strcpy(str[19],"::CorpusCallosum:ccSchzNorms:1156_2_128.buchar");
   strcpy(str[20],"::CorpusCallosum:ccSchzNorms:1157_2_128.buchar");
   strcpy(str[21],"::CorpusCallosum:ccSchzNorms:1158_2_128.buchar");
   strcpy(str[22],"::CorpusCallosum:ccSchzNorms:1159_2_128.buchar");
   strcpy(str[23],"::CorpusCallosum:ccSchzNorms:1182_3_128.buchar");
   strcpy(str[24],"::CorpusCallosum:ccSchzNorms:521_4_130.buchar");
   strcpy(str[25],"::CorpusCallosum:ccSchzNorms:522_4_128.buchar");
   strcpy(str[26],"::CorpusCallosum:ccSchzNorms:538_4_128.buchar");
   strcpy(str[27],"::CorpusCallosum:ccSchzNorms:546_4_128.buchar");
   strcpy(str[28],"::CorpusCallosum:ccSchzNorms:594_5_128.buchar");
   strcpy(str[29],"::CorpusCallosum:ccSchzNorms:595_4_128.buchar");
   strcpy(str[30],"::CorpusCallosum:ccSchzNorms:632_4_128.buchar");
   strcpy(str[31],"::CorpusCallosum:ccSchzNorms:634_4_128.buchar");
   strcpy(str[32],"::CorpusCallosum:ccSchzNorms:638_4_128.buchar");
   strcpy(str[33],"::CorpusCallosum:ccSchzNorms:649_4_128.buchar");
   strcpy(str[34],"::CorpusCallosum:ccSchzNorms:652_4_128.buchar");
   strcpy(str[35],"::CorpusCallosum:ccSchzNorms:653_4_128.buchar");
   strcpy(str[36],"::CorpusCallosum:ccSchzNorms:654_4_128.buchar");
   strcpy(str[37],"::CorpusCallosum:ccSchzNorms:673_4_128.buchar");
   strcpy(str[38],"::CorpusCallosum:ccSchzNorms:674_4_128.buchar");
   strcpy(str[39],"::CorpusCallosum:ccSchzNorms:681_4_128.buchar");
   strcpy(str[40],"::CorpusCallosum:ccSchzNorms:688_4_128.buchar");
   strcpy(str[41],"::CorpusCallosum:ccSchzNorms:794_4_128.buchar");
   strcpy(str[42],"::CorpusCallosum:ccSchzNorms:847_7_128.buchar");
   strcpy(str[43],"::CorpusCallosum:ccSchzNorms:851_6_128.buchar");
   strcpy(str[44],"::CorpusCallosum:ccSchzNorms:860_4_128.buchar");
   strcpy(str[45],"::CorpusCallosum:ccSchzNorms:871_4_128.buchar");
   strcpy(str[46],"::CorpusCallosum:ccSchzNorms:874_4_128.buchar");
   strcpy(str[47],"::CorpusCallosum:ccSchzNorms:880_4_128.buchar");

   strcpy(sstr[0],"::CorpusCallosum:ccSchzNorms:1019");
   strcpy(sstr[1],"::CorpusCallosum:ccSchzNorms:1039");
   strcpy(sstr[2],"::CorpusCallosum:ccSchzNorms:1042");
   strcpy(sstr[3],"::CorpusCallosum:ccSchzNorms:1044");
   strcpy(sstr[4],"::CorpusCallosum:ccSchzNorms:1045");
   strcpy(sstr[5],"::CorpusCallosum:ccSchzNorms:1046");
   strcpy(sstr[6],"::CorpusCallosum:ccSchzNorms:1050");
   strcpy(sstr[7],"::CorpusCallosum:ccSchzNorms:1069");
   strcpy(sstr[8],"::CorpusCallosum:ccSchzNorms:1077");
   strcpy(sstr[9],"::CorpusCallosum:ccSchzNorms:1110");
   strcpy(sstr[10],"::CorpusCallosum:ccSchzNorms:1112");
   strcpy(sstr[11],"::CorpusCallosum:ccSchzNorms:1120");
   strcpy(sstr[12],"::CorpusCallosum:ccSchzNorms:1134");
   strcpy(sstr[13],"::CorpusCallosum:ccSchzNorms:1150");
   strcpy(sstr[14],"::CorpusCallosum:ccSchzNorms:1151");
   strcpy(sstr[15],"::CorpusCallosum:ccSchzNorms:1152");
   strcpy(sstr[16],"::CorpusCallosum:ccSchzNorms:1153");
   strcpy(sstr[17],"::CorpusCallosum:ccSchzNorms:1154");
   strcpy(sstr[18],"::CorpusCallosum:ccSchzNorms:1155");
   strcpy(sstr[19],"::CorpusCallosum:ccSchzNorms:1156");
   strcpy(sstr[20],"::CorpusCallosum:ccSchzNorms:1157");
   strcpy(sstr[21],"::CorpusCallosum:ccSchzNorms:1158");
   strcpy(sstr[22],"::CorpusCallosum:ccSchzNorms:1159");
   strcpy(sstr[23],"::CorpusCallosum:ccSchzNorms:1182");
   strcpy(sstr[24],"::CorpusCallosum:ccSchzNorms:521");
   strcpy(sstr[25],"::CorpusCallosum:ccSchzNorms:522");
   strcpy(sstr[26],"::CorpusCallosum:ccSchzNorms:538");
   strcpy(sstr[27],"::CorpusCallosum:ccSchzNorms:546");
   strcpy(sstr[28],"::CorpusCallosum:ccSchzNorms:594");
   strcpy(sstr[29],"::CorpusCallosum:ccSchzNorms:595");
   strcpy(sstr[30],"::CorpusCallosum:ccSchzNorms:632");
   strcpy(sstr[31],"::CorpusCallosum:ccSchzNorms:634");
   strcpy(sstr[32],"::CorpusCallosum:ccSchzNorms:638");
   strcpy(sstr[33],"::CorpusCallosum:ccSchzNorms:649");
   strcpy(sstr[34],"::CorpusCallosum:ccSchzNorms:652");
   strcpy(sstr[35],"::CorpusCallosum:ccSchzNorms:653");
   strcpy(sstr[36],"::CorpusCallosum:ccSchzNorms:654");
   strcpy(sstr[37],"::CorpusCallosum:ccSchzNorms:673");
   strcpy(sstr[38],"::CorpusCallosum:ccSchzNorms:674");
   strcpy(sstr[39],"::CorpusCallosum:ccSchzNorms:681");
   strcpy(sstr[40],"::CorpusCallosum:ccSchzNorms:688");
   strcpy(sstr[41],"::CorpusCallosum:ccSchzNorms:794");
   strcpy(sstr[42],"::CorpusCallosum:ccSchzNorms:847");
   strcpy(sstr[43],"::CorpusCallosum:ccSchzNorms:851");
   strcpy(sstr[44],"::CorpusCallosum:ccSchzNorms:860");
   strcpy(sstr[45],"::CorpusCallosum:ccSchzNorms:871");
   strcpy(sstr[46],"::CorpusCallosum:ccSchzNorms:874");
   strcpy(sstr[47],"::CorpusCallosum:ccSchzNorms:880");
  */
  
 
 
 /* NORMALS */
 /* The first ten are females, the last ten males.*/
 
 /*
   strcpy(dir_name,"::CorpusCallosum:ccNorms:");

   strcpy(str[0],"::CorpusCallosum:ccNorms:100_23_127.buchar");
   strcpy(str[1],"::CorpusCallosum:ccNorms:16_3_134.buchar");
   strcpy(str[2],"::CorpusCallosum:ccNorms:17_3_132.buchar");
   strcpy(str[3],"::CorpusCallosum:ccNorms:191_3_133.buchar");
   strcpy(str[4],"::CorpusCallosum:ccNorms:202_3_131.buchar");
   strcpy(str[5],"::CorpusCallosum:ccNorms:205_3_120.buchar");
   strcpy(str[6],"::CorpusCallosum:ccNorms:1_24_130.buchar");
   strcpy(str[7],"::CorpusCallosum:ccNorms:2_4_131.buchar");
   strcpy(str[8],"::CorpusCallosum:ccNorms:4_8_128.buchar");
   strcpy(str[9],"::CorpusCallosum:ccNorms:5_8_131.buchar");

   strcpy(str[10],"::CorpusCallosum:ccNorms:110_3_127.buchar");
   strcpy(str[11],"::CorpusCallosum:ccNorms:111_2_122.buchar");
   strcpy(str[12],"::CorpusCallosum:ccNorms:112_2_121.buchar");
   strcpy(str[13],"::CorpusCallosum:ccNorms:11_3_131.buchar");
   strcpy(str[14],"::CorpusCallosum:ccNorms:12_3_137.buchar");
   strcpy(str[15],"::CorpusCallosum:ccNorms:13_3_127.buchar");
   strcpy(str[16],"::CorpusCallosum:ccNorms:15_3_137.buchar");
   strcpy(str[17],"::CorpusCallosum:ccNorms:6_10_131.buchar");
   strcpy(str[18],"::CorpusCallosum:ccNorms:7_8_132.buchar");
   strcpy(str[19],"::CorpusCallosum:ccNorms:8_4_133.buchar");

   strcpy(sstr[0],"::CorpusCallosum:ccNorms:100");
   strcpy(sstr[1],"::CorpusCallosum:ccNorms:16");
   strcpy(sstr[2],"::CorpusCallosum:ccNorms:17");
   strcpy(sstr[3],"::CorpusCallosum:ccNorms:191");
   strcpy(sstr[4],"::CorpusCallosum:ccNorms:202");
   strcpy(sstr[5],"::CorpusCallosum:ccNorms:205");
   strcpy(sstr[6],"::CorpusCallosum:ccNorms:1");
   strcpy(sstr[7],"::CorpusCallosum:ccNorms:2");
   strcpy(sstr[8],"::CorpusCallosum:ccNorms:4");
   strcpy(sstr[9],"::CorpusCallosum:ccNorms:5");

   strcpy(sstr[10],"::CorpusCallosum:ccNorms:110");
   strcpy(sstr[11],"::CorpusCallosum:ccNorms:111");
   strcpy(sstr[12],"::CorpusCallosum:ccNorms:112");
   strcpy(sstr[13],"::CorpusCallosum:ccNorms:11");
   strcpy(sstr[14],"::CorpusCallosum:ccNorms:12");
   strcpy(sstr[15],"::CorpusCallosum:ccNorms:13");
   strcpy(sstr[16],"::CorpusCallosum:ccNorms:15");
   strcpy(sstr[17],"::CorpusCallosum:ccNorms:6");
   strcpy(sstr[18],"::CorpusCallosum:ccNorms:7");
   strcpy(sstr[19],"::CorpusCallosum:ccNorms:8");
*/


/* KIDS */
/* The first 15 are females, the last 15 males. */

/*  strcpy(dir_name,"::CorpusCallosum:ccKids:"); 
  
   strcpy(str[0],"::CorpusCallosum:ccKids:132_5_123.buchar");
   strcpy(str[1],"::CorpusCallosum:ccKids:158_4_129.buchar");
   strcpy(str[2],"::CorpusCallosum:ccKids:161_5_126.buchar");
   strcpy(str[3],"::CorpusCallosum:ccKids:169_4_127.buchar");
   strcpy(str[4],"::CorpusCallosum:ccKids:267_6_135.buchar");
   strcpy(str[5],"::CorpusCallosum:ccKids:277_5_90.buchar");
   strcpy(str[6],"::CorpusCallosum:ccKids:278_5_124.buchar");
   strcpy(str[7],"::CorpusCallosum:ccKids:291_5_129.buchar");
   strcpy(str[8],"::CorpusCallosum:ccKids:292_6_127.buchar");
   strcpy(str[9],"::CorpusCallosum:ccKids:349_4_103.buchar");
   strcpy(str[10],"::CorpusCallosum:ccKids:376_4_128.buchar");
   strcpy(str[11],"::CorpusCallosum:ccKids:386_4_124.buchar");
   strcpy(str[12],"::CorpusCallosum:ccKids:392_4_131.buchar");
   strcpy(str[13],"::CorpusCallosum:ccKids:401_4_124.buchar");
   strcpy(str[14],"::CorpusCallosum:ccKids:9_12_132.buchar");

   strcpy(str[15],"::CorpusCallosum:ccKids:139_4_125.buchar");
   strcpy(str[16],"::CorpusCallosum:ccKids:140_4_125.buchar");
   strcpy(str[17],"::CorpusCallosum:ccKids:146_4_128.buchar");
   strcpy(str[18],"::CorpusCallosum:ccKids:159_4_130.buchar");
   strcpy(str[19],"::CorpusCallosum:ccKids:160_4_124.buchar");
   strcpy(str[20],"::CorpusCallosum:ccKids:295_4_129.buchar");
   strcpy(str[21],"::CorpusCallosum:ccKids:308_6_129.buchar");
   strcpy(str[22],"::CorpusCallosum:ccKids:309_7_127.buchar");
   strcpy(str[23],"::CorpusCallosum:ccKids:310_4_119.buchar");
   strcpy(str[24],"::CorpusCallosum:ccKids:312_4_135.buchar");
   strcpy(str[25],"::CorpusCallosum:ccKids:323_4_127.buchar");
   strcpy(str[26],"::CorpusCallosum:ccKids:330_5_125.buchar");
   strcpy(str[27],"::CorpusCallosum:ccKids:385_4_121.buchar");
   strcpy(str[28],"::CorpusCallosum:ccKids:387_4_130.buchar");
   strcpy(str[29],"::CorpusCallosum:ccKids:415_4_129.buchar");

   strcpy(sstr[0],"::CorpusCallosum:ccKids:132");
   strcpy(sstr[1],"::CorpusCallosum:ccKids:158");
   strcpy(sstr[2],"::CorpusCallosum:ccKids:161");
   strcpy(sstr[3],"::CorpusCallosum:ccKids:169");
   strcpy(sstr[4],"::CorpusCallosum:ccKids:267");
   strcpy(sstr[5],"::CorpusCallosum:ccKids:277");
   strcpy(sstr[6],"::CorpusCallosum:ccKids:278");
   strcpy(sstr[7],"::CorpusCallosum:ccKids:291");
   strcpy(sstr[8],"::CorpusCallosum:ccKids:292");
   strcpy(sstr[9],"::CorpusCallosum:ccKids:349");
   strcpy(sstr[10],"::CorpusCallosum:ccKids:376");
   strcpy(sstr[11],"::CorpusCallosum:ccKids:386");
   strcpy(sstr[12],"::CorpusCallosum:ccKids:392");
   strcpy(sstr[13],"::CorpusCallosum:ccKids:401");
   strcpy(sstr[14],"::CorpusCallosum:ccKids:9");


   strcpy(sstr[15],"::CorpusCallosum:ccKids:139");
   strcpy(sstr[16],"::CorpusCallosum:ccKids:140");
   strcpy(sstr[17],"::CorpusCallosum:ccKids:146");
   strcpy(sstr[18],"::CorpusCallosum:ccKids:159");
   strcpy(sstr[19],"::CorpusCallosum:ccKids:160");
   strcpy(sstr[20],"::CorpusCallosum:ccKids:295");
   strcpy(sstr[21],"::CorpusCallosum:ccKids:308");
   strcpy(sstr[22],"::CorpusCallosum:ccKids:309");
   strcpy(sstr[23],"::CorpusCallosum:ccKids:310");
   strcpy(sstr[24],"::CorpusCallosum:ccKids:312");
   strcpy(sstr[25],"::CorpusCallosum:ccKids:323");
   strcpy(sstr[26],"::CorpusCallosum:ccKids:330");
   strcpy(sstr[27],"::CorpusCallosum:ccKids:385");
   strcpy(sstr[28],"::CorpusCallosum:ccKids:387");
   strcpy(sstr[29],"::CorpusCallosum:ccKids:415");
*/

/* DRUGS */
/* first 7 belong to group I; the next 6 belong to
group II and the last belongs to group III  */

   strcpy(dir_name,"::CorpusCallosum:ccDrugs:");
   
   strcpy(str[0],"::CorpusCallosum:ccDrugs:2560_2");
   strcpy(str[1],"::CorpusCallosum:ccDrugs:2563_2");
   strcpy(str[2],"::CorpusCallosum:ccDrugs:2564_2");
   strcpy(str[3],"::CorpusCallosum:ccDrugs:2569_2");
   strcpy(str[4],"::CorpusCallosum:ccDrugs:2570_2");
   strcpy(str[5],"::CorpusCallosum:ccDrugs:2573_2");
   strcpy(str[6],"::CorpusCallosum:ccDrugs:2574_2");
   
   strcpy(str[7],"::CorpusCallosum:ccDrugs:2561_2");
   strcpy(str[8],"::CorpusCallosum:ccDrugs:2562_2");
   strcpy(str[9],"::CorpusCallosum:ccDrugs:2565_2");
   strcpy(str[10],"::CorpusCallosum:ccDrugs:2566_2");
   strcpy(str[11],"::CorpusCallosum:ccDrugs:2567_2");
   strcpy(str[12],"::CorpusCallosum:ccDrugs:2568_2");
   strcpy(str[13],"::CorpusCallosum:ccDrugs:2572_2");
   
   strcpy(sstr[0],"::CorpusCallosum:ccDrugs:2560");
   strcpy(sstr[1],"::CorpusCallosum:ccDrugs:2563");
   strcpy(sstr[2],"::CorpusCallosum:ccDrugs:2564");
   strcpy(sstr[3],"::CorpusCallosum:ccDrugs:2569");
   strcpy(sstr[4],"::CorpusCallosum:ccDrugs:2570");
   strcpy(sstr[5],"::CorpusCallosum:ccDrugs:2573");
   strcpy(sstr[6],"::CorpusCallosum:ccDrugs:2574");
   
   strcpy(sstr[7],"::CorpusCallosum:ccDrugs:2561");
   strcpy(sstr[8],"::CorpusCallosum:ccDrugs:2562");
   strcpy(sstr[9],"::CorpusCallosum:ccDrugs:2565");
   strcpy(sstr[10],"::CorpusCallosum:ccDrugs:2566");
   strcpy(sstr[11],"::CorpusCallosum:ccDrugs:2567");
   strcpy(sstr[12],"::CorpusCallosum:ccDrugs:2568");
   strcpy(sstr[13],"::CorpusCallosum:ccDrugs:2572");
  
   
   for (m=0;m<sample_count;m++) 
   {
      printf("m = %d\n",m);
	  strcpy(inname,str[m]);
	  strcpy(outname, sstr[m]);
	  strcat(outname, "dist"); 
   
   /* READ IN THE IMAGE */
	
	  printf("Opening %s\n",inname);
      if ((infile = fopen(inname,"r+")) == NULL)
      {
	    printf("Cannot open %s\n",inname);
	    exit(1);
      }   

	  n = fread(pic,sizeof(unsigned char),ydim*xdim,infile);
	  printf("%d bytes in image read\n",n);
	  fclose(infile);
	  
      calc_dist(); 
      output_results();
   }

   printf("\n\n done \n");
   systime = time(NULL);
   puts(ctime(&systime));
}