Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1358 1567 (N37W05 L= 82)  2451    536        1 (1278/1342-1471/1782) -50.6 E07-W17 N25-N53
      2  360 1636 (N14E70 L=  7)  3413    144        3 ( 310/1579- 421/1694) -81.9 E82-E61 N10-N18
      3  709 1527 (N15E33 L= 44)   437    104        1 ( 668/1496- 740/1553)  11.3 E35-E29 N13-N16
      4 1157 1497 (N25E06 L= 71)   428    118        1 (1104/1456-1204/1533) -14.0 E09-E01 N24-N26
      5  449 1443 (N04E47 L= 29)   340     55        0 ( 414/1428- 482/1460)  21.8 E50-E45 N02-N05
      6 1602 1305 (N27W28 L=105)   295     67        0 (1581/1277-1631/1338) -26.6 W26-W30 N26-N29
      7  392 1250 (S08E45 L= 32)   335     85        0 ( 377/1202- 404/1289) -68.2 E47-E43 S10-S05
      8 1684 1204 (N24W36 L=113)   571    138        0 (1634/1175-1720/1242)  29.7 W32-W39 N21-N27
      9  253 1188 (S14E56 L= 21)   327     57        2 ( 232/1162- 269/1218) -45.0 E58-E54 S15-S12
     10  939 1103 (S03E07 L= 70)   236     31        0 ( 918/1078- 958/1129) -56.3 E07-E05 S04-S01
     11  103 1041 (S25E77 L=360)  1063     76        1 (  80/1011- 142/1063)   7.1 E86-E68 S27-S23
     12  911  958 (S12E04 L= 72)   522    169        1 ( 832/ 918- 960/1013)  63.4 E08-E02 S15-S07
     13  496  985 (S21E31 L= 46)   265     76        0 ( 460/ 960- 524/1007)  -9.5 E34-E28 S21-S19
     14 1574  788 (S02W39 L=116)  1506    288        4 (1509/ 650-1604/ 915) -49.8 W34-W45 S09-N04
     15 1764  836 (N06W52 L=129)   590    103        0 (1739/ 786-1788/ 887) -71.6 W50-W53 N03-N09
     16 1798  764 (N03W59 L=135)  1113    179        2 (1757/ 696-1829/ 830)  90.0 W55-W62 S01-N07
  

    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