<div dir="ltr">Hello everyone, <div><br></div><div>I am trying to discern how I can accomplish the following at the command line with GDAL and could use some guidance. I have a GeoTIFF containing data from the National Land Cover Database and I would like to add another raster layer to the GeoTIFF that contains land management labels for the centroids of each pixel in the NLCD raster. I have a GeoJSON file that has polygons representing areas with a particular land management label. At the end of the day, I'd like to produce a GeoTIFF with two bands, one representing the land cover type and another representing the land management type. </div><div><br></div><div>I located <a href="https://gis.stackexchange.com/questions/365258/spatial-join-between-raster-and-shapefille-give-raster-the-shapefile-values-w">this example</a> that I believe is a Python-based solution to my problem but the code is producing an empty GeoTIFF. I assume if a solution path does exist, there's a way to do this at the command line which would be preferable. </div><div><br></div><div>If someone can point me in the right direction, I would be very grateful! I'm struggling to see how to get there from the GDAL docs. </div><div><br></div><div>Chris</div></div>