def ogr_georss_10():
try:
os.remove ('tmp/test32631.rss')
except:
pass
srs = osr.SpatialReference()
srs.ImportFromEPSG(32631)
ds = ogr.GetDriverByName('GeoRSS').CreateDataSource('tmp/test32631.rss')
gdal.PushErrorHandler( 'CPLQuietErrorHandler' )
try:
lyr = ds.CreateLayer('georss', srs = srs)
except:
lyr = None
gdal.PopErrorHandler()
if lyr is not None:
gdal.post_reason('should not have accepted EPSG:32631 with GEOM_DIALECT != GML')
return 'fail'
ds.Destroy()
try:
os.remove ('tmp/test32631.rss')
except:
pass
ds = ogr.GetDriverByName('GeoRSS').CreateDataSource('tmp/test32631.rss', options = [ 'GEOM_DIALECT=GML' ])
lyr = ds.CreateLayer('georss', srs = srs)
dst_feat = ogr.Feature( feature_def = lyr.GetLayerDefn() )
dst_feat.SetGeometry(ogr.CreateGeometryFromWkt('POINT (500000 4000000)'))
if lyr.CreateFeature( dst_feat ) != 0:
gdaltest.post_reason('CreateFeature failed.')
return 'fail'
dst_feat.Destroy()
ds.Destroy()
content = open('tmp/test32631.rss').read()
if content.find('<georss:where><gml:Point srsName="urn:ogc:def:crs:EPSG::32631"><gml:pos>500000 4000000') == -1:
print(('%s' % content))
return 'fail'
return 'success'
def pcraster_2():
if gdaltest.pcraster_drv is None:
return 'skip'
ds = gdal.Open( 'data/ldd.map' )
gt = ds.GetGeoTransform()
if gt[0] != 182140.0 or gt[1] != 10 or gt[2] != 0 \
or gt[3] != 327880.0 or gt[4] != 0 or gt[5] != -10:
gdal.post_reason( 'PCRaster geotransform wrong.' )
return 'fail'
band1 = ds.GetRasterBand(1)
if band1.GetNoDataValue() != 255:
gdal.post_reason( 'PCRaster NODATA value wrong or missing.' )
return 'fail'
return 'success'
def pcraster_2():
if gdaltest.pcraster_drv is None:
return 'skip'
ds = gdal.Open( 'data/ldd.map' )
gt = ds.GetGeoTransform()
if gt[0] != 182140.0 or gt[1] != 10 or gt[2] != 0 \
or gt[3] != 327880.0 or gt[4] != 0 or gt[5] != -10:
gdal.post_reason( 'PCRaster geotransform wrong.' )
return 'fail'
band1 = ds.GetRasterBand(1)
if band1.GetNoDataValue() != 255:
gdal.post_reason( 'PCRaster NODATA value wrong or missing.' )
return 'fail'
return 'success'
def lcp_2():
ds = gdal.Open("data/test_USGS_LFNM_Alb83.lcp")
if ds.RasterCount != 8:
gdal.post_reason("wrong number of bands")
return "fail"
metadata = [("LATITUDE", "48"), ("LINEAR_UNIT", "Meters"), ("DESCRIPTION", "")]
md = ds.GetMetadata()
for item in metadata:
if md[item[0]] != item[1]:
gdaltest.post_reason(
"wrong metadataitem for dataset. md['%s']='%s', expected '%s'" % (item[0], md[item[0]], item[1])
)
return "fail"
check_gt = (-1328145, 30, 0, 2961735, 0, -30)
new_gt = ds.GetGeoTransform()
for i in range(6):
if abs(new_gt[i] - check_gt[i]) > 1e-5:
print("")
print("old = ", check_gt)
print("new = ", new_gt)
gdaltest.post_reason("Geotransform differs.")
return "fail"
dataPerBand = [
(
28381,
[
("ELEVATION_UNIT", "0"),
("ELEVATION_UNIT_NAME", "Meters"),
("ELEVATION_MIN", "1064"),
("ELEVATION_MAX", "1492"),
("ELEVATION_NUM_CLASSES", "-1"),
("ELEVATION_FILE", "d:\\scratch\\dist\\79990093\\Output\\rastert_elevation_1.txt"),
],
),
(
25824,
[
("SLOPE_UNIT", "0"),
("SLOPE_UNIT_NAME", "Degrees"),
("SLOPE_MIN", "0"),
("SLOPE_MAX", "34"),
("SLOPE_NUM_CLASSES", "35"),
("SLOPE_FILE", "d:\\scratch\\dist\\79990093\\Output\\rastert_slope_1.txt"),
],
),
(
28413,
[
("ASPECT_UNIT", "2"),
("ASPECT_UNIT_NAME", "Azimuth degrees"),
("ASPECT_MIN", "0"),
("ASPECT_MAX", "357"),
("ASPECT_NUM_CLASSES", "-1"),
("ASPECT_FILE", "d:\\scratch\\dist\\79990093\\Output\\rastert_aspect_1.txt"),
],
),
(
19052,
[
("FUEL_MODEL_OPTION", "0"),
("FUEL_MODEL_OPTION_DESC", "no custom models AND no conversion file needed"),
("FUEL_MODEL_MIN", "1"),
("FUEL_MODEL_MAX", "10"),
("FUEL_MODEL_NUM_CLASSES", "5"),
("FUEL_MODEL_VALUES", "1,2,5,8,10"),
("FUEL_MODEL_FILE", "d:\\scratch\\dist\\79990093\\Output\\rastert_fuel1.txt"),
],
),
(
30164,
[
("CANOPY_COV_UNIT", "1"),
("CANOPY_COV_UNIT_NAME", "Percent"),
("CANOPY_COV_MIN", "0"),
("CANOPY_COV_MAX", "95"),
("CANOPY_COV_NUM_CLASSES", "10"),
("CANOPY_COV_FILE", "d:\\scratch\\dist\\79990093\\Output\\rastert_canopy1.txt"),
],
),
(
22316,
[
("CANOPY_HT_UNIT", "3"),
("CANOPY_HT_UNIT_NAME", "Meters x 10"),
("CANOPY_HT_MIN", "0"),
("CANOPY_HT_MAX", "375"),
("CANOPY_HT_NUM_CLASSES", "5"),
("CANOPY_HT_FILE", "d:\\scratch\\dist\\79990093\\Output\\rastert_height_1.txt"),
],
),
(
30575,
[
("CBH_UNIT", "3"),
("CBH_UNIT_NAME", "Meters x 10"),
("CBH_MIN", "0"),
("CBH_MAX", "100"),
("CBH_NUM_CLASSES", "33"),
("CBH_FILE", "d:\\scratch\\dist\\79990093\\Output\\rastert_base_1.txt"),
],
),
(
23304,
[
("CBD_UNIT", "3"),
("CBD_UNIT_NAME", "kg/m^3 x 100"),
("CBD_MIN", "0"),
("CBD_MAX", "21"),
("CBD_NUM_CLASSES", "20"),
("CBD_FILE", "d:\\scratch\\dist\\79990093\\Output\\rastert_density_1.txt"),
],
),
]
for i in range(8):
band = ds.GetRasterBand(i + 1)
if band.Checksum() != dataPerBand[i][0]:
gdaltest.post_reason(
"wrong checksum for band %d. Got %d, expected %d" % (i + 1, band.Checksum(), dataPerBand[i][0])
)
return "fail"
md = band.GetMetadata()
for item in dataPerBand[i][1]:
if md[item[0]] != item[1]:
gdaltest.post_reason(
"wrong metadataitem for band %d. md['%s']='%s', expected '%s'"
% (i + 1, item[0], md[item[0]], item[1])
)
return "fail"
ds = None
return "success"
def lcp_2():
ds = gdal.Open('data/test_USGS_LFNM_Alb83.lcp')
if ds.RasterCount != 8:
gdal.post_reason('wrong number of bands')
return 'fail'
metadata = [('LATITUDE', '48'), ('LINEAR_UNIT', 'Meters'),
('DESCRIPTION', '')]
md = ds.GetMetadata()
for item in metadata:
if md[item[0]] != item[1]:
gdaltest.post_reason(
'wrong metadataitem for dataset. md[\'%s\']=\'%s\', expected \'%s\''
% (item[0], md[item[0]], item[1]))
return 'fail'
check_gt = (-1328145, 30, 0, 2961735, 0, -30)
new_gt = ds.GetGeoTransform()
for i in range(6):
if abs(new_gt[i] - check_gt[i]) > 1e-5:
print('')
print('old = ', check_gt)
print('new = ', new_gt)
gdaltest.post_reason('Geotransform differs.')
return 'fail'
dataPerBand = [
(28381,
[('ELEVATION_UNIT', '0'), ('ELEVATION_UNIT_NAME', 'Meters'),
('ELEVATION_MIN', '1064'), ('ELEVATION_MAX', '1492'),
('ELEVATION_NUM_CLASSES', '-1'),
('ELEVATION_FILE',
'd:\\scratch\\dist\\79990093\\Output\\rastert_elevation_1.txt')]),
(25824, [('SLOPE_UNIT', '0'), ('SLOPE_UNIT_NAME', 'Degrees'),
('SLOPE_MIN', '0'), ('SLOPE_MAX', '34'),
('SLOPE_NUM_CLASSES', '35'),
('SLOPE_FILE',
'd:\\scratch\\dist\\79990093\\Output\\rastert_slope_1.txt')
]),
(28413, [('ASPECT_UNIT', '2'), ('ASPECT_UNIT_NAME', 'Azimuth degrees'),
('ASPECT_MIN', '0'), ('ASPECT_MAX', '357'),
('ASPECT_NUM_CLASSES', '-1'),
('ASPECT_FILE',
'd:\\scratch\\dist\\79990093\\Output\\rastert_aspect_1.txt')
]),
(19052, [('FUEL_MODEL_OPTION', '0'),
('FUEL_MODEL_OPTION_DESC',
'no custom models AND no conversion file needed'),
('FUEL_MODEL_MIN', '1'), ('FUEL_MODEL_MAX', '10'),
('FUEL_MODEL_NUM_CLASSES', '5'),
('FUEL_MODEL_VALUES', '1,2,5,8,10'),
('FUEL_MODEL_FILE',
'd:\\scratch\\dist\\79990093\\Output\\rastert_fuel1.txt')]),
(30164, [('CANOPY_COV_UNIT', '1'), ('CANOPY_COV_UNIT_NAME', 'Percent'),
('CANOPY_COV_MIN', '0'), ('CANOPY_COV_MAX', '95'),
('CANOPY_COV_NUM_CLASSES', '10'),
('CANOPY_COV_FILE',
'd:\\scratch\\dist\\79990093\\Output\\rastert_canopy1.txt')
]),
(22316,
[('CANOPY_HT_UNIT', '3'), ('CANOPY_HT_UNIT_NAME', 'Meters x 10'),
('CANOPY_HT_MIN', '0'), ('CANOPY_HT_MAX', '375'),
('CANOPY_HT_NUM_CLASSES', '5'),
('CANOPY_HT_FILE',
'd:\\scratch\\dist\\79990093\\Output\\rastert_height_1.txt')]),
(30575, [('CBH_UNIT', '3'), ('CBH_UNIT_NAME', 'Meters x 10'),
('CBH_MIN', '0'), ('CBH_MAX', '100'),
('CBH_NUM_CLASSES', '33'),
('CBH_FILE',
'd:\\scratch\\dist\\79990093\\Output\\rastert_base_1.txt')]),
(23304, [
('CBD_UNIT', '3'), ('CBD_UNIT_NAME', 'kg/m^3 x 100'),
('CBD_MIN', '0'), ('CBD_MAX', '21'), ('CBD_NUM_CLASSES', '20'),
('CBD_FILE',
'd:\\scratch\\dist\\79990093\\Output\\rastert_density_1.txt')
])
]
for i in range(8):
band = ds.GetRasterBand(i + 1)
if band.Checksum() != dataPerBand[i][0]:
gdaltest.post_reason(
'wrong checksum for band %d. Got %d, expected %d' %
(i + 1, band.Checksum(), dataPerBand[i][0]))
return 'fail'
md = band.GetMetadata()
for item in dataPerBand[i][1]:
if md[item[0]] != item[1]:
gdaltest.post_reason(
'wrong metadataitem for band %d. md[\'%s\']=\'%s\', expected \'%s\''
% (i + 1, item[0], md[item[0]], item[1]))
return 'fail'
ds = None
return 'success'
