Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1023 1938 (N64W54 L=  4)   861     76        0 (1010/1896-1039/1978)  10.7 W42-W79 N59-N66
      2  958 1897 (N63W38 L=348)   613     55        1 ( 930/1869- 980/1924)   0.0 W32-W46 N59-N65
      3 1309 1810 (N42W52 L=  2)  2025    351        2 (1160/1765-1465/1853) -15.9 W36-W64 N35-N46
      4  966 1663 (N43W16 L=325)   473     68        0 ( 941/1623- 996/1699)  33.7 W12-W19 N40-N45
      5  776 1596 (N44E01 L=309)   418     81        0 ( 728/1569- 821/1619)   7.1 E05-W03 N42-N45
      6 1742 1203 (S04W49 L=359)  2975    524        4 (1677/ 956-1785/1420)  58.0 W39-W56 S16-N08
      7  679 1255 (N25E15 L=294)   295     79        0 ( 640/1225- 710/1279) -36.9 E17-E13 N22-N26
      8  364 1183 (N27E40 L=270)   350    108        1 ( 350/1139- 383/1238)  63.4 E41-E37 N24-N30
      9  991 1175 (N13W02 L=311)   753    220        1 ( 902/1139-1083/1238) -59.0 E02-W06 N08-N18
     10  529 1180 (N24E27 L=283)   250     68        1 ( 500/1162- 565/1205) -36.9 E28-E24 N22-N25
     11  280 1106 (N24E48 L=261)  1626    343        2 ( 149/ 982- 382/1193) -33.7 E65-E40 N15-N32
     12  486 1034 (N16E32 L=278)   330    150        0 ( 418/1002- 551/1054)   0.0 E37-E27 N15-N17
     13  237  839 (N10E57 L=253)   780     94        2 ( 202/ 805- 276/ 885)  12.5 E61-E52 N07-N11
     14  392  798 (N05E44 L=266)   541    127        0 ( 331/ 775- 440/ 819) -26.6 E48-E40 N02-N06
     15  891  737 (S09E14 L=296)   346     80        0 ( 870/ 691- 915/ 779)  80.5 E14-E12 S12-S06
     16 1432  660 (S27W16 L=326)   332     33        0 (1410/ 632-1451/ 687)  45.0 W14-W17 S28-S25
     17  359  430 (S16E66 L=244)  5637    423        7 ( 259/ 247- 509/ 570)  87.5 E79-E60 S31-S06
     18 1494  410 (S47W18 L=327)   422     69        0 (1452/ 395-1524/ 430) -15.9 W13-W20 S49-S45
     19  901  331 (S35E28 L=281)   886    181        0 ( 832/ 271- 959/ 387)  14.9 E36-E21 S37-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