Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1168 1733 (N42W36 L=132)   549    156        1 (1105/1701-1234/1763)  -4.8 W29-W41 N40-N42
      2 1181 1495 (N27W24 L=120)   491    206        0 (1091/1458-1267/1543)  -4.1 W16-W30 N25-N27
      3 1572 1474 (N14W47 L=144)  1729    433        6 (1412/1394-1741/1539) -36.9 W38-W59 N04-N22
      4 1397 1384 (N14W32 L=128)   285     64        1 (1366/1366-1432/1409) -45.0 W30-W33 N12-N16
      5 1875 1132 (S14W61 L=157)  1898    338        5 (1837/1001-1905/1254)  47.3 W53-W69 S21-S07
      6 1238 1063 (N02W13 L=109)   277     94        1 (1187/1045-1279/1089) -45.0 W10-W14 N00-N04
      7  355  969 (N18E43 L= 53)   654    152        3 ( 310/ 924- 423/1017) -36.9 E46-E38 N14-N20
      8 1684  956 (S17W38 L=135)   667    135        2 (1638/ 905-1727/1003) -58.0 W35-W40 S20-S12
      9 1716  827 (S26W38 L=135)   298     40        1 (1699/ 802-1734/ 847)  14.0 W36-W40 S27-S24
     10 1368  769 (S18W13 L=109)   273    106        1 (1310/ 743-1418/ 802)   0.0 W08-W16 S19-S17
     11 1191  668 (S19E00 L= 96)   336    136        1 (1168/ 631-1225/ 706)  45.0 E01-W02 S21-S16
     12 1587  652 (S33W26 L=122)   249     52        1 (1566/ 629-1609/ 679)  56.3 W23-W27 S33-S30
     13 1349  520 (S33W06 L=102)   996    175        1 (1310/ 449-1378/ 607)  52.1 W01-W09 S37-S27
     14  791  411 (S24E33 L= 64)   977    231        3 ( 751/ 322- 843/ 496)  65.2 E35-E27 S31-S18
     15 1529  413 (S49W18 L=114)  1930    303        4 (1457/ 352-1600/ 485)  38.7 W11-W25 S54-S42
     16  526  326 (S24E56 L= 40)  1682    273        2 ( 499/ 261- 566/ 406)  47.7 E63-E52 S28-S17
     17 1438  367 (S49W09 L=105)   749    113        2 (1397/ 328-1475/ 406)  82.9 W05-W12 S52-S44
  

    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