[gdal-dev] gdal_rasterize parameters not working

[Clive Swan]

Hi

Thanks for the information.


 The Shp file (WGS84), the aim is to export to a Heatmap.

The Shp was exported to a Tif:

coastal_2000-te.tif
coastal_2000-tr.tif

The pixel output ranges from 0-3

This doesn't produce a Heatmap.


Any options you would recommend??


Thanks

Clive



gdalinfo "C:\Temp\PROJECTS\CCRT\ClimateChangeFME\Modelling\UK\Maps\coastal_2000-tr.tif"

>> 1Kb

Driver: GTiff/GeoTIFF

Files: C:\Temp\PROJECTS\CCRT\ClimateChangeFME\Modelling\UK\Maps\coastal_2000-tr.tif

Size is 1, 1

Coordinate System is:

GEOGCRS["WGS 84", ENSEMBLE["World Geodetic System 1984 ensemble",

        MEMBER["World Geodetic System 1984 (Transit)"], MEMBER["World Geodetic System 1984 (G730)"], MEMBER["World Geodetic System 1984 (G873)"], MEMBER["World Geodetic System 1984 (G1150)"],

        MEMBER["World Geodetic System 1984 (G1674)"], MEMBER["World Geodetic System 1984 (G1762)"], MEMBER["World Geodetic System 1984 (G2139)"], ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]], ENSEMBLEACCURACY[2.0]],

    PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], CS[ellipsoidal,2], AXIS["geodetic latitude (Lat)",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433]], AXIS["geodetic longitude (Lon)",east, ORDER[2], ANGLEUNIT["degree",0.0174532925199433]],

    USAGE[ SCOPE["Horizontal component of 3D system."], AREA["World."], BBOX[-90,-180,90,180]], ID["EPSG",4326]]

Data axis to CRS axis mapping: 2,1

Origin = (-8.263930151999944,60.841669081000077)

Pixel Size = (10.000000000000000,-10.000000000000000)

Metadata: AREA_OR_POINT=Area

Image Structure Metadata: INTERLEAVE=BAND

Corner Coordinates:

Upper Left  (  -8.2639302,  60.8416691) (  8d15'50.15"W, 60d50'30.01"N)

Lower Left  (  -8.2639302,  50.8416691) (  8d15'50.15"W, 50d50'30.01"N)

Upper Right (   1.7360698,  60.8416691) (  1d44' 9.85"E, 60d50'30.01"N)

Lower Right (   1.7360698,  50.8416691) (  1d44' 9.85"E, 50d50'30.01"N)

Center      (  -3.2639302,  55.8416691) (  3d15'50.15"W, 55d50'30.01"N)

Band 1 Block=1x1 Type=Float64, ColorInterp=Gray





________________________________
From: Rahkonen Jukka <jukka.rahkonen at maanmittauslaitos.fi>
Sent: Monday, July 3, 2023 4:22:41 pm
To: Clive Swan <cliveswan at gmail.com>; gdal-dev at lists.osgeo.org <gdal-dev at lists.osgeo.org>
Subject: Re: gdal_rasterize parameters not working

Hi,

What coordinate system and extent ogrinfo reports with
ogrinfo -al -so "C:\Temp\PROJECTS\CCRT\ClimateChangeFME\Modelling\UK\Shp\Coastal\Coastal_Flooding_2000.shp"
What coordinate system and extent gdalinfo reports with
gdalinfo "C:\Temp\PROJECTS\CCRT\ClimateChangeFME\Modelling\UK\Maps\coastal_2000.tif"

Did the coastal_2000.tif exist already or was it created by gdal_rasterize? Have you been thinking that the default datatype in the result is Float64 https://gdal.org/programs/gdal_rasterize.html#cmdoption-gdal_rasterize-ot and that the whole range  52 to 82 in the source data is very, very close to the black end of the target TIFF? Without LUT stretch the result looks certainly all black in image viewers.

-Jukka Rahkonen-

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20230703/9b5fcd53/attachment.htm>