[gdal-dev] gdal_rasterize parameters not working

[Scott]

gdal_grid is probably a better choice.


On 7/3/23 16:06, Clive Swan wrote:
> 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 <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-
> 
> 
> 
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev