Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1110 1790 (N43E22 L= 18)  2294    378        3 ( 993/1661-1242/1928) -21.3 E47-E06 N36-N53
      2  502 1725 (N22E62 L=338)  4801    279        9 ( 414/1631- 620/1815) -82.9 E72-E57 N14-N31
      3  682 1728 (N26E46 L=354)   538    169        1 ( 617/1661- 711/1776) -19.7 E54-E40 N21-N28
      4  546 1671 (N19E53 L=347)  1402    249        2 ( 448/1589- 628/1738)  72.9 E57-E49 N11-N24
      5  596 1567 (N13E42 L=358)   765    205        1 ( 534/1490- 668/1657)  80.5 E44-E40 N08-N20
      6 1698 1463 (N41W35 L= 75)   418     46        1 (1673/1435-1721/1493) -45.0 W32-W38 N39-N43
      7  444 1425 (N02E47 L=353)   399     64        0 ( 418/1397- 473/1455)  56.3 E48-E45 N00-N03
      8  363 1080 (S20E44 L=356)   407     76        2 ( 331/1050- 389/1118)  53.1 E46-E41 S21-S17
      9  297 1055 (S23E50 L=350)   785    174        1 ( 257/ 964- 324/1111) -58.0 E54-E45 S27-S19
     10  449 1018 (S22E36 L=  4)   354     89        0 ( 412/ 993- 480/1038)  21.8 E38-E33 S23-S20
     11 1232  249 (S44W51 L= 91)   427     84        0 (1196/ 225-1271/ 283)   0.0 W44-W56 S45-S43
  

    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