Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1222 1607 (N33W33 L=308)  2246    571        7 (1032/1527-1484/1670) -45.0 W21-W46 N20-N42
      2  723 1646 (N53E00 L=274)   426     53        0 ( 699/1623- 751/1669)   0.0 E03-W02 N51-N55
      3  406 1541 (N56E40 L=235)   979    117        2 ( 364/1493- 462/1568)  -7.1 E47-E31 N51-N58
      4 1794 1474 (N07W71 L=346)   924    148        1 (1761/1429-1821/1539)  24.4 W64-W77 N05-N10
      5 1569 1403 (N11W45 L=320)   518    127        1 (1519/1349-1620/1446) -56.3 W42-W46 N07-N14
      6  853 1332 (N28E02 L=273)   193     98        0 ( 808/1303- 893/1347) -45.0 E04-W00 N25-N29
      7 1348 1311 (N13W27 L=302)   299     89        1 (1311/1281-1381/1344) -78.7 W24-W28 N10-N15
      8  333 1195 (N33E44 L=231)   985    305        0 ( 229/1120- 459/1256) -28.4 E57-E33 N26-N39
      9 1815 1075 (S14W54 L=328)   371     35        0 (1801/1049-1830/1108) -90.0 W51-W55 S15-S12
     10  943 1012 (N08E05 L=270)   245    115        0 ( 900/ 989- 996/1051)   0.0 E07-W00 N07-N08
     11  203  928 (N20E61 L=213)   868    130        0 ( 140/ 899- 252/ 951) -10.6 E71-E55 N17-N21
     12  376  932 (N17E44 L=231)   350    114        0 ( 327/ 917- 432/ 951) -20.6 E47-E39 N14-N18
     13 1758  848 (S26W43 L=318)   543    158        1 (1702/ 773-1795/ 903)  29.1 W37-W47 S28-S22
     14 1540  666 (S29W22 L=297)  1160    221        2 (1454/ 599-1601/ 762)  22.6 W14-W27 S32-S24
     15 1799  626 (S43W50 L=325)  1589    146        3 (1738/ 585-1840/ 681)  15.3 W41-W57 S46-S38
     16 1677  613 (S38W34 L=309)   625     99        1 (1645/ 573-1702/ 662)  51.3 W30-W36 S40-S35
     17  788  430 (S22E32 L=242)   326     54        1 ( 771/ 404- 807/ 457)  76.0 E34-E30 S24-S20
     18 1221  302 (S45E11 L=263)  1416    257        2 (1103/ 239-1319/ 380)   5.4 E22-E01 S47-S41
  

    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