Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1305 1539 (N35W05 L= 77)  1273    242        2 (1244/1432-1355/1630) -42.3 E02-W10 N28-N40
      2  391 1490 (N09E56 L= 16)   774    132        1 ( 364/1441- 416/1549) -60.3 E58-E51 N06-N13
      3  690 1457 (N14E31 L= 41)   286     83        0 ( 668/1418- 722/1489)  78.7 E32-E30 N11-N16
      4  374 1342 (N00E50 L= 22)   304     36        2 ( 349/1321- 401/1363)   0.0 E52-E47 N00-N01
      5 1613 1163 (N19W34 L=106)  1126    245        3 (1561/1050-1669/1259)  81.9 W30-W37 N12-N26
      6 1834 1042 (N19W56 L=128)  1617    186        2 (1797/ 957-1867/1138) -52.1 W51-W61 N14-N23
      7 1320 1031 (N04W17 L= 89)   330     66        0 (1287/1007-1349/1059)  63.4 W15-W18 N02-N06
      8 1006  744 (S20W06 L= 78)   406    158        0 ( 944/ 716-1083/ 772)  40.6 W02-W09 S22-S16
      9  532  704 (S34E25 L= 47)   942    173        0 ( 460/ 674- 618/ 748)   4.4 E31-E18 S35-S31
     10 1599  555 (S15W52 L=124)   444    103        1 (1567/ 518-1640/ 592)  63.4 W49-W53 S17-S11
  

    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