Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1586 1667 (N24W65 L=275)   537     80        1 (1546/1645-1621/1689) -59.0 W62-W67 N21-N26
      2 1043 1545 (N35W18 L=229)  1765    519        6 ( 835/1502-1296/1591) -27.3 W03-W32 N25-N41
      3  438 1527 (N53E34 L=177)   320     75        0 ( 427/1490- 450/1566) -63.4 E35-E31 N50-N56
      4 1731 1470 (N09W63 L=273)   408     64        0 (1706/1448-1749/1507)   7.1 W59-W67 N07-N10
      5 1405 1336 (N12W31 L=242)   998    378        1 (1228/1256-1529/1391) -39.1 W20-W37 N04-N18
      6  904 1255 (N22E01 L=210)   753    223        1 ( 812/1222-1015/1294) -14.0 E06-W06 N20-N24
      7  392 1185 (N31E38 L=172)   319     56        1 ( 373/1147- 411/1212)  76.0 E39-E36 N28-N32
      8 1890 1175 (S12W66 L=277)   934    106        1 (1855/1144-1928/1205) -26.6 W62-W73 S14-S09
      9  892 1132 (N15E05 L=206)   380    205        0 ( 852/1066- 979/1172)   6.3 E08-W01 N12-N17
     10  337  945 (N18E47 L=164)   538    146        3 ( 278/ 928- 413/ 963) -14.0 E52-E40 N15-N19
     11 1673  860 (S22W36 L=246)   447    153        0 (1615/ 826-1723/ 918)   5.7 W30-W41 S23-S20
     12  224  742 (N08E64 L=146)  1463    129        4 ( 196/ 687- 260/ 803)  66.8 E68-E59 N04-N11
     13  310  707 (N05E56 L=155)   470     84        0 ( 278/ 672- 338/ 743) -59.0 E57-E53 N02-N07
     14  338  583 (S03E59 L=151)   380     81        1 ( 324/ 544- 357/ 626)  45.0 E61-E56 S05-N00
     15  922  520 (S20E20 L=190)   649    211        2 ( 901/ 423- 943/ 615)  57.5 E24-E17 S25-S13
     16 1785  584 (S46W49 L=260)  1077     92        1 (1741/ 559-1826/ 610)  -7.6 W42-W57 S48-S43
     17 1591  459 (S47W25 L=235)   747    155        0 (1521/ 422-1659/ 498)  -3.8 W17-W32 S49-S45
     18 1253  253 (S51E12 L=198)  1024    189        2 (1171/ 223-1347/ 296) -14.0 E19-E03 S53-S48
  

    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