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
