R Terra Package – Equivalent Function to isLonLat for Coordinate Checks

rastersfstarsterra

I am refactoring my code which extensively relied on raster package to now use the newer and better terra package. I was making use of the isLonLat() function that returns TRUE if the crs of a raster object is angular and FALSE if it's planar.

Is there a way to do this using terra instead?

I want to make my code future-proof and the reliance of this function on the now deprecated proj4string is not acceptable.

Best Answer

Well, you have terra::is.lonlat() https://www.rdocumentation.org/packages/terra/versions/1.5-17/topics/is.lonlat

Related Solutions

[GIS] Using spatial points to extract values from large raster in R

What happens if you run data.df <- sfSapply(list(veg.r), extract, y=data.sp) ? The question you link to is set up to extract data from a list of rasters, not just one.

Otherwise, how's performance if you rasterize your points layer and use it to extract values? Pseudocode:

presabs <- rasterize(points, vegraster, field = 'ID')
presabs[!is.na(values(presabs))] <- 0
xtrct <- presabs + vegraster
vals <- rasterToPoints(xtrct)

Then you could use this spatial join method to append the extracted data back to your original points. Suspect it might not be foolproof (multiple points in the same cell might cause issues) but maybe worth a shot.

Python Lidar – Handling NoData Values in Lidar Packages

Looks like that package was built with an old version of numpy. The newer versions did away with the np.float function.

You could downgrade your numpy version to 1.20 (you will get a deprecation warning rather than the error) or you could try changing the source code to use the built in Python float function such as: max_elev = float(np.max(dem[dem != no_data]))

See this post:

https://stackoverflow.com/questions/74844262/how-to-solve-error-numpy-has-no-attribute-float-in-python

A side note: I prefer to use PDAL for this type of lidar manipulation. There is a learning curve but it is a very powerful suite of tools. https://pypi.org/project/pdal/

Best Answer

Related Solutions

[GIS] Using spatial points to extract values from large raster in R

Python Lidar – Handling NoData Values in Lidar Packages

Related Question