landcover | Land Cover Mapping | Map library

I have the following problem. I want to build a model for landcover classification. My data are multitemporal Remote Sensing data with several bands. For training I created stratified randomly distributed points to extract spectral data at their positions. With these data a Random Forrest (Rpart) was trained using mlr3 package. For accuracy measurement a repeated spatial cross validation using mlr3spatiotempcv was performed. The resulting model of the training step is, after extraction, stored in an R Object of type rpart. In the terms field of this object are the variable names stored. These are all my used bands but also the spatial x and y coordinates. This brings problems when predicting new data. I used terra package and got an error the x and y layer are missing in my input data. Which kind of makes sense because they are stored in the terms field of the model. But from my understanding, the coordinates should not be a variable of the model. The coordinates are just used for spatial resampling and not for predicting. I "solved" this problem by removing x and y coordinates during the training process and perform just an ordinary non-spatial cross validation. After that I performed the prediction and it works perfectly.

So, my Question is, how can I train a model, using mlr3 package, with data containing coordinates, to perform spatial cross validation?, and then use this model to predict a new Raster.

I am trying to plot only certain values from a categorical land cover raster I am working with. I have loaded it in to R using the terra package and it plots fine. However, since the original data did not come with a legend, I am trying to find out which raster value corresponds to what on the map.

Similar to the answer provided here: How to subset a raster based on grid cell values

I have tried using the following line:

I'm trying to plot the dual facet plot of one bar graph and line graph.

But I can't get the line in my line graph without any error code.

I googled it before, and some people said they could add the aes(group=1) when drawing the line graph. But with the aes(group=1), it can't separate the data into groups as I ordered.

I used this code on other datasets, and they work well. I don't know why they don't work on this data set.

I'm quite new to R.

I try to follow the instructions of this script to draw two scales dummy facet plot in R (https://gist.github.com/sebastianrothbucher/de847063f32fdff02c83b75f59c36a7d).

But I met the problem when I tried to show the group of my data.

I hope I can get the graph like this:

But I only get the plot like this.

I attached my code here.

I am running a common RSF to evaluate landcover type usage by GPS collared deer. I have a factor column named Landcover_Type populated with numbers 1-5 which correspond to a specific landcover type (i.e., 2 = woody, 3 = herbaceous, 4 = water, etc.). How can I exclude water from being considered in my model?

The code below specifies the factor level to use as a reference and the model runs beautifully. However, I need to remove water from consideration.

I'm working with the Circumpolar Arctic Vegetation map. Stored as a SpatRaster with terra, the raster has 21 land cover classes.

I have 16 large datasets of landcover variables around routes. Example dataset "Trial1":

I'm trying to calculate landcover repartition for each US county. I have obtained NLCD for the Apache county using the FedData package (devtools version) and I'm using county shapefiles from the Census bureau.

The problem is that I get an area that is much larger than the official one and the one indicated in my shapefile, namely 51,000km^2 instead of 29,0000km^2 officially. There must be something I don't understand about the raster object but I'm a very confused after hours of websearching, any help appreciated.

The following describes the code used and the method used to calculate. The county data can be downloaded here: https://www2.census.gov/geo/tiger/TIGER2016/COUNTY/

The following code assumes the county shapefile is saved and unzipped.

1. Get and read the data

I am unable to load raster image into NetLogo 6.2.

Following the code: