走天涯徐小洋地理数据科学

其他

GSOD全球逐日气象站点数据介绍与下载

temperature最低温度,单位Fahrenheit华氏度MIN_ATTRIBUTES空值,自动观测到的最低温度*星号,从逐小时数据中提取得到的最低温度PRCP,Precipitation
2022年10月12日
其他

R语言+QGIS处理NOAA全球气象数据(二)气象站点数据下载

ISD数据的站点获取,接下来下载这个数据。在rnoaa包中,使用isd函数即可下载ISD数据,要指定三个参数:usaf,站点号,根据自己需要来选择wban,
2022年6月26日
其他

R语言+QGIS处理NOAA全球气象数据(一)气象站点信息提取

气象站点数据是地理相关研究非常重要的数据,由于国内数据共享程度低,获取困难,那么,没办法了吗?气象数据共享靠美帝啊。全球逐小时气象数据发现一个很方便的R包,可以获取NOAA数据,不光是这个ISD气象数据,其他也可以。R语言获取ISD站点信息R语言有一个rnoaa包,很神奇,可以直接下载NOAA数据。安装好这个包之后,加载包,就可以使用了。搜索中国及周边区域气象站为了方便查找中国区域的NOAA数据,我将中心点设置为经纬度103,36,搜索半径4000公里,搜索气象站点。代码如下:library(tidyverse)library(rnoaa)library(rstatatools)devtools::install_local("rstatatools.zip")#获取气象站x
2022年6月16日
其他

GIS基础之GIS属性表操作

GIS属性表操作对于数据量大的情况,QGIS由于是64位软件,且可以选择文本数据编码,对于大数据量的CSV数据处理我建议优先选择QGIS。下面以中国地面气候资料日值数据集(V3.0)逐日气温数据为例,介绍一下GIS属性数据的处理。下面以QGIS和ArcGIS这两种最常用的GIS软件平台分别进行介绍。文本数据导入数据为文本型,用空格分隔,根据数据文档我们得知13个字段分别为:区站号纬度(度、分)经度(度、分)观测场海拔高度(0.1m)年月日平均气温(0.1℃)日最高气温日最低气温平均气温质量控制码日最高气温质量控制码日最低气温质量控制码导入文本数据QGIS属性计算QGIS字段计算器使用的是Python语言,在这里我们需要使用python的函数进行计算。我们以下面几个例子进行介绍:经纬度计算原始数据中经纬度是以度、分的形式存储的,后两位是分,前面是度。由于中国疆域广阔,东西跨度大,所以经度的度有2位,也有3位,这个时候就不适合通过位数来进行分割。那么应该用什么方法计算呢?我是使用了整除100的方式进行划分:度分整除100,商整数部分作为经纬度的度度分整除100,余数即为分,由于GIS经纬度需要为小数点,余数处理60化为小数制度+余数/60即为小数的经纬度具体计算过程和代码如下,注意一下这个数据纬度字段在前(field2),经度字段在后(field3):经纬度度、分计算建立一个新的字段,双精度格式,输入python函数计算
2022年4月30日
其他

ArcGIS里面如何绘制箭头和矩形框?

有同学提了一个问题,“arcgis出区位图能直接在里面加→和矩形框吗?”当然是可以的。有两种方法:绘制矢量数据,建一个线要素和面要素,线要素做箭头,面要素做矩形框绘制图形,修改形状颜色和样式要想得到想要的样式,都是两步,创建图形,添加可视化效果(符号化)。第一种方法,矢量数据的创建和编辑前面讲过很多了,可以看以前的推文,不再多说,下面讲讲第二种。ArcGIS图形编辑命令在GIS中如何绘制胡焕庸线地图制图之ArcGIS制图技巧(鹰眼图、数据框联动、统一符号化效果)ArcGIS绘制图形矩形框绘制矩形框,本质上就是一个填充为透明的矩形,理解了这个,矩形框就很好画了。使用Drawing工具条上的矩形工具,画一个矩形矩形上面右击,Properties,调出属性面板,调整Fill
其他

【更新189篇】地理·生态·遥感科研必备技能合集

