This is what combined raster table looks like:
So Value of 1 corresponds to dry land before and after; 2 - was dry become flooded etc.
Alternatively one can use fact that sum of any two numbers in range (0,1,2,4,8..) gives unique combination. Those who used paper tapes for coding will understand this real quick:) So you can do:
You can us gdal_calc.py in a for loop. First read all the files name, then use:
gdal_calc.py -A input1.tif -B input2.tif --outfile=result.tif --calc="A+B"
And then in a loop:
gdal_calc.py -A inputX.tif -B result.tif --outfile=result.tif --calc="A+B"
Where inputX.tifis the current file you are looping over, so your loop should start at input3 since the ...
There is no direct way to calculate the minimum and maximum from the Raster Calculator in QGIS. The only workaround is to extract the minimum and maximum of the raster layer using Raster layer statistics tool from Processing Toolbox, and save it into a text editor:
Then use the above minimum and maximum values in the raster calculator from Processing ...
The "error" was that the min max value was just "estimate" and not actual, and as the max value of original raster was just above 0.1 (0.010298) the setting of symbology of the layer (it was estimated) was defined in way that it was not showing the last cathegory (the 4th one) thats why I was missing all the time the last cathegory.
I think you should use set the extent of the output in the Environments of the Raster Calculator to match your larger raster (Raster 1).
(here, you should add your Raster2 extent, otherwise use other options like Union of inputs that might work as well)
Then do a Conditioal expression like (not tested):
This should output a ...
Select several control points on your DEM to create the surface. Create a new point layer and digitize these points.
Obtain elevation data from the DEM at control points, by Add raster values to points tool in Processing Toolbox | SAGA | Vector <> raster
Start Polynomial regression tool in Processing Toolbox | SAGA | Geostatistics. Select ...
ArcGIS Help description of the error is the following:
This is an operating system error and can indicate different error
conditions. These errors are not specifically documented in the ArcGIS
Desktop Help system. However, in most cases, the error description
provided with the error code can be used to look up additional
information on the ...
cvr_rasters will be a list of only raster names with no path/workspace. When you change workspace after it is created the rasters will not be found. Try adding path to each raster:
arcpy.env.workspace = CVR
cvr_rasters = arcpy.ListRasters()
cvr_rasters = [os.path.join(CVR, r) for r in cvr_rasters]
Same goes for lss_rasters since you change ...
The arcpy.ListRasters() function will give you a list with the filenames of all the rasters in a given workspace. This way, you can specify the first GDB as your workspace and get a list of all the rasters in it. Then, you can repeat this process for the second GDB.
CVR = r"C:\Aquifer_recalc\Aquifer_recalc_NA\Treecover_NA\Treecover_Inputs.gdb"
LSS = r"C:\...
"Reclassify by table" in the Processing toolbox did the job.
There, you just need to define the range values and click the box "use no-data when no range matches value". Then you will get your mask layer.
Thanks to GBG for suggesting the post where I found the solution.