Skip to content

PCRaster LDD to GRIB error with large coordinates #494

Description

@gregNS

Source data was DEM in EPSG:8083 which has very large X coordinate since the false easting of the transverse mercator projection is 25500000 (twenty-five million, five hundred thousand).

DEM converted to PCRaster map okay.

LDDCREATE produced output with correct coordinates...
Layer properties reported extent:
25590742.2474000006914139,4956697.0966332228854299
25592536.0353999994695187,4957903.9539999999105930

Crayfish | Conversions | PCRaster LDD to GRIB output mesh file has incorrect X coordinates...
Layer properties reported extent:
-21474836.4699999988079071,4956697.4500000001862645 :
-21473043.4699999988079071,4957903.4500000001862645

Alternative testing:

  • coordinate transformation from EPGS:8083 to UTM Zone 20 North with the original DEM
  • convert DEM to PCRaster
  • PCRaster lddcreate
  • Crayfish conversion of PCRaster LDD to GRIB

Coordinate extent was produced as expected.

Version info:

QGIS version | 3.28.1-Firenze
Qt version | 5.15.3
Python version | 3.9.5
GDAL/OGR version | 3.6.0
PROJ version | 9.1.0
EPSG Registry database version | v10.074 (2022-08-01)
GEOS version | 3.11.1-CAPI-1.17.1
SQLite version | 3.39.4
PDAL version | 2.4.3
PostgreSQL client version | unknown
SpatiaLite version | 5.0.1
QWT version | 6.1.6
QScintilla2 version | 2.13.1
OS version | Windows 10 Version 2009

Active Python plugins
AppendFeaturesToLayer | 1.2.0
ClipMultipleLayers | 3.2.0
Coordtransform | 1.1
GarminCustomMap | 3.4.1
generate_swmm_inp | 0.24.1
merge_selected_features | 0.1
pcraster_tools | 0.2.0
qfieldsync | v4.2.0
raster_cutter | 0.6
slyr_community | 4.0.6
valuetool | 3.0.15
wbt_for_qgis | 1.0.5
db_manager | 0.1.20
grassprovider | 2.12.99
processing | 2.12.99
sagaprovider | 2.12.99
RasterAttributeTable | 1.2
crayfish | 3.6.0
MetaSearch | 0.3.6

Test.zip

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions