I am trying to analyze how sets of rasters relate to each other using some statistical techniques. As, I don't have much experience using the spatial statistics tools in ArcGIS I was exporting my rasters as Ascii files, and analyzing them using R (specifically the maptools
package, and readAsciiGrid()
). This has been functioning ok (but as the datasets have 90,000 points it is slow to run the analysis), but I don't know if I am recreating in R, existing functionality in ArcGIS.
For example, I want to perform regressions between each of these rasters using a few different transformations (logarithmic, exponential, etc). Can this be done within ArcGIS? A second broader question is if there are standard statistical methods for examining this type of data?
Each raster pair has matching data/no-data values and all parameters are identical, aside from the gridcell value.
Best Answer
I would stick to R. If speed is really a problem ( I doubt so 90.000 is not such a big number) you could try finding relationships between a subset of your data. Actually the first thing I would do is make a plot to look for obvious relationships.
Even if arcgis contains tools to compare rasters, R will always give you a lot more statistical tools.
Eg:
I should actually add that often it is more correct to work with a subset of your data then with your full dataset. In many cases grid cells are not independent from the surrounding data cells, which will result in overly optimistic p values for eg regression fits (you will find more info if you search on declustering).