- edited description
Classification Workflow fails with Raster + Vector Input
Classifcation Workflow cannot be used with the following data. While class names, -colors and strata sizes are found, it throws “Error: raster and reference grids do not match”.
Testdata is attached, raster and vector file are projected in EPSG:3035,, the vector layer needs to get a classified symbol renderer first.
Comments (4)
-
reporter -
reporter - edited description
-
The bug is specific to the given dataset. The raster and the vector have the same spatial reference WKT string, but when burning the points into a raster, GDAL slightly changes the WKT string.
So, the direct string comparison says that the WKT strings do not match.
Also ogr.SpatialReference.IsSame says that the spatial references do not match
As a quick and dirty fix, I simply disabled the consistency check … hope is doesn’t blow up later
# original wkt PROJCS["ETRS89-extended/LAEAEurope",GEOGCS["ETRS89",DATUM["European_Terrestrial_Reference_System_1989",SPHEROID["GRS1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4258"]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["latitude_of_center",52],PARAMETER["longitude_of_center",10],PARAMETER["false_easting",4321000],PARAMETER["false_northing",3210000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Northing",NORTH],AXIS["Easting",EAST],AUTHORITY["EPSG","3035"]] # wkt after burning PROJCS["ETRS89-extended/LAEAEurope",GEOGCS["ETRS89",DATUM["European_Terrestrial_Reference_System_1989",SPHEROID["GRS1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["latitude_of_center",52],PARAMETER["longitude_of_center",10],PARAMETER["false_easting",4321000],PARAMETER["false_northing",3210000],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH]
-
- changed status to resolved
resolves
#511→ <<cset 71de0c8b7f4c>>
- Log in to comment