Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1397 1665 (N35W39 L=140)  2127    553        1 (1173/1606-1639/1721)  -2.5 W19-W64 N31-N38
      2  405 1437 (N33E41 L= 61)   294     49        1 ( 380/1412- 426/1458) -36.9 E43-E38 N31-N34
      3  259 1291 (N25E54 L= 48)   485     85        3 ( 236/1262- 289/1320) -29.7 E56-E49 N23-N27
      4  237 1078 (N12E54 L= 47)  1531    181        6 ( 204/ 985- 260/1165) -79.7 E58-E51 N06-N17
      5 1643  846 (S18W39 L=140)  2530    678        3 (1509/ 576-1784/1139)  54.3 W28-W53 S34-S02
      6  158 1013 (N09E64 L= 37)  1676    164        5 ( 118/ 943- 208/1073)  41.2 E70-E57 N05-N12
      7 1841  808 (S23W59 L=161)   704    115        0 (1794/ 760-1882/ 852) -28.3 W53-W66 S26-S19
      8 1099  750 (S18W01 L=103)   427    105        0 (1061/ 711-1149/ 786)  31.0 E01-W04 S19-S15
      9  407  737 (S10E44 L= 57)   473     90        5 ( 365/ 706- 444/ 772)  15.9 E47-E40 S12-S08
     10 1127  602 (S27W01 L=103)   515    168        2 (1091/ 536-1154/ 693) -90.0 W00-W02 S31-S21
     11  601  617 (S20E33 L= 68)   711    265        2 ( 522/ 559- 701/ 659)  11.3 E40-E25 S22-S16
     12 1276  577 (S31W11 L=113)   656    171        0 (1181/ 523-1330/ 621)  18.4 W03-W15 S33-S28
  

    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