Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1114 1854 (N51E22 L= 82)   949     98        1 (1057/1835-1162/1871)  49.4 E26-E17 N47-N54
      2  688 1776 (N32E47 L= 57)   430     41        1 ( 664/1745- 708/1802) -29.7 E50-E43 N29-N34
      3 1459 1679 (N47W14 L=118)   553     84        0 (1429/1636-1486/1708)  90.0 W11-W16 N42-N49
      4  664 1645 (N22E40 L= 64)   505     91        0 ( 624/1612- 712/1667)  38.7 E41-E36 N19-N24
      5  975 1559 (N23E16 L= 88)   445    180        0 ( 928/1483-1013/1602) -36.9 E20-E11 N18-N25
      6  522 1460 (N08E42 L= 62)   353     70        1 ( 495/1430- 547/1498)  76.0 E42-E39 N05-N10
      7  394  645 (S44E40 L= 64)   344     84        0 ( 357/ 630- 431/ 668)   0.0 E44-E35 S44-S41
      8 1506  435 (S27W52 L=156)  1072    251        0 (1433/ 353-1614/ 491)  67.4 W49-W57 S33-S21
      9  930  246 (S59W24 L=129)   659    107        0 ( 898/ 199- 968/ 292) -17.4 W17-W33 S62-S55
     10 1128  155 (S58W57 L=161)   879     90        1 (1106/ 117-1163/ 199) -19.3 W47-W68 S61-S54
  

    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