File: process.py

package info (click to toggle)
arc-gui-clients 0.4.6-7
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 25,484 kB
  • sloc: cpp: 11,077; python: 181; perl: 171; pascal: 165; ansic: 12; makefile: 7
file content (17 lines) | stat: -rwxr-xr-x 331 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
#!/bin/env python

import numpy as np
import scipy
import pylab
import pymorph
import mahotas
from scipy import ndimage

dna = mahotas.imread('dna.jpeg')

dnaf = ndimage.gaussian_filter(dna, 16)
rmax = pymorph.regmax(dnaf)
pylab.imshow(pymorph.overlay(dna, rmax))
seeds,nr_nuclei = ndimage.label(rmax)
print nr_nuclei
pylab.show()