MODIS R

MODIS包也在github上,上面有常见问题模块,可以向开发者提问和提供bug.
1,在EARTHDATA注册用户。
到官方网站注册https://earthdata.nasa.gov/,比如你的用户名和密码分别是:用户名:james;密码:lileihanmeimei123

2,安装MODIS包
install.packages(“MODIS”)

3,下载数据流程
3.1 配置用户名和密码
第一个是earthdata的账户配置,在MODIS包中 EarthdataLogin函数来配置,假如你的earthdata账户和密码是下面的例子:

EarthdataLogin(usr = ‘james’, pwd = ‘lileihanmeimei123’)

在R中运行上面的代码后,会在电脑上生成一个不可见的密码文件。

3.2 配置下载和处理数据目录
lap = “D:/R_MODIS/MODIS_ARC” # ‘设置MODIS下载到本地的位置’

odp = file.path(lap, “PROCESSED”) # ‘设置MODIS处理数据的位置’ MODISoptions(localArcPath = lap, outDirPath = odp) # 配置文件

配置好后,运行 MODISoptions(),控制台出现下面的信息
STORAGE:
localArcPath : D:/R_MODIS/MODIS_ARC
outDirPath : D:/R_MODIS/MODIS_ARC/PROCESSED
DOWNLOAD: _______________
MODISserverOrder : LPDAAC, LAADS
dlmethod : auto
stubbornness : high
wait : 0.5
quiet : TRUE
PROCESSING: _______________
GDAL : 3.0.4
MRT : Enabled
pixelSize : asIn
outProj : asIn
resamplingType : NN
dataFormat : GTiff
cellchunk : 1
3.3 下载数据
(1)先看看可以下载哪些产品
运行 getProduct() 出现了详细的产品列表,目前用166种MODIS产品可以通过这个包下载。详细的介绍可以看官方网站:https://modis.gsfc.nasa.gov/data/dataprod/

https://github.com/MatMatt/MODIS
https://github.com/MatMatt/MODIS/tree/master/R
http://blog.sciencenet.cn/blog-437026-1248955.html
https://gisersqdai.top/2017/04/15/MODIS%E6%95%B0%E6%8D%AE%E7%9A%84%E7%AE%80%E4%BB%8B%E5%92%8C%E4%B8%8B%E8%BD%BD%EF%BC%88%E4%B8%89%EF%BC%89%E2%80%94%E2%80%94MODIS%E6%95%B0%E6%8D%AE%E4%B8%8B%E8%BD%BD%E6%96%B9%E5%BC%8F%EF%BC%88%E5%9F%BA%E4%BA%8EMODIS%20Web%20Service%EF%BC%89/

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是使用Matlab处理MODIS数据的示例代码: 1. 数据读取和可视化: ```matlab % 读取MODIS数据 filename = 'MOD09A1.A2008001.h21v09.006.2015181070345.hdf'; modis_data = hdfread(filename,'MOD_Grid_500m_Surface_Reflectance','Fields','sur_refl_b01'); % 显示MODIS数据 imagesc(modis_data); colormap('jet'); colorbar; ``` 2. 数据处理: ```matlab % 辐射校正 r = 0.293; % 定义地表反射率 L = modis_data / r; modis_data_calibrated = (pi * L) / 2.5; % 进行分类 threshold = 1000; % 定义分类阈值 modis_data_classified = modis_data_calibrated > threshold; % 显示分类结果 imagesc(modis_data_classified); colormap('gray'); colorbar; ``` 3. 数据分析和建模: ```matlab % 时间序列分析 % 读取多个MODIS文件 filenames = {'MOD09A1.A2008001.h21v09.006.2015181070345.hdf', ... 'MOD09A1.A2008002.h21v09.006.2015181070345.hdf', ... 'MOD09A1.A2008003.h21v09.006.2015181070345.hdf', ... 'MOD09A1.A2008004.h21v09.006.2015181070345.hdf'}; for i = 1:length(filenames) modis_data(:,:,i) = hdfread(filenames{i},'MOD_Grid_500m_Surface_Reflectance','Fields','sur_refl_b01'); end % 计算每个像素点的平均值 modis_data_mean = mean(modis_data, 3); % 绘制时间序列图 plot(modis_data_mean); xlabel('时间'); ylabel('地表反射率'); title('时间序列分析'); ``` 以上是一些基本的代码示例,实际上,Matlab有很多工具箱可以用来处理和分析MODIS数据,根据具体的需求选择相应的工具箱和函数,可以更加高效地处理和分析MODIS数据。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值