export raster with replaced values

[wxinl]

I used the rioxarray to read a dem tiff, then I tried to replace the raster values with a numpy array of the same shape and output the raster with replaced values using rio.to_raster. But when I read the new file, the values are all 0s. Image attached.

[snowman2]

@wxinl if you wouldn't mind making these changes, it would better enable us to help you:

• the input raster
• the code as text instead of an image for re-reproducibility

Also, what were the values contained by dem_rx right before writing to a raster?

UPDATE:

dem_rx.encoding

{'dtype': 'int16',
 '_FillValue': -32768.0,
 'grid_mapping': 'spatial_ref',
 'source': './ayt08_nasadem.tif',
 'rasterio_dtype': 'int16'}

Since you have a different dtype, you need to do:

dem_rx.encoding = {}

Then it should work as expected,

[wxinl]

slope.zip
Hi, thank you. The dem_rx contains the elevation of each pixel. A sample file and code are provided.

[snowman2]

Can you also provide what the values contained by dem_rx were right before writing to a raster?

[wxinl]

It's the slope calculated in the code.

[snowman2]

Updated the answer.