Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  833 1809 (N37E42 L=103)  2369    409        6 ( 767/1716- 887/1896) -33.0 E54-E34 N29-N42
      2  414 1695 (N18E67 L= 77)   691    100        1 ( 382/1673- 444/1714)  45.0 E69-E64 N16-N20
      3  430 1590 (N12E57 L= 87)  2379    397        7 ( 349/1506- 505/1697)  82.4 E64-E51 N05-N20
      4 1118 1438 (N20E06 L=138)   831    153        2 (1069/1384-1181/1481) -20.6 E09-E01 N16-N22
      5 1408  643 (S15W34 L=179)   943    190        0 (1293/ 597-1488/ 695)  48.0 W28-W37 S21-S10
      6 1722  533 (S10W64 L=209)   838    104        1 (1692/ 496-1747/ 574)  78.7 W62-W67 S12-S07
      7 1627  490 (S16W58 L=202)   898    195        3 (1570/ 455-1692/ 565)  68.2 W53-W62 S19-S09
      8 1555  426 (S22W56 L=200)   634    125        1 (1526/ 394-1585/ 463)  -6.3 W51-W60 S22-S19
      9 1027  373 (S43W25 L=169)  1326    360        1 ( 857/ 342-1205/ 411)  23.7 W12-W37 S49-S37
  

    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