Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  661 1404 (N33E14 L=178)  1101    194        3 ( 549/1336- 737/1446) -34.7 E22-E09 N28-N37
      2 1553 1329 (N08W39 L=232)   570    134        1 (1519/1292-1609/1360)  -8.1 W36-W43 N05-N09
      3  157 1061 (N23E61 L=132)  1159    197        3 ( 118/ 978- 183/1118) -66.0 E64-E57 N17-N26
      4  858  787 (S07E15 L=178)  3313    699        7 ( 677/ 574-1000/ 926) -10.0 E25-E07 S20-N03
      5 1263  873 (S11W10 L=203)   394    128        2 (1207/ 858-1340/ 891) -23.2 W06-W14 S12-S09
      6  171  809 (N08E66 L=126)   792    109        1 ( 152/ 744- 190/ 852) -90.0 E68-E64 N03-N10
      7  678  640 (S11E29 L=163)   324    106        0 ( 641/ 603- 708/ 688)  26.6 E32-E26 S12-S09
      8 1070  564 (S25E08 L=184)   500    137        1 (1006/ 529-1151/ 614)  11.3 E12-E02 S25-S23
      9  824  535 (S21E24 L=169)   434    117        0 ( 773/ 502- 878/ 567)   6.3 E28-E19 S21-S19
     10  515  308 (S29E57 L=135)   517     85        1 ( 492/ 275- 537/ 343)  71.6 E58-E55 S31-S25
  

    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