Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  590 1708 (N61E15 L= 67)   344     30        0 ( 569/1685- 607/1731)  26.6 E17-E12 N59-N63
      2  952 1320 (N23W03 L= 85)   899    177        2 ( 873/1289-1041/1347) -36.9 W00-W08 N20-N26
      3  872 1285 (N23E02 L= 79)   205     41        0 ( 842/1269- 898/1302) -33.7 E03-W00 N21-N24
      4 1817  917 (S23W50 L=132)  2318    536        4 (1710/ 695-1893/1143)  46.3 W37-W63 S33-S11
      5 1788 1103 (S11W51 L=132)   401     83        1 (1760/1079-1821/1131) -31.0 W48-W53 S13-S09
      6 1224 1098 (N04W13 L= 94)   242     92        0 (1181/1083-1266/1120) -21.8 W10-W15 N02-N05
      7 1767  989 (S17W46 L=128)   981    201        1 (1703/ 905-1839/1064)  14.9 W39-W54 S20-S14
      8  340  977 (N18E45 L= 37)  1949    418        3 ( 193/ 945- 561/1033) -15.7 E59-E27 N13-N23
      9  752  917 (N05E18 L= 64)   490    111        2 ( 721/ 876- 782/ 962)  45.0 E20-E15 N03-N07
     10  191  800 (N10E64 L= 17)  1388    201        2 ( 153/ 711- 232/ 870) -90.0 E68-E61 N04-N14
     11  588  845 (N05E29 L= 52)   350     84        0 ( 552/ 823- 637/ 870) -38.7 E31-E26 N02-N07
     12  543  657 (S05E38 L= 44)  1030    314        0 ( 419/ 587- 636/ 770)  13.4 E49-E28 S06-N00
     13 1578  499 (S43W24 L=106)   319     39        0 (1558/ 468-1598/ 522) -78.7 W22-W26 S45-S40
     14 1344  401 (S42W03 L= 84)   679    155        1 (1283/ 348-1412/ 447)  20.0 E02-W09 S44-S40
     15  472  351 (S22E61 L= 21)   432     49        2 ( 444/ 338- 501/ 369) -21.8 E63-E57 S22-S20
     16 1189  302 (S44E12 L= 69)   294     82        1 (1147/ 280-1221/ 327)   7.1 E16-E08 S45-S42
     17  931  209 (S43E37 L= 45)   942    158        2 ( 858/ 184-1013/ 233)  -8.1 E44-E28 S45-S40
  

    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