Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  552 1127 (N23E27 L=318)   579    198        2 ( 467/1094- 646/1155) -24.8 E33-E20 N20-N26
      2 1552  988 (S08W30 L= 16)   443    130        0 (1514/ 941-1611/1013)   7.1 W27-W35 S09-S06
      3  860  873 (N02E13 L=333)   708    206        4 ( 838/ 769- 887/ 965)  68.2 E15-E10 S03-N07
      4 1693  814 (S23W38 L= 23)   866    251        1 (1644/ 691-1718/ 940)  60.3 W32-W41 S29-S15
      5 1395  839 (S12W17 L=  3)   723    183        0 (1333/ 756-1427/ 925)  48.4 W11-W20 S16-S07
      6 1366  592 (S27W10 L=356)  1167    341        1 (1205/ 470-1488/ 723)  20.0 E02-W20 S31-S22
      7 1723  664 (S34W40 L= 25)   571     86        0 (1682/ 632-1763/ 699)   8.1 W35-W43 S35-S32
      8 1018  489 (S25E14 L=331)   280     71        0 ( 989/ 464-1054/ 519)  18.4 E16-E10 S25-S23
      9  965  276 (S39E27 L=318)   490    166        1 ( 939/ 228- 995/ 339)  38.7 E31-E21 S43-S34
  

    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