Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  923 1765 (N58W08 L= 55)   857    198        1 ( 840/1729-1011/1791)  -3.0 E01-W18 N55-N60
      2  768 1733 (N59E11 L= 37)   307     43        0 ( 747/1712- 796/1747)   0.0 E13-E07 N56-N60
      3 1188 1485 (N32W19 L= 66)   382    200        0 (1091/1431-1244/1510) -37.9 W12-W21 N27-N34
      4  242 1332 (N34E62 L=345)  3625    414        9 ( 166/1141- 337/1486)  67.8 E67-E53 N21-N43
      5  996 1417 (N31W04 L= 52)   424    170        0 ( 947/1380-1076/1470) -41.2 W01-W09 N27-N34
      6 1698 1309 (N12W51 L= 98)   828    189        2 (1662/1259-1762/1386)  82.9 W48-W56 N08-N16
      7  367 1188 (N24E44 L=  3)   498    111        3 ( 330/1144- 395/1237) -45.0 E48-E41 N20-N26
      8  614 1096 (N16E25 L= 22)   654    285        2 ( 533/1020- 732/1192) -45.0 E30-E17 N09-N22
      9 1845  914 (S15W58 L=106)  1349    197        5 (1817/ 820-1880/ 995) -71.6 W55-W63 S21-S09
     10  237  908 (N07E58 L=349)   417     91        0 ( 211/ 868- 256/ 949) -68.2 E60-E55 N04-N09
     11 1409  796 (S13W21 L= 68)   430    125        2 (1368/ 760-1463/ 838)  23.2 W17-W24 S14-S11
     12  547  634 (S12E37 L= 11)   405    129        0 ( 523/ 581- 564/ 692)  82.9 E39-E35 S15-S07
     13  270  600 (S12E64 L=343)   620     81        1 ( 255/ 561- 289/ 641)  51.3 E67-E61 S14-S09
     14  355  605 (S12E54 L=353)   454     77        1 ( 325/ 579- 394/ 629) -36.9 E57-E51 S14-S10
     15  493  535 (S18E44 L=  3)   452    112        0 ( 445/ 504- 529/ 578)  24.0 E49-E40 S19-S15
     16 1400  394 (S42W18 L= 66)   420     62        1 (1371/ 370-1440/ 412)   0.0 W15-W22 S43-S39
     17 1227  304 (S46W03 L= 50)   870    112        2 (1178/ 272-1285/ 332)   6.3 E01-W08 S48-S42
     18 1034  215 (S51E17 L= 30)  1123    142        2 ( 990/ 168-1074/ 275)  45.0 E22-E11 S56-S45
  

    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