Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1117 1789 (N56E08 L=167)   494     41        1 (1092/1766-1142/1808) -21.8 E10-E05 N53-N58
      2 1032 1779 (N53E16 L=160)   363     57        1 (1010/1749-1056/1799)  90.0 E18-E13 N50-N55
      3 1231 1769 (N56W05 L=180)   494     75        0 (1199/1751-1267/1787)  20.6 W01-W09 N54-N58
      4 1301 1760 (N57W13 L=189)   408     75        0 (1280/1723-1325/1780) -71.6 W10-W16 N53-N59
      5  971 1718 (N46E18 L=157)   340     82        1 ( 954/1676- 988/1753) -54.5 E20-E14 N42-N49
      6 1411 1730 (N56W26 L=202)   468     57        1 (1379/1707-1445/1746)   0.0 W22-W30 N54-N58
      7  672 1582 (N31E35 L=140)   364     71        1 ( 646/1550- 693/1615) -78.7 E37-E33 N28-N33
      8 1325 1469 (N33W15 L=190)   258     56        1 (1294/1454-1358/1485)   0.0 W12-W17 N31-N34
      9  278 1446 (N16E66 L=109)   948     93        3 ( 247/1407- 314/1482)  68.2 E70-E63 N13-N18
     10 1717 1336 (N29W49 L=224)   547    132        0 (1655/1299-1768/1364)  21.0 W41-W55 N25-N31
     11  611  990 (S07E25 L=150)   350     65        2 ( 587/ 952- 635/1022) -36.9 E27-E23 S09-S04
     12 1336  485 (S29W31 L=206)   271     36        1 (1322/ 462-1352/ 516) -56.3 W29-W32 S30-S27
  

    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