Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  879 1876 (N61E10 L=258)  6507    955       18 ( 704/1752-1088/1973)  26.6 E36-W36 N48-N78
      2  891 1915 (N67E10 L=258)   733    129        1 ( 866/1886- 914/1937)  41.2 E16-E06 N62-N71
      3  660 1816 (N56E35 L=234)  2231    344        7 ( 590/1769- 729/1874) -14.0 E49-E23 N50-N64
      4  582 1758 (N51E40 L=229)  1681    300        5 ( 528/1692- 630/1809) -34.4 E55-E34 N44-N57
      5  450 1724 (N49E56 L=212)  1065    182        2 ( 415/1700- 498/1747)  10.3 E65-E49 N47-N51
      6  503 1708 (N47E46 L=223)   727    132        2 ( 461/1687- 546/1729)   6.3 E51-E40 N45-N49
      7 1032  341 (S48E04 L=264)   289     26        1 (1007/ 326-1056/ 360) -11.3 E06-E01 S48-S45
      8 1324  213 (S62W32 L=300)  1109    219        1 (1291/ 163-1361/ 260) -50.2 W25-W41 S68-S56
      9  921  219 (S58E19 L=249)   381    101        3 ( 895/ 198- 953/ 249)  24.0 E22-E13 S60-S55
  

    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