Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  851 1820 (N62W24 L=317)   843    112        0 ( 795/1795- 912/1838) -24.4 W17-W30 N58-N64
      2  698 1763 (N64W03 L=296)   473     59        0 ( 678/1725- 719/1795)  40.6 W00-W07 N60-N66
      3 1174 1731 (N42W38 L=331)   837    194        2 (1083/1676-1262/1770)  -6.0 W28-W47 N39-N43
      4  375 1598 (N61E47 L=247)  1202    179        1 ( 302/1549- 454/1649)   6.6 E61-E33 N57-N63
      5 1474 1603 (N24W48 L=342)   402     96        0 (1443/1568-1514/1627)   0.0 W44-W51 N22-N25
      6 1837 1266 (S05W61 L=355)   395     24        0 (1820/1246-1854/1286) -71.6 W59-W63 S06-S03
      7  515 1167 (N27E28 L=266)  1079    301        1 ( 380/1114- 637/1221) -30.1 E38-E19 N20-N32
      8 1796  935 (S21W48 L=341)  3471    438        4 (1694/ 780-1887/1123)  27.5 W36-W62 S27-S13
      9 1351  884 (S10W14 L=308)   568    126        0 (1334/ 824-1370/ 957)  76.0 W12-W16 S13-S05
     10  308  869 (N14E50 L=243)   501    153        2 ( 253/ 823- 359/ 903) -36.9 E54-E46 N10-N16
     11  393  846 (N11E43 L=250)   257      5        0 ( 376/ 828- 412/ 866)  26.6 E44-E41 N09-N12
     12  384  684 (N02E49 L=245)   291     59        0 ( 364/ 658- 410/ 702)  21.8 E51-E46 N00-N03
     13 1444  578 (S31W13 L=307)   843    180        1 (1370/ 535-1519/ 638)   4.8 W07-W19 S33-S29
     14  595  571 (S09E38 L=256)   305     60        1 ( 567/ 543- 620/ 597)   0.0 E40-E35 S10-S07
     15 1180  398 (S35E09 L=284)  1556    284        1 (1063/ 342-1313/ 476)   3.0 E18-W01 S37-S32
     16  600  392 (S20E47 L=247)   354     80        2 ( 565/ 361- 629/ 414) -63.4 E48-E44 S22-S18
     17  807  165 (S42E53 L=241)  2703    208        6 ( 730/ 115- 891/ 231)  -2.6 E73-E39 S46-S38
  

    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