Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  422 1587 (N38E52 L=209)   548    138        1 ( 401/1554- 447/1624) -29.7 E55-E48 N35-N40
      2 1234 1492 (N34W15 L=276)   261     38        1 (1209/1473-1263/1511) -33.7 W13-W17 N32-N34
      3  965 1278 (N19E04 L=257)   214     47        1 ( 937/1261- 995/1291) -14.0 E06-E02 N17-N19
      4  428 1258 (N16E42 L=219)   559    102        2 ( 392/1221- 460/1288) -15.9 E45-E38 N13-N18
      5 1121 1180 (N13W06 L=267)   332    172        1 (1063/1134-1171/1261)  53.1 W02-W09 N10-N18
      6  224 1189 (N11E61 L=200)   786    185        1 ( 178/1139- 255/1229) -32.0 E67-E56 N07-N13
      7  954 1183 (N13E05 L=256)   276    120        1 ( 924/1129-1010/1229) -50.2 E06-E01 N09-N15
      8  407  954 (S03E41 L=220)   374     89        2 ( 371/ 923- 443/ 991) -36.9 E44-E38 S04-N00
      9 1024  860 (S07W00 L=261)   253    120        0 ( 990/ 823-1077/ 879)  14.0 E01-W03 S09-S05
     10 1423  743 (S14W26 L=287)   536     88        1 (1398/ 703-1455/ 787)  59.0 W24-W28 S16-S11
     11  322  709 (S19E51 L=210)   717    158        1 ( 281/ 645- 370/ 771) -53.1 E54-E47 S22-S14
     12 1036  641 (S21W01 L=262)   632    218        3 (1008/ 526-1081/ 711)  76.0 W00-W04 S28-S16
     13  288  639 (S24E56 L=205)   606    110        0 ( 261/ 616- 317/ 671) -36.9 E58-E53 S25-S21
     14  731  564 (S27E20 L=241)   352    148        0 ( 711/ 506- 773/ 613)  66.8 E21-E16 S30-S23
     15 1289  555 (S27W19 L=280)   581    110        2 (1225/ 532-1336/ 588)   7.1 W14-W22 S28-S24
     16  417  428 (S39E51 L=210)   487    110        1 ( 406/ 396- 433/ 458)  51.3 E53-E48 S41-S36
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y