Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1310 1738 (N40W49 L=241)   537     96        0 (1264/1724-1360/1751) -39.8 W46-W53 N37-N42
      2 1616 1662 (N25W69 L=260)   654    109        1 (1582/1612-1653/1701)  60.3 W66-W71 N20-N27
      3  992 1659 (N46W20 L=211)   272     55        0 ( 959/1648-1028/1669) -18.4 W17-W23 N44-N47
      4 1551 1587 (N23W55 L=246)   624    128        0 (1494/1548-1592/1627) -74.1 W52-W57 N20-N27
      5 1488 1460 (N19W42 L=234)  1529    460        0 (1299/1332-1649/1540) -51.3 W32-W48 N06-N27
      6 1857 1084 (S14W59 L=251)   426     37        2 (1836/1064-1879/1109)  63.4 W56-W62 S14-S12
      7  146  963 (N21E71 L=120)   861    125        1 ( 116/ 912- 189/1012) -26.6 E77-E64 N17-N23
      8 1902  915 (S26W64 L=256)  1254     93        2 (1871/ 869-1933/ 952)  16.7 W58-W71 S27-S23
      9 1292  882 (S08W12 L=203)   401    157        0 (1271/ 814-1319/ 946)  76.0 W10-W14 S12-S03
     10 1633  528 (S42W31 L=222)  1443    226        1 (1539/ 461-1704/ 600)  10.0 W21-W38 S45-S38
     11  450  425 (S16E58 L=133)   779    145        2 ( 424/ 364- 490/ 492)  68.2 E60-E54 S21-S11
     12  724  256 (S34E48 L=143)   418     50        0 ( 699/ 241- 750/ 277) -45.0 E50-E45 S35-S31
  

    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