Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  313 1549 (N11E64 L=115)  1925    305        4 ( 274/1474- 399/1634)  56.3 E69-E57 N06-N17
      2  461 1424 (N05E45 L=135)   551    137        2 ( 408/1382- 520/1461)  26.6 E49-E41 N02-N08
      3  524 1391 (N04E39 L=141)   249     63        0 ( 491/1365- 564/1413)  36.9 E40-E36 N02-N06
      4 1064 1299 (N09E04 L=176)   538    159        0 (1042/1217-1098/1370) -58.0 E05-W00 N05-N13
      5  185 1127 (S16E61 L=118)  2990    484        7 ( 139/1031- 239/1248) -77.0 E68-E57 S21-S08
      6  136 1166 (S14E69 L=111)  1070    149        3 ( 116/1126- 167/1203)  76.0 E72-E64 S15-S11
      7  164 1014 (S23E62 L=118)   410    144        0 ( 136/ 992- 197/1039)  31.0 E64-E57 S23-S20
      8 1513  984 (N02W29 L=209)   232     69        1 (1482/ 970-1547/1001)   0.0 W27-W31 N01-N02
      9  466  889 (S26E33 L=147)  1482    334        3 ( 363/ 772- 570/1000)  50.5 E40-E26 S34-S17
     10 1452  820 (S09W30 L=209)   323    116        4 (1412/ 796-1496/ 857)  -8.1 W26-W33 S09-S07
     11 1441  674 (S17W34 L=213)  1272    341        2 (1280/ 580-1560/ 785)  57.1 W26-W39 S25-S08
     12 1585  702 (S11W43 L=222)   569    145        1 (1547/ 649-1615/ 743) -51.3 W40-W46 S14-S08
     13 1563  608 (S17W45 L=224)   300     89        3 (1541/ 585-1592/ 629) -11.3 W42-W47 S18-S15
     14 1029  527 (S36W13 L=193)   506    144        0 ( 956/ 506-1111/ 550)  15.3 W08-W19 S37-S34
     15 1575  507 (S22W51 L=230)   429     64        0 (1539/ 487-1611/ 529)  36.9 W48-W52 S23-S20
     16 1314  486 (S31W34 L=213)   244    102        1 (1287/ 449-1343/ 521) -15.9 W30-W37 S32-S29
     17 1678  439 (S22W64 L=243)  2082    366        1 (1565/ 384-1790/ 493)  35.8 W56-W75 S28-S15
  

    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