Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  778 1904 (N43E61 L= 18)   620     52        0 ( 759/1879- 795/1927) -18.4 E69-E54 N40-N46
      2  712 1656 (N22E39 L= 39)  1626    341        4 ( 652/1543- 816/1782) -24.8 E53-E27 N16-N29
      3 1061 1646 (N31E17 L= 62)   535    158        1 ( 995/1584-1110/1692) -20.0 E22-E11 N27-N32
      4  156 1026 (S26E67 L= 11)  1182    123        5 ( 112/ 985- 206/1054)  10.6 E75-E59 S28-S24
      5 1684 1007 (N11W40 L=119)   255     59        1 (1658/ 992-1723/1026)  21.8 W38-W43 N10-N12
      6 1728  648 (S06W58 L=136)   577     70        2 (1700/ 616-1752/ 679) -21.8 W55-W61 S07-S04
      7 1487  455 (S24W50 L=128)   447    100        0 (1443/ 421-1518/ 498) -81.9 W48-W51 S27-S20
      8 1640  466 (S18W62 L=140)   531     88        0 (1605/ 445-1682/ 482)  29.7 W57-W66 S19-S15
      9 1581  368 (S25W66 L=144)  2054    313        3 (1427/ 329-1677/ 420)  43.0 W55-W75 S32-S18
     10 1350  402 (S32W44 L=122)   388    142        0 (1283/ 393-1416/ 417)  40.6 W39-W47 S34-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