Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1555 1698 (N43W44 L=357)   533     49        0 (1523/1671-1583/1721) -36.9 W41-W46 N41-N45
      2 1674 1544 (N32W48 L=  1)  1068    169        1 (1602/1482-1732/1609) -41.2 W42-W52 N27-N36
      3  898 1554 (N29E10 L=303)   257     71        0 ( 858/1536- 943/1567)  -9.5 E12-E06 N28-N30
      4  998 1382 (N18E02 L=311)   339     78        0 ( 975/1338-1020/1423) -68.2 E03-W00 N15-N20
      5  393 1289 (N12E43 L=271)   499    147        1 ( 333/1232- 457/1334)  36.9 E47-E37 N08-N14
      6  300 1247 (N10E50 L=263)   970    167        2 ( 274/1169- 336/1329)  84.3 E52-E47 N04-N14
      7  230 1284 (N12E58 L=255)   396    122        1 ( 199/1254- 272/1326)  33.7 E60-E53 N10-N14
      8  842 1266 (N10E12 L=302)   599    181        0 ( 753/1232- 931/1316) -20.0 E17-E06 N08-N13
      9 1858 1213 (N11W60 L= 14)   430     89        0 (1824/1194-1886/1231)  -8.1 W56-W63 N10-N11
     10 1526 1069 (N00W31 L=344)   347     65        1 (1493/1039-1555/1095)  33.7 W29-W33 S01-N01
     11 1760  984 (S04W50 L=  3)   293     76        0 (1744/ 954-1777/1020)  63.4 W48-W51 S05-S01
     12 1499  855 (S13W31 L=344)  1109    410        3 (1389/ 679-1614/ 996)  59.0 W24-W37 S23-S03
     13  236  846 (S14E56 L=257)  1786    441        9 ( 159/ 771- 338/ 930) -28.3 E64-E46 S18-S09
     14  142  819 (S15E68 L=245)   977    286        4 ( 105/ 789- 193/ 840)   0.0 E74-E61 S17-S13
     15  775  662 (S26E16 L=297)  1153    270        2 ( 657/ 603- 837/ 762)  34.7 E24-E11 S30-S19
     16 1711  655 (S24W51 L=  4)   362     60        0 (1682/ 618-1734/ 687)  63.4 W48-W52 S25-S21
     17 1144  628 (S28W09 L=322)   481    151        1 (1108/ 588-1181/ 672)  71.6 W06-W11 S30-S24
     18 1456  522 (S34W33 L=346)   468    101        0 (1404/ 497-1508/ 548)  29.7 W29-W36 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