wgs84 | A zero-dependency Go package for coordinate transformations | Map library

Write the copy and the deletion of the old keys

You could use package ggspatial

Based on dput, you only have one column in your data. It is good practice to examine your data! If you look at str(tracks.utm@data) the only column is TripID and you are attempting to specify the third column.

If you look at the functions help you will see the xy argument that you are passing specifies:

An object inheriting the class SpatialPoints containing the x and y relocations of the animal. If xy inherits the class SpatialPointsDataFrame, it should contain only one column (factor) corresponding to the identity of the animals for each relocation.

This means that TripID should be the only column, be a factor and correspond to the individual animal id(s). Since it is expected that there is only a single column it does not need to be specified, only the sp object. However, apparently the need for a single column does not seem to be the case and you can have multiple columns but need to specify the column containing the unique animal ids, thus the column bracket index.

To ensure that your animal ids are a factor (per help), I would recommend coercing the appropriate column (eg., TripID) to a factor.

The data you are looking to extract is JSON data. It is however not completely valid JSON. You'll have to clean it up a little bit. You can use sites like jsonlint.com to validate JSON data.

The issues with the JSON data are (1) at the beginning, you need to remove the part that says "Highcharts.maps["countries/in/custom/in-all-disputed"] =" up until the first curly bracket {; (2) three lines have a "comment" with two slashes followed by four digits, like this "// 0000". These need to be removed before ingesting the JSON data (I see "// 8440" twice and "// 1227" once).

After doing this you can extract the data from the txt file in R with the jsonlite package.

What version of terra are you using? I ask because I believe that the current version will give a warning but proceed. Also, you could of course aggregate with 100 instead of 113 so that the origins are not changed.

They easiest approach may be to make a virtual raster like this:

v <- vrt(files, "my.vrt")

and then proceed with

r <- rast("my.vrt")

and perhaps

a <- aggregate(r, 100)

Or whatever it is you would like to do.

I have followed a slightly different approach, that is reading and projecting the shapefile with sf, convert it to Spatial* and create an owin. See reprex:

Have a look to the extract function

In your case something like

Since you're not sharing a dataset I will try to explain the problem with mtcars.

The problem you have is that the geom_text function in the following

This is but tricky to debug because you mix different packages, and you do not show(object). Anyway, here is a terra approach, showing where it would be useful to see the objects metadata; and a plot that shows the raster and vector data together.