Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1218 1693 (N40W27 L=115)   333     68        0 (1198/1660-1241/1719)  45.0 W24-W29 N37-N41
      2 1103 1694 (N42W18 L=107)   609    137        0 (1050/1657-1172/1714)  10.3 W13-W24 N39-N42
      3  736 1651 (N45E12 L= 76)   408     52        3 ( 702/1629- 765/1672)   0.0 E15-E09 N42-N46
      4 1466 1559 (N26W39 L=128)   298     65        0 (1440/1535-1489/1584)  21.8 W36-W41 N25-N27
      5 1585 1403 (N15W43 L=131)   531     84        1 (1538/1378-1628/1425) -26.6 W39-W45 N13-N16
      6  159 1065 (N14E63 L= 26)   927    243        0 ( 110/1037- 212/1101) -21.0 E69-E56 N11-N16
      7 1760  741 (S27W49 L=137)   531     81        1 (1735/ 704-1775/ 789)  90.0 W46-W50 S29-S24
  

    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