Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  663 1738 (N26E48 L= 31)   580     90        2 ( 645/1701- 682/1783) -40.6 E53-E45 N23-N29
      2 1059 1645 (N30E17 L= 62)   548    150        0 (1004/1581-1104/1694) -24.4 E22-E11 N27-N32
      3  780 1571 (N18E30 L= 49)   275     65        0 ( 757/1548- 809/1600) -21.8 E32-E27 N16-N19
      4  154 1029 (S26E67 L= 11)  1244    102        6 ( 111/ 993- 203/1055)   0.0 E76-E59 S28-S24
      5 1726  648 (S06W57 L=136)   482     45        2 (1703/ 628-1750/ 672)  56.3 W54-W60 S07-S04
      6 1484  458 (S24W49 L=128)   446     94        0 (1443/ 427-1513/ 499)  90.0 W48-W50 S26-S20
      7 1648  469 (S17W62 L=141)   471     80        0 (1618/ 448-1681/ 484)  45.0 W60-W66 S19-S15
      8 1619  377 (S23W69 L=148)  2127    149        3 (1577/ 327-1681/ 422)  60.3 W64-W77 S25-S18
      9 1344  403 (S32W43 L=122)   354    110        0 (1287/ 394-1400/ 412)  33.7 W39-W47 S33-S29
     10 1386  298 (S36W56 L=135)   347     79        0 (1359/ 271-1415/ 332) -68.2 W54-W58 S38-S33
  

    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