Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  983 1680 (N31E23 L=328)   377    115        0 ( 929/1641-1032/1697)   0.0 E27-E18 N29-N31
      2 1494 1619 (N45W13 L=  5)   824    129        1 (1454/1582-1526/1656)  58.0 W10-W16 N41-N49
      3 1560 1596 (N46W20 L= 11)   474     83        1 (1542/1561-1586/1628)  90.0 W18-W22 N42-N48
      4 1723 1378 (N36W37 L= 28)  3674    503        6 (1601/1190-1799/1538) -41.8 W24-W45 N25-N44
      5  411 1478 (N04E52 L=299)   548    143        1 ( 371/1435- 439/1506) -20.6 E57-E49 N01-N05
      6 1076 1244 (N07E03 L=348)   244     77        1 (1059/1203-1094/1288) -59.0 E04-E01 N04-N09
      7 1571 1154 (N16W28 L= 20)   486    123        0 (1524/1090-1599/1217) -32.0 W23-W31 N13-N18
      8 1488 1000 (N05W27 L= 18)   787    166        1 (1431/ 951-1561/1063) -11.3 W22-W32 N03-N07
      9 1879  974 (N17W59 L= 51)   589     56        0 (1863/ 945-1894/1001) -45.0 W56-W61 N15-N19
     10 1734  929 (N09W46 L= 38)  1153    259        1 (1644/ 900-1816/ 948)  21.8 W39-W54 N05-N12
     11 1352  814 (S09W24 L= 16)   220     89        0 (1334/ 764-1364/ 860) -59.0 W22-W25 S11-S06
     12 1251  765 (S14W20 L= 12)   339    131        0 (1184/ 749-1318/ 787)  26.6 W16-W24 S16-S12
     13  752  738 (S30E08 L=343)   269    112        0 ( 704/ 708- 789/ 782)  60.3 E10-E06 S33-S26
     14 1724  724 (S02W53 L= 45)   649    106        1 (1675/ 684-1761/ 764)  56.3 W50-W55 S05-N01
     15 1093  673 (S24W14 L=  6)   582    120        0 (1021/ 646-1146/ 700)  20.6 W09-W17 S25-S21
     16 1238  562 (S26W27 L= 19)   543    139        1 (1192/ 510-1288/ 614) -16.7 W22-W32 S27-S23
     17 1508  558 (S17W44 L= 36)   294     60        0 (1475/ 532-1544/ 583)  68.2 W42-W45 S19-S14
  

    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