Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1633 1667 (N21W65 L= 47)   917    177        2 (1595/1638-1674/1689) -29.7 W61-W70 N19-N23
      2 1846 1326 (S03W67 L= 48)  1456    197        3 (1815/1267-1878/1399)  66.8 W62-W70 S06-N01
      3 1738 1168 (S07W48 L= 30)   348    105        0 (1709/1117-1754/1216)  80.5 W46-W50 S09-S03
      4  728 1105 (N18E15 L=327)   345    135        1 ( 670/1080- 787/1148)  -7.1 E18-E10 N16-N19
      5  329  816 (N11E50 L=292)  2791    823        6 ( 122/ 769- 566/ 866) -15.0 E74-E32 N03-N15
      6 1128  613 (S20E05 L=336)   322    100        1 (1101/ 558-1151/ 656)  45.0 E08-E03 S22-S17
      7  682  355 (S24E43 L=299)  1941    408        3 ( 637/ 230- 719/ 463)  60.9 E50-E36 S34-S16
      8 1031  268 (S41E26 L=315)   371     47        0 (1013/ 239-1053/ 293)  45.0 E28-E23 S43-S39
  

    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