Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  751 1843 (N36E53 L=168)   475     92        1 ( 715/1819- 775/1862)   0.0 E58-E47 N34-N38
      2  809 1775 (N33E41 L=180)   378     58        1 ( 789/1749- 829/1802) -53.1 E43-E38 N30-N34
      3  858 1654 (N25E30 L=191)   752    225        1 ( 823/1544- 897/1748) -42.9 E37-E23 N18-N31
      4  513 1637 (N16E53 L=168)   887    141        1 ( 478/1590- 549/1687) -24.4 E59-E48 N12-N19
      5  646 1613 (N17E42 L=180)   304     67        0 ( 618/1589- 672/1630) -14.0 E43-E38 N16-N18
      6  448 1588 (N11E56 L=165)   532    118        0 ( 418/1553- 503/1628) -14.0 E61-E49 N10-N13
      7  718 1538 (N14E33 L=188)   285     81        0 ( 687/1504- 741/1580)  90.0 E34-E32 N11-N17
      8  622 1502 (N10E38 L=183)   416    144        0 ( 604/1433- 651/1563) -77.5 E39-E35 N05-N14
      9 1426 1463 (N30W12 L=233)   347     76        0 (1396/1426-1448/1498) -36.9 W08-W13 N28-N32
     10  694 1470 (N10E32 L=189)   294     80        1 ( 663/1440- 735/1496) -18.4 E34-E28 N08-N10
     11 1769 1422 (N42W44 L=265)   462     68        1 (1730/1400-1800/1437)  11.3 W38-W48 N40-N43
     12 1751 1295 (N32W41 L=262)   841    216        0 (1727/1208-1781/1405) -65.6 W38-W45 N26-N38
     13  797 1266 (N01E19 L=202)   285     96        1 ( 760/1243- 835/1294)   0.0 E21-E16 N00-N01
     14  658  841 (S27E17 L=204)   564    151        1 ( 591/ 808- 722/ 875)   5.7 E22-E12 S27-S25
     15 1239  503 (S29W31 L=252)   382    159        0 (1170/ 485-1311/ 523)  18.4 W26-W35 S30-S27
     16 1651  386 (S20W72 L=293)  1049     90        0 (1623/ 349-1686/ 421) -71.6 W68-W75 S23-S17
  

    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