Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  440 1676 (N27E59 L=336)  1226    166        4 ( 372/1618- 484/1735)  90.0 E63-E55 N22-N32
      2  502 1618 (N23E48 L=347)   471     73        0 ( 474/1584- 527/1652)  90.0 E50-E46 N20-N25
      3  517 1559 (N19E44 L=351)   324     44        2 ( 487/1541- 539/1579)  33.7 E45-E42 N18-N20
      4 1844 1142 (N16W56 L= 91)   760     99        1 (1819/1094-1865/1188) -45.0 W53-W58 N13-N18
      5  848  825 (S21E08 L= 27)   554    183        0 ( 767/ 794- 935/ 864)  21.8 E12-E02 S22-S17
      6  336  700 (S35E48 L=347)   289     43        0 ( 309/ 686- 367/ 711)   8.1 E51-E44 S36-S34
      7  288  628 (S40E57 L=338)   774     89        0 ( 274/ 578- 309/ 677)  41.2 E61-E51 S43-S36
      8 1391  501 (S31W36 L= 71)   422    120        0 (1360/ 445-1433/ 551)  90.0 W34-W38 S34-S26
      9  562  360 (S58E33 L=  2)   724     71        3 ( 527/ 335- 589/ 388)  20.0 E39-E28 S60-S54
  

    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