예제 #1
0
def main():
    import mrc
    im1 = mrc.read(
        '/ami/data00/leginon/07sep26cal/rawdata/07sep26cal_00006m.mrc')
    mag1 = 14500
    im2 = mrc.read(
        '/ami/data00/leginon/07sep26cal/rawdata/07sep26cal_00007m.mrc')
    mag2 = 11500
    '''
        im1 = mrc.read('/ami/data00/leginon/07sep26cal/rawdata/07sep26cal_00002m.mrc')
        mag1 = 50000
        im2 = mrc.read('/ami/data00/leginon/07sep26cal/rawdata/07sep26cal_00003m.mrc')
        mag2 = 29000
        '''
    '''
        im1 = mrc.read('/ami/data00/leginon/07sep25jim/rawdata/07sep25jim_00005hl3.mrc')
        mag1 = 7800
        im2 = mrc.read('/ami/data00/leginon/07sep25jim/rawdata/07sep25jim_00005hl2.mrc')
        mag2 = 6500
        '''
    '''
        im1 = mrc.read('/ami/data00/leginon/07sep25jim/rawdata/07sep25jim_00004fa.mrc')
        mag1 = 50000
        im2 = mrc.read('/ami/data00/leginon/07sep25jim/rawdata/07sep25jim_00004fa2.mrc')
        mag2 = 29000
        '''
    '''
        im1 = mrc.read('../images/07apr03a/07apr03a_00003gr_00005sq_v01_00006hl.mrc')
        mag1 = 5000
        im2 = mrc.read('../images/07apr03a/07apr03a_00003gr_00005sq_v01.mrc')
        mag2 = 800
        '''
    #im1 = mrc.read('../images/07sep17jim4/07sep17jim4_00001a.mrc')
    #im2 = mrc.read('../images/07sep17jim4/07sep17jim4_00003a.mrc')
    #im1 = mrc.read('im00.mrc')
    #im2 = mrc.read('im01.mrc')

    scale = float(mag2) / mag1
    scalestart = scale - 0.02
    scaleend = scale + 0.02
    scaleinc = 0.005
    anglestart = -3
    angleend = 3
    angleinc = 0.25
    bin = 4
    binim1 = imagefun.bin(im1, bin)
    binim2 = imagefun.bin(im2, bin)
    rotation, scale = findRotationScale(binim1, binim2, anglestart, angleend,
                                        angleinc, scalestart, scaleend,
                                        scaleinc)

    shift = findShift(im1, im2, scale, rotation)
    print 'BEST SHIFT', shift
예제 #2
0
파일: msc.py 프로젝트: nramm/maskiton
def main():
    import mrc

    im1 = mrc.read("/ami/data00/leginon/07sep26cal/rawdata/07sep26cal_00006m.mrc")
    mag1 = 14500
    im2 = mrc.read("/ami/data00/leginon/07sep26cal/rawdata/07sep26cal_00007m.mrc")
    mag2 = 11500
    """
	im1 = mrc.read('/ami/data00/leginon/07sep26cal/rawdata/07sep26cal_00002m.mrc')
	mag1 = 50000
	im2 = mrc.read('/ami/data00/leginon/07sep26cal/rawdata/07sep26cal_00003m.mrc')
	mag2 = 29000
	"""
    """
	im1 = mrc.read('/ami/data00/leginon/07sep25jim/rawdata/07sep25jim_00005hl3.mrc')
	mag1 = 7800
	im2 = mrc.read('/ami/data00/leginon/07sep25jim/rawdata/07sep25jim_00005hl2.mrc')
	mag2 = 6500
	"""
    """
	im1 = mrc.read('/ami/data00/leginon/07sep25jim/rawdata/07sep25jim_00004fa.mrc')
	mag1 = 50000
	im2 = mrc.read('/ami/data00/leginon/07sep25jim/rawdata/07sep25jim_00004fa2.mrc')
	mag2 = 29000
	"""
    """
	im1 = mrc.read('../images/07apr03a/07apr03a_00003gr_00005sq_v01_00006hl.mrc')
	mag1 = 5000
	im2 = mrc.read('../images/07apr03a/07apr03a_00003gr_00005sq_v01.mrc')
	mag2 = 800
	"""
    # im1 = mrc.read('../images/07sep17jim4/07sep17jim4_00001a.mrc')
    # im2 = mrc.read('../images/07sep17jim4/07sep17jim4_00003a.mrc')
    # im1 = mrc.read('im00.mrc')
    # im2 = mrc.read('im01.mrc')

    scale = float(mag2) / mag1
    scalestart = scale - 0.02
    scaleend = scale + 0.02
    scaleinc = 0.005
    anglestart = -3
    angleend = 3
    angleinc = 0.25
    bin = 4
    binim1 = imagefun.bin(im1, bin)
    binim2 = imagefun.bin(im2, bin)
    rotation, scale = findRotationScale(binim1, binim2, anglestart, angleend, angleinc, scalestart, scaleend, scaleinc)

    shift = findShift(im1, im2, scale, rotation)
    print "BEST SHIFT", shift
예제 #3
0
def readMrc():
        if len(sys.argv) <= 1:
                return numpy.zeros((1,1,1))
        mrcfile = sys.argv[1]
        if not os.path.isfile(mrcfile):
                return numpy.zeros((1,1,1))
        import mrc
        voldata = mrc.read(mrcfile)
        array = numpy.where(voldata > 0.5, 1.0, 0.0)
        return array