Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1061 1659 (N32E17 L=140)  1591    378        1 ( 990/1537-1153/1720)  86.4 E21-E09 N21-N37
      2  532 1683 (N20E54 L=102)   961    124        3 ( 479/1647- 579/1715)  56.3 E57-E50 N17-N23
      3 1397 1638 (N41W07 L=163)  1130    268        2 (1294/1535-1494/1705) -14.7 E03-W16 N36-N44
      4  441 1618 (N14E58 L= 99)   812    132        2 ( 425/1555- 462/1679) -41.6 E62-E53 N10-N18
      5  583 1627 (N18E46 L=110)   808    186        1 ( 496/1604- 656/1648)  24.4 E52-E40 N14-N20
      6  585 1539 (N12E41 L=115)   342    110        1 ( 566/1486- 615/1579) -66.8 E42-E38 N08-N15
      7 1070 1141 (N00E00 L=156)   321    112        1 (1048/1080-1099/1186) -53.1 E02-W02 S02-N02
      8 1803  979 (N13W52 L=208)   591    121        1 (1770/ 938-1835/1027) -23.2 W47-W55 N11-N14
      9 1835  866 (N08W59 L=215)   649     62        3 (1808/ 832-1864/ 898)  45.0 W56-W61 N05-N10
     10 1057  774 (S20W09 L=165)   335    147        1 (1015/ 713-1080/ 829) -35.5 W04-W11 S23-S18
     11  403  758 (S38E37 L=119)   797    167        1 ( 334/ 732- 458/ 812)  30.3 E44-E32 S40-S32
     12 1839  719 (N01W66 L=222)   506     48        0 (1815/ 700-1857/ 745) -63.4 W62-W69 N00-N02
     13 1734  648 (S07W58 L=214)   670     99        2 (1694/ 620-1764/ 676)   6.3 W53-W62 S08-S04
     14 1410  639 (S17W34 L=191)   351    123        0 (1370/ 582-1435/ 670) -33.7 W30-W37 S20-S15
     15 1463  543 (S21W42 L=199)   324     90        0 (1416/ 529-1502/ 557)  45.0 W39-W45 S23-S18
     16 1517  370 (S28W59 L=215)  1256    237        1 (1411/ 342-1625/ 411)  31.0 W49-W66 S31-S21
     17  757  268 (S64W10 L=166)  1462    147        3 ( 708/ 225- 806/ 327) -19.4 W01-W19 S68-S59
  

    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