公众号技术文章合集目录,目前已整理以下8个部分,共计189篇文章GIS基础GIS制图与空间分析地理数据分享与处理方法景观生态R语言/Python地理数据综合处理遥感与影像信息提取河口海岸各种软件工具和技巧总结文章附带练习数据、资料本来是开放获取的,被某树、某研修、某谷抄吐了,以后只给培训班会员提供了。等哪天无节操号们不搬运抄袭了我再继续公开共享数据,感谢大家的支持!推荐使用浏览器或者微信自带搜索功能进行页面搜索,快速定位问题!使用微信搜索功能,输入关键词快速定位问题点击蓝色标题即可跳转阅读对应文章1
其他

rgee尝鲜Landsat9,以及Landsat8对比

LT9mosaic$clip(ee_roi)Map$centerObject(ee_roi)Map$addLayer(LT9_ROI,
其他

能过审的论文插图获取方法

目前国内期刊投稿对地图管理的越来越严格,很多期刊要么要求有审图号,要么要求以带审图号的国家标准地图作为底图。以一个5A级景区全国分布图为例,介绍一下如何将专题要素弄到标准地图上的方法。标准地图下载既然要在标准地图上叠加专题要素,那么首先就要拿到标准地图。标准地图可以去标准地图服务系统下载。标准地图服务系统网址:http://bzdt.ch.mnr.gov.cn/index.html为了方便放置专题要素,建议选择比较简单的标准地图。在这里我选择了一个审图号为GS(2019)4344的中国主要河流湖泊分布图作为底图。选择EPS格式的矢量图下载。下载EPS矢量图专题要素绘制标准地图底图拿到了,我们接下来要绘制专题要素,在这里是一些经纬度的5A级景区数据。5A级景区展点使用Add
其他

含视频讲解 | rgee获取全国高程、坡度、坡向数据

高程、坡度、坡向是经常用的基础地理数据,要想获取全国的高精度高程、坡度、坡向数据往往需要花费大量的时间进行下载、镶嵌、投影转换和计算,在这给大家分享一个基于rgee的全国高程、坡度、坡向数据获取方法,能够大大加快速度。视频号直播讲解如下:还请大家多多关注一下这个没啥人的视频号吧!ALOS
其他

rgee:R语言+Google Earth Engine遥感大数据分析神器

前面已经给大家分享了不少rgee相关的内容,在这里再给大家系统总结一下,顺便分享一下我在学习rgee过程中的一些心得体会。为什么我认为R语言+Google
其他

菌菌:移动窗口法原理是什么?然后我给他找了一个栅格计算的原理教程

的值,这是一个全局操作。例如,欧几里得距离就是全局运算的一个例子。通过计算距源的最近距离,它在栅格中全局应用该函数。左边1为起始位置,右边为欧几里德距离值,sqrt(2)=1.4,
其他

1:100万基础地理数据库怎么看?怎么用?能干嘛?

上个月就开始写的推文,由于生娃一直没写完,前天终于出院了,把这篇“难产”的推文写完了,今天是大年初一,分享给大家。前面给大家讲了1:100万公众版基础地理信息数据(2021)的获取,接下来讲讲这个数据如何使用。如何看懂里面数据代表的含义呢?2021版国家基础地理数据库批量下载与拼接方法相关标准规范要想看懂里面数据的含义,必须了解1:100万公众版基础地理信息数据(2021)是根据什么制作的,本数据涉及的部分标准规范如下:GB/T
其他

Google Earth Engine大数据应用:荟萃分析和系统性综述

Geoinformation(JAG)》是GEE应用文章发文量最高的3个期刊。GEE相关文章>5篇的期刊研究领域349篇论文经过汇总整理划分为11个研究领域:作物制图(Crop
自由知乎 自由微博
其他

Sentinel-5P数据介绍与处理

Sentinel-5P数据介绍与处理Sentinel-5P是欧空局于2017年10月13日发射的一颗全球大气污染监测卫星。卫星搭载了对流层观测仪(Tropospheric
其他

在GIS中如何绘制胡焕庸线

805–824.专题要素到标准地图的整合——以中国5A级景区分布图为例R语言降水量数据处理QGIS绘制等降水量线图全国1:100W地理数据库的制作与合成
其他

2021版国家基础地理数据库批量下载与拼接方法

2021版国家基础地理数据库批量下载与拼接方法最近全国地理信息资源目录服务系统又放了个大招,公开了2021版全国基础地理数据库。本套数据现势性2019年,快去下载看看吧!数据获取网址:https://www.webmap.cn/main.do?method=index数据获取网址从网站上检索数据,然后全选,加入购物车,注意需要账号注册,根据指引操作即可。全选加入购物车数据成果组织数据使用前一定要看明白数据说明,另外对基础地理数据库里面的字段和代码不了解的一定要去查询相应的国标,可以看:【资源分享】如何查找靠谱的国标,全文免费看!全文免费看!全文免费看!全国1:100万公众版地形数据(2021)内容包括水系、居民地及设施、交通、管线、境界与政区、地貌与土质、植被、地名及注记9个数据集。数据分层的命名采用四个字符,第一个字符代表数据分类,第二三个字符是数据内容的缩写,第四个字符代表几何类型。图层名称和类型介绍数据批量下载默认是一个个下载,挨个点77下好麻烦啊。有没有批量下载的方法呢?网页上提供了77个下载按钮可以使用Google
其他

文献计量三大定律与Bibliometrix

Environment(RSE)》https://www.sciencedirect.com/journal/remote-sensing-of-environment《Remote
其他

Landsat Collection2 Surface Reflectance产品介绍、数据获取、去云

日获取通过考虑气溶胶散射和薄云等大气效应,地表反射率数据提高了同一区域多期影像间的可比性,这有助于检测和表征地球表面变化。地表反射率数据由满足
其他

全国分省、市、县NPP面板数据是如何制作的

ee$ImageCollection("MODIS/006/MOD17A3HGF")$select("Npp")#数据预览npp_composite
其他

2022新年快乐!迟来的2021年总结,R语言分析公众号文章阅读

2022新年快乐!迟来的2021年总结,R语言分析公众号文章阅读最近实在是太忙了,都没怎么更新原创文章。2021公众号总结我也来看看这一年公众号的运营成果。用户关注情况从2021年1月1日的10471个用户,到2021年12月31日我写这个文章的时候的28671的关注用户,平均每天增长50个用户。本来想用R统计一下,结果发现公众号后台导出的EXCEL有问题,直接打开会报错,没办法,EXCEL里面手动粘贴一下,结果如下:累计用户增长从性别看,男生较多,女生略少。男女比例从年龄看,绝大多数粉丝年龄低于35岁,也有一部分忘年交。青年居多从地区分布看,北京高居榜首,其次广东、江苏、湖北,陕西。粉丝已经遍布全国,有133个香港同胞,11个澳门同胞,还有4个台湾同胞。北京高居榜首感谢近5K铁粉的支持,还希望大家能给我加个星标,第一时间收到我的推送,如果可以的话还请多多转发、分享支持一下!常读用户占比18%单篇阅读接下来分析一下单篇阅读量library(tidyverse)library(readxl)userdata
其他

【数据分享】2001-2020年中国分省市县净初级生产力(NPP)面板数据工具变量

【数据分享】2001-2020年中国分省市县净初级生产力(NPP)面板数据工具变量公众号后台很多同学留言NPP,在这给大家分享一套2001-2020年中国分省市县净初级生产力(NPP)面板数据。数据介绍使用MOD17A3HGF.006数据,根据2020版中国行政区划数据分省市县进行空间统计,统计方法为均值。空间分辨率:500m时间分辨率:年统计方法:均值时间范围:2001-2020数据包含:NPPSheng.csv
其他

QGIS和R语言栅格数据重分类

QGIS和R语言栅格数据重分类重分类(Reclassify)是栅格数据处理中常用的操作,可以用于将土地覆被地类合并,或者将连续的像元值合并为离散的值(NDVI时间序列分析之Sen+MK分析全过程梳理)等等。常用的GIS软件ArcGIS、QGIS都可以实现栅格数据的重分类。ArcGIS重分类网上已经很多了,不再赘述,这里主要讲讲QGIS和R语言里面的重分类。相比ArcGIS有什么特点呢?QGIS优势QGIS相比ArcGIS有这更高的栅格处理效率QGIS免费QGIS栅格显示不出错,多波段的栅格ArcGIS显示经常出问题栅格数据的一些概念和常用处理方法,ArcGIS栅格数据显示错误问题如何解决R语言优势下面我们分别来看在QGIS和R语言中如何进行栅格数据的重分类。可以并行计算:对比R语言Raster包和Terra包栅格计算rgee可以调用Google
其他

​QGIS+R语言分析不同地区土地覆被占比变化

QGIS+R语言分析不同地区土地覆被变化完全使用开源软件给大家做一期土地覆被变化的实验!实验目的研究多个地区不同年份的土地覆被类型占比变化情况。结果图如下:1990-2020年京津冀土地覆被占比变化实验数据4期土地覆被数据分区矢量数据,在这里使用了京津冀分市的行政区矢量土地覆被重分类原始的土地覆被类别比较多,二级类,先合并为一级类。土地覆被对应代码表一般栅格数据重分类可以使用重分类工具来做,但是这个土地覆被二级类转一级类如果用普通的重分类工具来做这个重分类要用的重分类对照表格实在是有些冗长。那么有没有什么优雅简便的方法呢?可以使用栅格计算的方法!通过观察我们发现,一级类的代码也就是栅格值,其实是二级类除以10取整得到的。由于QGIS栅格计算器可用的函数太少了,我决定用R语言,terra包并行计算来搞定它!代码如下:library(terra)L2C2020
其他

使用rgee和QGIS合成非洲ALOS DSM

DSM(数字地表模型)数据GEE下载方法我选择了这种方法,不过对于这种洲一级别的30米分辨率数据来说,GEE也不是很方便,主要是下载过程容易中断,也需要下载完之后手动拼接数据。GEE下载ALOS
其他

一些GISer必备的电脑知识

一些GISer必备的电脑知识最近见识了一些对电脑基础知识严重匮乏的同学,作为一个GISer,经常要进行数据处理和编程,怎么能没有电脑常识呢?下面总结一些常用必备技能。对于计算机同学来说我这可能讲的不一定准确,就是我个人的一些理解,但是用在GIS和日常电脑操作上应该是没问题的。Windows系统相关32位与64位计算机是二进制的,目前主流的电脑系统是64位系统,可以使用4GB以上内存,而32位的系统和软件只能利用4GB以下的内存空间。64位的电脑可以使用4G以上的内存从任务管理器可以查看软件是32位还是64位32位的软件在64位任务管理器中,进程名称后面会有一个括号(32位)QGIS是64位软件ArcMap32位软件由于ArcMap是32位的,这就决定了ArcMap无法充分利用4GB以上内存,限制了电脑性能的发挥,所以,处理栅格数据,推荐QGIS32位还是64位的软件电脑任务管理器里面能看出来什么?任务管理器是什么???任务管理器任务管理器,是windows电脑管理进程的工具,简单来说,就是电脑里一个进程就是你脑子里的“思路”,当你思路混乱的时候,就是类似电脑进程太多了,内存要爆掉了,或者你的大脑宕机了,就是出现了任务管理器里面的进程“未响应”。如何查看任务管理器呢?最简单的方法,在开始菜单栏上右击,任务管理器在开始菜单栏上右击,任务管理器当电脑程序卡死,就可以从任务管理器中结束进程,进程结束后,程序停止执行,需要注意的是,如果没保存,你程序正在执行的东西会丢失。结束任务,强制把某个进程断掉在R语言并行计算中,往往需要说明核心数是多少,使用任务管理器就可以查看电脑的核心数。核心数指的是内核的数量!核心说是内核的数量,4核后缀名前面fragstats景观格局指数批量计算当中涉及过后缀名的修改,有人居然不会,我在这再讲一次,这应该是初中或者小学计算机课就上过的啊以最常用的win10系统为例:GlobalBedrock.tif
其他

栅格数据分区统计转面板数据,单波段、多波段数据都可以找到合适的方法

Statistics)可以用于统计矢量区域范围内的像元值,对该区域内的所有像元进行统计,计算本区域像元的数量和(和上面的Zonal
其他

地理学数学方法SPSS与R语言应用

2)单击选项,选中均值、标准差、最大最小值、方差。再单击继续,确定。如下图所示:结果:描述统计量N极小值极大值均值标准差方差销售额40254937.855.95935.515有效的
其他

科研软件工具汇总

公众号中提供了一些常用软件下载方式,关键词回复不太方便,整理成一个文档吧。下载这些软件工具之前,建议先看看这个:【工具分享】如何快速下载网络资源1
其他

随机森林R语言回归学习笔记和一个失败的试验记录

2001)。可以处理小n大p问题,高阶相互作用,相关的预测变量等。随机森林可以进行分类或回归分析,得到变量的重要性、模型的误差等指标,并可以进行预测。随机森林算法发展Breiman
其他

如何消除栅格镶嵌后的零星NODATA像元

如何消除栅格镶嵌后的零星NODATA像元一套分省的栅格数据,拼接后发现有零星的NODATA像元。白色部分为NODATA的像元怎么办呢?ArcGIS消除NODATA像元从网上搜索了一下,有人说ArcGIS可以通过滑动窗口的方法消除NODATA像元,本次试验数据为土地覆被数据,对此我使用了“MAJORITY”的方法进行插值消除NODATA。具体代码如下:Con(IsNull("L2_2020.tif"),
其他

R语言月均温栅格转年均温栅格

R语言月均温转年均温前面给大家推荐过彭守璋老师发布的1km分辨率全国的月度气候数据。本套数据是逐月的,对于年际变化研究来说使用不太方便,需要把月度数据合称为年度数据,在这举个栗子,教大家如何进行年度数据的合成。【数据分享】中国1km分辨率逐月降水、最低最高平均气温数据1901-2020上图数据下载地址:http://www.geodata.cn/data/index.html?ownername=%E5%BD%AD%E5%AE%88%E7%92%8B月均温转年均温本套数据多数为3年一个nc文件,我采用的处理思路如下:使用terra包读取nc文件提取每年12个月的数据对12个月的数据进行合成,月均温转年均温采用求均值的方法,栅格计算app函数对每年的数据进行波段合成c导出处理好的数据
其他

​R语言获取和处理MCD12Q1土地覆被数据

R语言处理MCD12Q1土地覆被数据MCD12Q1数据介绍MCD12Q1是基于MODIS的一套年度土地覆被数据,其中的IGBP分类较为常用。目前MCD12Q1数据为2001年至2019年,每年一期。具体IGBP分类如下:IGBP分类1常绿针叶林,树冠高于2m,覆盖度>60%2常绿阔叶林3落叶针叶林4落叶阔叶林5混合林,针叶阔叶在40-60%之间,树冠高>2米,覆盖度>60%6茂盛的灌木,覆盖度>60%7稀疏的灌木,覆盖度10-60%8有林草原,具有草本和其它林下系统的土地,树覆盖度30-60%,树冠高大于2m9稀树草原,具有草本和其它林下系统的土地,树木覆盖度10-30%,树冠高大于2m10草地,草高10%12农田,60%以上耕地13城市和建成区,至少30%的不透水面,包含建筑材料、沥青路面、车辆等14农田/自然植被混合区域,较小块的耕地,有40-60%的部分是树木和草本植物15永久性积雪,60%以上区域,一年10个月以上被积雪覆盖16荒地,60%以上区域无植被覆盖,植被覆盖度低于10%17水域,60%以上的水域MCD12Q1数据获取常规方法从网上下载MODIS数据,然后使用HEG或者MRT等软件进行坐标系转换和拼接。具体可以看以前的两篇推文,不再多说,太麻烦了,反正我是不想用了,尤其是MODIS正弦投影的转换,太麻烦,HEG或者MRT软件都很难用。wget批量下载遥感影像MODIS
其他

免费!超大空间!史上最强文献管理工具Zotero+TeraCLOUD网盘多端同步!

al点击宏按钮创建一个宏,将下面代码复制到宏编辑页面中,然后运行宏,即可实现英文文献引用中的“等”自动替换为“et
其他

R语言栅格时间序列回归分析——整体和逐像元计算,并行计算,快到飞起!

R语言栅格时间序列回归分析——整体和逐像元计算,并行计算,快到飞起!前面给大家分享了2001-2020年1km分辨率中国降水量数据,有同学问如何做回归分析,接下来以一元线性回归举例,从整体上和逐像元两种方法来计算。使用的R语言程序包terra包,栅格计算,支持并行计算tidyverse包,R语言数据处理可视化必备,内含ggplot2ggpmisc包,作为ggplot2的补充,时间序列图中标注回归方程用整体回归为了减小点计算量,我把数据进行了裁剪,保留京津冀地区。降水量为年降水量,单位为0.1mm。京津冀降水分布,0.1mm对整个区域的计算和时间序列分析较为简单,只需对每个栅格进行统计,计算整体均值,生成时间序列,使用ggplot2包进行制图,还使用了ggpmisc包进行了统计结果的标注。library(terra)library(tidyverse)library(ggpmisc)pre_cn
其他

已知某点经纬度,如何快速查看该点周围卫星影像

已知某点经纬度,如何快速查看该点周围卫星影像如题,一个小技能,使用QGIS很容易实现。问:107°36′44.57028″,38°10′27.11964″
其他

专题要素到标准地图的整合——以中国5A级景区分布图为例

目前国内期刊投稿对地图管理的越来越严格,很多期刊要么要求有审图号,要么要求以带审图号的国家标准地图作为底图。以一个5A级景区全国分布图为例,介绍一下如何将专题要素弄到标准地图上的方法。标准地图下载既然要在标准地图上叠加专题要素,那么首先就要拿到标准地图。标准地图可以去标准地图服务系统下载。标准地图服务系统网址:http://bzdt.ch.mnr.gov.cn/index.html为了方便放置专题要素,建议选择比较简单的标准地图。在这里我选择了一个审图号为GS(2019)4344的中国主要河流湖泊分布图作为底图。选择EPS格式的矢量图下载。下载EPS矢量图专题要素绘制标准地图底图拿到了,我们接下来要绘制专题要素,在这里是一些经纬度的5A级景区数据。5A级景区展点使用Add
其他

R语言栅格时间序列回归分析——整体和逐像元计算,并行计算,快到飞起!

R语言栅格时间序列回归分析——整体和逐像元计算,并行计算,快到飞起!前面给大家分享了2001-2020年1km分辨率中国降水量数据,有同学问如何做回归分析,接下来以一元线性回归举例,从整体上和逐像元两种方法来计算。使用的R语言程序包terra包,栅格计算,支持并行计算tidyverse包,R语言数据处理可视化必备,内含ggplot2ggpmisc包,作为ggplot2的补充,时间序列图中标注回归方程用整体回归为了减小点计算量,我把数据进行了裁剪,保留京津冀地区。降水量为年降水量,单位为0.1mm。京津冀降水分布,0.1mm对整个区域的计算和时间序列分析较为简单,只需对每个栅格进行统计,计算整体均值,生成时间序列,使用ggplot2包进行制图,还使用了ggpmisc包进行了统计结果的标注。library(terra)library(tidyverse)library(ggpmisc)pre_cn
其他

GIS基础之矢量数据,概念、编辑、转换、空间分析,一文讲透!

GIS基础之矢量数据去年的基础课主要是从图形角度出发,从角度,距离、面积方面给大家进行的讲解,这次我决定换个视角,从数据角度给大家讲讲GIS的基础。矢量和栅格是最常用的两类数据格式,本次先给大家介绍矢量数据。1.
其他

R语言降水量数据处理QGIS绘制等降水量线图

R语言降水量数据处理QGIS绘制等降水量线图降水量数据是地理研究中经常用到的。下面介绍一下如何获取降水量数据并绘制等降水量线图。降水量数据获取和处理降水量数据获取前面给大家分享过多套降水量数据,具体可以进入公众号,点击发现宝藏,数据共享来查看。本次使用的是下面的数据:【数据共享】全球2.5分分辨率最低、最高气温,累积降水栅格数据集在geodata公开发布下载其它降水数据可以查阅:【数据分享】中国1km分辨率逐月降水、最低最高平均气温数据1901-2019如何找到靠谱的科研数据?免费的国家数据中心为你服务!降水量数据处理前面的数据是年累积降水量,为了进行气候分区我们需要进行多年平均,计算多年平均的降水量进而提取等降水量线。那么我们就需要给多年累积降水量做一个均值。在这里我使用了R语言terra包进行多年降水量的均值计算。开始计算前先把降水量数据放到同一个文件夹里library(terra)f
其他

R语言Raster包和Terra包栅格读写、计算和一些使用经验分享

R语言Raster包和Terra包栅格读写、计算和一些使用经验分享栅格读写Raster包和Terra包在读写方面语法略有差异。个人感觉Terra包的读写更方便一些。大多数时候我推荐使用Terra包,除非需要使用for循环或者借助其他包进行并行计算的时候使用Raster包,多数情况尽量使用Terra包进行栅格计算操作。Terra读取NC文件Terra包可以直接读取NetCDF文件大多数栅格文件都可以直接使用rast函数读取,Raster包读取的栅格也可以通过rast函数调用。vect读取矢量文件,可以用于裁剪栅格library(terra)JJJ
其他

GIS基础教程之坐标系(2021版)

GIS基础教程之坐标系(2021版)坐标系在我看来是GIS课程中最重要也是最难以掌握的点,结合以前讲过的GIS坐标系相关内容和同学们的反馈,我决定再次更新《GIS基础教程之坐标系》。个人精力有限,如有什么纰漏还请留言交流指正。本教程从以下几个方面入手:坐标系的基本概念地理坐标系投影坐标系如何选择坐标系根据研究区域大小根据研究目的(等角?等面积?等距离?其他)推荐一个在线坐标系选择网站GIS坐标系几种情况数据的坐标系GIS显示坐标系动态投影GIS软件中如何定义坐标系ArcGIS中坐标系的定义QGIS中坐标系的定义坐标系的转换矢量数据的转换栅格数据的转换本次更新主要是以思维导图的形式,给大家缕清楚各种坐标系和转换方法之间的关系,帮助大家建立起坐标系的知识体系。有了知识体系,才能以不变应万变。1.
其他

在GIS中如何绘制胡焕庸线

805–824.专题要素到标准地图的整合——以中国5A级景区分布图为例R语言降水量数据处理QGIS绘制等降水量线图全国1:100W地理数据库的制作与合成
其他

栅格数据处理答疑

栅格数据处理答疑针对上次课学员提出的栅格处理问题做一个答疑。学员甲提问从栅格中提取的数据写为excel之后是否可以写回栅格?两幅栅格影像之间怎么做运算,eg用降水与气温求饱和水汽压差?多波段栅格如何分层提取,多个单波段数据如何写为一个栅格?答:从栅格中提取的数据写为excel之后是否可以写回栅格?从栅格提取数据有以下几种情况:提取栅格值到点:在ArcGIS中可以使用Extract
其他

QGIS栅格图例设置详解——如何制作连续的和间断的图例

QGIS栅格图例设置详解——栅格如何统一色带?如何制作连续的和间断的图例?前面给大家介绍过QGIS绘制等降水量线图,今天在这个的基础上讲讲QGIS栅格的图例。QGIS3.18版本开始支持连续的色带图例,以前的老版本是不支持的,只能是间断的色块图例。在本文中给大家介绍一下如何给多波段或者多幅栅格应用统一的图例,以及连续或间断图例的设置。京津冀春夏秋冬降水分布,使用了统一的间断色块图例栅格统一色带要想做出上面的京津冀春夏秋冬降水分布图,有以下几个要点:栅格统一色带QGIS中分幅图制作,四个分幅图制作方法参考下面的推文:推荐一个野外考察神器!+QGIS野外考察路径制图在QGIS中色带统一很容易先根据值域设置好色带,以这个降水量数据为例,要选择夏季的,夏季降水量大,色带能覆盖全部值域如下图,保存色带改到其他季节(波段),应用前面保存的色带文件通过引用色带文件可以保持色带的统一保存和打开引用色带连续图例和间断图例QGIS3.18及以后的新版本是支持连续图例的,且默认的栅格渐变图例是连续图例。那么连续和间断图例在哪里切换呢?这个是在图层属性中控制的。Use
其他

R语言栅格处理、遥感应用总结,如何删除栅格处理中大量的垃圾文件

R语言栅格处理、遥感应用总结,如何删除栅格处理中大量的垃圾文件前面的推文中给大家介绍过不少R语言栅格计算相关内容,在这给大家做一个总结。另外今天有同学提问:想问一下,怎么可以清理R语言运行占用的C盘内存呀!感觉我用R语言计算栅格数据,就看着C盘占用内存越来越多。下面对R栅格计算相关推文做一个总结,还有答疑:R语言栅格计算总结点击下面蓝色标题链接,即可查看对应推文:优势R语言栅格计算优势:R语言是开源软件,免费R语言安装部署基础R语言支持并行计算对比R语言Raster包和Terra包栅格计算R语言栅格数据并行计算,HDF批量转TIFFR语言拥有大量的统计学算法,可以方便的调用R语言并行计算Sen+MK趋势分析R语言NDVI和降水量逐像元相关性分析BFAST分析MODIS时间序列数据NDVI时间序列分析原理与实现(CV和Sen+MK趋势分析)R语言NDVI变异系数计算与结果分析NDVI时间序列逐像元计算Hurst指数R语言GIMMS
其他

R语言栅格计算之月降水量栅格合成季节降水量

月降水量计算季节降水量以月度数据计算季节数据是很常用的操作,在这里介绍一下如何用R语言由月度栅格数据合成季节栅格数据。以降水量数据为例做个介绍。目的有月度数据,合成季节数据。根据气象学上对季节的划分,3–5月为春季,6–8月为夏季,9–11月为秋季,12月至翌年2月为冬季。操作过程数据选取根据气象学上对季节的划分,涉及两个年份,这个时候我们需要选出第一年3-12月的数据,还有第二年1、2月的数据。使用terra包进行栅格的操作rast读取栅格c挑选组成新的栅格writeRaster输出新的栅格数据library(terra)JJJ
其他

Fragstats景观指数计算常见问题

Fragstats常见问题解答常见问题总结没有选择投影坐标系,景观指数计算一定要投影坐标系文件目录必须全英文!不得有特殊符号!“-、*&……%¥#@汉字”这些统统不能有!Fragstats是32位软件,你就算是电脑有64GB内存它也顶多能用不到4GB,太大的数据计算不了!所以要用网格法切分栅格。csv、fbt是文件后缀名,别问我,去重修小学、中学、大学计算机,谢谢!景观指数不能一股脑全选中进行计算,不同的景观指数需要输入不同的参数,报错往往是你输入参数不够或者错误,具体去看详细帮助,Fragstats官网有帮助!官方详细帮助在这http://www.umass.edu/landeco/research/fragstats/documents/fragstats.help.4.2.pdf我觉得这篇文章可以封笔了。关于如何提问欢迎大家在交流群中讨论问题,我有时间的时候会给大家解答。但是请大家提问的时候,不要发送手机拍屏幕的照片,很难看清!提问不要用手机拍屏幕,一定要电脑微信截屏!电脑微信截屏!电脑微信截屏!最烦的就是拿照片来提问,看不清!看不清!看不清!常见问题Row
其他

GIS和R语言空间数据分析可视化技术

我们提供GIS+R语言的全套课程,涵盖论文数据处理、分析、可视化全过程的培训指导!课程体系涵盖基础课,即使是GIS和R语言的小白也可以无忧入门!SVIP会员尊享资深老师一对一答疑,不用担心学习过程中遇到问题无人指导!会员可以参与每周的直播课,所有直播课均可在平台上查看录播,直播课可以随时提问,现场答疑!地理数据科学培训班全新升级,更好的学习体验,更丰富的课程体系,快来报名学习吧!手机端地址:https://geodatasci.duanshu.com/#/PC端地址:https://geodatasci-pc.duanshu.com你可以从培训班学习到什么各种地理计算的
其他

使用rgee进行Climate Data Record NDVI月度最大值合成

NDVI月度最大值合成前面的推文给大家介绍了1982-2020年逐月NDVI数据产品的获取:【数据共享】1982–2020年中国5