Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1264 1694 (N36W38 L=140)   619    168        0 (1195/1669-1356/1718) -37.9 W33-W43 N31-N39
      2 1701 1635 (N19W71 L=173)  1219    111        3 (1680/1578-1717/1689)  20.6 W65-W84 N15-N21
      3 1341 1584 (N27W36 L=138)   626    135        1 (1279/1531-1391/1617)  10.3 W30-W41 N25-N28
      4 1109 1414 (N23W15 L=117)   558    245        0 (1015/1324-1183/1483) -69.0 W11-W17 N16-N29
      5  912 1280 (N21E00 L=101)  1102    369        3 ( 829/1195-1068/1386) -63.4 E03-W07 N12-N28
      6 1762 1174 (S07W50 L=152)  2197    395        2 (1708/ 977-1801/1315)  62.1 W42-W55 S16-N01
      7  463 1269 (N31E30 L= 71)   904    197        2 ( 386/1227- 559/1302)  -8.1 E37-E23 N28-N32
      8  729 1200 (N21E13 L= 89)  1251    322        2 ( 589/1085- 825/1262) -49.4 E21-E08 N12-N27
      9  508 1090 (N19E30 L= 72)   261    115        0 ( 454/1079- 565/1104) -26.6 E33-E25 N17-N21
     10 1855  865 (S27W55 L=156)  2422    227        3 (1749/ 806-1945/ 927)  -2.2 W42-W70 S30-S24
     11 1571  855 (S19W29 L=130)  1547    329        2 (1450/ 739-1698/ 925)   8.1 W18-W39 S23-S16
     12  399  798 (N05E44 L= 58)   565     98        1 ( 366/ 749- 432/ 842)  45.0 E47-E40 N02-N07
     13  782  729 (S07E21 L= 81)   422    130        0 ( 739/ 684- 805/ 782)  51.3 E24-E18 S08-S03
     14  302  655 (S01E57 L= 44)   635     71        1 ( 274/ 628- 330/ 687) -33.7 E60-E54 S02-N00
     15  403  559 (S08E52 L= 49)  3936    300       17 ( 317/ 446- 474/ 685)  56.3 E63-E43 S15-N00
     16  365  640 (S03E52 L= 50)   414     82        0 ( 340/ 608- 384/ 675) -76.0 E53-E49 S04-N00
     17  479  461 (S16E51 L= 50)   564     90        2 ( 445/ 428- 515/ 499)  12.5 E55-E46 S18-S14
     18  597  335 (S26E50 L= 51)  1111    208        0 ( 499/ 314- 686/ 365) -19.7 E57-E42 S28-S22
  

    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