Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  875 1767 (N60W12 L=272)   371    103        0 ( 825/1738- 905/1799)   8.7 W04-W17 N58-N62
      2 1177 1637 (N41W29 L=288)   359     89        0 (1158/1591-1209/1672)  26.6 W24-W33 N37-N43
      3 1317 1646 (N38W40 L=300)   743    207        4 (1229/1604-1406/1672) -38.7 W33-W44 N32-N41
      4  502 1637 (N59E38 L=222)   535     74        1 ( 463/1617- 534/1661)  12.5 E44-E33 N56-N60
      5 1593  750 (S22W31 L=291)   273     80        1 (1563/ 725-1616/ 789)  21.8 W28-W33 S23-S20
      6 1709  646 (S33W42 L=302)   369     67        1 (1688/ 627-1729/ 670)  33.7 W40-W44 S34-S30
      7  332  546 (S11E62 L=198)   699    104        0 ( 283/ 519- 370/ 581)  11.3 E70-E55 S12-S09
      8  441  416 (S21E59 L=201)  1029    105        4 ( 388/ 387- 490/ 451) -16.7 E65-E54 S23-S18
  

    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