Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1254 1825 (N57E15 L=353)  1484    177        2 (1174/1784-1333/1853)   3.0 E23-E04 N53-N60
      2  624 1727 (N27E52 L=316)  2496    272        7 ( 586/1607- 671/1842) -27.3 E75-E44 N19-N35
      3 1672 1622 (N57W41 L= 49)   866     74        1 (1650/1577-1690/1660)  40.2 W36-W53 N52-N63
      4  654 1569 (N18E39 L=329)   359     77        1 ( 622/1548- 696/1594)   0.0 E42-E35 N16-N19
      5 1893  906 (N14W67 L= 75)   720    118        1 (1873/ 879-1917/ 950) -20.6 W63-W72 N12-N16
      6  528  837 (S28E26 L=341)   481     79        1 ( 513/ 791- 551/ 889) -71.6 E28-E24 S30-S24
      7  906  768 (S22W00 L=  8)   259     90        0 ( 889/ 731- 925/ 801) -76.0 E01-W02 S23-S19
      8  682  727 (S30E13 L=355)   238     99        0 ( 634/ 712- 728/ 746)  39.8 E16-E10 S32-S27
      9  739  635 (S34E07 L=  1)   225     58        0 ( 708/ 620- 766/ 653)  11.3 E09-E04 S35-S33
     10 1370  586 (S19W34 L= 42)   397    125        0 (1307/ 571-1430/ 604)  35.5 W30-W37 S21-S16
     11  634  326 (S62E08 L=  0)   570     57        0 ( 620/ 292- 649/ 367) -82.9 E10-E05 S66-S58
     12 1417  277 (S35W61 L= 69)  1080    109        2 (1364/ 256-1466/ 312)  39.8 W55-W66 S38-S32
  

    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