def test_gdalbuildvrt_2():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
if test_cli_utilities.get_gdaltindex_path() is None:
return "skip"
try:
os.remove("tmp/tileindex.shp")
except:
pass
try:
os.remove("tmp/tileindex.dbf")
except:
pass
try:
os.remove("tmp/tileindex.shx")
except:
pass
try:
os.remove("tmp/mosaic.vrt")
except:
pass
gdaltest.runexternal(
test_cli_utilities.get_gdaltindex_path()
+ " tmp/tileindex.shp tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif"
)
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + " tmp/mosaic.vrt tmp/tileindex.shp")
return test_gdalbuildvrt_check()
def vrtmask_4():
import test_cli_utilities
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
gtiff_drv = gdal.GetDriverByName("GTiff")
md = gtiff_drv.GetMetadata()
if md["DMD_CREATIONOPTIONLIST"].find("JPEG") == -1:
return "skip"
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() + " tmp/vrtmask_4.vrt ../gcore/data/ycbcr_with_mask.tif"
)
src_ds = gdal.Open("../gcore/data/ycbcr_with_mask.tif")
expected_msk_cs = src_ds.GetRasterBand(1).GetMaskBand().Checksum()
src_ds = None
ds = gdal.Open("tmp/vrtmask_4.vrt")
msk_cs = ds.GetRasterBand(1).GetMaskBand().Checksum()
ds = None
os.remove("tmp/vrtmask_4.vrt")
if msk_cs != expected_msk_cs:
gdaltest.post_reason("did not get expected mask band checksum")
print(msk_cs)
print(expected_msk_cs)
return "fail"
return "success"
def test_gdalbuildvrt_6():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' -separate tmp/stacked.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif')
ds = gdal.Open('tmp/stacked.vrt')
if ds.GetProjectionRef().find('WGS 84') == -1:
gdaltest.post_reason('Expected WGS 84\nGot : %s' % (ds.GetProjectionRef()))
return 'fail'
gt = ds.GetGeoTransform()
expected_gt = [2, 0.1, 0, 49, 0, -0.1]
for i in range(6):
if abs(gt[i] - expected_gt[i] > 1e-5):
gdaltest.post_reason('Expected : %s\nGot : %s' % (expected_gt, gt))
return 'fail'
if ds.RasterXSize != 20 or ds.RasterYSize != 20:
gdaltest.post_reason('Wrong raster dimensions : %d x %d' % (ds.RasterXSize, ds.RasterYSize))
return 'fail'
if ds.RasterCount != 4:
gdaltest.post_reason('Wrong raster count : %d ' % (ds.RasterCount))
return 'fail'
if ds.GetRasterBand(1).Checksum() != 0:
gdaltest.post_reason('Wrong checksum')
return 'fail'
return 'success'
def test_gdalbuildvrt_16():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
(out, err) = gdaltest.runexternal_out_and_err(
test_cli_utilities.get_gdalbuildvrt_path() +
' /non_existing_dir/non_existing_subdir/out.vrt ../gcore/data/byte.tif'
)
if 'TRAVIS_BRANCH' in os.environ:
val = os.environ['TRAVIS_BRANCH']
else:
val = ''
if val.find('mingw') < 0:
if err.find('ERROR ret code = 1') < 0:
gdaltest.post_reason('fail')
print(out)
print(err)
return 'fail'
else:
# We don't get the error code on Travis mingw
if err.find('ERROR') < 0:
gdaltest.post_reason('fail')
print(out)
print(err)
return 'fail'
return 'success'
def test_gdalbuildvrt_11():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
out_ds = gdal.GetDriverByName('GTiff').Create(
'tmp/test_gdalbuildvrt_11_1.tif', 10, 10, 1)
out_ds.GetRasterBand(1).Fill(255)
cs1 = out_ds.GetRasterBand(1).Checksum()
out_ds = None
out_ds = gdal.GetDriverByName('GTiff').Create(
'tmp/test_gdalbuildvrt_11_2.tif', 10, 10, 1)
out_ds.GetRasterBand(1).Fill(127)
cs2 = out_ds.GetRasterBand(1).Checksum()
out_ds = None
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() +
' -separate tmp/gdalbuildvrt11.vrt tmp/test_gdalbuildvrt_11_1.tif tmp/test_gdalbuildvrt_11_2.tif'
)
ds = gdal.Open('tmp/gdalbuildvrt11.vrt')
assert ds.GetRasterBand(1).Checksum() == cs1, 'Wrong checksum'
assert ds.GetRasterBand(2).Checksum() == cs2, 'Wrong checksum'
ds = None
def test_gdalbuildvrt_7():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
out_ds = gdal.GetDriverByName("GTiff").Create("tmp/vrtnull1.tif", 20, 10, 3, gdal.GDT_UInt16)
out_ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
srs = osr.SpatialReference()
srs.SetFromUserInput("EPSG:4326")
out_ds.SetProjection(srs.ExportToWkt())
out_ds.GetRasterBand(1).SetRasterColorInterpretation(gdal.GCI_RedBand)
out_ds.GetRasterBand(2).SetRasterColorInterpretation(gdal.GCI_GreenBand)
out_ds.GetRasterBand(3).SetRasterColorInterpretation(gdal.GCI_BlueBand)
out_ds.GetRasterBand(1).SetNoDataValue(256)
try:
ff = "\xff".encode("latin1")
except:
ff = "\xff"
out_ds.GetRasterBand(1).WriteRaster(0, 0, 10, 10, ff, buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(2).WriteRaster(0, 0, 10, 10, "\x00", buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(3).WriteRaster(0, 0, 10, 10, "\x00", buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds = None
out_ds = gdal.GetDriverByName("GTiff").Create("tmp/vrtnull2.tif", 20, 10, 3, gdal.GDT_UInt16)
out_ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
srs = osr.SpatialReference()
srs.SetFromUserInput("EPSG:4326")
out_ds.SetProjection(srs.ExportToWkt())
out_ds.GetRasterBand(1).SetRasterColorInterpretation(gdal.GCI_RedBand)
out_ds.GetRasterBand(2).SetRasterColorInterpretation(gdal.GCI_GreenBand)
out_ds.GetRasterBand(3).SetRasterColorInterpretation(gdal.GCI_BlueBand)
out_ds.GetRasterBand(1).SetNoDataValue(256)
out_ds.GetRasterBand(1).WriteRaster(10, 0, 10, 10, "\x00", buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(2).WriteRaster(10, 0, 10, 10, ff, buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(3).WriteRaster(10, 0, 10, 10, "\x00", buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds = None
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() + " tmp/gdalbuildvrt7.vrt tmp/vrtnull1.tif tmp/vrtnull2.tif"
)
ds = gdal.Open("tmp/gdalbuildvrt7.vrt")
if ds.GetRasterBand(1).Checksum() != 1217:
gdaltest.post_reason("Wrong checksum")
return "fail"
if ds.GetRasterBand(2).Checksum() != 1218:
gdaltest.post_reason("Wrong checksum")
return "fail"
if ds.GetRasterBand(3).Checksum() != 0:
gdaltest.post_reason("Wrong checksum")
return "fail"
ds = None
return "success"
def vrtmask_4():
import test_cli_utilities
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
gtiff_drv = gdal.GetDriverByName('GTiff')
md = gtiff_drv.GetMetadata()
if md['DMD_CREATIONOPTIONLIST'].find('JPEG') == -1:
return 'skip'
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() +
' tmp/vrtmask_4.vrt ../gcore/data/ycbcr_with_mask.tif')
src_ds = gdal.Open('../gcore/data/ycbcr_with_mask.tif')
expected_msk_cs = src_ds.GetRasterBand(1).GetMaskBand().Checksum()
src_ds = None
ds = gdal.Open('tmp/vrtmask_4.vrt')
msk_cs = ds.GetRasterBand(1).GetMaskBand().Checksum()
ds = None
os.remove('tmp/vrtmask_4.vrt')
if msk_cs != expected_msk_cs:
gdaltest.post_reason('did not get expected mask band checksum')
print(msk_cs)
print(expected_msk_cs)
return 'fail'
return 'success'
def vrtmask_4():
import test_cli_utilities
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
gtiff_drv = gdal.GetDriverByName('GTiff')
md = gtiff_drv.GetMetadata()
if md['DMD_CREATIONOPTIONLIST'].find('JPEG') == -1:
return 'skip'
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/vrtmask_4.vrt ../gcore/data/ycbcr_with_mask.tif')
src_ds = gdal.Open('../gcore/data/ycbcr_with_mask.tif')
expected_msk_cs = src_ds.GetRasterBand(1).GetMaskBand().Checksum()
src_ds = None
ds = gdal.Open('tmp/vrtmask_4.vrt')
msk_cs = ds.GetRasterBand(1).GetMaskBand().Checksum()
ds = None
os.remove('tmp/vrtmask_4.vrt')
if msk_cs != expected_msk_cs:
gdaltest.post_reason('did not get expected mask band checksum')
print(msk_cs)
print(expected_msk_cs)
return 'fail'
return 'success'
def test_vrtmask_4():
import test_cli_utilities
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
gtiff_drv = gdal.GetDriverByName('GTiff')
md = gtiff_drv.GetMetadata()
if md['DMD_CREATIONOPTIONLIST'].find('JPEG') == -1:
pytest.skip()
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() +
' tmp/vrtmask_4.vrt ../gcore/data/ycbcr_with_mask.tif')
src_ds = gdal.Open('../gcore/data/ycbcr_with_mask.tif')
expected_msk_cs = src_ds.GetRasterBand(1).GetMaskBand().Checksum()
src_ds = None
ds = gdal.Open('tmp/vrtmask_4.vrt')
msk_cs = ds.GetRasterBand(1).GetMaskBand().Checksum()
ds = None
os.remove('tmp/vrtmask_4.vrt')
assert msk_cs == expected_msk_cs, 'did not get expected mask band checksum'
def test_gdalbuildvrt_11():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
out_ds = gdal.GetDriverByName('GTiff').Create('tmp/test_gdalbuildvrt_11_1.tif', 10, 10, 1)
out_ds.GetRasterBand(1).Fill(255)
cs1 = out_ds.GetRasterBand(1).Checksum()
out_ds = None
out_ds = gdal.GetDriverByName('GTiff').Create('tmp/test_gdalbuildvrt_11_2.tif', 10, 10, 1)
out_ds.GetRasterBand(1).Fill(127)
cs2 = out_ds.GetRasterBand(1).Checksum()
out_ds = None
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' -separate tmp/gdalbuildvrt11.vrt tmp/test_gdalbuildvrt_11_1.tif tmp/test_gdalbuildvrt_11_2.tif')
ds = gdal.Open('tmp/gdalbuildvrt11.vrt')
if ds.GetRasterBand(1).Checksum() != cs1:
print(ds.GetRasterBand(1).Checksum())
print(cs1)
gdaltest.post_reason('Wrong checksum')
return 'fail'
if ds.GetRasterBand(2).Checksum() != cs2:
print(ds.GetRasterBand(2).Checksum())
print(cs2)
gdaltest.post_reason('Wrong checksum')
return 'fail'
ds = None
return 'success'
def test_gdalbuildvrt_12():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
(out, err) = gdaltest.runexternal_out_and_err(
test_cli_utilities.get_gdalbuildvrt_path() +
' -tap tmp/gdalbuildvrt12.vrt ../gcore/data/byte.tif',
check_memleak=False)
if err.find('-tap option cannot be used without using -tr') == -1:
gdaltest.post_reason('expected error')
return 'fail'
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() +
' -tr 100 50 -tap tmp/gdalbuildvrt12.vrt ../gcore/data/byte.tif')
ds = gdal.Open('tmp/gdalbuildvrt12.vrt')
gt = ds.GetGeoTransform()
expected_gt = [440700.0, 100.0, 0.0, 3751350.0, 0.0, -50.0]
for i in range(6):
if abs(gt[i] - expected_gt[i] > 1e-5):
gdaltest.post_reason('Expected : %s\nGot : %s' % (expected_gt, gt))
return 'fail'
if ds.RasterXSize != 13 or ds.RasterYSize != 25:
gdaltest.post_reason('Wrong raster dimensions : %d x %d' %
(ds.RasterXSize, ds.RasterYSize))
return 'fail'
return 'success'
def test_gdalbuildvrt_14():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
if test_cli_utilities.get_gdal_translate_path() is None:
return 'skip'
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() +
' tmp/test_gdalbuildvrt_14.vrt ../gcore/data/byte.tif -r cubic -tr 30 30'
)
gdaltest.runexternal(
test_cli_utilities.get_gdal_translate_path() +
' -of VRT ../gcore/data/byte.tif tmp/test_gdalbuildvrt_14_ref.vrt -r cubic -outsize 40 40'
)
ds = gdal.Open('tmp/test_gdalbuildvrt_14.vrt')
ds_ref = gdal.Open('tmp/test_gdalbuildvrt_14_ref.vrt')
cs = ds.GetRasterBand(1).Checksum()
cs_ref = ds_ref.GetRasterBand(1).Checksum()
ds = None
ds_ref = None
if cs != cs_ref:
gdaltest.post_reason('fail')
print(cs)
print(cs_ref)
return 'fail'
return 'success'
def test_gdalbuildvrt_11():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
out_ds = gdal.GetDriverByName('GTiff').Create('tmp/test_gdalbuildvrt_11_1.tif', 10, 10, 1)
out_ds.GetRasterBand(1).Fill(255)
cs1 = out_ds.GetRasterBand(1).Checksum()
out_ds = None
out_ds = gdal.GetDriverByName('GTiff').Create('tmp/test_gdalbuildvrt_11_2.tif', 10, 10, 1)
out_ds.GetRasterBand(1).Fill(127)
cs2 = out_ds.GetRasterBand(1).Checksum()
out_ds = None
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' -separate tmp/gdalbuildvrt11.vrt tmp/test_gdalbuildvrt_11_1.tif tmp/test_gdalbuildvrt_11_2.tif')
ds = gdal.Open('tmp/gdalbuildvrt11.vrt')
if ds.GetRasterBand(1).Checksum() != cs1:
print(ds.GetRasterBand(1).Checksum())
print(cs1)
gdaltest.post_reason('Wrong checksum')
return 'fail'
if ds.GetRasterBand(2).Checksum() != cs2:
print(ds.GetRasterBand(2).Checksum())
print(cs2)
gdaltest.post_reason('Wrong checksum')
return 'fail'
ds = None
return 'success'
def test_gdalbuildvrt_16():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
(out, err) = gdaltest.runexternal_out_and_err(test_cli_utilities.get_gdalbuildvrt_path() + ' /non_existing_dir/non_existing_subdir/out.vrt ../gcore/data/byte.tif')
if 'TRAVIS_BRANCH' in os.environ:
val = os.environ['TRAVIS_BRANCH']
else:
val = ''
if val.find('mingw') < 0:
if err.find('ERROR ret code = 1') < 0:
gdaltest.post_reason('fail')
print(out)
print(err)
return 'fail'
else:
# We don't get the error code on Travis mingw
if err.find('ERROR') < 0:
gdaltest.post_reason('fail')
print(out)
print(err)
return 'fail'
return 'success'
def test_gdalbuildvrt_2():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
if test_cli_utilities.get_gdaltindex_path() is None:
return 'skip'
try:
os.remove('tmp/tileindex.shp')
except:
pass
try:
os.remove('tmp/tileindex.dbf')
except:
pass
try:
os.remove('tmp/tileindex.shx')
except:
pass
try:
os.remove('tmp/mosaic.vrt')
except:
pass
gdaltest.runexternal(
test_cli_utilities.get_gdaltindex_path() +
' tmp/tileindex.shp tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif'
)
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() +
' tmp/mosaic.vrt tmp/tileindex.shp')
return test_gdalbuildvrt_check()
def test_gdalbuildvrt_2():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
if test_cli_utilities.get_gdaltindex_path() is None:
pytest.skip()
try:
os.remove('tmp/tileindex.shp')
except OSError:
pass
try:
os.remove('tmp/tileindex.dbf')
except OSError:
pass
try:
os.remove('tmp/tileindex.shx')
except OSError:
pass
try:
os.remove('tmp/mosaic.vrt')
except OSError:
pass
gdaltest.runexternal(test_cli_utilities.get_gdaltindex_path() + ' tmp/tileindex.shp tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif')
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/mosaic.vrt tmp/tileindex.shp')
return gdalbuildvrt_check()
def test_gdalbuildvrt_14():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
if test_cli_utilities.get_gdal_translate_path() is None:
return "skip"
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path()
+ " tmp/test_gdalbuildvrt_14.vrt ../gcore/data/byte.tif -r cubic -tr 30 30"
)
gdaltest.runexternal(
test_cli_utilities.get_gdal_translate_path()
+ " -of VRT ../gcore/data/byte.tif tmp/test_gdalbuildvrt_14_ref.vrt -r cubic -outsize 40 40"
)
ds = gdal.Open("tmp/test_gdalbuildvrt_14.vrt")
ds_ref = gdal.Open("tmp/test_gdalbuildvrt_14_ref.vrt")
cs = ds.GetRasterBand(1).Checksum()
cs_ref = ds_ref.GetRasterBand(1).Checksum()
ds = None
ds_ref = None
if cs != cs_ref:
gdaltest.post_reason("fail")
print(cs)
print(cs_ref)
return "fail"
return "success"
def test_gdalbuildvrt_10():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
out_ds = gdal.GetDriverByName('GTiff').Create('tmp/test_gdalbuildvrt_10_1.tif', 10, 10, 1, gdal.GDT_Byte, options=['NBITS=1', 'PHOTOMETRIC=MINISWHITE'])
out_ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
srs = osr.SpatialReference()
srs.SetFromUserInput('EPSG:4326')
out_ds.SetProjection(srs.ExportToWkt())
out_ds.GetRasterBand(1).WriteRaster(1, 1, 3, 3, '\x01', buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds = None
out_ds = gdal.GetDriverByName('GTiff').Create('tmp/test_gdalbuildvrt_10_2.tif', 10, 10, 1, gdal.GDT_Byte, options=['NBITS=1', 'PHOTOMETRIC=MINISWHITE'])
out_ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
srs = osr.SpatialReference()
srs.SetFromUserInput('EPSG:4326')
out_ds.SetProjection(srs.ExportToWkt())
out_ds.GetRasterBand(1).WriteRaster(6, 6, 3, 3, '\x01', buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds = None
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' -srcnodata 0 tmp/gdalbuildvrt10.vrt tmp/test_gdalbuildvrt_10_1.tif tmp/test_gdalbuildvrt_10_2.tif')
ds = gdal.Open('tmp/gdalbuildvrt10.vrt')
assert ds.GetRasterBand(1).Checksum() == 18, 'Wrong checksum'
ds = None
def test_gdalbuildvrt_12():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
(out, err) = gdaltest.runexternal_out_and_err(test_cli_utilities.get_gdalbuildvrt_path() + ' -tap tmp/gdalbuildvrt12.vrt ../gcore/data/byte.tif',
check_memleak = False)
if err.find('-tap option cannot be used without using -tr') == -1:
gdaltest.post_reason('expected error')
return 'fail'
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' -tr 100 50 -tap tmp/gdalbuildvrt12.vrt ../gcore/data/byte.tif')
ds = gdal.Open('tmp/gdalbuildvrt12.vrt')
gt = ds.GetGeoTransform()
expected_gt = [ 440700.0, 100.0, 0.0, 3751350.0, 0.0, -50.0 ]
for i in range(6):
if abs(gt[i] - expected_gt[i] > 1e-5):
gdaltest.post_reason('Expected : %s\nGot : %s' % (expected_gt, gt) )
return 'fail'
if ds.RasterXSize != 13 or ds.RasterYSize != 25:
gdaltest.post_reason('Wrong raster dimensions : %d x %d' % (ds.RasterXSize, ds.RasterYSize) )
return 'fail'
return 'success'
def test_gdalbuildvrt_12():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
(_, err) = gdaltest.runexternal_out_and_err(
test_cli_utilities.get_gdalbuildvrt_path() +
' -tap tmp/gdalbuildvrt12.vrt ../gcore/data/byte.tif',
check_memleak=False)
assert err.find('-tap option cannot be used without using -tr') != -1, \
'expected error'
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() +
' -tr 100 50 -tap tmp/gdalbuildvrt12.vrt ../gcore/data/byte.tif')
ds = gdal.Open('tmp/gdalbuildvrt12.vrt')
gt = ds.GetGeoTransform()
expected_gt = [440700.0, 100.0, 0.0, 3751350.0, 0.0, -50.0]
for i in range(6):
assert not abs(gt[i] - expected_gt[i] > 1e-5), \
('Expected : %s\nGot : %s' % (expected_gt, gt))
assert ds.RasterXSize == 13 and ds.RasterYSize == 25, \
('Wrong raster dimensions : %d x %d' % (ds.RasterXSize, ds.RasterYSize))
def test_gdalbuildvrt_12():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
(out, err) = gdaltest.runexternal_out_and_err(
test_cli_utilities.get_gdalbuildvrt_path() + " -tap tmp/gdalbuildvrt12.vrt ../gcore/data/byte.tif",
check_memleak=False,
)
if err.find("-tap option cannot be used without using -tr") == -1:
gdaltest.post_reason("expected error")
return "fail"
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() + " -tr 100 50 -tap tmp/gdalbuildvrt12.vrt ../gcore/data/byte.tif"
)
ds = gdal.Open("tmp/gdalbuildvrt12.vrt")
gt = ds.GetGeoTransform()
expected_gt = [440700.0, 100.0, 0.0, 3751350.0, 0.0, -50.0]
for i in range(6):
if abs(gt[i] - expected_gt[i] > 1e-5):
gdaltest.post_reason("Expected : %s\nGot : %s" % (expected_gt, gt))
return "fail"
if ds.RasterXSize != 13 or ds.RasterYSize != 25:
gdaltest.post_reason("Wrong raster dimensions : %d x %d" % (ds.RasterXSize, ds.RasterYSize))
return "fail"
return "success"
def test_gdalbuildvrt_9():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path()
+ " -te 1 46 5 50 tmp/mosaic2.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif"
)
ds = gdal.Open("tmp/mosaic2.vrt")
gt = ds.GetGeoTransform()
expected_gt = [1, 0.1, 0, 50, 0, -0.1]
for i in range(6):
if abs(gt[i] - expected_gt[i] > 1e-5):
gdaltest.post_reason("Expected : %s\nGot : %s" % (expected_gt, gt))
return "fail"
if ds.RasterXSize != 40 or ds.RasterYSize != 40:
gdaltest.post_reason("Wrong raster dimensions : %d x %d" % (ds.RasterXSize, ds.RasterYSize))
return "fail"
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + " -te 2 47 4 49 tmp/mosaic.vrt tmp/mosaic2.vrt")
return test_gdalbuildvrt_check()
def test_gdalbuildvrt_6():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() +
' -separate tmp/stacked.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif'
)
ds = gdal.Open('tmp/stacked.vrt')
assert ds.GetProjectionRef().find('WGS 84') != -1, \
('Expected WGS 84\nGot : %s' % (ds.GetProjectionRef()))
gt = ds.GetGeoTransform()
expected_gt = [2, 0.1, 0, 49, 0, -0.1]
for i in range(6):
assert not abs(gt[i] - expected_gt[i] > 1e-5), \
('Expected : %s\nGot : %s' % (expected_gt, gt))
assert ds.RasterXSize == 20 and ds.RasterYSize == 20, \
('Wrong raster dimensions : %d x %d' % (ds.RasterXSize, ds.RasterYSize))
assert ds.RasterCount == 4, ('Wrong raster count : %d ' % (ds.RasterCount))
assert ds.GetRasterBand(1).Checksum() == 0, 'Wrong checksum'
def test_gdalbuildvrt_6():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path()
+ " -separate tmp/stacked.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif"
)
ds = gdal.Open("tmp/stacked.vrt")
if ds.GetProjectionRef().find("WGS 84") == -1:
gdaltest.post_reason("Expected WGS 84\nGot : %s" % (ds.GetProjectionRef()))
return "fail"
gt = ds.GetGeoTransform()
expected_gt = [2, 0.1, 0, 49, 0, -0.1]
for i in range(6):
if abs(gt[i] - expected_gt[i] > 1e-5):
gdaltest.post_reason("Expected : %s\nGot : %s" % (expected_gt, gt))
return "fail"
if ds.RasterXSize != 20 or ds.RasterYSize != 20:
gdaltest.post_reason("Wrong raster dimensions : %d x %d" % (ds.RasterXSize, ds.RasterYSize))
return "fail"
if ds.RasterCount != 4:
gdaltest.post_reason("Wrong raster count : %d " % (ds.RasterCount))
return "fail"
if ds.GetRasterBand(1).Checksum() != 0:
gdaltest.post_reason("Wrong checksum")
return "fail"
return "success"
def test_gdalbuildvrt_9():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() +
' -te 1 46 5 50 tmp/mosaic2.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif'
)
ds = gdal.Open('tmp/mosaic2.vrt')
gt = ds.GetGeoTransform()
expected_gt = [1, 0.1, 0, 50, 0, -0.1]
for i in range(6):
if abs(gt[i] - expected_gt[i] > 1e-5):
gdaltest.post_reason('Expected : %s\nGot : %s' % (expected_gt, gt))
return 'fail'
if ds.RasterXSize != 40 or ds.RasterYSize != 40:
gdaltest.post_reason('Wrong raster dimensions : %d x %d' %
(ds.RasterXSize, ds.RasterYSize))
return 'fail'
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() +
' -te 2 47 4 49 tmp/mosaic.vrt tmp/mosaic2.vrt')
return test_gdalbuildvrt_check()
def test_gdalbuildvrt_7():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
out_ds = gdal.GetDriverByName('GTiff').Create('tmp/vrtnull1.tif', 20, 10, 3, gdal.GDT_UInt16)
out_ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
srs = osr.SpatialReference()
srs.SetFromUserInput('EPSG:4326')
out_ds.SetProjection(srs.ExportToWkt())
out_ds.GetRasterBand(1).SetRasterColorInterpretation(gdal.GCI_RedBand)
out_ds.GetRasterBand(2).SetRasterColorInterpretation(gdal.GCI_GreenBand)
out_ds.GetRasterBand(3).SetRasterColorInterpretation(gdal.GCI_BlueBand)
out_ds.GetRasterBand(1).SetNoDataValue(256)
try:
ff = '\xff'.encode('latin1')
except UnicodeDecodeError:
ff = '\xff'
out_ds.GetRasterBand(1).WriteRaster(0, 0, 10, 10, ff, buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(2).WriteRaster(0, 0, 10, 10, '\x00', buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(3).WriteRaster(0, 0, 10, 10, '\x00', buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds = None
out_ds = gdal.GetDriverByName('GTiff').Create('tmp/vrtnull2.tif', 20, 10, 3, gdal.GDT_UInt16)
out_ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
srs = osr.SpatialReference()
srs.SetFromUserInput('EPSG:4326')
out_ds.SetProjection(srs.ExportToWkt())
out_ds.GetRasterBand(1).SetRasterColorInterpretation(gdal.GCI_RedBand)
out_ds.GetRasterBand(2).SetRasterColorInterpretation(gdal.GCI_GreenBand)
out_ds.GetRasterBand(3).SetRasterColorInterpretation(gdal.GCI_BlueBand)
out_ds.GetRasterBand(1).SetNoDataValue(256)
out_ds.GetRasterBand(1).WriteRaster(10, 0, 10, 10, '\x00', buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(2).WriteRaster(10, 0, 10, 10, ff, buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(3).WriteRaster(10, 0, 10, 10, '\x00', buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds = None
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/gdalbuildvrt7.vrt tmp/vrtnull1.tif tmp/vrtnull2.tif')
ds = gdal.Open('tmp/gdalbuildvrt7.vrt')
if ds.GetRasterBand(1).Checksum() != 1217:
gdaltest.post_reason('Wrong checksum')
return 'fail'
if ds.GetRasterBand(2).Checksum() != 1218:
gdaltest.post_reason('Wrong checksum')
return 'fail'
if ds.GetRasterBand(3).Checksum() != 0:
gdaltest.post_reason('Wrong checksum')
return 'fail'
ds = None
return 'success'
def test_gdalbuildvrt_7():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
out_ds = gdal.GetDriverByName('GTiff').Create('tmp/vrtnull1.tif', 20, 10, 3, gdal.GDT_UInt16)
out_ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
srs = osr.SpatialReference()
srs.SetFromUserInput('EPSG:4326')
out_ds.SetProjection(srs.ExportToWkt())
out_ds.GetRasterBand(1).SetRasterColorInterpretation(gdal.GCI_RedBand)
out_ds.GetRasterBand(2).SetRasterColorInterpretation(gdal.GCI_GreenBand)
out_ds.GetRasterBand(3).SetRasterColorInterpretation(gdal.GCI_BlueBand)
out_ds.GetRasterBand(1).SetNoDataValue(256)
try:
ff = '\xff'.encode('latin1')
except:
ff = '\xff'
out_ds.GetRasterBand(1).WriteRaster(0, 0, 10, 10, ff, buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(2).WriteRaster(0, 0, 10, 10, '\x00', buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(3).WriteRaster(0, 0, 10, 10, '\x00', buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds = None
out_ds = gdal.GetDriverByName('GTiff').Create('tmp/vrtnull2.tif', 20, 10, 3, gdal.GDT_UInt16)
out_ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
srs = osr.SpatialReference()
srs.SetFromUserInput('EPSG:4326')
out_ds.SetProjection(srs.ExportToWkt())
out_ds.GetRasterBand(1).SetRasterColorInterpretation(gdal.GCI_RedBand)
out_ds.GetRasterBand(2).SetRasterColorInterpretation(gdal.GCI_GreenBand)
out_ds.GetRasterBand(3).SetRasterColorInterpretation(gdal.GCI_BlueBand)
out_ds.GetRasterBand(1).SetNoDataValue(256)
out_ds.GetRasterBand(1).WriteRaster(10, 0, 10, 10, '\x00', buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(2).WriteRaster(10, 0, 10, 10, ff, buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds.GetRasterBand(3).WriteRaster(10, 0, 10, 10, '\x00', buf_type=gdal.GDT_Byte, buf_xsize=1, buf_ysize=1)
out_ds = None
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/gdalbuildvrt7.vrt tmp/vrtnull1.tif tmp/vrtnull2.tif')
ds = gdal.Open('tmp/gdalbuildvrt7.vrt')
if ds.GetRasterBand(1).Checksum() != 1217:
gdaltest.post_reason('Wrong checksum')
return 'fail'
if ds.GetRasterBand(2).Checksum() != 1218:
gdaltest.post_reason('Wrong checksum')
return 'fail'
if ds.GetRasterBand(3).Checksum() != 0:
gdaltest.post_reason('Wrong checksum')
return 'fail'
ds = None
return 'success'
def test_gdalbuildvrt_13():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/gdalbuildvrt13.vrt ../gcore/data/byte.tif -a_srs EPSG:4326')
ds = gdal.Open('tmp/gdalbuildvrt13.vrt')
assert ds.GetProjectionRef().find('4326') != -1
ds = None
def test_gdalbuildvrt_3():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
open('tmp/filelist.txt', 'wt').write('tmp/gdalbuildvrt1.tif\ntmp/gdalbuildvrt2.tif\ntmp/gdalbuildvrt3.tif\ntmp/gdalbuildvrt4.tif')
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' -input_file_list tmp/filelist.txt tmp/mosaic.vrt')
return gdalbuildvrt_check()
def test_gdalbuildvrt_16():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
(out, err) = gdaltest.runexternal_out_and_err(test_cli_utilities.get_gdalbuildvrt_path() + ' /non_existing_dir/non_existing_subdir/out.vrt ../gcore/data/byte.tif')
if not gdaltest.is_travis_branch('mingw'):
assert 'ERROR ret code = 1' in err, out
else:
# We don't get the error code on Travis mingw
assert 'ERROR' in err, out
def test_gdalbuildvrt_15():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/test_gdalbuildvrt_15.vrt ../gcore/data/byte.tif -b 1')
ds = gdal.Open('tmp/test_gdalbuildvrt_15.vrt')
cs = ds.GetRasterBand(1).Checksum()
ds = None
assert cs == 4672
def test_gdalbuildvrt_13():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/gdalbuildvrt13.vrt ../gcore/data/byte.tif -a_srs EPSG:4326')
ds = gdal.Open('tmp/gdalbuildvrt13.vrt')
if ds.GetProjectionRef().find('4326') == -1:
return 'fail'
ds = None
return 'success'
def test_gdalbuildvrt_13():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/gdalbuildvrt13.vrt ../gcore/data/byte.tif -a_srs EPSG:4326')
ds = gdal.Open('tmp/gdalbuildvrt13.vrt')
if ds.GetProjectionRef().find('4326') == -1:
return 'fail'
ds = None
return 'success'
def test_gdalbuildvrt_3():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
open("tmp/filelist.txt", "wt").write(
"tmp/gdalbuildvrt1.tif\ntmp/gdalbuildvrt2.tif\ntmp/gdalbuildvrt3.tif\ntmp/gdalbuildvrt4.tif"
)
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() + " -input_file_list tmp/filelist.txt tmp/mosaic.vrt"
)
return test_gdalbuildvrt_check()
def test_gdalbuildvrt_13():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() + " tmp/gdalbuildvrt13.vrt ../gcore/data/byte.tif -a_srs EPSG:4326"
)
ds = gdal.Open("tmp/gdalbuildvrt13.vrt")
if ds.GetProjectionRef().find("4326") == -1:
return "fail"
ds = None
return "success"
def test_gdalbuildvrt_4():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
drv = gdal.GetDriverByName('GTiff')
wkt = 'GEOGCS[\"WGS 72\",DATUM[\"WGS_1972\"]]'
ds = drv.Create('tmp/gdalbuildvrt5.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([47, 0.1, 0, 2, 0, -0.1])
ds = None
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/mosaic.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif tmp/gdalbuildvrt5.tif')
return gdalbuildvrt_check()
def test_gdalbuildvrt_15():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/test_gdalbuildvrt_15.vrt ../gcore/data/byte.tif -b 1')
ds = gdal.Open('tmp/test_gdalbuildvrt_15.vrt')
cs = ds.GetRasterBand(1).Checksum()
ds = None
if cs != 4672:
gdaltest.post_reason('fail')
print(cs)
return 'fail'
return 'success'
def test_gdalbuildvrt_cleanup():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
ogr.GetDriverByName('ESRI Shapefile').DeleteDataSource('tmp/tileindex.shp')
gdal.GetDriverByName('VRT').Delete('tmp/mosaic.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/mosaic2.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/stacked.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/gdalbuildvrt7.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/gdalbuildvrt10.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/gdalbuildvrt11.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/gdalbuildvrt12.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/gdalbuildvrt13.vrt')
drv = gdal.GetDriverByName('GTiff')
drv.Delete('tmp/gdalbuildvrt1.tif')
drv.Delete('tmp/gdalbuildvrt2.tif')
drv.Delete('tmp/gdalbuildvrt3.tif')
drv.Delete('tmp/gdalbuildvrt4.tif')
drv.Delete('tmp/gdalbuildvrt5.tif')
drv.Delete('tmp/vrtnull1.tif')
drv.Delete('tmp/vrtnull2.tif')
drv.Delete('tmp/test_gdalbuildvrt_10_1.tif')
drv.Delete('tmp/test_gdalbuildvrt_10_2.tif')
drv.Delete('tmp/test_gdalbuildvrt_11_1.tif')
drv.Delete('tmp/test_gdalbuildvrt_11_2.tif')
try:
os.remove('tmp/filelist.txt')
except:
pass
return 'success'
def test_gdalbuildvrt_cleanup():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
ogr.GetDriverByName('ESRI Shapefile').DeleteDataSource('tmp/tileindex.shp')
gdal.GetDriverByName('VRT').Delete('tmp/mosaic.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/mosaic2.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/stacked.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/gdalbuildvrt7.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/gdalbuildvrt10.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/gdalbuildvrt11.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/gdalbuildvrt12.vrt')
gdal.GetDriverByName('VRT').Delete('tmp/gdalbuildvrt13.vrt')
drv = gdal.GetDriverByName('GTiff')
drv.Delete('tmp/gdalbuildvrt1.tif')
drv.Delete('tmp/gdalbuildvrt2.tif')
drv.Delete('tmp/gdalbuildvrt3.tif')
drv.Delete('tmp/gdalbuildvrt4.tif')
drv.Delete('tmp/gdalbuildvrt5.tif')
drv.Delete('tmp/vrtnull1.tif')
drv.Delete('tmp/vrtnull2.tif')
drv.Delete('tmp/test_gdalbuildvrt_10_1.tif')
drv.Delete('tmp/test_gdalbuildvrt_10_2.tif')
drv.Delete('tmp/test_gdalbuildvrt_11_1.tif')
drv.Delete('tmp/test_gdalbuildvrt_11_2.tif')
try:
os.remove('tmp/filelist.txt')
except:
pass
return 'success'
def test_gdalbuildvrt_15():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/test_gdalbuildvrt_15.vrt ../gcore/data/byte.tif -b 1')
ds = gdal.Open('tmp/test_gdalbuildvrt_15.vrt')
cs = ds.GetRasterBand(1).Checksum()
ds = None
if cs != 4672:
gdaltest.post_reason('fail')
print(cs)
return 'fail'
return 'success'
def test_gdalbuildvrt_5():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
drv = gdal.GetDriverByName('GTiff')
srs = osr.SpatialReference()
srs.SetWellKnownGeogCS('WGS84')
wkt = srs.ExportToWkt()
ds = drv.Create('tmp/gdalbuildvrt5.tif', 10, 10, 2)
ds.SetProjection(wkt)
ds.SetGeoTransform([47, 0.1, 0, 2, 0, -0.1])
ds = None
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/mosaic.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif tmp/gdalbuildvrt5.tif')
return gdalbuildvrt_check()
def test_gdalbuildvrt_14():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
if test_cli_utilities.get_gdal_translate_path() is None:
pytest.skip()
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/test_gdalbuildvrt_14.vrt ../gcore/data/byte.tif -r cubic -tr 30 30')
gdaltest.runexternal(test_cli_utilities.get_gdal_translate_path() + ' -of VRT ../gcore/data/byte.tif tmp/test_gdalbuildvrt_14_ref.vrt -r cubic -outsize 40 40')
ds = gdal.Open('tmp/test_gdalbuildvrt_14.vrt')
ds_ref = gdal.Open('tmp/test_gdalbuildvrt_14_ref.vrt')
cs = ds.GetRasterBand(1).Checksum()
cs_ref = ds_ref.GetRasterBand(1).Checksum()
ds = None
ds_ref = None
assert cs == cs_ref
def test_gdalbuildvrt_4():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
drv = gdal.GetDriverByName("GTiff")
wkt = 'GEOGCS["WGS 72",DATUM["WGS_1972"]]'
ds = drv.Create("tmp/gdalbuildvrt5.tif", 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([47, 0.1, 0, 2, 0, -0.1])
ds = None
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path()
+ " tmp/mosaic.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif tmp/gdalbuildvrt5.tif"
)
return test_gdalbuildvrt_check()
def test_gdalbuildvrt_1():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return "skip"
drv = gdal.GetDriverByName("GTiff")
srs = osr.SpatialReference()
srs.SetWellKnownGeogCS("WGS84")
wkt = srs.ExportToWkt()
ds = drv.Create("tmp/gdalbuildvrt1.tif", 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
ds.GetRasterBand(1).Fill(0)
ds = None
ds = drv.Create("tmp/gdalbuildvrt2.tif", 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([3, 0.1, 0, 49, 0, -0.1])
ds.GetRasterBand(1).Fill(63)
ds = None
ds = drv.Create("tmp/gdalbuildvrt3.tif", 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([2, 0.1, 0, 48, 0, -0.1])
ds.GetRasterBand(1).Fill(127)
ds = None
ds = drv.Create("tmp/gdalbuildvrt4.tif", 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([3, 0.1, 0, 48, 0, -0.1])
ds.GetRasterBand(1).Fill(255)
ds = None
(out, err) = gdaltest.runexternal_out_and_err(
test_cli_utilities.get_gdalbuildvrt_path()
+ " tmp/mosaic.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif"
)
if not (err is None or err == ""):
gdaltest.post_reason("got error/warning")
print(err)
return "fail"
return test_gdalbuildvrt_check()
def test_gdalbuildvrt_1():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
drv = gdal.GetDriverByName('GTiff')
srs = osr.SpatialReference()
srs.SetWellKnownGeogCS('WGS84')
wkt = srs.ExportToWkt()
ds = drv.Create('tmp/gdalbuildvrt1.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
ds.GetRasterBand(1).Fill(0)
ds = None
ds = drv.Create('tmp/gdalbuildvrt2.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([3, 0.1, 0, 49, 0, -0.1])
ds.GetRasterBand(1).Fill(63)
ds = None
ds = drv.Create('tmp/gdalbuildvrt3.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([2, 0.1, 0, 48, 0, -0.1])
ds.GetRasterBand(1).Fill(127)
ds = None
ds = drv.Create('tmp/gdalbuildvrt4.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([3, 0.1, 0, 48, 0, -0.1])
ds.GetRasterBand(1).Fill(255)
ds = None
(out, err) = gdaltest.runexternal_out_and_err(
test_cli_utilities.get_gdalbuildvrt_path() +
' tmp/mosaic.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif'
)
if not (err is None or err == ''):
gdaltest.post_reason('got error/warning')
print(err)
return 'fail'
return test_gdalbuildvrt_check()
def test_gdalbuildvrt_4():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
drv = gdal.GetDriverByName('GTiff')
wkt = """GEOGCS["WGS 72",
DATUM["WGS_1972",
SPHEROID["WGS 72",6378135,298.26],
TOWGS84[0,0,4.5,0,0,0.554,0.2263]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]]"""
ds = drv.Create('tmp/gdalbuildvrt5.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([47, 0.1, 0, 2, 0, -0.1])
ds = None
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' tmp/mosaic.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif tmp/gdalbuildvrt5.tif')
return gdalbuildvrt_check()
def test_gdalbuildvrt_9():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' -te 1 46 5 50 tmp/mosaic2.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif')
ds = gdal.Open('tmp/mosaic2.vrt')
gt = ds.GetGeoTransform()
expected_gt = [1, 0.1, 0, 50, 0, -0.1]
for i in range(6):
assert not abs(gt[i] - expected_gt[i] > 1e-5), \
('Expected : %s\nGot : %s' % (expected_gt, gt))
assert ds.RasterXSize == 40 and ds.RasterYSize == 40, \
('Wrong raster dimensions : %d x %d' % (ds.RasterXSize, ds.RasterYSize))
gdaltest.runexternal(test_cli_utilities.get_gdalbuildvrt_path() + ' -te 2 47 4 49 tmp/mosaic.vrt tmp/mosaic2.vrt')
return gdalbuildvrt_check()
def test_gdalbuildvrt_1():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
return 'skip'
drv = gdal.GetDriverByName('GTiff')
srs = osr.SpatialReference()
srs.SetWellKnownGeogCS('WGS84')
wkt = srs.ExportToWkt()
ds = drv.Create('tmp/gdalbuildvrt1.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([2, 0.1, 0, 49, 0, -0.1])
ds.GetRasterBand(1).Fill(0)
ds = None
ds = drv.Create('tmp/gdalbuildvrt2.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([3, 0.1, 0, 49, 0, -0.1])
ds.GetRasterBand(1).Fill(63)
ds = None
ds = drv.Create('tmp/gdalbuildvrt3.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([2, 0.1, 0, 48, 0, -0.1])
ds.GetRasterBand(1).Fill(127)
ds = None
ds = drv.Create('tmp/gdalbuildvrt4.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([3, 0.1, 0, 48, 0, -0.1])
ds.GetRasterBand(1).Fill(255)
ds = None
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() +
' tmp/mosaic.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif'
)
return test_gdalbuildvrt_check()
def test_gdalbuildvrt_4():
if test_cli_utilities.get_gdalbuildvrt_path() is None:
pytest.skip()
drv = gdal.GetDriverByName('GTiff')
wkt = """GEOGCS["WGS 72",
DATUM["WGS_1972",
SPHEROID["WGS 72",6378135,298.26],
TOWGS84[0,0,4.5,0,0,0.554,0.2263]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]]"""
ds = drv.Create('tmp/gdalbuildvrt5.tif', 10, 10, 1)
ds.SetProjection(wkt)
ds.SetGeoTransform([47, 0.1, 0, 2, 0, -0.1])
ds = None
gdaltest.runexternal(
test_cli_utilities.get_gdalbuildvrt_path() +
' tmp/mosaic.vrt tmp/gdalbuildvrt1.tif tmp/gdalbuildvrt2.tif tmp/gdalbuildvrt3.tif tmp/gdalbuildvrt4.tif tmp/gdalbuildvrt5.tif'
)
return gdalbuildvrt_check()
