Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1215 1786 (N47E15 L=107)   409     51        2 (1191/1757-1238/1813) -38.7 E18-E13 N44-N49
      2  650 1719 (N24E48 L= 74)   436     76        1 ( 634/1685- 667/1763) -45.0 E52-E44 N21-N27
      3  508 1642 (N16E54 L= 68)  1813    297        3 ( 418/1531- 598/1736)  86.6 E59-E50 N07-N24
      4 1537 1686 (N53W17 L=139)   462     46        0 (1507/1668-1569/1705)  41.2 W13-W21 N50-N57
      5  622 1505 (N10E38 L= 84)   511     77        2 ( 588/1462- 656/1543) -51.3 E40-E34 N07-N12
      6 1335 1519 (N31W04 L=126)   233     32        0 (1315/1494-1355/1540)  90.0 W02-W05 N28-N32
      7 1823 1210 (N29W50 L=172)   441    103        0 (1781/1172-1861/1244)  -5.2 W44-W55 N28-N30
      8 1608  485 (S17W57 L=179)   331     62        0 (1581/ 458-1634/ 511)  90.0 W56-W58 S19-S15
      9 1187  429 (S35W32 L=154)   640    141        0 (1154/ 362-1219/ 493) -26.6 W26-W38 S38-S31
     10 1356  399 (S31W45 L=167)   324     68        1 (1330/ 372-1389/ 421)  90.0 W43-W46 S32-S28
  

    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