Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1431 1774 (N55W54 L=163)  1257    100        2 (1401/1733-1474/1816)  11.7 W45-W74 N51-N57
      2  523 1679 (N49E51 L= 58)   478     69        0 ( 495/1655- 549/1706)  68.2 E53-E48 N46-N51
      3  725 1434 (N31E21 L= 89)  2292    561        1 ( 525/1286- 928/1563) -26.6 E41-E05 N20-N39
      4 1729 1484 (N30W63 L=172)   511     47        0 (1713/1456-1744/1514)  29.7 W60-W67 N27-N31
      5 1267 1265 (N18W17 L=126)   528    215        1 (1234/1184-1336/1331)  56.3 W14-W21 N13-N22
      6  332 1143 (N12E49 L= 60)   333     90        0 ( 311/1107- 363/1167)  45.0 E50-E46 N09-N13
      7  972 1008 (N03E03 L=106)   533    133        3 ( 920/ 966-1030/1052)  33.7 E06-W00 N00-N05
      8  836  670 (S18E13 L= 96)  1614    407        3 ( 715/ 497- 910/ 815)  56.3 E23-E08 S29-S08
      9 1472  721 (S16W29 L=138)   260     83        0 (1448/ 684-1496/ 758)  63.4 W27-W30 S18-S14
     10 1561  522 (S31W40 L=150)  4879    711        8 (1289/ 363-1767/ 736)  37.1 W19-W55 S42-S17
     11 1122  447 (S34W06 L=115)  1657    325        3 ( 984/ 353-1252/ 535)  34.4 E05-W14 S41-S28
  

    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