Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1254 1810 (N52E14 L=220)   498     63        0 (1233/1780-1283/1833)  90.0 E15-E12 N48-N54
      2  719 1737 (N27E44 L=189)   523     89        0 ( 691/1689- 741/1780) -60.3 E47-E42 N23-N30
      3 1429 1722 (N51W05 L=239)   694     88        0 (1379/1693-1480/1748)  16.7 W00-W10 N49-N53
      4  413 1660 (N16E65 L=169)  1236    181        1 ( 381/1575- 447/1741) -45.0 E72-E58 N09-N21
      5  387 1393 (S01E50 L=183)   423    104        0 ( 347/1348- 418/1419)  -7.1 E54-E46 S03-N00
      6  239 1130 (S19E56 L=177)  1262    139        3 ( 180/1092- 311/1164)  21.0 E62-E49 S22-S16
      7  844 1091 (S07E12 L=222)   361    161        0 ( 811/1027- 907/1140)  76.0 E12-E09 S11-S03
      8 1525 1082 (N11W27 L=260)   294     76        0 (1496/1043-1550/1110) -21.8 W24-W29 N09-N12
      9 1221  970 (S04W12 L=246)   503    147        1 (1189/ 905-1265/1036) -45.0 W09-W16 S06-N00
     10  604  868 (S26E21 L=213)   552    168        0 ( 571/ 794- 661/ 921)  77.5 E22-E18 S30-S21
     11 1219  726 (S17W20 L=253)   217     73        0 (1190/ 710-1252/ 753)   0.0 W17-W21 S17-S15
     12 1429  685 (S13W33 L=267)   355    115        0 (1368/ 668-1494/ 701)  26.6 W29-W37 S14-S10
     13  540  528 (S50E20 L=214)   291     49        0 ( 520/ 496- 560/ 557) -45.0 E22-E17 S51-S47
     14  685  473 (S48E05 L=228)   761    147        2 ( 608/ 450- 764/ 502)  12.1 E12-W02 S50-S45
     15 1358  437 (S28W42 L=275)   460     83        0 (1318/ 421-1406/ 462)  18.4 W37-W44 S29-S25
     16 1534  333 (S27W64 L=298)  1509    154        2 (1486/ 289-1574/ 381) -74.1 W60-W70 S30-S23
     17 1172  218 (S48W50 L=283)  1066    152        2 (1118/ 172-1226/ 264)  -2.9 W40-W61 S49-S45
  

    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