![](https://img-blog.csdnimg.cn/20201014180756738.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
R语言
文章平均质量分 65
撼沧
From zero to hero.
展开
-
R语言 | 儒略日与日期互相转换 the conversion between Julian day and year-month-day
基于package modiscloudlibrary(modiscloud)简介:R tools for processing Level 2 Cloud Mask products from MODISyearmonthday_to_julianday:Get the julian day for a year/month/day dateyear=2011; month=06; day=15yearmonthday_to_julianday(year, month, day)原创 2021-10-23 22:39:55 · 891 阅读 · 0 评论 -
基于R语言进行栅格数据统计及Raster包简介
基于R语言进行栅格数据统计及Raster包简介写在前面:日常数据处理中常需要对单个栅格图层进行统计,包括计算其均值、最值和总和等等。使用R语言可方面快速地实现该需求,主要借助的package为Raster包。library(raster)x = raster('data/data.tif')cellStats(x,stat='mean')x表示待计算的栅格;stat表示需要统计的内容,包括sum, mean, min, max, sd(标准差), skew(偏度) 和 rms(均方根);原创 2021-05-25 15:41:58 · 11284 阅读 · 0 评论 -
基于R语言读取NC数据并提取单个波段(netcdf转为tiff或其他格式)
写在前面netcdf是气象、生态、地理领域非常常用的一种数据格式,读写的方法很多,arcgis、arcpy、python、matlab等等均有相应的方法。本文基于R语言实现netcdf的数据读取,提取数据中的单个波段(某一时间的数据),展示和写出。注:本文仅展示基础的代码,多个波段或多个数据的读写操作加上相应的循环即可。具体代码加载包library(raster)library(rasterVis)library(ncdf4)library(lattice)打开nc数据,并展示其中原创 2020-12-26 00:33:43 · 10456 阅读 · 2 评论 -
基于R语言或Matlab直接读取压缩包中的数据
写在前面以gosif网站上的tiff数据为例。该数据是地学上常用的一种数据,希望读取该网站上的gosif数据:http://data.globalecology.unh.edu/data/GOSIF_v2/8day/但数据全是压缩包格式的,解压后内存占用又会变大。所以需要直接从压缩包中读取数据。基于R语言library(R.utils)library(stringr)library(raster)library(rasterVis)# download and read the dat原创 2020-12-23 12:47:01 · 3170 阅读 · 0 评论 -
基于R语言进行热浪指数的计算_HWMid
写在前面当前对于高温热浪事件指标的定义缺乏统一标准,不同学者有不同的算法。Russo(2015)提出了一种每日热浪幅度指数(HWMId,the heat wave magnitude index daily)计算方法,适合对全球的热浪强度进行计算和对比分析。该方法也被wiki百科对热浪指数的介绍中采用(wiki)。该方法作者提供了相关的R语言包和函数,本文对此进行原理和代码的介绍。放上wiki上的一张图(Zampieri, M., 2016)1. 原理介绍HWMId定义为每年最大热浪事件的幅原创 2020-12-26 22:57:58 · 3482 阅读 · 11 评论 -
【R语言】使用r语言进行单个或多个栅格数据间的计算
本文介绍如何使用r语言进行单个或多个栅格数据间的计算,类似于ArcGIS中的栅格计算器。单个栅格数据:library(raster)calc(x, fun, filename='', na.rm, forcefun=FALSE, forceapply=FALSE, ...)元素:x:Raster* objectfun:functionfilename:character. Output filename (optional)Examples# NOT RUN {r <- rast原创 2020-12-07 10:08:31 · 7615 阅读 · 0 评论 -
【R语言】栅格数据快速制图-levelplot函数(2)
rasterVis包中的levelplot函数除了出单张图,也可以快速随多张地图进行制作。单张地图制图参考:【R语言】栅格数据快速制图-levelplot函数此处主要使用raster包中的stack函数将多幅图合为影像集。#栅格数据包括.tiff、.nc、.flt等常见格式library(raster)library(rasterVis)library(lattice)#设定工作空间setwd('E:\\hancang')#读取文件夹中的全部文件files = list.files(原创 2020-12-07 09:53:21 · 5308 阅读 · 5 评论 -
【R语言】栅格数据快速制图-levelplot函数
使用R语言中的levelplot函数可以快速方便地对栅格数据进行制图,非常实用,下面做简单地讲解实例。栅格制图# 加载packagelibrary(raster)library(rasterVis)library(lattice)基础展示:根据坐标轴展示的栅格图像;两个边缘分布展示为行列的平均值;颜色图例;levelplot(rprob)对图像的栅格值进行对数变换(zscaleLog=TRUE),并添加水平等值线(contour=TRUE)。levelplot(rprob,翻译 2020-12-06 23:05:47 · 6497 阅读 · 2 评论 -
【R语言】根据经纬度从tif图像中提取兴趣点
library(stringr)library(raster)library(sp)#set work spacegetwd()setwd("G:/data")#初始设置,提取点lai <- matrix(nr=50,nc=2)point<-cbind(128.5,45.9) #lon,lati <- 1n <- 1repeat { #读取栅...原创 2020-04-03 10:45:28 · 3761 阅读 · 2 评论