Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1369 1839 (N41W68 L=355)   643     71        1 (1342/1813-1403/1870)  -3.0 W61-W80 N40-N42
      2 1569 1582 (N20W56 L=342)  1723    410        6 (1406/1469-1693/1717) -77.2 W52-W61 N10-N32
      3  972 1554 (N38W14 L=301)  1319    316        5 ( 861/1512-1085/1588) -21.8 W06-W21 N33-N42
      4  340 1519 (N56E50 L=237)  1226    138        2 ( 294/1478- 407/1543)   6.7 E61-E38 N53-N58
      5 1280 1481 (N24W30 L=316)   476    137        1 (1221/1449-1343/1513) -54.5 W27-W32 N20-N27
      6 1716 1413 (N07W57 L=344)   999    147        2 (1684/1352-1747/1485)  28.6 W51-W62 N04-N10
      7 1404 1359 (N14W32 L=319)   358     53        1 (1375/1339-1436/1387) -71.6 W29-W33 N12-N15
      8 1170 1232 (N14W14 L=301)   384    102        3 (1125/1200-1220/1264) -53.1 W11-W16 N11-N16
      9  206 1177 (N35E60 L=227)   834    123        1 ( 161/1135- 251/1218) -21.8 E68-E53 N31-N37
     10  780 1136 (N19E11 L=276)   372     97        1 ( 734/1097- 816/1172) -63.4 E13-E09 N16-N22
     11 1717 1075 (S11W44 L=331)   339     46        0 (1699/1044-1735/1101)  56.3 W42-W45 S12-S09
     12 1863  975 (S22W57 L=344)   442     96        0 (1827/ 958-1895/1000) -11.3 W52-W62 S23-S19
     13  767  927 (N08E17 L=270)   396    152        1 ( 713/ 903- 829/ 977)   7.1 E20-E12 N06-N08
     14 1614  718 (S29W29 L=316)  1038    314        2 (1553/ 594-1705/ 841)  34.5 W22-W38 S35-S23
     15 1388  589 (S29W10 L=297)  1012    222        1 (1311/ 524-1450/ 682)  26.6 W03-W15 S32-S24
     16 1718  585 (S42W39 L=326)  1349    118        1 (1683/ 534-1749/ 643)  66.8 W34-W42 S46-S37
     17  690  532 (S13E34 L=253)   306     81        2 ( 659/ 513- 729/ 555) -63.4 E36-E31 S15-S11
     18 1079  236 (S45E26 L=261)  1609    248        3 ( 963/ 180-1173/ 290)   4.8 E39-E15 S47-S42
  

    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