def numpy_rw_17(): if gdaltest.numpy_drv is None: return 'skip' import numpy from osgeo import gdal_array # Disabled by default array = numpy.empty([1, 1], numpy.uint8) with gdaltest.error_handler(): ds = gdal.Open(gdal_array.GetArrayFilename(array)) if ds is not None: gdaltest.post_reason('failure') return 'fail' gdal.SetConfigOption('GDAL_ARRAY_OPEN_BY_FILENAME', 'TRUE') ds = gdal.Open(gdal_array.GetArrayFilename(array)) gdal.SetConfigOption('GDAL_ARRAY_OPEN_BY_FILENAME', None) if ds is None: gdaltest.post_reason('failure') return 'fail' # Invalid value with gdaltest.error_handler(): ds = gdal.Open('NUMPY:::invalid') if ds is not None: gdaltest.post_reason('failure') return 'fail' return 'success'
def copy_geo(array, prototype=None, xoffset=0, yoffset=0): """Copy geotransfrom from prototype dataset to array but account for x, y offset of clipped array.""" ds = gdal.Open(gdal_array.GetArrayFilename(array)) prototype = gdal.Open(prototype) gdal_array.CopyDatasetInfo(prototype, ds, xoff=xoffset, yoff=yoffset) return ds
def test_numpy_rw_17(): # Disabled by default array = numpy.empty([1, 1], numpy.uint8) with gdaltest.error_handler(): ds = gdal.Open(gdal_array.GetArrayFilename(array)) assert ds is None gdal.SetConfigOption('GDAL_ARRAY_OPEN_BY_FILENAME', 'TRUE') ds = gdal.Open(gdal_array.GetArrayFilename(array)) gdal.SetConfigOption('GDAL_ARRAY_OPEN_BY_FILENAME', None) assert ds is not None # Invalid value with gdaltest.error_handler(): ds = gdal.Open('NUMPY:::invalid') assert ds is None
def Copy_Geo(array, prototype=None, xoffset=0, yoffset=0): #从原数据集拷贝数组 ds = gdal.Open(gdal_array.GetArrayFilename(array)) prototype = gdal.Open(prototype) gdal_array.CopyDatasetInfo(prototype, ds, xoff=xoffset, yoff=yoffset) return ds