Error message that occurs processing hazard map

Community
1

Hi
Joining an exposure dataset and some hazard maps, I am getting the error message

Failed to validate model for execution
  - Failed to validate `join` step (exposures_join_hazards) on line 25 for execution
    - Problems found with 'bookmark' riskscape function
      - Problems found with 'mapfilepath.tif' bookmark
        - Failed to open 'mapfilepath.tif' - GEOTIFF Module Error Report
No code "EPSG:2193" from authority "EPSG" found for object of type "interface org.geotools.api.referencing.crs.EngineeringCRS".
ModelPixelScaleTag: [8.0,8.0,0.0]
ModelTiePointTag: (1 tie points)
TP #0: [0.0,0.0,0.0] -> [1778060.0,5582308.0,0.0]
ModelTransformationTag: NOT AVAILABLE
GeoKey #1: Key = 1024, Value = 1
GeoKey #2: Key = 3072, Value = 2193
GeoKey #3: Key = 1025, Value = 1
GeoKey #4: Key = 2049, Value = NZGD2000
GeoKey #5: Key = 1026, Value = NZGD2000 / New Zealand Transverse Mercator 2000
GeoKey #6: Key = 3076, Value = 9001
GeoKey #7: Key = 2054, Value = 9102

        - org.geotools.api.data.DataSourceException: GEOTIFF Module Error Report
No code "EPSG:2193" from authority "EPSG" found for object of type "interface org.geotools.api.referencing.crs.EngineeringCRS".
ModelPixelScaleTag: [8.0,8.0,0.0]
ModelTiePointTag: (1 tie points)
TP #0: [0.0,0.0,0.0] -> [1778060.0,5582308.0,0.0]
ModelTransformationTag: NOT AVAILABLE
GeoKey #1: Key = 1024, Value = 1
GeoKey #2: Key = 3072, Value = 2193
GeoKey #3: Key = 1025, Value = 1
GeoKey #4: Key = 2049, Value = NZGD2000
GeoKey #5: Key = 1026, Value = NZGD2000 / New Zealand Transverse Mercator 2000
GeoKey #6: Key = 3076, Value = 9001
GeoKey #7: Key = 2054, Value = 9102

The hazard map runs without problems with a different exposure dataset and it was generated as part of a map series that all work fine.

Do you have an idea what this means and if this might have to do with the exposure dataset? (That one also runs without problem with other hazard maps of the series.)
In future, it would be nice if RiskScape could add at which exposure feature an error occur.

I have tried explicitly adding crs-name = EPSG:2193 to the hazard map bookmark but this does not seem to make a difference.
Looking forward to your ideas.

Thanks,
Luisa

2

Hi Luisa,

This looks quite odd, but it looks as though there’s something wrong with the projection that’s defined for that particular tiff - mapfilepath.tif.

Do you have gdal installed? If you do, can you run gdalinfo mapfilepath.tif and paste the output in here, and then do the same for one of the rasters that isn’t being a problem. This might show that there’s a difference (and be a bit of a “smoking gun”).

Assuming you do have gdal installed, you can try and fix this using gdal_edit.py to overwrite the projection that’s defined - you might want to make a copy of the original first.

Something like:

  • mv mapfilepath.tif backup.tif
  • gdal_edit.py -a_srs EPSG:2193 mapfilepath.tif

Try your model again and see if that improves things.

Cheers,
Nick

3

Hi Nick,
Comparing with a map that works we get this output

diff <(gdalinfo mapfilepath.tif) <(gdalinfo workingmap.tif)
2,3c2,3
< Files: mapfilepath.tif
< Size is 10496, 10752
---
> Files: workingmap.tif
> Size is 7296, 11392
44c44
< Origin = (1778060.000000000000000,5582308.000000000000000)
---
> Origin = (1749916.000000000000000,5649612.000000000000000)
51,58c51,58
< Upper Left  ( 1778060.000, 5582308.000) (175d 4'57.34"E, 39d53'26.20"S)
< Lower Left  ( 1778060.000, 5496292.000) (175d 6'23.41"E, 40d39'54.15"S)
< Upper Right ( 1862028.000, 5582308.000) (176d 3'48.87"E, 39d52' 7.80"S)
< Lower Right ( 1862028.000, 5496292.000) (176d 5'55.35"E, 40d38'33.58"S)
< Center      ( 1820044.000, 5539300.000) (175d35'16.22"E, 40d16' 4.27"S)
< Band 1 Block=10496x1 Type=Float32, ColorInterp=Gray
<   Min=0.030 Max=50.975
<   Minimum=0.030, Maximum=50.975, Mean=15.310, StdDev=18.762
---
> Upper Left  ( 1749916.000, 5649612.000) (174d44'18.09"E, 39d17'23.63"S)
> Lower Left  ( 1749916.000, 5558476.000) (174d45'32.69"E, 40d 6'38.43"S)
> Upper Right ( 1808284.000, 5649612.000) (175d24'52.87"E, 39d16'40.19"S)
> Lower Right ( 1808284.000, 5558476.000) (175d26'36.43"E, 40d 5'53.71"S)
> Center      ( 1779100.000, 5604044.000) (175d 5'19.75"E, 39d41'40.84"S)
> Band 1 Block=7296x1 Type=Float32, ColorInterp=Gray
>   Min=0.030 Max=93.208
>   Minimum=0.030, Maximum=93.208, Mean=22.696, StdDev=15.607
61,65c61,65
<     STATISTICS_MAXIMUM=50.975048065186
<     STATISTICS_MEAN=15.309991057795
<     STATISTICS_MINIMUM=0.030000040307641
<     STATISTICS_STDDEV=18.762306021418
<     STATISTICS_VALID_PERCENT=17.25
---
>     STATISTICS_MAXIMUM=93.208351135254
>     STATISTICS_MEAN=22.696151747388
>     STATISTICS_MINIMUM=0.030000016093254
>     STATISTICS_STDDEV=15.60651094595
>     STATISTICS_VALID_PERCENT=15.4

Both have the projcrs

Coordinate System is:
PROJCRS["NZGD2000 / New Zealand Transverse Mercator 2000",
    BASEGEOGCRS["NZGD2000",
        DATUM["New Zealand Geodetic Datum 2000",
            ELLIPSOID["GRS 1980",6378137,298.257222101,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4167]],
    CONVERSION["New Zealand Transverse Mercator 2000",
        METHOD["Transverse Mercator",
            ID["EPSG",9807]],
        PARAMETER["Latitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",173,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["Scale factor at natural origin",0.9996,
            SCALEUNIT["unity",1],
            ID["EPSG",8805]],
        PARAMETER["False easting",1600000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",10000000,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["northing (N)",north,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["easting (E)",east,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    USAGE[
        SCOPE["Engineering survey, topographic mapping."],
        AREA["New Zealand - North Island, South Island, Stewart Island - onshore."],
        BBOX[-47.33,166.37,-34.1,178.63]],
    ID["EPSG",2193]]