Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1407 1685 (N37W50 L=299)  1402    357        3 (1272/1643-1546/1716) -15.6 W36-W61 N31-N42
      2  563 1651 (N59E29 L=219)   384     49        2 ( 536/1629- 586/1672)  21.8 E32-E25 N56-N60
      3  343 1362 (N40E48 L=200)   602    112        0 ( 300/1336- 384/1381)   0.0 E53-E43 N37-N41
      4  650  906 (N07E25 L=223)   224     65        2 ( 631/ 876- 672/ 935)  56.3 E26-E23 N05-N08
      5  511  880 (N08E35 L=213)   223     70        0 ( 478/ 853- 542/ 900)  11.3 E37-E32 N06-N08
      6 1750  836 (S21W46 L=294)   421     93        1 (1727/ 795-1776/ 875)  68.2 W43-W48 S23-S18
      7  457  560 (S11E48 L=200)   555     83        1 ( 417/ 535- 499/ 590)   7.1 E52-E44 S12-S09
      8  590  444 (S21E43 L=205)   376    105        2 ( 536/ 427- 640/ 458) -15.9 E47-E38 S21-S19
  

    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