Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  671 1574 (N50E09 L=191)   267     45        1 ( 646/1556- 698/1599) -11.3 E11-E06 N48-N51
      2 1539 1528 (N20W49 L=249)   300     71        0 (1502/1517-1578/1541) -31.0 W46-W52 N18-N21
      3 1530 1392 (N13W41 L=241)   445    109        2 (1500/1357-1581/1440)  12.5 W38-W47 N11-N14
      4 1578 1337 (N09W42 L=242)   415     99        1 (1547/1297-1613/1380)  26.6 W38-W45 N07-N11
      5  723 1278 (N29E12 L=188)   325    115        1 ( 676/1261- 783/1297) -33.7 E14-E08 N26-N31
      6 1019 1079 (N10W01 L=201)   315     87        1 ( 978/1049-1054/1105) -68.2 W00-W02 N07-N12
      7  743 1075 (N17E15 L=185)   337    146        1 ( 688/1041- 794/1103) -39.8 E18-E12 N13-N18
      8 1622  898 (S16W33 L=232)   294     58        1 (1594/ 871-1641/ 931)  63.4 W30-W34 S18-S14
      9 1232  708 (S16W04 L=204)   900    318        5 (1141/ 612-1335/ 795)  13.0 E02-W11 S21-S12
     10 1869  708 (S38W59 L=259)   677     58        2 (1841/ 685-1902/ 730)  -4.4 W54-W68 S40-S36
     11  324  602 (S02E59 L=141)   928    144        2 ( 253/ 561- 379/ 633) -26.6 E65-E54 S05-N00
     12 1067  574 (S20E09 L=191)   221     87        0 (1037/ 537-1084/ 606)  78.7 E10-E08 S22-S17
     13  373  462 (S12E62 L=138)  1799    159        3 ( 343/ 395- 404/ 544)  48.0 E71-E56 S16-S06
     14 1599  446 (S47W26 L=226)  1260    177        2 (1530/ 412-1673/ 497)   0.0 W19-W34 S49-S44
     15  738  323 (S28E41 L=159)  2506    309        6 ( 706/ 151- 781/ 491)  40.9 E60-E29 S43-S17
     16  463  423 (S16E56 L=144)   908    120        3 ( 434/ 374- 499/ 472)  74.1 E59-E52 S19-S12
     17  509  281 (S27E66 L=134)  1402    180        2 ( 426/ 232- 570/ 330) -59.0 E71-E61 S31-S21
  

    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