Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  898 1754 (N46E27 L=146)   391     76        2 ( 864/1738- 930/1778)  21.8 E30-E24 N43-N48
      2 1035 1677 (N42E13 L=160)   519     81        2 ( 994/1657-1077/1702)   0.0 E16-E09 N40-N44
      3  952 1633 (N37E17 L=156)   266     43        1 ( 935/1612- 974/1656) -26.6 E19-E15 N35-N38
      4  309 1540 (N20E66 L=107)   982    178        1 ( 290/1496- 326/1589) -63.4 E68-E62 N17-N23
      5  287 1447 (N14E63 L=110)  4916    587       12 ( 204/1309- 390/1537)  58.0 E75-E54 N04-N20
      6  430 1347 (N10E45 L=128)   762    189        2 ( 398/1304- 488/1394) -18.4 E48-E39 N07-N12
      7  918 1290 (N14E11 L=162)   229     67        1 ( 889/1258- 940/1320) -36.9 E13-E09 N12-N15
      8  885 1099 (N02E09 L=164)   328     91        2 ( 850/1075- 927/1130) -26.6 E11-E06 N00-N03
      9  428 1005 (S11E38 L=135)  1192    218        3 ( 399/ 903- 479/1113) -63.4 E41-E33 S16-S04
     10 1095 1012 (N00W04 L=177)   303    103        1 (1044/ 994-1149/1038)  26.6 W01-W07 S01-N02
     11 1043  894 (S08W03 L=176)   216     89        0 (1020/ 841-1062/ 937)  81.9 W02-W03 S11-S04
     12 1441  405 (S32W43 L=216)  1179    322        0 (1320/ 375-1565/ 451)   8.5 W32-W52 S34-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