Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1180 1822 (N52W24 L= 64)   439     88        1 (1158/1790-1210/1848)  35.5 W20-W27 N49-N54
      2 1658 1409 (N18W47 L= 86)  1027    211        2 (1577/1349-1767/1462) -25.0 W40-W55 N13-N21
      3  776 1372 (N21E13 L= 26)   356     84        1 ( 741/1348- 827/1397) -21.8 E15-E10 N19-N22
      4  631 1169 (N10E23 L= 17)   284     75        1 ( 614/1127- 649/1204)  63.4 E24-E21 N07-N11
      5  369 1065 (N05E42 L=357)   528    159        1 ( 344/1002- 402/1154) -77.5 E44-E40 N01-N10
      6  264 1037 (N05E52 L=348)  1621    359        3 ( 166/ 900- 305/1137) -51.8 E62-E48 S03-N11
      7  126  979 (N03E69 L=331)  5705    657       12 (  74/ 860- 175/1090) -69.0 E80-E62 S04-N09
      8 1759 1019 (S06W50 L= 89)   293     50        0 (1744/ 987-1777/1050)  63.4 W48-W51 S08-S04
      9 1848  789 (S21W59 L= 99)   525    143        3 (1821/ 752-1892/ 817)  23.2 W56-W64 S22-S19
     10  272  589 (S22E60 L=339)   970    160        2 ( 220/ 549- 303/ 641) -50.2 E64-E57 S24-S18
     11  805  365 (S42E24 L= 16)   694    164        0 ( 732/ 343- 879/ 403)   8.1 E30-E16 S44-S40
  

    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