|
Retested on two small regions, one with some null cells and one with no null cells, using both the GRASS r.covar module and the R layerStats function, and this time results were almost identical. In addition, r.covar outputs "N" the number of cells used, and I saw that in the case of null cells indeed the value of N was smaller. So that answers my previous post :-)
Now I redo the process on the full rasters to see what caused the different in results.
Regards, Micha
On 11/19/18 12:41 PM, Markus Neteler
wrote:
On Sun, Nov 18, 2018 at 10:32 PM Micha Silver <tsvi...@gmail.com> wrote:I am preparing a correlation matrix for 7 raster layers. The results using the r.covar module are different from the R layerStats function. I suspect this is due to handling of null cells. The R function has a parameter to remove NA cells, but the GRASS module, I think, just loops over all cells, including no value.Can anyone confirm that GRASS does not deal with null cells, and that this would cause the difference in correlation results?Here how r.covar treats NULL cells: https://trac.osgeo.org/grass/browser/grass/trunk/raster/r.covar/main.c#L93 Markus -- Micha Silver Ben Gurion Univ. Sde Boker, Remote Sensing Lab cell: +972-523-665918 |
_______________________________________________ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user
