Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  902 1858 (N61W06 L=158)   614     84        1 ( 865/1831- 938/1885)  15.3 W00-W11 N58-N64
      2 1011 1540 (N32W06 L=158)   518    160        0 ( 951/1486-1071/1591) -53.1 W02-W09 N27-N35
      3 1862 1309 (N07W67 L=219)   828     81        3 (1837/1279-1896/1354) -33.7 W63-W72 N04-N09
      4 1771 1284 (N06W55 L=207)   475     77        0 (1752/1247-1787/1325)  38.7 W52-W58 N04-N08
      5  379 1246 (N21E41 L=111)  1429    226        4 ( 290/1192- 479/1290)   3.4 E49-E32 N17-N23
      6 1579  886 (S15W34 L=185)  1532    376        1 (1533/ 712-1612/1046)  81.5 W31-W36 S25-S05
      7 1872  912 (S17W62 L=214)  4312    350       11 (1800/ 799-1930/1045) -71.6 W53-W72 S24-S09
      8 1468  771 (S21W25 L=177)   246     93        0 (1441/ 732-1490/ 816)  45.0 W23-W27 S23-S18
      9 1711  602 (S35W49 L=200)   800    114        0 (1665/ 568-1765/ 640)  16.7 W44-W54 S37-S33
     10  305  565 (S19E62 L= 90)   490     58        2 ( 283/ 541- 330/ 589) -36.9 E63-E58 S20-S17
     11 1511  485 (S41W31 L=182)   410     97        0 (1457/ 470-1555/ 503)  -5.7 W25-W35 S42-S39
  

    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