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 <[email protected]>
Sent: Monday, July 3, 2023 4:22:41 pm
To: Clive Swan <[email protected]>; [email protected] 
<[email protected]>
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-

_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to