Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1495 1790 (N52W45 L=359)   578    148        1 (1482/1754-1507/1827)  53.1 W41-W50 N47-N56
      2  988 1765 (N46E05 L=309)   309     56        3 ( 964/1746-1010/1785) -45.0 E07-E03 N44-N48
      3 1645 1567 (N33W46 L=  1)  2058    301        4 (1519/1443-1751/1712) -62.2 W40-W53 N25-N44
      4  878 1553 (N29E11 L=303)   257     82        0 ( 835/1536- 926/1567)  -8.1 E14-E07 N28-N30
      5  976 1383 (N18E04 L=310)   461    103        0 ( 953/1337-1001/1431) -71.6 E05-E02 N15-N21
      6  385 1295 (N12E43 L=271)   468    163        2 ( 318/1242- 447/1337)  29.1 E48-E38 N09-N14
      7  285 1246 (N10E51 L=263)   954    203        6 ( 250/1169- 323/1328)  78.7 E54-E48 N04-N14
      8  822 1263 (N10E13 L=301)   650    162        0 ( 740/1229- 909/1309) -20.0 E18-E07 N08-N13
      9 1841 1218 (N11W58 L= 13)   447    117        2 (1802/1197-1874/1240) -14.0 W54-W62 N10-N12
     10 1581 1066 (N00W35 L=349)   980    366        2 (1468/1030-1740/1093)   2.9 W27-W47 S01-N02
     11 1449  808 (S16W28 L=342)  1173    276        3 (1366/ 681-1513/ 939)  66.4 W23-W30 S23-S07
     12 1714  890 (S09W47 L=  1)   324     73        1 (1693/ 858-1734/ 921)  63.4 W45-W48 S11-S07
     13  249  834 (S15E55 L=259)  1157    342        6 ( 169/ 796- 329/ 867)  -7.1 E63-E47 S17-S13
     14  238  732 (S21E58 L=257)   507    140        3 ( 211/ 701- 275/ 766) -31.0 E60-E54 S23-S18
     15  759  659 (S26E17 L=297)  1139    241        2 ( 652/ 604- 822/ 752)  36.9 E25-E12 S30-S20
     16 1127  619 (S28W08 L=322)   745    138        2 (1080/ 559-1164/ 669)  74.1 W04-W09 S32-S25
     17 1594  556 (S31W43 L=358)   446     58        1 (1572/ 521-1611/ 597) -68.2 W41-W45 S32-S27
     18 1433  517 (S34W32 L=346)   567    122        0 (1364/ 485-1490/ 548)  24.0 W26-W35 S36-S31
  

    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