Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1721 1446 (N25W57 L=233)   627     69        0 (1689/1401-1749/1478) -68.2 W54-W58 N22-N27
      2  925 1278 (N21E05 L=171)   203     79        0 ( 884/1262- 969/1292) -11.3 E07-E02 N19-N21
      3 1269 1088 (N07W16 L=192)   312    126        0 (1250/1025-1289/1154) -82.9 W14-W16 N03-N11
      4  661  831 (S06E24 L=152)   380    111        2 ( 644/ 779- 676/ 888)  74.1 E25-E22 S09-S02
      5 1721  768 (S16W48 L=225)   452     76        0 (1706/ 731-1734/ 815)  90.0 W46-W49 S18-S13
      6  545  608 (S20E35 L=141)   260     89        0 ( 524/ 580- 569/ 637)  38.7 E37-E32 S22-S18
      7  747  581 (S22E21 L=155)   615    178        0 ( 671/ 545- 830/ 608)  14.0 E26-E14 S24-S20
      8 1134  317 (S45W05 L=181)   288     70        0 (1097/ 307-1174/ 326)  -8.1 W01-W08 S46-S44
  

    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