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 <[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 <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
_______________________________________________
gdal-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
