- 博客(191)
- 资源 (16)
- 收藏
- 关注
RSS订阅
原创 mapbox+eCharts路线大数据渲染
eCharts可以实现大量数据迅速渲染，笔者在项目中，测试包含25w余拐点的总计7w余条线渲染总共用时1秒左右，性能表现相当优秀，显示效果也相当不错。eCharts天生支持mapbox，就像其天生支持百度地图一样，而选择前者的原因是其地图相当漂亮，而且支持相当程度的用户自定义地图。效果图：1，添加相关JS引用包含mapbox底图引用以及eCharts引用 <s...
2018-10-27 16:19:21
8262
7
原创 Web端人脸识别登录
人脸识别技术在当下已经十分成熟，但主要在移动端应用上较为普及，而在Web端并不多见。本文介绍在Web端人脸识别的简单实现。Web端人脸识别主要有三个技术思路：1.前端的人脸识别，例如使用Tensorflow.js，2.后台人脸识别，有很多开源或者免费的SDK可以使用，3.前后端结合，即结合以上两种方法，虽然系统复杂度提高，但对于系统的安全性，以及减轻服务器负担都有很大提升。效果图：...
2018-08-04 20:29:37
29426
60
原创 增强现实之开源AR库——AR.js
AR.js是一个web端的AR库，它完全开源免费，获得了很高的热度。我们要实现的效果如下：首先去github下载AR.js库：AR.js下载建议顺带看下作者给出的介绍。介绍里给出一个示例，我们在此示例的源码进行分析并尝试修改示例中的三维模型。解压缩后目录如下：示例存three.js目录的example目录下。只有在服务器环境内才可以访问。我们可以将AR.js拷贝至Apa...
2018-04-26 09:12:39
25246
94
原创 ArcGIS 安装插件后缺失模块 Missing
删除重装ArcGIS和插件ArcSWAT后问题仍然存在。在社区也没有找到相应答案。事件起因是ArcGIS 10.6安装 ArcSWAT后显示缺失一个模块，原本应该是HRU的模块显示为Missing，2.弹出的窗口中找到ArcSWAT。1.点击 自定义——自定义模式。在重置后问题得到解决。
2023-02-27 15:48:15
1164
1
原创 1960-2020年中国气温线性趋势变化的时空差异及影响因素
说明：文章为笔者原创随笔，引用请标明出处。摘 要 全球变暖显著影响，然而，我们对气温变化的季节和空间差异了解有限。分析了中国1960-2020年气温变化的时空差异和影响因素。首先，分析了全国和5个不同气候区域年平均气温的变化趋势，结果表明1990前后全国和不同气候区域的年平均气温均出现显著突变。1990年后全国和不同区域年平均气温变暖加快，不同季节变化不同：春、夏季增温加快，冬季增温减缓，而秋季全国气温趋势变化不大。随后，利用地理探测器定量比较了两个阶段气温变化与下垫面影响因子的关系。植被类型、纬度地带
2021-10-21 10:38:06
9485
1
原创 ggplot2 一页多图+共享图例
结果图：ggpubr包提供了方便的一页多图方法，并且提供了方便的共享图例方法。因为应用中，一个页面上多幅单一类型的图一般只需一个图例即可。参照以下代码：library(ggpubr)dsamp <- diamonds[sample(nrow(diamonds), 1000), ]p1 <- qplot(carat, price, data = dsamp, colour = clarity)p2 <- qplot(cut, price, data = d.
2021-08-19 20:40:41
4501
原创 Pearson 相关分析 理论与应用
1 前提假设使用Pearson相关分析时，需要考虑满足5个假设。1.两个变量都是连续变量。2.两个连续变量应当是配对的，即来源于同一个个体。3.两个连续变量之间存在线性关系，通常做散点图检验该假设。4.两个变量均没有明显的异常值。Pearson相关系数易受异常值影响。5.两个变量符合双变量正态分布。2 验证前提假设2.1两个连续变量之间存在线性关系Pearson要求两个变量之间存在线性关系。要确定是否存在线性关系，需要查看两个变量的散点图。如果散点图大致呈一条...
2021-08-09 10:27:55
19315
原创 数据分享及分析方法——热门气象数据集
1、NASA的 Land Surface Temperature 基于Modis卫星获取地表反射率，进而转化为地表温度。
2021-06-15 18:55:04
1025
原创 R语言提取单个cru格点气象数据
inputpath <- 'data/cru_ts4.05.1901.2020.tmp.dat.nc'library(ncdf4) # package for netcdf manipulationlibrary(raster) # package for raster manipulationlibrary(rgdal) # package for geospatial analysislibrary(ggplot2) # package for plottingnc_data <.
2021-06-09 15:55:54
817
原创 R语言数据基础离散化（分箱）之——手动定义间断点
R语言本身提供手动定义间断点进行离散化的函数 cut参数如下：cut(x,breaks,labels=NULL, include.lowest=FALSE,right=TRUE,dig.lab=3, ordered_result=FALSE,...)其示例代码如下：...
2021-06-08 22:22:13
1039
原创 R语言树轮数据可视化二——以contour图展示
一般而言，树轮数据多以折线图进行展示，清楚展示每条序列的变化，但序列一多通常难以比较。笔者尝试使用contour图进行展示，可以清楚的展示高低生长期。使用 filled.contour(x, y, z, levels, col) 函数对一个采样点的不同树芯宽度进行可视化，效果如下：其中横轴表示时间（年），纵轴表示不同树轮序列的宽度。可以看出，每条序列基本都存在先高生长后低生长的趋势。另外，不同序列的敏感性是不同的，需要标准化后进行进一步分析。全部代码如下：library(dplR)
2021-05-24 09:23:18
2111
1
原创 R语言提取单点的cru格点数据
以下以提取降水的cru格点数据为例，温度的数据类似。inputpath <- 'data/cru_ts4.05.1901.2020.pre.dat.nc' #文件存储路径library(ncdf4) # package for netcdf manipulationnc_data <- nc_open(inputpath)print(nc_data)# Save the print(nc) dump to a text file{ sink('ncdata.txt')
2021-05-23 22:39:43
917
3
原创 Google Trends 的玩法1——时间序列分析
Google Trends是谷歌公司的公共网络设施,它基于谷歌搜索,它显示整个世界各地区的一个特定搜索项搜索量。我们可以查看2004年来逐月的某关键词搜索量，以及其以国家/城市为单位的空间分布特征。下面以 “柿子” 为例，其搜索量的时间序列有明显的季节周期变化：其每到 10-11月份 会达到峰值，其空间分布特征如下：由于我们是输入中搜索，其热点区域都是中文为主的城市，当我们换为英文“Persimmon”：时间序列结果类似，而其空间分布特征发生了变化：其他季节性的水果.
2021-05-12 10:34:13
582
原创 r语言实现小波分析
实现效果左图为原始数据，右图为小波分析加显著性检测结果所需包：WaveletComp一、安装并加载WaveletComp包install.packages('WaveletComp')加载：library('WaveletComp')二、创造数据并执行小波分析及模拟series.length = 6*128*24x1 = periodic.series(start.period = 1*24, length = series.length)x2 = p..
2021-04-10 19:21:17
5817
3
翻译 R语言netcdf格式数据读取与处理
在本教程中，我们将打开一些存储在netCDF文件中的地理空间数据。我们将选择感兴趣的变量和时间范围，并将数据导出到GeoTIFF，以便在R或其他地理空间软件中继续分析。（译者注：新手看可能不知道每个函数的用法，可以先根据文档仿写，记住即可，想要深入研究每个函数的用法，使用RStudio打开包帮助阅读即可）本文读取netcdf并可视乎数据的结果：NetCDF数据集示例首先，我们需要一些数据。作为一个例子，我们将使用北极地区植被绿色度趋势的一些数据。可通过以下链接从ORNL DAAC获取此数据：
2021-04-10 16:46:50
9199
2
原创 气象数据处理——数据表重排列（长表宽表转换）
一、长表转换为宽表好多气象数据（如气象站点、格点等）是用长表存储的，其格式如下：长表格式常用于数据的存储，尤其对于多维数据的存储十分方便，然而并不利于数据分析。R语言tidyr包中的pivot_wider函数 可以将长表变为宽表，适用于将多个变量保存到了一列的情况。如上图示例，我们需要将表中数据转换成：每年（列）不同站点（行）的观测值，我们使用以下代码：d.tempwider <-d.monthtemp1 %>%pivot_wider(names_from = stat
2021-04-08 18:17:48
596
原创 R语言计算多变量两两之间的相关关系及显著性
当数据存储为data.frame格式时，我们通常先将其转为可以使用 PerformanceAnalytics 包之间求相关，并将结果可视化。chart.Correlation(dda, histogram=TRUE, pch=19)其中，dda 为data.frame格式的数据，要求数据格式为 numeric ,如果不是，要先转化...
2021-04-08 11:48:34
7322
1
原创 R语言读取.csv数据报错
示例：Error in nchar(x[is_na], type = "width") : invalid multibyte string, element 1出现读取数据报错，很多时候由数据读取默认的编码格式 与 数据本身的编码格式不同 引起的表格中，如果含有中文或中文符号，一般默认为 GBK 编码，而一边编程语言读取时 默认为 UTF-8 编码此时设置 读取数据编码为 GBK 即可：示例：read_csv('data/treeringTP.csv',locale=lo
2021-04-07 20:37:41
10857
原创 ArcGIS将太平洋调至地图中间（改变地图中心经度）
很多时候，我们加载数据后结果大西洋在地图中间（下图左），而我们想展示太平洋全貌（下图右），类似这种要求，即需要调整起始经线位置。1、新建图层坐标系一般地理坐标系投影左端点为-180°，我们需新建图层的坐标系，并指定新的坐标原点：在数据视图上 右键—数据框属性—坐标系，然后双击当前坐标，弹出地图坐标系属性窗口，如下图修改坐标名称，如WGS-from180,最重要的一步，修改本初子午线，经度设置为180°，设置如下图所示：点击应用，返回数据框属性页面。此时地图并未变化，因为...
2021-04-06 19:13:03
5977
6
原创 R语言将第一行作为列名
colnames(ClimateData)=ClimateData[1,]格式为：colnames(表名)=表名[1,]将第一列作为行名 类似：rownames(表名)=表名[,1]
2021-01-20 23:57:15
27150
1
原创 R语言DataFrame(数据帧)数据格式改为数值
使用as.numeric()即可将一行数据强制更改为数值，参加数值运算。e.g. 更改前："2","2","2","2","2","2","2","2","2","2","2","2","2","2","2"更改后：2，2，2，2，2，2，2，2，2，2，2，2，2，2，2...
2021-01-20 23:01:02
5500
原创 R语言求取DataFrame(数据帧)多行/列平均值
在Excel 使用average先求取两格子的平均值，然后可以进一步拉两个格子求取整行/整列的平均值。这在R语言中使用rowMeans求取多行平均值，或colMeans求取多列平均值示例如下：#create data framedata <- data.frame(var1 = c(0, NA, 2, 2, 5), var2 = c(5, 5, 7, 8, 9), var3 = c(2, 7, 9, 9, 7)
2021-01-20 22:14:58
27063
原创 R语言在矢量地图上绘制分级设色散点图
实现效果：R语言中ggplot2包提供绘制地图、散点图的方法，是实现在矢量地图上绘制分级设色散点图核心包绘制多边形geom_polygon（data,aes,fill, colour） 绘制点 geom_point .....0、需要用到的包library(maptools) # 读取shp数据常用，可以将shp数据读取为SpatialPolygonsDataFrame 格式，为DataFrame（数据帧）子类，也称为空间多边形数据帧library(ggplot2) #绘图核...
2021-01-20 20:36:51
1909
原创 统计中的“不相关”与“线性无关”
以上思维导图，看完即可理解。下述是文字介绍。这二者是统计新手与老手都很容易混淆的两个概念，以下辨明一下：两变量“不相关”不相关是指二者互相独立，没有相关关系。注如森林里每棵树的树叶个数与村子里每个村民的体重...二者显然无关两变量“线性无关”线性这里单指直线，线性无关是指:没有一次函数关系。如果两变量是3次方的相关关系，这里也叫线性无关。其图示如下：‘可以看出，如果两变量不相关，则必有线性无关。反之则不一定。...
2020-12-19 10:43:13
5059
1
原创 R 语言实现 Superposed Epoch Analysis 时序叠加分析
一、安装并引用 “burnr” 包Superposed Epoch Analysis可以使用 “burnr” 包 来实现，首先安装 “burnr” 包：命令行输入：install.packages('burnr')其次将 burnr 包添加到引用：library(burnr)二、准备输入数据Superposed Epoch Analysis 的输入数据格式如下：1、 准备时间序列数据放入EXCEL,格式如下：之后，读入R，并修改行名称为年份：dat...
2020-12-17 11:26:23
2191
8
原创 R语言实现EMD分析
R语言实现EMD分析非常简单，使用EMD包即可。数据以frame存储时，只需要设置一列作为emd()函数的参数即可：remd=emd(all$runoff)remd即输出的所有IMF，同样以frame存储。可以直接导出为csv.
2020-12-13 18:04:44
2067
原创 mapbox加载turf.js创建的格网
turf.js是mapbox官方推荐的进行空间分析的包，其优点在于强大且全面的空间分析功能，且支持geojson格式的输入，因此与Mapbox GL JS集成效果很好。我们此处介绍Mapbox GL JS加载turf.js实时创建格网（渔网）的方法。其效果如下：turf.js能够指定格网大小、坐标等信息实时生成，其效率非常高，也可以在前端根据用户需求实时的生成格网了。生成结果为GeoJson格式，Mapbox根据其生成的结果直接加载。步骤如下：1、turf.js生成格网（渔网）tur
2020-10-29 18:38:29
1626
1
原创 mapbox获取图层中要素的属性值
Mapbox GL JS提供了强大的查询检索功能，我们可以不需要写后台程序实现对于图层中要素属性值的查询检索功能。例如：我们希望用一个表格或图表展示某个图层中所有或部分要素的某个属性值，我们可以用mapbox的map.querySourceFeatures()方法搞定。该方法可以查询出某一图层的所有要素（包括要素的属性值），如果你想筛选部分要素，只需要在参数中设置Filter即可。使用方法如下：var relatedFeatures = map.querySourceFeatures('c
2020-10-23 13:20:42
4992
原创 ArcGIS Editor for OSM安装出错
ArcGIS Editor for OSM目前不支持ArcGIS10.2.0版本（只支持ArcGIS10.2.1和ArcGIS10.2.2）。如果安装过程中出现如下所示错误：可尝试升级ArcGIS版本，或者请求插件作者大大开发一个ArcGIS10.2.0版本插件。
2020-10-22 10:27:16
1078
1
原创 ArcGIS加载与编辑OpenStreetMap(.osm)文档【最新避坑】
ArcGIS Editor for OSM下载ArcGIS Editor for OSM目前支持ArcMap 10..0至10.3共四个版本。使用ArcGIS Editor for OSM可以在ArcMap中加载.osm文件，应用符号系统，将数据贡献回OSM并从OSM数据创建网络数据集。ArcGIS Editor for OSM下载地址：http://www.arcgis.com/home/item.html?id=16970017f81349548d0a9eead0ebba39从中选择自己
2020-10-22 10:22:39
1788
原创 迈克菲添加白名单
迈克菲是一款非常强力的杀毒软件，我们日常使用许多程序，包括破解、梯子等可能被迈克菲实时检测并阻止运行。而迈克菲提供白名单功能，可以为我们使用一些敏感的软件提供方便。下面是具体操作步骤：1、在主页 打开实时扫描窗口点击红色方框的三个小点。。。2、打开已排除的文件这里面即为白名单的程序内容3、点击添加文件，选择主程序即可这里我们以pr为例子，添加之后，迈克菲扫描时便会跳过该程序啦。...
2020-10-08 11:30:06
15927
2
原创 免费人文数据分享网站（更新中）
Rdatasets大量有关自然及人文环境的统计数据首页：https://vincentarelbundock.github.io/Rdatasets/GitHub地址：https://github.com/vincentarelbundock/Rdatasetsfivethirtyeight大量有关政治及人文信息的数据首页：https://data.fivethirt...
2020-10-05 10:35:17
404
原创 Matplotlib绘制半圆形
Matplotlab绘制图形使用matplotlib.patches如果我们希望绘制一个圆，我们引入其中的Circle，而如果我们打算绘制一个部分圆形（半圆或者1/4圆），我们采用Wedge其参数如下所示：Wedge(center,radius,theta1,theta2)center代表其中心点，即圆心位置radius代表其半径theta1,theta2 代表其起始于结束角度，以x轴正向为0度，逆时针旋转以下为绘制半圆的程序...
2020-09-26 17:13:10
2975
1
原创 Excel实现一个基础的蒙特卡洛模拟
蒙特卡洛模拟貌似名字非常学术，其实其原理非常简单。蒙特卡洛是一个著名赌场的名字，将其命名为蒙特卡洛模拟是为了形象表示其方法含义。如果放在今天命名的话，也可以将其称为拉斯维加斯模拟或者澳门模拟。简而言之，蒙塔卡洛模拟通过将问题转化为概率统计问题进行求解。比如求圆周率的数值，有很多很多方法。蒙特卡洛模拟的方法是，以下图为例子进行说明：向图中x,y范围为[0,1]的正方形区域随机投点，统计落入1/4圆内的次数。在投了足够多的次数以后，计算落入圆内的次数/总次数，即为落入圆内的概率。这个概率值即接近..
2020-09-26 13:53:44
9508
1
原创 如何防止破解软件失效？
在电脑日常使用中，作为个人用户，我们常常不可避免的使用一些破解软件。而一些软件在联网后通常会自己将软件升级而使破解程序失效，为了阻止这种情况发生，常用的一种方法是：用防火墙阻止应用联网。具体操作过程如下：一、打开高级安全Windows 防火墙在底部状态栏搜索框中输入 高级安全Windows Defenser防火墙 打开。二、设置出站规则1、点击左侧出站规则进行设置，继而点击右侧新建规则：2、弹出窗口中选择程序，点击下一步3、选择程序位置此处以chrome浏览器为例..
2020-09-20 19:18:49
1576
原创 Excel 简单线性回归图表制作
Excel内嵌了回归分析与制图功能，用Excel进行回归分析的制图十分方便快捷。1、数据输入以一下数据为例：一共两列数据，一列为自变量X，一列为因变量Y2、插入散点图选中两列数据后，点击 插入|散点图；选择X、Y参数，在表格中插入散点图，如下图所示：3、插入拟合后的线及R方点击图中的点，右键选择 添加趋势线；在面板中勾选 显示公式，显示R平方值，如图：结果如下图所示：...
2020-09-19 18:51:39
15640
相关与回归分析基础教程
2020-12-20
ArcGIS简化面工具Bug修复工具
2019-05-14
python空间分析指南
2019-03-31
Echarts3Layer 修改
2019-01-11
Echarts3Layer修改
2019-01-11
ArcGIS水文分析教程
2018-12-19
eCharts格网统计图
2018-10-27
web端人脸识别登录源码
2018-08-24
jsapi_vsdoc12_v38.js
2018-03-16
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人