Classification Workflow fails with Raster + Vector Input

Issue #511 resolved
Benjamin Jakimow created an issue

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)

  1. Andreas Janz

    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]
    

  2. Log in to comment