Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1275 1762 (N55W30 L=167)   814    139        1 (1213/1730-1335/1790)   9.5 W22-W40 N52-N57
      2  393 1467 (N32E52 L= 85)  1764    318        2 ( 312/1357- 506/1574) -23.6 E69-E37 N24-N38
      3 1148 1552 (N38W11 L=148)   464    138        0 (1101/1519-1218/1574) -26.6 W07-W15 N35-N39
      4 1333 1504 (N34W24 L=161)   321     59        0 (1297/1488-1367/1518)   9.5 W21-W27 N33-N35
      5  852 1247 (N18E11 L=126)   580    186        0 ( 808/1169- 919/1302)  53.1 E13-E06 N13-N21
      6  148 1225 (N15E76 L= 61)  1849     88        2 ( 124/1180- 184/1261)  40.6 E83-E69 N11-N17
      7  556 1026 (N04E30 L=107)   967    246        0 ( 482/ 942- 657/1100)  33.7 E35-E23 S01-N08
      8  488  866 (S06E36 L=101)   299     50        1 ( 461/ 849- 518/ 883)   0.0 E37-E33 S06-S04
      9  431  622 (S22E44 L= 93)  2184    385        3 ( 333/ 445- 499/ 739)  36.3 E69-E36 S36-S14
     10 1178  481 (S32W10 L=147)  4039    710        8 ( 913/ 322-1422/ 701)  36.9 E10-W26 S44-S17
     11  822  503 (S30E15 L=122)  1583    494        1 ( 674/ 375-1018/ 661)  33.7 E30-W00 S40-S18
  

    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