Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  172 1125 (S13E63 L=259)   754    168        1 ( 133/1086- 216/1151)  26.6 E67-E58 S15-S10
      2 1777 1056 (N11W48 L= 11)   700    205        0 (1728/ 953-1795/1128) -52.1 W43-W51 N05-N14
      3  341  983 (S20E44 L=278)  2102    683        0 ( 158/ 834- 571/1080)  25.2 E61-E27 S30-S13
      4 1910  865 (N05W69 L= 32)  1335    231        2 (1873/ 820-1945/ 898)  36.9 W64-W75 N01-N07
      5  836  725 (S28E06 L=317)   224     48        0 ( 809/ 706- 864/ 749)  45.0 E07-E04 S29-S25
      6  367  639 (S41E43 L=280)   909    231        3 ( 318/ 577- 417/ 705)  42.3 E48-E37 S46-S36
      7 1504  642 (S19W39 L=  1)   353    124        0 (1443/ 626-1554/ 677)  35.5 W34-W41 S20-S15
      8 1786  621 (S12W64 L= 26)  1316    342        2 (1708/ 564-1841/ 660)   6.7 W54-W71 S14-S09
      9 1568  396 (S31W55 L= 18)  3214    489        5 (1432/ 357-1739/ 438)  28.4 W45-W69 S35-S22
  

    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