Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  618 1790 (N54E38 L= 38)   613     95        0 ( 589/1767- 647/1809)  15.9 E43-E33 N51-N55
      2 1689 1563 (N29W57 L=134)  1796    222        3 (1594/1491-1776/1637) -39.8 W49-W65 N24-N35
      3 1039 1509 (N28W03 L= 80)   328    138        0 (1006/1461-1103/1542)  39.8 W00-W07 N24-N29
      4 1500 1471 (N24W35 L=112)   484    203        0 (1438/1420-1531/1530) -60.3 W31-W37 N20-N27
      5  278 1286 (N18E52 L= 24)   478     83        2 ( 237/1271- 321/1301)  -6.3 E57-E48 N16-N18
      6  183 1232 (N15E63 L= 14)   768    112        1 ( 149/1181- 213/1275)  35.5 E66-E59 N12-N17
      7 1493  946 (S09W29 L=106)   608    233        0 (1414/ 865-1577/1018)  33.7 W23-W35 S13-S04
      8 1425  687 (S24W25 L=102)   384    124        1 (1406/ 629-1448/ 749) -82.9 W24-W26 S28-S20
      9 1754  653 (S27W55 L=132)   828    131        2 (1698/ 630-1820/ 673)   8.1 W49-W63 S28-S26
     10 1755  508 (S37W65 L=142)  1102    180        0 (1717/ 434-1782/ 580) -55.0 W58-W70 S42-S32
     11  291  525 (S29E66 L= 11)   594     65        1 ( 258/ 508- 328/ 540)   4.8 E72-E59 S29-S27
     12  689  288 (S50E38 L= 39)   498     54        1 ( 658/ 263- 712/ 322)  45.0 E42-E34 S51-S46
  

    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