Dear community,

I am exploring the functionality of raster package by following the online 
textbook Geocomputation with R, 1st ed. (I know the newer version of the 
textbook is online, I just want to finish what I left). The problem I 
encountered is associated with the example from Section 5.4.2, raster 
extraction:
https://bookdown.org/robinlovelace/geocompr/geometric-operations.html#raster-extraction
When I want to combine it with st_sample() from Section 5.2.3, i.e. extract the 
pixel values from say 15 randomly chosen locations, I met an error.
Here comes the script:

library(sf)
library(raster)

filename=system.file('raster/srtm.tif', package='spDataLarge')
srtm=raster(filename)
filename=system.file('vector/zion.gpkg', package='spDataLarge')
zion=st_read(filename)
zion=st_transform(zion, projection(srtm))

data('zion_points', package='spDataLarge')
extract(srtm, zion_points) # textbook example, work as expected

sampled_sites=st_sample(zion, 15)
extract(srtm, sampled_sites) # error occurred.

The error message says "����: ����extract����ǩ��x = "RasterLayer", y = 
"sfc_POINT"���Ҳ����̳з���"
There is Chinese text embedded because my OS is in Chinese. Translated into 
English, it reads "ERROR: function 'extract' labelling ��x = "RasterLayer", y = 
"sfc_POINT"�� cannot find the inherited methods".

My question: what's the reason of this problem, and how to fix it?

Thank you in advance!

Ҷ�� YE, Xiang
THINKING SPATIALLY<http://www.linkedin.com/in/spatialyexiang>.
Ph.D. in Spatial Statistics

        [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to