走天涯徐小洋地理数据科学-CSDN博客

原创 NDVI时间序列分析之Sen+MK分析全过程梳理

Sen斜率估计用于计算趋势值，通常与MK非参数检验结合使用，即先计算Sen趋势值，然后使用MK方法判断趋势显著性。Theil-Sen Median方法又被称为Sen斜率估计，是一种稳健的非参数统计的趋势计算方法。该方法计算效率高，对于测量误差和离群数据不敏感，常被用于长时间序列数据的趋势分析中。β=median(xj−xij−i),∀j>i\beta=median(\frac{x_j-x_i}{j-i}),∀j>i β=median(j−ixj−xi),∀j>i式中：xjx_jxj 和 xix_

2022-07-04 15:40:25 17299 6

原创 rgee尝鲜Landsat9，以及Landsat8对比

rgee尝鲜Landsat9，以及Landsat8对比看到吴秋生老师发布的推文，Landsat 9数据已经加入Google Earth Engine，我也迫不及待想要试一试。吴秋生老师 GEE（谷歌地球引擎）已经开始收录Landsat 9数据试验内容使用rgee查看Landsat9数据交互式地图划定感兴趣区ROILandsat8和Landsat9数据对比使用rgee查看Landsat9数据这部分代码基本上是根据吴秋生老师的推文和我以前写的Landsat的文章中的代码改的：Landsa

2022-03-18 22:09:43 2261

原创矢量科研数据质量检查方法

矢量科研数据质量检查方法矢量数据是地理科研工作当中经常接触的一类数据，那么如何判断自己制作的矢量数据是否合格呢？建议从以下几点进行判断：数据完整性检查位置精度检查，查看矢量数据空间范围是否正确，空间位置定位精度是否符合要求；图形完整性检查，查看矢量数据图形是否完整，查矢量数据图形是否有明显不合逻辑的空缺、飞地、孔洞等属性精度检查，检查数据属性表命名是否规范，和数据字典进行核对，是否有属性数据缺失或错误拓扑检查，通过建立拓扑规则，查看矢量数据是否存在不合理的拓扑错误下面详细逐一解释：1 数

2021-07-06 20:49:51 1900

原创喜大普奔！rgee能用了！R语言也可以使用Google Earth Engine了！

喜大普奔！rgee能用了！R语言也可以使用Google Earth Engine了！去年的时候曾经写过一个推文，关于rgee，但是由于rgee的不知道是BUG还是啥问题，认证一直不成功，没有办法愉快的使用，今天我收到了来自rgee的作者之一的邮件，告诉我了这个BUG已经有了解决方法，我再次试验了rgee，成功运行！rgee安装我使用了R4.0.3版本，建议大家使用最新的R语言版本，R4.0版本我目前已经用了几个月，感觉和3差别不大，以前的一些包也都可以使用，强烈大家更新R4.0最新版本。安装RGE

2021-07-06 20:47:25 2483 4

原创对比R语言Raster包和Terra包栅格计算

对比R语言Raster包和Terra包栅格计算Terra包可以认为是Raster的改进版，据说使用C语言重写了函数，提高了效率，同时支持并行计算，对于大数据量的栅格计算很有优势，我接下来就给大家展示一下我的试验。这次实验数据是MODIS250m的NDVI数据，数据量比较大，有2272×4299个像元，20个波段。为了保证计算公平性，我都是从数据读取开始计时，统一计算变异系数cv，然后绘图，得到运行时间。Raster计算CV#Raster计算CVt1=proc.time()#程序体flname

2021-03-04 16:59:41 2761 1

原创 GIS和R语言教学

数据时代来临，各个学科研究面临的数据处理需求越来越多。为满足广大硕士博士和科研人员地理数据处理分析需求，我们开设了这个地理数据科学培训班。为了满足高校教师和学生的学习需求，我们特别开设了高校教师1+2大礼包，只需900元（可开发票），即可获得一个SVIP卡+两个VIP卡，方便老师学生同时参加学习（详情请查看我们是如何收费的部分介绍，或咨询地理喵）。地理数据科学培训班全新升级，更好的学习体验，更丰富的课程体系，快来报名学习吧！手机端地址：https://geodatasci.duanshu.com

2021-02-01 16:17:54 1210

原创 Sentinel-5P数据介绍与处理

Sentinel-5P数据介绍与处理Sentinel-5P是欧空局于2017年10月13日发射的一颗全球大气污染监测卫星。卫星搭载了对流层观测仪（Tropospheric Monitoring Instrument，TROPOMI），可以有效的观测全球各地大气中痕量气体组分，包括NO2、O3、SO2、 HCHO、CH4和CO等重要的与人类活动密切相关的指标，加强了对气溶胶和云的观测。本文由三大主要部分组成：TROPOMI数据介绍数据下载数据处理与可视化Panoply可视化S5Processor包

2021-01-31 19:44:12 6862 10

原创 Sentinel-5P数据介绍与预处理

Sentinel-5P是欧空局于2017年10月13日发射的一颗全球大气污染监测卫星。卫星搭载了对流层观测仪（Tropospheric Monitoring Instrument，TROPOMI），可以有效的观测全球各地大气中痕量气体组分，包括NO2、O3、SO2、 HCHO、CH4和CO等重要的与人类活动密切相关的指标，加强了对气溶胶和云的观测。01TROPOMI数据介绍—TROPOMI是目前世界上技术最先进、空间分辨率最高的大气监测光谱仪。成像幅宽达2600km，...

2020-06-14 11:11:37 9215 11

原创 R语言经纬度度分秒转小数点

经纬度有两种常见的形式，度分秒形式和小数点形式。在GIS软件中，利用经纬度生成坐标点一般都需要将度分秒形式的经纬度转换为小数点形式，如何方便的进行转换呢？EXCEL确实可以，但是每次都要修改很多公式，而且经常不知道为啥报错，今天我终于用R语言写出来了一个小程序，可以比较方便的将EXCEL中的度分秒经纬度转化为小数点，在这里记录分享一下。EXCEL转换方式可以参阅往期文章：GIS基础教程之坐标系R语言安装部署请参阅：R语言安装部署基础下面是今天拿到的经纬度度分秒数据，表头E, N东经北纬，用.

2020-06-06 10:53:19 2136 1

原创 R语言安装部署基础

“R语言是一个开源的统计学程序环境，可以很方便的用于数学统计、机器学习、制图分析等一系列的应用。在这里讲一下如何进行R语言的安装部署。在这里介绍一下R语言和R studio的安装部署，为有兴趣学习R语言的同学迈出第一步。”01R语言环境搭建—首先下载安装R语言程序包。R语言的安装程序可以从R项目网站获取。网址如下：https://www.r-project.org/网站打开后如下图所示，目前最新的R语言版本是3.6.1，下载R程序包，需要...

2020-06-05 21:11:01 643