Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  446 1665 (N36E54 L=357)  2062    160        3 ( 389/1593- 495/1719) -19.7 E72-E44 N29-N40
      2  906 1681 (N38E13 L= 38)   403     98        0 ( 866/1661- 954/1711) -14.0 E17-E09 N36-N40
      3 1640 1427 (N25W42 L= 94)   609     96        0 (1591/1391-1683/1458) -23.2 W37-W45 N22-N26
      4  614 1387 (N16E28 L= 23)   639    124        1 ( 581/1325- 657/1441)  74.1 E30-E26 N12-N19
      5 1344 1399 (N21W19 L= 70)   321    100        0 (1294/1365-1384/1434) -26.6 W15-W21 N18-N22
      6  782 1357 (N15E17 L= 35)   440    129        1 ( 721/1317- 833/1381) -15.9 E20-E13 N12-N16
      7  185 1145 (N01E62 L=349)  1741    254        4 ( 155/1044- 220/1265) -85.6 E68-E57 S05-N08
      8  586 1179 (N03E28 L= 24)   413     99        0 ( 563/1137- 619/1228) -71.6 E29-E25 N00-N06
      9  305  846 (S17E50 L=  1)   582    137        1 ( 262/ 797- 342/ 891)  40.6 E55-E46 S20-S14
     10  553  866 (S16E30 L= 22)   558    216        0 ( 450/ 842- 642/ 887)   8.1 E37-E23 S17-S14
     11 1377  812 (S15W23 L= 75)   218     59        0 (1361/ 783-1393/ 844) -71.6 W22-W24 S16-S12
     12 1487  745 (S18W32 L= 84)   301     61        0 (1457/ 726-1516/ 765)  45.0 W30-W34 S19-S16
     13 1542  705 (S20W37 L= 89)   681    185        1 (1481/ 646-1609/ 757)  48.8 W32-W41 S24-S16
     14 1275  706 (S22W18 L= 70)   341    128        0 (1224/ 683-1343/ 739) -12.5 W14-W23 S23-S20
     15 1581  481 (S34W49 L=101)   835     94        2 (1547/ 452-1617/ 520)  90.0 W44-W53 S36-S31
     16 1216  424 (S42W20 L= 71)   318     64        1 (1185/ 403-1255/ 449)  11.3 W17-W22 S43-S39
  

    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