Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  573 1798 (N45E55 L=265)  1340    106        2 ( 518/1769- 619/1821)  -9.5 E72-E46 N42-N47
      2  809 1732 (N40E25 L=294)   326     91        0 ( 774/1706- 855/1750)  18.4 E27-E21 N37-N41
      3 1654 1557 (N37W46 L=  6)  1514    141        1 (1594/1500-1708/1614) -26.6 W40-W50 N32-N41
      4  898 1345 (N13E11 L=309)   739    306        3 ( 848/1248- 979/1484) -49.8 E15-E04 N07-N22
      5 1264 1053 (S01W14 L=334)   831    214        1 (1206/ 966-1345/1116)  48.8 W11-W18 S06-N03
      6 1684 1009 (N01W43 L=  3)   489     78        2 (1641/ 986-1734/1029)   8.1 W39-W47 N00-N02
      7 1580  970 (S03W36 L=356)   276     88        0 (1552/ 926-1603/1005)  68.2 W34-W37 S05-N00
      8  881  558 (S35E06 L=314)  1451    362        2 ( 736/ 484- 999/ 691)  43.3 E16-W01 S41-S25
      9  317  603 (S36E56 L=264)   808    156        0 ( 305/ 526- 333/ 680)  57.5 E60-E51 S41-S30
     10  493  613 (S35E37 L=283)   408     76        2 ( 452/ 598- 532/ 629)  -7.1 E40-E32 S35-S33
  

    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