Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  753 1837 (N36E50 L=105)  2134    367        6 ( 716/1768- 796/1906) -34.5 E59-E43 N30-N42
      2  340 1593 (N10E68 L= 87)  1719    226        1 ( 304/1538- 376/1642)  90.0 E72-E64 N06-N13
      3  969 1504 (N20E17 L=137)   757    158        2 ( 920/1451-1036/1549) -16.7 E21-E11 N16-N22
      4  316  771 (S39E43 L=111)   372     72        0 ( 292/ 751- 336/ 800)  51.3 E45-E41 S41-S36
      5 1839  771 (N06W63 L=218)   506    181        0 (1794/ 757-1878/ 781)   9.5 W57-W69 N04-N07
      6 1270  706 (S16W23 L=178)   687    205        1 (1159/ 659-1367/ 767)  56.3 W18-W27 S21-S09
      7 1414  519 (S22W40 L=195)   643    294        1 (1333/ 458-1491/ 619)  -9.5 W30-W48 S22-S18
      8 1665  589 (S10W55 L=209)   468     85        0 (1639/ 548-1694/ 617)  63.4 W53-W56 S11-S07
      9  922  427 (S44W15 L=169)   667    158        0 ( 847/ 398-1008/ 462)  13.0 W08-W21 S45-S41
  

    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