Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  813 1794 (N35E42 L=299)   538     66        0 ( 794/1761- 833/1828) -56.3 E44-E38 N31-N37
      2  623 1743 (N26E52 L=289)   844    189        3 ( 542/1703- 708/1775)  48.4 E57-E45 N21-N30
      3  714 1712 (N26E42 L=299)   321     53        1 ( 699/1679- 734/1742) -59.0 E43-E40 N23-N28
      4 1013 1524 (N21E15 L=326)   729    234        0 ( 919/1454-1100/1552)   4.8 E19-E07 N18-N23
      5  784 1434 (N10E25 L=316)   607    243        1 ( 749/1333- 819/1535) -45.0 E30-E19 N04-N15
      6  607 1363 (N02E33 L=308)   294     69        1 ( 588/1333- 636/1397) -36.9 E35-E30 N00-N03
      7 1711 1328 (N31W37 L= 18)  1265    195        3 (1645/1256-1772/1393) -17.1 W30-W44 N27-N33
      8 1430 1160 (N11W19 L=  0)   631    171        0 (1409/1076-1463/1236) -78.7 W17-W21 N06-N16
      9  461 1167 (S12E38 L=303)   469    148        0 ( 402/1148- 533/1189)  21.8 E42-E32 S14-S10
     10 1626 1071 (N12W34 L= 15)   767    232        0 (1574/ 974-1680/1166) -33.7 W28-W40 N08-N16
     11  207 1103 (S21E61 L=280)   686    112        0 ( 160/1056- 235/1148)  36.9 E66-E57 S23-S17
     12  654  966 (S20E20 L=321)   390    103        0 ( 604/ 948- 703/1000)  45.0 E23-E17 S21-S16
     13  432  904 (S28E35 L=305)   428    119        0 ( 374/ 883- 490/ 933)  21.8 E40-E30 S30-S25
     14  601  727 (S35E19 L=322)   823    192        2 ( 532/ 685- 675/ 766)  38.7 E24-E14 S38-S30
     15 1736  531 (S12W67 L= 48)   567    113        0 (1689/ 494-1782/ 567)  54.5 W64-W70 S15-S08
     16 1356  485 (S27W39 L= 20)  1274    227        3 (1263/ 441-1457/ 543)  54.0 W34-W42 S32-S20
     17 1521  319 (S30W66 L= 47)   638    120        1 (1468/ 292-1582/ 344)  49.4 W61-W70 S33-S26
     18 1128  328 (S44W35 L= 16)   398     56        0 (1096/ 313-1165/ 342)  18.4 W31-W37 S45-S42
  

    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