Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1311 1643 (N45W21 L=262)   461     87        0 (1271/1617-1356/1664)  24.0 W17-W26 N42-N46
      2 1577 1602 (N43W48 L=289)  1885    262        2 (1451/1563-1687/1633)   0.0 W33-W64 N39-N44
      3  242 1444 (N24E70 L=171)   660     72        1 ( 217/1414- 257/1485) -51.3 E73-E66 N21-N26
      4  545 1375 (N22E35 L=206)   242     57        1 ( 521/1357- 577/1398)  26.6 E36-E32 N20-N23
      5 1587 1132 (N11W37 L=278)  2171    403        4 (1511/ 956-1638/1297) -76.6 W31-W40 N00-N21
      6  687 1228 (N13E22 L=219)   268     76        2 ( 668/1189- 704/1265)  90.0 E23-E21 N10-N15
      7  275 1179 (N08E54 L=187)   908    126        1 ( 253/1134- 316/1246) -50.2 E57-E49 N05-N11
      8  835 1197 (N12E13 L=228)   292    121        1 ( 801/1152- 881/1246) -50.2 E14-E09 N09-N15
      9 1054 1112 (N08W01 L=242)   238     68        1 (1033/1078-1079/1148) -63.4 W00-W03 N06-N10
     10  882  997 (N00E09 L=232)   804    216        1 ( 829/ 906- 928/1085)  90.0 E11-E05 S05-N05
     11  469  863 (S10E36 L=205)   826    240        0 ( 424/ 774- 542/ 963) -57.5 E39-E30 S15-S04
     12 1113  699 (S17W07 L=248)  1815    612        1 (1013/ 483-1213/ 940)  73.8 W01-W11 S32-S01
     13 1663  775 (S11W45 L=286)   297     81        0 (1648/ 735-1689/ 810)  78.7 W43-W46 S13-S08
     14 1756  689 (S16W56 L=297)  1408    174        4 (1711/ 626-1803/ 748) -45.0 W51-W62 S20-S12
     15 1672  645 (S19W49 L=290)   485     98        1 (1638/ 607-1710/ 697) -33.7 W44-W53 S21-S15
     16 1359  277 (S48W38 L=279)   371     27        1 (1332/ 262-1385/ 291)   9.5 W34-W40 S49-S46
     17  803  243 (S55E18 L=223)  1895    184        5 ( 717/ 216- 902/ 282)   5.2 E29-E07 S58-S51
     18 1063  259 (S51W09 L=250)   462     61        1 (1027/ 237-1093/ 273)   0.0 W05-W11 S53-S49
  

    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