Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  684 1526 (N43E18 L=156)   315     89        0 ( 652/1498- 712/1559) -50.2 E21-E15 N40-N46
      2 1618 1417 (N20W48 L=223)  3102    458        4 (1455/1349-1767/1485) -32.5 W37-W60 N12-N26
      3 1149 1269 (N19W12 L=186)   416    164        0 (1078/1227-1214/1290)   0.0 W06-W15 N16-N20
      4 1527 1222 (N10W35 L=210)   340    122        0 (1468/1185-1576/1254) -35.5 W31-W38 N06-N12
      5 1878 1006 (S10W63 L=238)  4055    657       10 (1796/ 863-1943/1134)  53.7 W54-W73 S19-S01
      6  540 1084 (N16E30 L=144)   258     26        0 ( 517/1067- 562/1101)  14.0 E32-E28 N14-N16
      7 1759 1054 (S05W51 L=225)   427    127        2 (1723/1012-1793/1089)  14.0 W46-W54 S06-S02
      8 1769  980 (S09W51 L=225)   624    175        3 (1734/ 945-1811/1006)  -8.1 W47-W54 S11-S07
      9 1693  870 (S14W43 L=217)   373     95        0 (1682/ 817-1706/ 918)  80.5 W41-W43 S17-S11
     10 1037  850 (S04E02 L=173)   626     69        3 ( 979/ 817-1117/ 876)   0.0 E05-W03 S05-S02
     11 1121  723 (S13W02 L=176)  1025    225        4 (1079/ 613-1157/ 829)  63.4 E02-W05 S20-S07
     12 1032  508 (S26E08 L=167)   278     79        1 ( 992/ 477-1049/ 542)  36.9 E10-E05 S28-S24
     13  346  451 (S23E63 L=111)   848    197        2 ( 320/ 425- 372/ 478)  18.4 E68-E59 S25-S21
  

    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