读取shp文件
R包：sf
数据来源：.shp
library(sf)
#st_read()函数可以读取.shp文件
data<-st_read("Income_schooling.shp")
注意：文件路径上不要出现汉字和空格，否则读出的空间数据会出现乱码
查看空间数据data的前几条记录
head(data,n=4)
# Simple feature collection with 4 features and 5 fields
# Geometry type: MULTIPOLYGON
# Dimension: XY
# Bounding box: xmin: 379071.8 ymin: 4936182 xmax: 596500.1 ymax: 5255569
# Projected CRS: NAD83 / UTM zone 19N
# NAME Income NoSchool NoSchoolSE IncomeSE geometry
# 1 Aroostook 21024 0.01338720 0.00140696 250.909 MULTIPOLYGON (((513821.1 51...
# 2 Somerset 21025 0.00521153 0.00115002 390.909 MULTIPOLYGON (((379071.8 50...
# 3 Piscataquis 21292 0.00633830 0.00212896 724.242 MULTIPOLYGON (((445039.5 51...
# 4 Penobscot 23307 0.00684534 0.00102545 242.424 MULTIPOLYGON (((472271.3 49...
geometry为多边形的顶点信息
进行可视化表达：
#得到各个属性的可视化表达，若属性值为连续变量，则会自动 调用冷暖色调来做数值区分；如果是类别变量，利用不同的颜色区分
plot(data)
数据框创建空间数据类
有时候我们没有.shp文件，只有点的信息数据。这时候可以利用st_as_sf函数进行创建空间sf类。
df <- data.frame(lon = c(-68.783, -69.6458, -69.7653),
lat = c(44.8109, 44.5521, 44.3235),
Name= c("Bangor", "Waterville", "Augusta"))
#coords=设置坐标,crs=设置坐标系，4326代表WGS84坐标系
point<-st_as_sf(df, coords = c("lon", "lat"), crs = 4326)
point
读取栅格文件
R包：raster
library(raster)
img<-raster("elev.img")
sf类转换为其他空间数据类
转换函数as()
#data为sp类
data.sp<-as(data,"Spatial")
class(data.sp)
#[1] "SpatialPolygonsDataFrame"
#将sf多边形类转换为owin类
#需要maptools包
library(maptools)
data.owin<-as(data,"owin")
#将sf点数据转换为ppp类，ppp类和spatstat包关联
#需要先把sf类转换为sp类，然后将sp类转为ppp类
#需要maptools包
library(maptools)
data.sp<-as(data,"Spatial")
data.ppp<-as(data.sp,"ppp")
如果sp类转为ppp类时，出现如下错误：
Error in as.ppp.SpatialPointsDataFrame(point) : Only projected coordinates may be converted to spatstat class objects
ppp类与spatstat包相关联，该包是在笛卡尔投影坐标系下进行运算的。错误消息是：地理坐标系不能在该包下运算，因此需将地理坐标系转换为投影坐标系。
data.utm<-st_transform(data,32619)#st_transform()表示坐标系统转换，32619代表epsg代码，utm投影坐标系
data.sp<-as(data.utm,"Spatial")
data.ppp<-as(data.sp,"ppp")
sp类转换为sf类
st<-st_as_sf(sp)
空间数据文件导出
st_write(data, "shapefile_out.shp", driver="ESRI Shapefile")