def test_ingr_10():
tst = gdaltest.GDALTest('INGR', 'ingr/frmt24.cit', 1, 23035)
return tst.testOpen()
def test_gtiff_create(filename, checksum, testfunction):
ut = gdaltest.GDALTest('GTiff', filename, 1, checksum)
getattr(ut, testfunction)()
def byn_2():
tst = gdaltest.GDALTest('BYN', 'cgg2013ai08_reduced.byn', 1, 64764)
return tst.testCreateCopy(new_filename='tmp/byn_test_2.byn')
def test_iris_1():
tst = gdaltest.GDALTest('IRIS', 'fakeiris.dat', 1, 65532)
return tst.testOpen()
def test_ignfheightasciigrid_7():
tst = gdaltest.GDALTest('IGNFHeightASCIIGrid',
'ignfheightasciigrid_ar1_noprec.mnt', 1, 21)
gt = (-152.125, 0.25, 0.0, -16.375, 0.0, -0.25)
return tst.testOpen(check_gt=gt, check_prj='WGS84')
def blx_2():
prj = 'WGS84'
gt = [ 20.0004166, 0.0008333, 0.0, 50.0004166, 0.0, -0.0008333 ]
tst = gdaltest.GDALTest( 'BLX', 's4103.xlb', 1, 47024 )
return tst.testOpen( check_prj = prj, check_gt = gt )
def blx_4():
tst = gdaltest.GDALTest( 'BLX', 's4103.blx', 1, 47024, options = [ 'BIGENDIAN=YES' ] )
return tst.testCreateCopy( check_gt = 1, check_srs = 1 )
def test_ingr_6():
tst = gdaltest.GDALTest('INGR', 'ingr/frmt02.cot', 1, 26968)
return tst.testCreate()
os.symlink('../data/byte.vrt', 'tmp/byte.vrt')
os.symlink('../tmp/byte.vrt', 'tmp/other_byte.vrt')
ds = gdal.Open('tmp/other_byte.vrt')
os.remove('tmp/other_byte.vrt')
os.remove('tmp/byte.vrt')
if ds is None:
return 'fail'
return 'success'
for item in init_list:
ut = gdaltest.GDALTest( 'VRT', item[0], item[1], item[2] )
if ut is None:
print( 'VRT tests skipped' )
sys.exit()
gdaltest_list.append( (ut.testOpen, item[0]) )
gdaltest_list.append( vrt_read_1 )
gdaltest_list.append( vrt_read_2 )
gdaltest_list.append( vrt_read_3 )
gdaltest_list.append( vrt_read_4 )
gdaltest_list.append( vrt_read_5 )
gdaltest_list.append( vrt_read_6 )
gdaltest_list.append( vrt_read_7 )
gdaltest_list.append( vrt_read_8 )
gdaltest_list.append( vrt_read_9 )
gdaltest_list.append( vrt_read_10 )
def test_ingr_2():
tst = gdaltest.GDALTest('INGR', 'ingr/uint32.cot', 1, 4672)
return tst.testOpen()
def test_ingr_4():
tst = gdaltest.GDALTest('INGR', 'ingr/frmt02.cot', 1, 26968)
return tst.testOpen()
def test_ingr_1():
tst = gdaltest.GDALTest('INGR', 'ingr/8bit_rgb.cot', 2, 4855)
return tst.testOpen()
def test_ingr_16():
tst = gdaltest.GDALTest('INGR', 'ingr/frmt09t.cot', 1, 3178)
return tst.testOpen()
def test_ingr_12():
tst = gdaltest.GDALTest('INGR', 'ingr/frmt28.cot', 2, 45616)
return tst.testOpen()
def isce_4():
tst = gdaltest.GDALTest('isce', 'isce.slc', 1, 350)
return tst.testCreateCopy(check_gt=0, new_filename='isce.tst.slc', vsimem=1)
def test_bt_4():
tst = gdaltest.GDALTest('BT', 'float32.tif', 1, 4672)
return tst.testCreate(out_bands=1)
def test_jdem_1():
tst = gdaltest.GDALTest('JDEM', 'fakejdem.mem', 1, 15)
return tst.testOpen()
def test_bt_5():
tst = gdaltest.GDALTest('BT', 'float32.tif', 1, 4672)
return tst.testSetProjection()
def blx_3():
tst = gdaltest.GDALTest( 'BLX', 's4103.xlb', 1, 47024 )
return tst.testCreateCopy( check_gt = 1, check_srs = 1 )
def test_bt_6():
tst = gdaltest.GDALTest('BT', 'float32.tif', 1, 4672)
return tst.testSetGeoTransform()
def fujibas_1():
tst = gdaltest.GDALTest('FujiBAS', 'fakefujibas.pcb', 1, 1)
return tst.testOpen()
def xyz_1():
tst = gdaltest.GDALTest( 'XYZ', 'byte.tif', 1, 4672 )
return tst.testCreateCopy( vsimem = 1, check_gt = ( -67.00041667, 0.00083333, 0.0, 50.000416667, 0.0, -0.00083333 ) )
gdaltest.post_reason("Did not get expected block size")
return 'fail'
cs = ds.GetRasterBand(1).Checksum()
if cs != 20976:
gdaltest.post_reason('did not get expected checksum')
print(cs)
return 'fail'
ds = None
return 'success'
for item in init_list:
ut = gdaltest.GDALTest('HDF4Image', item[0], item[1], item[2])
if ut is None:
print('HDF4 tests skipped')
sys.exit()
gdaltest_list.append((ut.testOpen, item[0]))
gdaltest_list.append(hdf4_read_online_1)
gdaltest_list.append(hdf4_read_online_2)
gdaltest_list.append(hdf4_read_online_3)
gdaltest_list.append(hdf4_read_online_4)
gdaltest_list.append(hdf4_read_online_5)
gdaltest_list.append(hdf4_read_online_6)
gdaltest_list.append(hdf4_read_online_7)
gdaltest_list.append(hdf4_read_online_8)
gdaltest_list.append(hdf4_read_online_9)
gdaltest_list.append(hdf4_read_online_10)
def xyz_2():
tst = gdaltest.GDALTest( 'XYZ', 'float.img', 1, 23529, options = ['COLUMN_SEPARATOR=,', 'ADD_HEADER_LINE=YES'] )
return tst.testCreateCopy()
def test_gtiff_set(testfunction):
ut = gdaltest.GDALTest('GTiff', 'byte.tif', 1, 4672)
getattr(ut, testfunction)()
def test_mrsid_1():
gdaltest.mrsid_drv = gdal.GetDriverByName('MrSID')
if gdaltest.mrsid_drv is None:
pytest.skip()
tst = gdaltest.GDALTest('MrSID', 'mercator.sid', 1, None)
gt = (-15436.385771224039, 60.0, 0.0, 3321987.8617962394, 0.0, -60.0)
#
# Old, internally generated.
#
prj = """PROJCS["MER E000|",
GEOGCS["NAD27",
DATUM["North_American_Datum_1927",
SPHEROID["Clarke 1866",6378206.4,294.9786982138982,
AUTHORITY["EPSG","7008"]],
AUTHORITY["EPSG","6267"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4267"]],
PROJECTION["Mercator_1SP"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]"""
#
# MrSID SDK getWKT() method.
#
prj = """PROJCS["MER E000|",
GEOGCS["NAD27",
DATUM["North_American_Datum_1927",
SPHEROID["Clarke 1866",6378206.4,294.9786982139006,
AUTHORITY["EPSG","7008"]],
AUTHORITY["EPSG","6267"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4267"]],
PROJECTION["Mercator_1SP"],
PARAMETER["latitude_of_origin",1],
PARAMETER["central_meridian",1],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",1],
PARAMETER["false_northing",1],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]"""
#
# MrSID SDK getWKT() method - DSDK 8 and newer?
#
prj = """PROJCS["MER E000|",
GEOGCS["NAD27",
DATUM["North_American_Datum_1927",
SPHEROID["Clarke 1866",6378206.4,294.9786982139006,
AUTHORITY["EPSG","7008"]],
AUTHORITY["EPSG","6267"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4267"]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]"""
tst.testOpen(check_gt=gt,
check_stat=(0.0, 255.0, 103.319, 55.153),
check_approx_stat=(2.0, 243.0, 103.131, 43.978))
ds = gdal.Open('data/mercator.sid')
got_prj = ds.GetProjectionRef()
ds = None
if prj.find('North_American_Datum_1927') == -1 or \
prj.find('Mercator_1SP') == -1:
print(got_prj)
pytest.fail('did not get expected projection')
if got_prj != prj:
print('Warning: did not get exactly expected projection. Got %s' % got_prj)
def byn_1():
tst = gdaltest.GDALTest('BYN', 'cgg2013ai08_reduced.byn', 1, 64764)
return tst.testOpen()
def jpeg_15():
tst = gdaltest.GDALTest('JPEG', 'albania.jpg', 2, 17016)
return tst.testCreateCopy(vsimem=1, interrupt_during_copy=True)
def gtx_1():
tst = gdaltest.GDALTest( 'GTX', 'hydroc1.gtx', 1, 64183 )
gt = (276.725, 0.05, 0.0, 42.775, 0.0, -0.05)
return tst.testOpen( check_gt = gt, check_prj = 'WGS84' )
def test_ingr_9():
tst = gdaltest.GDALTest('INGR', 'ingr/frmt10.cot', 1, 47031)
return tst.testOpen()
