Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  767 1760 (N53E30 L=159)   503     84        1 ( 716/1740- 801/1778)  14.0 E35-E26 N51-N54
      2 1589 1680 (N48W59 L=248)   598     43        0 (1560/1659-1616/1701)   8.1 W51-W65 N46-N49
      3 1227 1516 (N35W14 L=203)   750    248        1 (1148/1448-1314/1568) -32.5 W07-W19 N30-N38
      4  305 1524 (N31E67 L=123)   671     80        0 ( 274/1506- 348/1545)  -3.6 E75-E59 N30-N32
      5  469 1474 (N29E44 L=145)  1111    213        3 ( 394/1442- 575/1509)  -6.0 E52-E33 N26-N31
      6 1274 1321 (N22W16 L=205)   328     92        1 (1258/1273-1292/1366)  71.6 W14-W17 N18-N24
      7  290 1168 (N09E53 L=137)   807    173        0 ( 265/1097- 321/1254) -51.3 E56-E48 N04-N14
      8 1649 1143 (N11W42 L=232)   484    101        0 (1602/1105-1685/1178) -33.7 W38-W44 N08-N13
      9  608 1070 (N04E26 L=163)   258     89        0 ( 583/1033- 653/1098)  45.0 E27-E23 N02-N06
     10 1705  989 (N02W46 L=236)  1090    153        1 (1674/ 917-1728/1069)  90.0 W43-W48 S02-N06
     11 1850  807 (S10W64 L=254)   790    106        3 (1826/ 756-1877/ 850)  63.4 W61-W67 S13-S07
     12 1590  752 (S13W39 L=228)   361     53        1 (1566/ 726-1611/ 778)  33.7 W37-W40 S14-S11
     13  270  696 (S21E57 L=132)   702    122        0 ( 220/ 668- 313/ 736)  19.7 E65-E51 S22-S17
     14 1762  585 (S25W62 L=251)  2840    183        3 (1682/ 519-1814/ 663)  48.4 W54-W70 S29-S19
     15  876  580 (S25E09 L=181)   884    199        1 ( 813/ 544- 963/ 636) -21.8 E13-E02 S27-S21
     16  675  588 (S26E23 L=166)   462    151        0 ( 611/ 559- 733/ 630)  26.6 E28-E18 S27-S22
     17 1283  218 (S55W33 L=222)  2579    282        8 (1164/ 184-1400/ 261)   6.1 W19-W48 S59-S50
  

    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