Пример #1
0
    print 'FAIL: dials_regression not configured'
    raise

  # The XDS values
  xds_sigma_d = 0.060
  xds_sigma_m = 0.154

  # Read the reflections file
  filename = os.path.join(dials_regression,
      'centroid_test_data', 'reflections.h5')
  #filename = os.path.join('/home/upc86896/Data/X1_strong_M1S1_1_', 'reflections.h5')
  handle = NexusFile(filename, 'r')

  # Get the reflection list
  print 'Reading reflections.'
  reflections = handle.get_reflections()
  print 'Read {0} reflections.'.format(len(reflections))

  # Read images
  template = os.path.join(dials_regression,
      'centroid_test_data', 'centroid_*.cbf')
  #template = os.path.join('/home/upc86896/Data/X1_strong_M1S1_1_', 'X1_strong_M1S1_1_*.cbf')
  filenames = glob(template)

  # Load the sweep
  print 'Loading sweep'
  sweep = SweepFactory.sweep(filenames)[0:10]
  print 'Loaded sweep of {0} images.'.format(len(sweep))

  # Select the strong pixels to use in the divergence calculation
  print 'Select the strong pixels from the images.'
Пример #2
0
        print("FAIL: dials_regression not configured")
        raise

    # The XDS values
    xds_sigma_d = 0.060
    xds_sigma_m = 0.154

    # Read the reflections file
    filename = os.path.join(dials_regression, "centroid_test_data",
                            "reflections.h5")
    # filename = os.path.join('/home/upc86896/Data/X1_strong_M1S1_1_', 'reflections.h5')
    handle = NexusFile(filename, "r")

    # Get the reflection list
    print("Reading reflections.")
    predicted = handle.get_reflections()
    print("Read {0} reflections.".format(len(predicted)))

    # Read images
    template = os.path.join(dials_regression, "centroid_test_data",
                            "centroid_*.cbf")
    # template = os.path.join('/home/upc86896/Data/X1_strong_M1S1_1_', 'X1_strong_M1S1_1_*.cbf')
    filenames = glob(template)

    # Load the sweep
    print("Loading sweep")
    sweep = SweepFactory.sweep(filenames)
    print("Loaded sweep of {0} images.".format(len(sweep)))

    from dials.algorithms.peak_finding.spot_finder import SpotFinder
Пример #3
0
        print("FAIL: dials_regression not configured")
        raise

    # The XDS values
    xds_sigma_d = 0.060
    xds_sigma_m = 0.154

    # Read the reflections file
    filename = os.path.join(dials_regression, "centroid_test_data",
                            "reflections.h5")
    # filename = os.path.join('/home/upc86896/Data/X1_strong_M1S1_1_', 'reflections.h5')
    handle = NexusFile(filename, "r")

    # Get the reflection list
    print("Reading reflections.")
    reflections = handle.get_reflections()
    print("Read {0} reflections.".format(len(reflections)))

    # Read images
    template = os.path.join(dials_regression, "centroid_test_data",
                            "centroid_*.cbf")
    # template = os.path.join('/home/upc86896/Data/X1_strong_M1S1_1_', 'X1_strong_M1S1_1_*.cbf')
    filenames = glob(template)

    # Load the sweep
    print("Loading sweep")
    sweep = SweepFactory.sweep(filenames)
    print("Loaded sweep of {0} images.".format(len(sweep)))

    # Select the strong pixels to use in the divergence calculation
    print("Select the strong pixels from the images.")
Пример #4
0
    print 'FAIL: dials_regression not configured'
    raise

  # The XDS values
  xds_sigma_d = 0.060
  xds_sigma_m = 0.154

  # Read the reflections file
  filename = os.path.join(dials_regression,
      'centroid_test_data', 'reflections.h5')
  #filename = os.path.join('/home/upc86896/Data/X1_strong_M1S1_1_', 'reflections.h5')
  handle = NexusFile(filename, 'r')

  # Get the reflection list
  print 'Reading reflections.'
  predicted = handle.get_reflections()
  print 'Read {0} reflections.'.format(len(predicted))

  # Read images
  template = os.path.join(dials_regression,
      'centroid_test_data', 'centroid_*.cbf')
  #template = os.path.join('/home/upc86896/Data/X1_strong_M1S1_1_', 'X1_strong_M1S1_1_*.cbf')
  filenames = glob(template)

  # Load the sweep
  print 'Loading sweep'
  sweep = SweepFactory.sweep(filenames)
  print 'Loaded sweep of {0} images.'.format(len(sweep))

  from dials.algorithms.peak_finding.spot_finder import SpotFinder
  sweep.reader().set_max_cache(1)