Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1148 1738 (N44W04 L= 63)   274     48        0 (1123/1719-1174/1753)  36.9 W02-W06 N42-N45
      2 1344 1668 (N40W20 L= 79)   464     66        0 (1318/1642-1377/1691)  38.7 W17-W22 N37-N42
      3  902 1619 (N32E13 L= 46)   879    159        0 ( 847/1546- 952/1673) -48.4 E17-E08 N27-N36
      4 1444 1625 (N38W27 L= 86)   612     85        1 (1412/1580-1473/1654)  71.6 W25-W30 N34-N40
      5 1662 1508 (N31W44 L=103)  1842    279        3 (1558/1438-1783/1563) -11.9 W35-W55 N27-N34
      6  781 1308 (N11E17 L= 42)   269     48        0 ( 744/1287- 806/1327) -14.0 E19-E15 N09-N12
      7  213 1260 (N07E60 L=359)   490     87        0 ( 186/1233- 238/1287) -31.0 E63-E57 N05-N08
      8  436 1196 (N03E39 L= 20)   253     51        0 ( 404/1174- 466/1220)  31.0 E41-E36 N01-N04
      9  834  690 (S26E10 L= 49)   291     78        0 ( 810/ 657- 858/ 736)  59.0 E11-E08 S28-S23
     10  398  573 (S35E45 L= 14)   845    164        0 ( 333/ 544- 473/ 609) -15.3 E50-E39 S37-S33
  

    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