Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  503 1690 (N59E35 L=356)  1041    112        1 ( 459/1658- 541/1723)   4.8 E43-E27 N54-N61
      2 1399 1673 (N33W45 L= 75)  1797    320        3 (1263/1614-1521/1723) -38.7 W34-W51 N26-N39
      3  458 1644 (N55E39 L=351)   991     85        1 ( 426/1592- 486/1682) -42.0 E46-E34 N50-N59
      4 1730 1586 (N20W71 L=101)   992     85        1 (1684/1562-1771/1602)  -3.2 W64-W82 N18-N20
      5  762 1535 (N38E07 L= 24)  1253    253        1 ( 667/1459- 865/1581) -40.2 E13-W00 N31-N43
      6  846 1322 (N23E05 L= 25)   385     67        2 ( 811/1297- 884/1347)  11.3 E07-E02 N21-N24
      7  955 1223 (N15E00 L= 30)   504    127        0 ( 898/1191-1015/1264) -60.3 E02-W02 N11-N18
      8  278 1167 (N24E50 L=341)   974    187        2 ( 234/1115- 327/1236) -45.0 E55-E44 N20-N29
      9  378 1116 (N19E41 L=350)   762    113        3 ( 339/1062- 422/1168)  63.4 E44-E36 N16-N22
     10 1781 1061 (S11W51 L= 81)  1093    207        0 (1755/ 957-1806/1151)  51.3 W47-W55 S16-S06
     11  294 1025 (N15E49 L=341)  1137    195        5 ( 239/ 951- 355/1080) -41.6 E54-E44 N10-N18
     12  142 1001 (N16E68 L=323)  1177    108        2 ( 114/ 959- 191/1049) -31.0 E72-E60 N13-N19
     13 1609  839 (S21W34 L= 64)   460    128        1 (1558/ 800-1676/ 875) -37.9 W30-W39 S24-S17
     14  342  710 (S04E52 L=338)   318     45        0 ( 323/ 677- 357/ 736)  71.6 E53-E50 S05-S02
     15  841  554 (S22E22 L=  9)   375    164        0 ( 800/ 490- 876/ 627)  37.9 E26-E17 S25-S18
     16 1115  558 (S28E04 L= 26)   308    103        0 (1077/ 517-1142/ 595)  31.0 E07-E02 S29-S24
     17  378  505 (S17E59 L=332)  2255    260       10 ( 315/ 421- 425/ 592)  57.5 E65-E50 S21-S10
     18  579  297 (S34E56 L=334)   445     57        2 ( 551/ 275- 610/ 319)   0.0 E59-E50 S35-S32
  

    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