Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  994 1682 (N49W18 L=350)   348     66        0 ( 964/1660-1021/1702) -53.1 W15-W19 N46-N51
      2 1194 1624 (N39W29 L=  2)  1316    189        2 (1125/1585-1274/1678) -63.4 W24-W33 N34-N44
      3 1668 1543 (N21W62 L= 35)   490     81        0 (1636/1520-1687/1575)  18.4 W59-W66 N19-N22
      4 1513 1492 (N22W45 L= 17)   632     98        2 (1483/1440-1539/1539)  35.5 W42-W49 N19-N24
      5 1686 1383 (N12W54 L= 26)   659    107        0 (1646/1336-1720/1426)  10.3 W48-W59 N10-N13
      6  478 1327 (N36E33 L=299)  1091    230        2 ( 361/1287- 567/1359) -14.0 E45-E25 N32-N38
      7  745 1246 (N25E13 L=319)   882    240        3 ( 643/1169- 825/1296) -39.3 E19-E08 N19-N30
      8  704 1070 (N16E18 L=314)   498    213        0 ( 649/1011- 762/1108) -45.0 E22-E15 N11-N18
      9 1751 1022 (S10W48 L= 20)   332     69        1 (1722/ 997-1772/1042) -36.9 W45-W49 S11-S08
     10  158  948 (N16E69 L=263)  2151    170        3 ( 125/ 863- 211/1025) -73.3 E74-E60 N10-N20
     11  470  937 (N12E37 L=295)   523    155        3 ( 421/ 887- 530/ 988) -53.1 E40-E33 N07-N15
     12 1730  903 (S17W44 L= 16)  1132    223        4 (1672/ 848-1805/ 955)  14.0 W39-W52 S20-S12
     13 1271  845 (S08W11 L=343)   433    197        0 (1202/ 790-1331/ 911)  21.8 W05-W15 S10-S05
     14  869  748 (S06E15 L=318)   495    171        0 ( 815/ 711- 928/ 780)  -8.1 E18-E11 S08-S03
     15  608  496 (S17E39 L=293)   258     35        0 ( 584/ 476- 634/ 516) -33.7 E40-E37 S18-S15
     16  659  340 (S29E44 L=288)   471     93        0 ( 616/ 327- 716/ 354) -15.9 E47-E39 S30-S27
  

    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