利用R中趋势包的Pettitt.test (x)函数检测栅格叠加的变化点
Pettitt.test (x)函数是R语言中趋势包中的一个函数，用于检测栅格叠加的变化点。该函数基于Pettitt检验方法，可以帮助我们确定时间序列数据中是否存在突变点。
栅格叠加是一种常见的数据处理方法，用于将多个栅格数据叠加在一起，以便进行分析和可视化。栅格数据通常用于表示地理空间上的属性，如地表温度、降雨量等。
Pettitt.test (x)函数的输入参数x是一个时间序列数据，可以是栅格数据或其他类型的时间序列数据。该函数会计算出时间序列数据中的突变点，并返回突变点的位置和显著性水平。
该函数的应用场景包括但不限于以下几个方面：
- 气象学：用于检测气象数据中的突变点，如温度、降雨量等的突变。
- 地理信息系统（GIS）：用于检测栅格数据中的突变点，如地表温度、土地利用等的突变。
- 经济学：用于检测经济数据中的突变点，如股票价格、GDP等的突变。
- 环境监测：用于检测环境数据中的突变点，如水质、空气质量等的突变。
腾讯云提供了一系列与云计算相关的产品，可以帮助用户进行数据处理和分析。以下是一些推荐的腾讯云产品和产品介绍链接地址，可以用于支持栅格叠加的变化点检测：
- 腾讯云云服务器（CVM）：提供弹性计算能力，用于处理大规模数据计算任务。产品介绍链接：https://cloud.tencent.com/product/cvm
- 腾讯云云数据库MySQL版：提供高性能、可扩展的数据库服务，用于存储和管理时间序列数据。产品介绍链接：https://cloud.tencent.com/product/cdb_mysql
- 腾讯云人工智能平台（AI Lab）：提供丰富的人工智能算法和工具，可用于数据分析和模型训练。产品介绍链接：https://cloud.tencent.com/product/ai
- 腾讯云物联网平台（IoT Hub）：提供物联网设备接入和管理的能力，可用于接收和处理传感器数据。产品介绍链接：https://cloud.tencent.com/product/iothub
通过以上腾讯云产品的组合，用户可以构建一个完整的云计算环境，支持栅格叠加的变化点检测任务。同时，腾讯云还提供了丰富的文档和技术支持，帮助用户快速上手和解决问题。
相关·内容
我希望这成为R.中各种时间序列突破/变化/干扰检测方法的标志，我的问题是用下面的每一个包描述方法的动机和不同之处。也就是说，什么时候使用一种方法比另一种方法更有意义，异同等等。
所涉一揽子方案：
(示例)
(示例)
(链接包括简单示例)
(教程)
也许(?)在较小程度上：和
我希望有针对性的答案。也许每种方法都有一个段落。很容易跨越时间序列，但这可能是以滥用/违反假设为代价的。有些资源为ML监督/非监督技术提供了指导方针。我(当然还有其他人)会欣赏一些关于时间序列分析领域的指南/指点。
浏览 0提问于2015-03-23得票数 10
回答已采纳
2回答
我需要一些关于转换点如何在时间序列中工作的指导。我正在尝试使用R来检测一些转换点，以及称为“转换点”()的包。
对于如何检测方差(cpt.var)和均值(cpt.mean)的变化，有一些选项，但我要寻找的是时间序列变化趋势。
也许我搞不懂真正的转换点是什么，但是有什么方法可以得到这些信息吗？
我正在显示使用cpt.var()函数的结果，并添加了一些箭头，显示了我想要实现的目标。
有办法做到这一点吗？我想应该有点像拐点。
我希望你能对此有任何指点。
谢谢你，乔恩
编辑
我尝试过使用diff()的方法，但没有正确地检测到更改：
我使用的数据如下：
[1] 10.695 10.715
浏览 5提问于2021-05-29得票数 2
回答已采纳
2回答
我有大量的数据要放入一个数据库，这是一个GPS点的列表。
我想遍历这个数据库并创建一个“热点”表，其中在一定大小的区域(无论是正方形区域还是圆形区域--我不需要精确地说)中有大量的数据库点。
有没有人能推荐一些对我有帮助的现有算法？
提前感谢！
r3mo
浏览 0提问于2011-01-26得票数 1
回答已采纳
1回答
在使用prophet库时获取以下错误：
[<-中的错误(*tmp*，m$history$t >= m$changepoints.ti，i，value = 1)：下标超出界限
代码：m <- prophet(data)我从csv文件中加载的数据。
我的数据集如下所示：
ds y
1 2017-05-23 08:07:00 21.16641
2 2017-05-23 08:07:10 16.79345
3 2017-05-23 08:07:20 16.40846
4 2017-05-23 08:07:30 16.246
浏览 3提问于2017-05-24得票数 3
回答已采纳
2回答
我正在使用R，我有每周的数据(总共660个数据)，我想使用来自季节包的X-13阿里玛座位来季节性地调整我的数据。我将数据存储在ts对象中：
library(lubridate)
x <- ts(data, freq=365.25/7, start=decimal_date(ymd("2004-02-01")))
library(seasonal)
x_sa <- seas(x)
但是，我得到了错误：
Error: X-13 run failed
Errors:
- Seasonal period too large. See Section 2.7 of the
浏览 1提问于2018-02-01得票数 3
回答已采纳
2回答
R中的MCMC变点模型
、、、、
我想运行MCMC线性高斯多变点模型来检测连续值的时间序列向量的变点。
在这样做的时候，我正在考虑使用MCMCregressChange函数，但我在这里有几个问题：
(1)如何获得这些模型的对数边际似然？
(2) MCMCregressChange函数和MCMCresidualBreakAnalysis函数有什么区别？
R脚本如下所示。如果你能帮我解决这个问题，我将非常高兴。
library(MCMCpack)
set.seed(1234)
n <- 100
x1 <- runif(n, min = 0, max = 1)
x2 <- runif(n, min = 1, ma
浏览 20提问于2020-07-31得票数 0
1回答
我感兴趣的是每月数据趋势中的多个平均变化年份，即平均变化显著的年份。我遇到了一个R包，并将我的问题描述为
library(changepoint)
df <- read.csv("test.csv")
tdf <- ts(df$x, frequency = 12, start = c(1991, 1))
d <- decompose(tdf)
tr <- as.vector(d$trend)
tr2 <- tr[13:432]
cp <- cpt.mean(tr2,penalty="SIC",pen.value=0.05,
浏览 1提问于2015-11-18得票数 0
1回答
在R中检测业务销售数据中的异常值(这意味着可能存在趋势和季节性)需要哪些步骤？
我已经学习了ACF，PACF，残差，ARIMA模型(基本上是时间序列分析和建模)。我可以使用这些知识来帮助我识别异常值吗？
还有没有可能要求R指出哪个数据点是异常值？
非常感谢。
浏览 0提问于2014-07-15得票数 3
回答已采纳
我有的行政边界数据。使用这些数据绘制细分非常容易。例如，如果我想在佛罗里达州的县上叠加文件mydata中包含的地理编码数据，我会下载美国的"level 2“，然后使用以下代码：
load("/home/anindya/Desktop/DELETE/USA_adm2.RData")
temp = as.data.frame(gadm) *Code for eyeballing structure of data
florida = gadm[344:410,]
plot(florida); points(mydata$longitude, mydata$latitude)
浏览 1提问于2012-02-09得票数 1
回答已采纳
当我浏览异常检测器API文档时，我发现有两个API(或“模式”)：/last和/entire。医生说他们是流模式和批处理模式。但是，我不认为这条消息非常清晰，这两种API/模式似乎具有非常相似的功能。我有一些来自工厂层传感器的IoT数据；我可以对数据进行预处理，以确保它符合API的要求；我可以在我的应用程序中使用C#编写代码。有人能帮我详细说明如何为我的场景选择更好的API吗？
我在Azure笔记本中尝试过这两种APIs。
浏览 2提问于2019-04-19得票数 1
回答已采纳
我有一个shapefile，它是智利的国境线。如果可能的话，我想在这个区域内随机生成1000个多边形。我想知道在R中做这件事最好的方法是什么？我想也许可以，使用'shapefiles‘包加载shapefile，然后使用'spatstat’中提供的东西创建一个点数据集。然后我可以在这些点周围创建圆形多边形。这在很大程度上是假设的，我不知道是否可能。
我在R方面有一些经验，但之前还没有真正尝试过任何空间分析。
如有任何想法和建议，我们将不胜感激。
浏览 5提问于2012-04-03得票数 3
回答已采纳
1回答
gis面地图叠加相交操作
、、、
矢量数据格式的二值地图叠加运算有许多算法，它们以两层地图作为输出，并生成结果层，即叠加层。我想知道是否有任何算法可以同时使用两个以上的层，比如三个层，并产生覆盖结果？
浏览 0提问于2011-07-18得票数 5
回答已采纳
4回答
是否有一个全面的开放源码包(最好在python或R中)可以用于时间序列中的异常检测？
在scikit中有一个支持向量机软件包--学习，但它不是针对时间序列数据的。我正在寻找更复杂的包，例如，使用贝叶斯网络进行异常检测。
浏览 0提问于2018-05-24得票数 24
回答已采纳
1回答
我有几个图层的堆叠栅格数据集，但是，我想计算每个单元格的总和，用于不同的图层选择，并最终生成一个新的图层，谁有一些好的建议，使用计算或叠加或一些其他栅格计算在R？我可以通过循环进行计算，但当我有很多层时，它会消耗很多次，也会使用很多存储空间，我的脚本如下所示：
## library(raster)
make_calc <- function(rr, start, end) {
rr <- as.array(rr)
start <- as.array(start)
end <- as.array(end)
dms <- dim(raster)
t
浏览 1提问于2018-04-12得票数 0
我正在尝试在R中构建四个时间序列模型
1)趋势恒定的HoltWinters，
2)具有线性趋势的HoltWinters，
3)具有恒定趋势的分步自回归。
2)具有线性趋势的分步自回归
在SAS中，我可以使用PROC Forecast并指定方法和趋势选项来执行此操作。
你能帮我做这件事吗？谢谢。
浏览 3提问于2013-03-04得票数 0
回答已采纳
1回答
我在时间序列异常检测方面非常新，所以我的问题对你们中的一些人来说是显而易见的。今天，我使用lstm和聚类技术来检测时间序列上的异常，但这些方法不能识别出随着时间的推移而变慢的异常(我认为这称为趋势)，即机器温度在一个月内缓慢上升(lstm将了解这一趋势，并在没有任何特殊误差的情况下预测增长)。有这样一种方法来检测这种花椰菜吗？
浏览 3提问于2017-09-09得票数 0
回答已采纳
我正在开发一个iOS应用程序，它使用Polyline选项在源到目的地之间绘制地图上的路径。
这个很好用。
我需要的是，如果用户从源开始，而不是选择向目的地绘制的路由，而是选择另一条路径，则应该重新绘制地图上的路径，以获取用户的当前路径。
我正在使用Directions绘制路径。
我怎么能这么做？
提前谢谢。
浏览 1提问于2018-11-30得票数 0
回答已采纳
3回答
R可以用于GIS吗？
、、、、
我想要创建一些GIS地块，我想知道R是否可以用于这个项目。下面是一些我想要制作的情节的例子，在概念上和我想做的相似：
美国的温度图(或等高线图)，颜色(或高度)由州GDP决定。因此，国家边界将在所产生的地块中产生不连续性。
美国的一种温度图，用高度作为数据。在这种情况下，生成的地块应该在不同的州边界之间平稳地变化。
以上两幅图之和(加一些比例)。
我刚开始学习R，想知道它是否是适合这种工作的工具。从美国央行( of ggplot2 )的数据来看，似乎有可能将数据叠加到美国身上。但是要获得数据以尊重州的边界可能是非常困难的。
有什么建议吗？
浏览 6提问于2011-04-22得票数 30
回答已采纳
因此，我正在尝试开发一个非平凡的Angular2应用程序。我有些头疼，因为UI有时没有更新以反映应用程序的状态。我已经搜索了很多关于如何解决这个问题的方法，我也找到了一些我尝试遵循的建议，但是我并没有完全理解它应该是这样的方式或原因。在这里，我想问一些关于我收集到的一些事情的问题。请参阅此参考代码：
Person.ts
export class Person {
name: string;
constructor(name: string) {
this.name = name;
}
}
DataService.ts
@Injectable()
export class D
浏览 2提问于2016-06-16得票数 0
我已经从包含R中某一景观的高度值的光栅中创建了以下过渡对象：
wd <- "C:/Users/LG/Dropbox/Random Walk"
setwd(wd)
library(gdistance)
library(raster)
r <- raster("altitude.tif")
alt
浏览 0提问于2016-04-11得票数 0
1回答
我目前正在实施一个预测模型。整个过程建立在回归的基础上。第一种模型是基于Keras的LSTM网络，第二种是基于sklearn的随机森林模型。
对于这两个模型，我使用相同的输入值和目标变量。LSTM网络的预报更接近实际观测值(较好的rmse和r2) (可能是由于记忆/回溯)，但随机森林模型提供了更好的高度和深度。随机森林模型似乎发现了比LSTM模型更好的依赖关系。
LSTM：
📷
随机-森林：
📷
LSTM模型理解更好的依赖关系的最佳方法是什么？我应该使用来自LSTM + RF的集成或叠加方法吗？说得通吗？
浏览 0提问于2018-05-20得票数 0
我有一个我认为是一个棘手的问题，所以我期待着听到一些选择--这是我最喜欢的工作示例：
cellID X Y Area AvgGFP DeviationGFP AvgRFP DeviationsRFP Slice GUI.ID
1 1 18.20775 26.309859 568 5.389085 7.803248 12.13028 5.569880 0 1
2 2 39.78755 9.505495 546 5.260073 6.638375 17.44505
浏览 22提问于2013-10-02得票数 0
1回答
如何将趋势线添加到图表？
、、、
对于外行来说，这是一条趋势线：
我如何在像MPAndroidChart这样的安卓图表库中做到这一点？我已经设置了必要的代码来生成所需的折线图，但我不知道如何在其上叠加趋势线。
浏览 3提问于2017-03-19得票数 1
4回答
我在一家大型金融服务公司担任分析工作。我们每天都要对很少有意义地改变的指标做大量的报告。从这份每日报告中，我们的管理层必须提取昨天重要的内容，以及随着时间的推移已经形成/正在发展的重要趋势。
我想把这个转变为每日异常报告和每周趋势报告的模式。
其特点可包括：
用户报告整合(所以每天只有一封电子邮件)
基于与以往绩效差异级别的报告排序(首先请参阅最重要的内容)
HTML电子邮件支持(对我的观众来说，非常重要)
Web接口，允许更改首选项，包括LDAP支持(使管理更容易)
在报表级别取消订阅功能
我想知道的是：
我可能遇到的实际问题是什么？
显示新报告的最佳方法是什么？
我应该如何定义“异常”？我
浏览 0提问于2014-10-25得票数 6
我正在使用Socket.IO服务器运行Node.JS，它通常使用大约400 MB的内存，因为有大量的数据被缓存以发送到客户端。然而，几个小时后，它突然开始在大约40分钟内增加到1.4GB的使用量。有人告诉我要用堆转储来检查是否有内存泄漏。
问题是堆转储结果仅为317 MB，其中没有任何东西看起来很平常，所以我一直在调试。我还使用nodetime运行了它，其中显示V8堆的使用率约为400 MB，但V8堆的总大小为1.4GB。
如何找出剩余的1GB来自何处？
浏览 3提问于2014-03-28得票数 5
1回答
我正在尝试使用海拔光栅包获取高程光栅，但我一直在“在选择函数'spTransform'：NA的方法时计算参数'CRSobj‘时出错。”
它最初在我的笔记本上使用相同的数据，但现在在我的台式机上，我一直收到这个错误。什么地方可能出了问题？
浏览 11提问于2022-03-14得票数 0
2回答
我正在使用rgdal包在R中编写一个图像分类脚本。所讨论的光栅是一个具有28个通道的PCIDSK文件:一个训练数据通道、一个验证数据通道和26个光谱数据通道。目标是填充一个数据帧，该数据帧包含训练数据通道中非零的每个像素的值，以及26个频带中的相关谱值。
在Python/Numpy中，我可以轻松地将整个图像的所有波段导入到多维数组中，然而，由于内存限制，R中唯一的选择似乎是逐块导入此数据块，这非常慢：
library(rgdal)
raster = "image.pix"
training_band = 2
validation_band = 1
BlockWidth = 5
浏览 2提问于2010-11-16得票数 2
回答已采纳
我使用Leaflet插件"Leaflet.ImageOverlay.Rotated.js“来使用它的L.imageOverlay.rotated(...)为了在法线贴图顶部的不同位置叠加特定的贴图片段。
它通过获取一张图像并让我告诉它左上角、右上角和左下角的坐标来确定如何正确地旋转、倾斜和拉伸/挤压它。
我花了很长时间才用手算出这些坐标。出于这个原因，我正在寻找某种“地理定位模式”，可能是由这个扩展启用的，它可以让我简单地在地图上点击三次，告诉它这些点的位置。这对开发人员来说很简单，而且会有很大帮助。这是一件显而易见的事情，我强烈怀疑它已经实现并准备好了。
有没有这样的“模式”？如果不
浏览 0提问于2018-04-08得票数 1
有人知道从R的ESRI v10文件地理数据库中获取(栅格)数据的方法吗？
ESRI提供了一个 (用于Linux和Windows)，所以我猜原则上R包应该可以检索(和写入)数据到Geodatabase。不过，我找不到任何能够做到这一点的软件包。
浏览 1提问于2012-12-11得票数 5
2回答
我正在尝试从一个巨大的销售-交易数据集中(超过100万个观察值)中发现异常，该数据集中有数千个独立的客户。同一客户可以在同一日期多次购买。数据集包含随机和季节性交易的混合。下面是我的数据的一个虚拟样本：
Date CustomerID TransactionType CompanyAccountNum Amount
01.01.19 1 Sales 111xxx 100
01.01.19 1 Credit 111
浏览 34提问于2020-05-18得票数 1
5回答
我需要帮助分解我的每月数据，其中有季节性，但它不工作，因为NA值没有被删除。可能还有另外一个问题。请看我的数据和错误如下。
ts.monthly<-ts(monthly$rBC.median, frequency=12, start=c(2006, 4))
ts.monthly
Jan Feb Mar Apr May Jun
2006 5.1656479 6.2847959 19.4833690
200
浏览 3提问于2014-07-11得票数 10
回答已采纳
我使用DoCheck检测来自ngFor的输入属性中的重要更改。正常工作。但是，当我做了一些更改(从数组中删除一个由ngFor迭代的项)时，我的子组件不会看到任何更改，也不会重新呈现自己。
也许我做错了什么？
ngDoCheck() {
console.log('zz');
var changes = this._differ.diff(this.schedulerData);
console.log(changes);
if(changes) {
changes.forEachAddedItem(r => {
浏览 3提问于2016-07-08得票数 0
回答已采纳
我正在寻找一些教程，我如何使画布或SVG (不确定哪一个)的图像在一个等轴测地图的RTS游戏的目的。图像应该只能在可见部分上单击，否则用户应该单击背景图像。它可以通过自动检测PNG的透明区域或从颜色叠加来完成。我不知道哪一个是最好的。
谢谢。
浏览 3提问于2013-03-06得票数 0
我很难分析从下载的MODIS数据(Terra和Aqua)。
数据在一定程度上跨越了一段时间，所以我将不同的层堆叠在一个文件中。
library(raster)
setwd("~/NDVI/")
NDVI <- list.files(full.names = TRUE, pattern = ".tif$")
sNDVI <- stack(NDVI)
这里没有问题，但是当我看每一层的观测数时，它们在层间是不同的。这使得层间的进一步分析变得非常困难。
layers <- 0; layers <- as.data.frame(layers)
浏览 2提问于2018-08-06得票数 0
1回答
我有一个由5个栅格组成的光栅堆栈(这里有一些例子)，我想从这个堆栈中绘制两个最大和最小值的图--使用max/ min (a)非常简单。
r <- raster(nrows=10,ncols=100)
r[] <- rnorm(1000)
s <- stack(r,r,r,r,r)
a<-max(s, na.rm=T)
plot(a)
然而，我想做的是情节，但应用光栅特定的颜色。例如，如果单元格中的最大值来自堆栈中的栅格3，则图中的单元格将被着色为红色，在下一个单元格中，最大值来自光栅5，则为绿色等等。对于这样的光栅叠加图，是否有应用光栅特定颜色的方法？
浏览 2提问于2017-02-26得票数 1
回答已采纳
3回答
颤振图像悬停覆盖效应
、、、、
我正在做一个颤振的网页项目，我想要对我的图像的叠加效果，每当光标盘旋在图像上，几个按钮应该出现，而图像的不透明度被降低(或模糊)，同时徘徊。
我试过、InkWell、和InkWell，但似乎没有什么效果。
这基本上就是我想要达到的目标：。
浏览 4提问于2021-03-26得票数 2
回答已采纳
1回答
我试图设置一个GLMM，但在用glmer表示R中的分层数据结构时遇到了问题。我的数据有以下结构：
y (dependent variable);
visit (L1: nested independent variable);
year (L2: nested independent variable);
site (L3: independent variable);
如果visit嵌套在year中，则嵌套在site中(例如，每年可能有多达5次访问，其中每个站点最多可能有10年的数据)。
我感兴趣的是year对y的影响，同时也希望将site视为一种随机效应(即拦截量和坡度随站点的不同而变
浏览 5提问于2014-07-09得票数 2
回答已采纳
我有一些高质量的LiDAR数据。我需要动态处理它，并根据我的输入和当时的高程更新一个字段(而不是高程字段)。一旦完成此操作，我需要在Web地图(必应地图，谷歌地图等)上显示此栅格数据。尽可能的实时。这是我到目前为止所能做的:我已经想出了一些使用ArcGIS Python Lib的脚本，它们处理我的激光雷达剪辑，添加一个字段，根据输入更新该字段。所以现在我基本上每5分钟就把栅格上的所有点分成3类。如果我在arcGIS中查看这个栅格，我可以看到叠加在底图上的漂亮的3色栅格。
但是现在我需要在Web地图上渲染这些光栅剪辑。我的第一个想法是对这些分类剪辑进行多边形化，并在JavaScript地图上显示
浏览 3提问于2011-08-17得票数 2
回答已采纳
1回答
我正在使用QGIS，并试图在栅格层中每个像素的中心生成均匀的空间点。我能够生成均匀分布的点，但无法使它们足够靠近或位于每个像素的中心。每个像素应该是1.85m，并且在UTM中的投影CRS上。我没有访问ARIGIS的权限。在QGIS中有没有办法做到这一点？甚至是在R中？我尝试生成一个网格以及生成点，但默认距离是0.0001，当我试图变小时，它会自动恢复。考虑到我在UTM，我也不确定如何解释这个默认值。
浏览 6提问于2018-01-31得票数 1
我有两个尺寸相同的二进制文件(栅格)：第一个表示2个数据之间的相关性，第二个表示具有10个类别的土地覆盖图。我想取基于土地覆盖类别的相关性的平均值。因此，最终我们将得到一个与土地覆盖图相同的地图，但具有相关性的平均值，而不是类别编号。
下面是两个文件：
1- to read the first file correlation map:
conne <- file("C:\\corr.bin","rb")
corr<- readBin(conne, numeric(), size=4, n=1440*720, signed=TRUE
浏览 0提问于2013-02-23得票数 2
回答已采纳
1回答
如何求三线性插值函数的导数？
、、、、
我是python的新用户。我有一个3D规则网格数据作为h5文件格式。我可以使用RegularGridInterpolator对我的数据进行插值(三线性插值)。但是，我不知道如何从我的插值函数中求导数。
(我的问题与类似，但似乎没有合适的解决方案，所以我想我应该再问一遍。)
实际上，我正在尝试对我的h5文件使用与这个问题相同的代码()。我的h5文件可以从这个链接()下载。
期待专家的好意帮助。
我的代码：
import numpy as np
import h5py
import matplotlib.pyplot as plt
from scipy.interpolate import
浏览 0提问于2019-02-04得票数 0
我正在开发一个使用谷歌脸检测API的应用程序，我已经使用了示例项目，我需要的是，我想添加基于人脸方向变化的覆盖图像(掩码)，例如:如果脸被旋转到右侧或左侧，我想通过使用坐标values.How来做this?.Can anyone help.Thanks来更新叠加图像。
浏览 3提问于2017-05-24得票数 2
回答已采纳
2回答
我在R中使用非常大的栅格，我有一个代码，我以前已经能够用相同的大栅格运行过，但是当我试图再次运行我的代码以获得可重复性时，我得到了相同的错误“在栅格IO期间失败”。
例如，在代码的这一部分中，我试图将一个大的栅格(~ 5 GB)分成较小的块，这样我就可以在R中加快处理时间。
large_raster <- stack("D:/large_raster.tif")
tileno <- 2 # split the large raster into 2 tiles
cl <- makeCluster(3)
registerDoSNOW(cl)
foreach(i
浏览 67提问于2020-05-29得票数 0
回答已采纳
问题背景：
用于分析糖尿病患者血糖读数的软件通常有一个叫做“每日模式”的图表。
此图表类型是一张线图，显示X轴上的一天的时间，Y轴上的血糖读数，并在指定的日期范围内显示这一点。例如，我可以看到一张我女儿的读数图表，比如上周--每天的读数都叠加在一张图表上(这是关键--将每天的血糖读数叠加在同一天时间的图表上)。这是一天中观察趋势的最好方法。此外，如果有一条趋势线显示一天中某段时间内的平均数，那就更好了。
问题：
有人知道有一个JavaScript图表库可以让我构建这样一个图表吗？我想知道Flot是这么做的，还是Google可视化API呢？
提前感谢
浏览 1提问于2010-03-01得票数 3
回答已采纳
1回答
对于时间序列数据，自动检测趋势或绘制趋势线(上升趋势、下降趋势、无趋势)的好算法是什么？如果你能向我介绍任何好的研究论文或在python，R或Matlab的好图书馆，非常感谢。
理想情况下，该算法的输出将有4列：
from_time
to_time
趋势(向上/向下/无趋势/未知)
probability_of_trend或degree_of_trend
非常感谢你抽出时间。
浏览 10提问于2017-10-03得票数 6
我正在尝试复制这种效果：唯一的区别是我希望粒子以这样一种方式定位，即它们类似于地球-如果你愿意的话，它是一个“有纹理的脸”。
浏览他们的代码，这是他们用来设置粒子组的代码：
var map = THREE.ImageUtils.loadTexture('/admin/wp-content/themes/hys/assets/img/particle.png');
var m = new THREE.ParticleSystemMaterial({
color: 0x000000,
size: 1.5,
map: map,
浏览 2提问于2015-02-03得票数 0
1回答
在R中按属性提取栅格
、、、
我一直在尝试创建一个新的栅格对象，它只包含现有栅格中的几个值。我使用的是这里找到的类栅格：。
class : RasterLayer dimensions : 14902, 19161, 285537222 (nrow, ncol, ncell)
resolution : 0.002349, 0.002349 (x, y)
extent : 110, 155.0092, -45.0048, -9.999999 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs
浏览 0提问于2012-03-03得票数 0
回答已采纳
1回答
我有一个栅格，其中每个属性字段表示火灾风险模型中的不同输入。我正在试着测试我的模型预测火场位置的能力。我用多边形表示两年的火圈(没有用来训练原始模型的年份)。
我想在基于着火多边形的栅格中创建一个二项式存在/缺失字段，其中存在表示由着火多边形覆盖的单元，而不存在表示在研究期间未燃烧的区域。我想包括完全和部分被火多边形覆盖的单元作为存在。在ArcMap或R中有没有好的方法可以做到这一点呢？
我考虑过栅格化我的多边形，但这看起来很笨拙。如果可能的话，我希望避免丢失分辨率。
浏览 5提问于2020-11-14得票数 0
1回答
在R中，我试图创建一个合唱图。我建立了一个企业数据库，其中一些是连锁企业(如McDonalds)的一部分，另一些则是独立的。我想计算出地图上每点30公里内有多少家企业，但把不同的连锁店位置当作一个企业。
例如，如果您的要点是：
Chick-Fil-A20km距5公里，距塔科贝尔10公里，距KFC15公里，距KFC25公里，距<code>F 213</code>
颜色显示在30公里内有4个快餐店。
我很高兴使用任何R包，但我大多熟悉tmap和ggplot2地图。
在这个阶段，我能想到的最好的方法是为每个链创建多边形，并将它们叠加成透明的相同颜色的层。我不认为这将是非常有效
浏览 1提问于2020-12-07得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
即时通信 IM
对象存储活动推荐
腾讯云开发者
