撼沧

基于package modiscloud library(modiscloud) 简介： R tools for processing Level 2 Cloud Mask products from MODIS yearmonthday_to_julianday: Get the julian day for a year/month/day date year=2011; month=06; day=15 yearmonthday_to_julianday(year, month, day)

基于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（均方根）；

写在前面 netcdf是气象、生态、地理领域非常常用的一种数据格式，读写的方法很多，arcgis、arcpy、python、matlab等等均有相应的方法。 本文基于R语言实现netcdf的数据读取，提取数据中的单个波段（某一时间的数据），展示和写出。 注：本文仅展示基础的代码，多个波段或多个数据的读写操作加上相应的循环即可。 具体代码 加载包 library(raster) library(rasterVis) library(ncdf4) library(lattice) 打开nc数据，并展示其中

写在前面 以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

写在前面 当前对于高温热浪事件指标的定义缺乏统一标准，不同学者有不同的算法。Russo（2015）提出了一种每日热浪幅度指数（HWMId，the heat wave magnitude index daily）计算方法，适合对全球的热浪强度进行计算和对比分析。该方法也被wiki百科对热浪指数的介绍中采用（wiki）。 该方法作者提供了相关的R语言包和函数，本文对此进行原理和代码的介绍。 放上wiki上的一张图（Zampieri, M., 2016） 1. 原理介绍 HWMId定义为每年最大热浪事件的幅

本文介绍如何使用r语言进行单个或多个栅格数据间的计算，类似于ArcGIS中的栅格计算器。 单个栅格数据： library(raster) calc(x, fun, filename='', na.rm, forcefun=FALSE, forceapply=FALSE, ...) 元素： x：Raster* object fun：function filename：character. Output filename (optional) Examples # NOT RUN { r <- rast

rasterVis包中的levelplot函数除了出单张图，也可以快速随多张地图进行制作。 单张地图制图参考：【R语言】栅格数据快速制图-levelplot函数 此处主要使用raster包中的stack函数将多幅图合为影像集。 #栅格数据包括.tiff、.nc、.flt等常见格式 library(raster) library(rasterVis) library(lattice) #设定工作空间 setwd('E:\\hancang') #读取文件夹中的全部文件 files = list.files(

使用R语言中的levelplot函数可以快速方便地对栅格数据进行制图，非常实用，下面做简单地讲解实例。 栅格制图 # 加载package library(raster) library(rasterVis) library(lattice) 基础展示: 根据坐标轴展示的栅格图像； 两个边缘分布展示为行列的平均值； 颜色图例； levelplot(rprob) 对图像的栅格值进行对数变换(zscaleLog=TRUE)，并添加水平等值线(contour=TRUE)。 levelplot(rprob,

library(stringr) library(raster) library(sp) #set work space getwd() setwd("G:/data") #初始设置，提取点 lai <- matrix(nr=50,nc=2) point<-cbind(128.5,45.9) #lon,lat i <- 1 n <- 1 repeat { #读取栅...