r语言 siar 代码_R语言实践代码汇总

持续更新一波最近遇到的R语言实践过程中的问题。竖着生成一列data frame并把生产的值传给另一个data frame的一列

ID

label

label1

ID

label

label2

data.label

data.label

x$label = data.label$label

2. 给data frame重新命名

names(x)

3. 生成n个重复的sequence

rep(1, times=40)

4. 生成n个不重复按照步数迭代的sequence

seq(1, n,1) #迭代步数为1

5. 去除plot坐标轴的默认坐标

plot(x, s1,xaxt='n')#去除x做标

6. 给x做标对应位置添加标记

axis(1,at=seq(0.00, 0.28, 0.04),labels=seq(0.00, 0.28, 0.04))

7. 数据框同一乘以10

x$x1 = x$x1*10

8. VIF检测多重共线性

library(car)

vif(model)

说明:

得到各个系数的方差膨胀因子,一般认为,当0=100,多重共线性非常严重。

9. 用函数kappa()检验多重共线性

df

df_cor=cor(df)

kappa(df_cor, exact=T)

说明:

当κ<100κ<100, 说明共线性程度小;

当1

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
要在R中利用MODIS数据和像元二分模型计算植被覆盖度,可以按照以下步骤进行: 1. 下载MODIS数据,可以使用MODIS数据下载工具,如MODIS Reprojection Tool (MRT)或MODIS Download Script。下载的数据应该是以HDF格式存储的。 2. 安装并加载raster和MODIS包,以便在R中进行数据处理和分析。 3. 利用raster包中的brick()函数将下载的HDF格式数据转换为raster格式,并加载到R中。例如: ```R library(raster) library(MODIS) #设置下载路径和文件名 modis_path <- "path/to/modis/data" modis_file <- "MOD13Q1.A2001001.h19v04.006.2015130075102.hdf" #下载MODIS数据 runGdal("MOD13Q1", begin='2001-01-01', end='2001-01-01', tileH='h19v04', outDir=modis_path, outProj='+proj=longlat +datum=WGS84', verbose=T) #将HDF格式文件转换为raster格式 modis_raster <- brick(file.path(modis_path, modis_file)) ``` 4. 对MODIS数据进行预处理,包括云掩膜、地表温度计算等。可以使用MODIS包中的函数,如MODIS_process()和MODIS_LST()。例如: ```R #云掩膜 cloud_mask <- MODIS_process(modis_raster, SDSstring='250.0', scale_factor = 0.0001) #计算地表温度 LST <- MODIS_LST(cloud_mask, emissivity = 0.95, method = 'MOD11A2', scale_factor = 0.02) ``` 5. 利用像元二分模型计算植被覆盖度。可以使用像元二分模型的R包,如MixSIARSIAR。例如,可以使用MixSIAR包中的mixsiar()函数进行计算: ```R library(MixSIAR) #将MODIS数据转换为MixSIAR格式 LST_mixsiar <- rasterToMixSIAR(LST) #计算植被覆盖度 veg_cover <- mixsiar(LST_mixsiar) ``` 6. 可以将计算结果可视化,以便更好地理解和分析。可以使用raster包中的plot()函数或ggplot2包中的ggplot()函数进行可视化。例如: ```R #使用raster包进行可视化 plot(veg_cover) #使用ggplot2包进行可视化 library(ggplot2) ggplot(data.frame(veg_cover), aes(x = x, y = y, fill = value)) + geom_raster() + scale_fill_gradient(low = "white", high = "darkgreen") + theme_void() ``` 以上是利用MODIS数据和像元二分模型在R中计算植被覆盖度的基本步骤。具体的数据处理和分析步骤可能会因数据类型和分析目的而异。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值