I have tried searching on here for my issues, read many threads, but to no avail.
Essentially, I am trying to convert a line feature to a raster so that I can run proximity on it – because I want to use it in the raster calculator to find areas that are both X elevation and X distance from a road. I have buffered the line (road) by 200m, and want to convert the buffered layer to a raster.
Whichever way I run the converter, I get the same non-existent layer with min/max 1.79769e+ 308 values.
Any suggestions?
v3.4.4
Processing algorithm…
Algorithm 'Rasterize (vector to raster)' starting…
Input parameters:
{ 'BURN' : 1, 'DATA_TYPE' : 5, 'EXTENT' : '-16.95314021128497,-16.0864476950926,27.984934352626247,28.611619710488423 [EPSG:4326]', 'FIELD' : None, 'HEIGHT' : 1000, 'INIT' : None, 'INPUT' : 'C:/Users/Ben/Documents/University/EnvGIS/qgis/buffer_road.gpkg|layername=buffer_road', 'INVERT' : False, 'NODATA' : 0, 'OPTIONS' : '', 'OUTPUT' : 'C:/Users/Ben/Documents/University/EnvGIS/qgis/raster_road_w_buffer.tif', 'UNITS' : 0, 'WIDTH' : 1000 }
GDAL command:
gdal_rasterize -l buffer_road -burn 1.0 -ts 1000.0 1000.0 -a_nodata 0.0 -te -16.95314021128497 27.984934352626247 -16.0864476950926 28.611619710488423 -ot Float32 -of GTiff C:/Users/Ben/Documents/University/EnvGIS/qgis/buffer_road.gpkg C:/Users/Ben/Documents/University/EnvGIS/qgis/raster_road_w_buffer.tif
GDAL command output:
0...10...20...30...40...50...60...70...80...90...100 - done.
Execution completed in 0.30 seconds
Results:
{'OUTPUT': <QgsProcessingOutputLayerDefinition {'sink':C:/Users/Ben/Documents/University/EnvGIS/qgis/raster_road_w_buffer.tif, 'createOptions': {'fileEncoding': 'System'}}>}
Loading resulting layers
Algorithm 'Rasterize (vector to raster)' finished
Raster information log
Processing algorithm…
Algorithm 'Raster information' starting…
Input parameters:
{ 'INPUT' : 'C:/Users/Ben/Documents/University/EnvGIS/qgis/raster_road_w_buffer.tif', 'MIN_MAX' : False, 'NOGCP' : False, 'NO_METADATA' : False, 'OUTPUT' : 'C:/Users/Ben/AppData/Local/Temp/processing_068aa7c8382d44f4974580e307eb6f9f/630ed4178e7946888b29335bca0f80b6/OUTPUT.html', 'STATS' : True }
GDAL command:
gdalinfo -stats C:/Users/Ben/Documents/University/EnvGIS/qgis/raster_road_w_buffer.tif
GDAL command output:
ERROR 1: C:/Users/Ben/Documents/University/EnvGIS/qgis/raster_road_w_buffer.tif, band 1: Failed to compute statistics, no valid pixels found in sampling.
Driver: GTiff/GeoTIFF
Files: C:/Users/Ben/Documents/University/EnvGIS/qgis/raster_road_w_buffer.tif
Size is 1000, 1000
Coordinate System is:
PROJCS["WGS 84 / UTM zone 28N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-15],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AXIS["Easting",EAST],
AXIS["Northing",NORTH],
AUTHORITY["EPSG","32628"]]
Origin = (-16.953140211284971,28.611619710488423)
Pixel Size = (0.000866692516192,-0.000626685357862)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -16.9531402, 28.6116197) ( 19d29'20.02"W, 0d 0' 0.93"N)
Lower Left ( -16.9531402, 27.9849344) ( 19d29'20.02"W, 0d 0' 0.91"N)
Upper Right ( -16.0864477, 28.6116197) ( 19d29'20.00"W, 0d 0' 0.93"N)
Lower Right ( -16.0864477, 27.9849344) ( 19d29'20.00"W, 0d 0' 0.91"N)
Center ( -16.5197940, 28.2982770) ( 19d29'20.01"W, 0d 0' 0.92"N)
Band 1 Block=1000x2 Type=Float32, ColorInterp=Gray
NoData Value=0
Metadata:
STATISTICS_VALID_PERCENT=0
Execution completed in 0.16 seconds
Results:
{'OUTPUT': 'C:/Users/Ben/AppData/Local/Temp/processing_068aa7c8382d44f4974580e307eb6f9f/630ed4178e7946888b29335bca0f80b6/OUTPUT.html'}
Loading resulting layers
Algorithm 'Raster information' finished
HTML output has been generated by this algorithm.
Open the results dialog to check it.
Best Answer
The problem is in the definition of the extents of the output raster.
When using an input layer with a Coordinate Reference System, but defining the extension of the output raster through a layer with a different CRS, the Rasterize tool creates a raster with the CRS of the input layer, but with the coordinates of its extents referred to the other system.
As a result, nothing that was in the input layer is found in those coordinates and the output raster contains only nodata values.
The solution is to use both layers defined under the same CRS (reprojecting some of them), when it is intended to Raterize a layer with the extents defined by another layer.