Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1650 1604 (N22W66 L=  3)  1239    179        1 (1579/1581-1731/1632) -29.7 W60-W76 N17-N25
      2  722 1598 (N51E06 L=290)   305     48        1 ( 706/1568- 740/1629)  68.2 E08-E04 N48-N53
      3 1053 1531 (N36W16 L=313)   513    123        1 ( 992/1499-1116/1547) -12.5 W12-W21 N34-N37
      4 1508 1451 (N19W43 L=340)   639    139        1 (1440/1417-1551/1472) -50.2 W39-W47 N15-N21
      5 1225 1357 (N21W21 L=317)   247     91        1 (1210/1305-1236/1394)  51.3 W18-W22 N18-N23
      6  568 1121 (N23E27 L=270)   937    339        2 ( 412/1081- 711/1151) -23.2 E38-E17 N17-N27
      7 1746 1090 (S08W48 L=345)   665    116        1 (1720/1039-1773/1150)  39.8 W45-W52 S10-S05
      8  386 1025 (N20E42 L=254)   744    141        2 ( 337/ 975- 422/1099)  81.9 E46-E38 N17-N24
      9 1007 1075 (N10W00 L=297)   258     91        0 ( 961/1063-1051/1089) -14.0 E02-W02 N09-N11
     10 1363 1044 (N00W20 L=316)   197     50        0 (1347/1015-1377/1067)  45.0 W18-W21 S01-N01
     11  173 1003 (N22E66 L=230)   940    113        1 ( 137/ 960- 208/1045) -28.6 E72-E61 N18-N24
     12 1850  951 (S20W57 L=353)  1123    194        0 (1782/ 915-1895/ 994) -19.7 W49-W63 S22-S16
     13 1303  931 (S05W14 L=310)   313     82        1 (1270/ 899-1333/ 961) -68.2 W12-W15 S07-S02
     14  434  894 (N12E40 L=256)   420    156        0 ( 361/ 867- 497/ 909) -21.8 E46-E36 N09-N13
     15 1364  650 (S23W12 L=308)  1217    367        0 (1270/ 529-1471/ 798)  32.3 W02-W21 S29-S16
     16  767  485 (S19E30 L=267)   274     72        0 ( 749/ 459- 789/ 524) -90.0 E30-E28 S21-S16
     17 1152  292 (S44E13 L=283)  2990    394        7 (1004/ 173-1252/ 432)  28.7 E32-E01 S51-S34
     18  959  164 (S50E40 L=257)   597     78        0 ( 928/ 150- 996/ 184) -45.0 E43-E35 S52-S47
  

    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