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_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_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_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()