Why is terra's focal function changing raster location?
I am using the focal function in the terra package to fill holes in a DEM raster. However, after running focal(m, w=19, fun=mean, na.policy="only", na.rm=TRUE) my DEM is moved ~11cm relative to it's original position (see images below). This problem does not occur when I simply read in the raster using rast() and then write it using writeRaster(), but it happens when I use focal() regardless of whether I save the output using writeRaster() or using the filename = argument within focal().
I have also noticed that the larger the window size used, the greater the movement of the output raster. A window size of 9 produced ~5cm error. Raster pixels are ~1.2cm in size so focal() appears to be causing a movement of ~1/2 the width of the window used. As far as I can tell the movement is straight up/north.
Any idea why this might be occuring and how it could be avoided?
This the correctly aligned raster:
And this is after running focal (w = 19):
That is a bug that, I believe, has been fixed in the development version of terra available at https://github.com/rspatial/terra (and that is also the best place to report bugs). Can you please try it?
- Rcpp Rf_warningcall compiler warnings
- Modify the name of factor variables in lm function(summary function)
- Extreme value analysis and quantile estimation using log Pearson type 3 (Pearson III) distribution - R vs Python
- How to hide NAs when using xlsx::saveWorkbook?
- How do I retrieve a simple numeric value from a named numeric vector in R?
- Matching pair-wise columns from left to right across rows in one dataframe to another dataframe and adding new columns with matching values
- Income to outcome flow chart in Sankey plotly R
- color mapping in geom_conn_bundle not showing correctly
- Print R package startup message AFTER automatic package conflict messages instead of before
- Summing a set of R dataframe rows (column-wise), while retaining the first n columns
- Added variable / partial regression plots for groups in an interaction?
- how to make a topoplot in R with coordinates variable distribution
- List of all functions in base R?
- Plotting multiple plots for different initial conditions in one graph
- Printing repetitively on the same line in R
- Generating UI/Server based on initial selection
- Subset dataframe based on pickerInput
- How to let user pick the data in R-shiny?
- Couldn't show my simple bar charts separately on Shiny R dashboardBody
- How to programmatically filter contents of a second shiny app displayed via iframe
- How to select specific interesting groups for the boxplot in R Shiny app?
- Crosstable and Plot grouping with reactive values
- Is there a way to make multiple Shiny picker inputs where the selections must be disjoint?
- Delay/avoid duplication of shiny server side functions until after credentials
- Predictions only returns value "1"
- How to display a busy indicator in a shiny app?
- Append doesn't work when writing to CSV in R
- Changing the start date of a gantt chart in DiagrammeR
- Check for installed packages before running install.packages()
- Compare two columns element-wise