Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  580 1705 (N58E37 L= 63)  1781    328        1 ( 409/1677- 715/1733)   0.0 E70-E18 N55-N60
      2  488 1365 (N33E35 L= 65)  1109    340        1 ( 400/1232- 555/1462) -46.7 E46-E29 N23-N40
      3 1780 1363 (N14W62 L=162)   458     71        0 (1755/1340-1811/1391) -56.3 W59-W65 N12-N15
      4  371 1270 (N28E45 L= 55)   358    112        0 ( 347/1229- 411/1321) -45.0 E48-E40 N24-N31
      5 1891  793 (S23W68 L=168)  1659    176        1 (1844/ 743-1923/ 886)  45.0 W60-W74 S26-S17
      6 1536  810 (S14W30 L=130)   465     89        2 (1509/ 763-1567/ 852)  59.0 W27-W32 S16-S11
      7 1437  550 (S30W22 L=122)  2312    447        2 (1242/ 461-1616/ 642)  13.6 W07-W36 S33-S25
      8  537  603 (S15E38 L= 62)   720    167        1 ( 463/ 569- 615/ 640)   8.7 E44-E31 S16-S12
      9 1091  357 (S39E05 L= 95)   453    166        0 (1020/ 321-1156/ 389)  15.3 E11-W00 S41-S37
     10  685  204 (S48E52 L= 48)   923     69        1 ( 654/ 184- 727/ 224)   7.6 E61-E45 S50-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