- 博客(131)
- 收藏
- 关注
RSS订阅
原创 探索Google Earth Engine：利用MODIS数据和R语言进行2000-2021年遥感生态指数（RSEI）的时空趋势分析
本文教会你在GEE上面实现RSEI计算和数据后处理
2024-11-14 22:25:31
568
原创 “傻瓜式”软件Panoply教程—一文搞定全球降雨量分布图与风场图
本期带你走进气象绘图软件Panoply的应用，一期带你搞定降雨量分布与风场图
2024-03-20 18:44:05
9697
2
原创 Grapher教程—重建长江中下游降雨量时间变化序列
我们发现拟合曲线的效果还是比较好的，把近150年降雨量变化趋势大致拟合出来了，从中我们发现，经历若干年，降雨量便存在峰和谷的交替，当然，此时的峰谷数量还是很多的。数据导入后，接下来我们就要生成时间序列图了，全选A、B两列，在标题栏找到“Home”并点击，找到新图表建立栏目“New Graph”，选择基础图形“Basic”，下拉箭头，选择第一种“Line Plot”，点击即可生成降雨量随着年份的线性时间序列变化图表。不错呦，接下来，我们想知道，近150年来，长江中下游降雨量的变化趋势到底怎么样？
2024-03-13 22:20:51
2977
4
原创 看这里！地理人必备工具箱：15个让你效率翻倍的神器网站
从坐标转换到三维建模，从遥感处理到数据爬取，这份清单将彻底革新你的工作流。今天的分享就到这里，很久没有好好写文章，大家有想看的内容也可以私信梧桐君哦。三维神器：一键下载任意区域DEM数据生成STL模型，也就是说不需要自己制作，直接下载就行。特色：支持WGS84/GCJ02/BD09等加密坐标系互转，一键批量转换CSV文件。隐藏福利：包含罕见的文化要素数据（宗教分布、语言分区）
2025-03-25 16:11:37
906
原创 解惑Python：一文解决osgeo库安装失败问题
总结一下，要安装osgeo模块，关键在GDAL依赖库的安装，如果能够顺利安装GDAL依赖库，因其是osego的模块，因此osgeo的部分模块也能够得到自动安装，此时再运行代码，发现代码并不报错，能够正常运行了。第一步，我们要安装imagecodecs库，它支持多种常见的图像格式，当然也包含我们的目标tif格式，为方便后续的图像处理和安装osgeo模块，先安装此库。目前GDAL已经更新到了3.10版本，这里主要看你的Python版本是上面，小编用的是Python3.13，电脑是64位的，因此，我们选择。
2025-02-17 23:44:37
1977
2
原创 寻找长春的酒吧——获取任意城市的POI（AOI），不用再为收集不到POI数据发愁……
来到新的城市游玩，第一件事找喝酒的地方。在这个被冰雪覆盖的季节里，长春的每一条街道，每一座建筑，都似乎在诉说着它独特的故事。城市的灯火阑珊处，隐藏着无数的温馨与感动……R语言中的osmdata 包提供有关全球各种空间属性和对象的空间数据。我是加拿大一枝黄花，我们下期见。我们今天想获取的是长春市的酒吧（bar），它属于"amenity"类别。现在长春市的边界数据拿到了，这个范围的所有酒吧别想跑了。我们在QGIS中打开看一下数据是否正常，看起来不错，名字都挺不错的。第一步需要做的是获取长春市的经纬度边界信息。
2024-12-31 17:11:29
268
原创 使用R语言绘制交互地图
在现代地理信息系统（GIS）应用中，交互地图成为了数据展示的重要工具。相比传统的静态地图，交互地图不仅能够更生动地呈现空间数据，还能增强用户的参与感和数据探索性。本文将介绍如何使用R语言绘制交互地图，并展示如何利用R的leaflet包创建具有交互功能的地图。随着地理信息数据的增多，传统的静态地图无法满足复杂数据展示的需求。交互地图允许用户通过操作地图实现数据的动态查看和分析，提供了更多的自定义选项和功能，使得地图更加具有探索性和应用价值。
2024-12-31 17:05:30
822
原创 使用R语言绘制标准的中国地图和世界地图
在日常的学习和生活中，有时我们常常需要制作带有国界线的地图。这个时候绘制标准的国家地图就显得很重要。目前国家标准地图服务系统向全社会公布的标准中国地图数据，是最权威的地图数据。总的来说，这个R包对于绘制带有国界线的标准地图非常方便，大家可以自己去尝试一下。以上就是今天的内容，如果觉得有用，请多多支持转发！今天介绍的R包“ggmapcn”，就是基于最新公布的地图数据制作的（审图号为：GS（2024）0650号）ggmapcn是基于ggplot框架开发的，语法简洁方便，用来绘制中国地图很方便。
2024-12-31 16:59:49
1377
原创 看这里！各类地理数据免费下载教程与链接
这个网站是我认为下载DEM高程数据最方便的网站，自从发现它后，再也没有去其他地方下载过高程数据，选择你感兴趣地区，直接点击就可以下载。在各种研究和科学问题中，我们需要各种各样的数据，虽然现在有很多朋友会分享数据，但是有的时候学会自己区官网找数据下载也是我们必备的技能。该网站也同样涵盖了大量数据，种类很多，尺度也很大，部分数据可以免费下载，像气温，降水之类的，总的来说还是比较方便。这个站可以下载全球任意地区的人口数据，也可以自由选择年份，非常方便，界面非常好看帅气，而且可免费下载。
2024-12-31 16:50:06
1332
原创 GEE必须会教程——一文教会你下载火灾信息数据
这几年频频登上热搜的澳洲山火，加州山火以及国内四川地区的山火等等，造成的较大生态系统损失和生命健康安全，令人惋惜。今天我们要用到的火灾信息数据是NASA提供的“FIRMS: Fire Information for Resource Management System——资源管理系统中的火灾信息”，该数据集包含栅格化形式的LANCE火灾探测产品。2019年，澳大利亚发生震惊世界的山火，据报道，自9月6起，该山火燃烧时间长达7个月之久，给当地的景观生态造成了巨大的破坏。如果出现平局，则保留较亮的那一个。
2024-12-04 13:54:01
688
原创 GEE Download Data——10mESA土地覆盖数据
高精度的土地覆盖数据类型不多，欧洲航天局（ESA）的WorldCover 10 m 2020产品基于Sentinel-1和Sentinel-2数据，以10 m分辨率提供了2020年的全球土地覆盖图。WorldCover产品包含11个土地覆盖类别，是在欧空局WorldCover项目的框架内生成的，该项目是欧洲航天局第五个地球观测包络计划（EOEP-5）的一部分。GEE下载数据第三弹，今天我们来看看土地覆盖数据的下载。了解完了数据之后，我们来看看怎么下载数据。
2024-12-04 13:00:02
229
原创 GEE Download Data——气温数据的下载
气温数据我们要用到的是MODIS数据产品，MOD11A2 V6.1 产品提供 1200 x 1200 公里网格内 8 天平均陆地表面温度 (LST)。MOD11A2 中的每个像素值都是该 8 天内收集的所有相应 MOD11A1 LST 像素的简单平均值。选择 8 天的合成周期是因为该周期是 Terra 和 Aqua 平台的精确地面轨迹重复周期的两倍。这个数据产品的波段比较多，有白天的气温，也有晚上的气温，还有一些其他的数据波段。和降水数据一样，我们也是按照最大值合成的逐年气温数据。GEE数据下载第二弹！
2024-12-03 17:13:56
367
原创 GEE Download Data——降雨数据的下载
该数据的空间分辨率为5km，在时间范围上，为1981至今，当我们研究区域较大时，选用该数据不失为一种可行的选择。在进行长时间序列分析时，我们往往需要的年度数据，此处我们选择最大值合成，合成2013-2023年的年度数据进行导出。感兴趣的朋友们也可以选择中值合成，均值合成等方式进行数据的合成，对比一下，取得最后的结果。上述代码定义了研究的结果如上，我们点击RUN ALL便可以下载合成的年度降水数据了。降雨数据想必大家在研究中用的比较多，今天我们在GEE上面看看如何实现降雨数据的使用。
2024-12-03 17:00:57
267
原创 使用R语言优雅的获取任意区域的POI，道路，河流等数据
我们可以开心的可视化了：osm 数据支持“点”，“线”，“”面多种的数据类型，像医院我们就可以有面数据和点数据。今天使用R语言调用的就是OSM的各类数据，R语言中的osmdata 包提供有关全球各种空间属性和对象的空间数据。是地图上任何非地理意义的有意义的点，如商店，酒吧，加油站，医院，车站等。我们今天想使用的是某地区的医院，道路，河流等数据，它们分别属于"amenity"，"highway"，"waterway"。目前这些数据的获取都有一定的难度，有的下载需要付费，使用代码爬取又需要编写复杂的代码。
2024-12-02 17:10:47
911
原创 使用R语言进行美国失业率时空分析（包括绘图）
数据当中"UE60", "UE70", "UE80", "UE90"字段分别表示1960，1970，1980，1990年的美国各区县的失业率。这次的示例数据非常简单，只有一个shp格式的美国区县矢量数据，我们在QGIS中打开数据查看一下它的属性表。我们还可以进一步的做分析，通过计算四个年份的失业率的均值，对失业率的数据进行空间自相关分析。1、引入必要的包，将本地的数据导入，查看数据是否正确。4、热点分析，进行可视化。
2024-11-29 00:10:43
672
原创 使用R的数据包快速获取、调用各种地理数据
rnaturalearth 包使 Natural Earth 数据可用。自然地球特征包括 1：10m、1：50m 和 1：1.1 亿比例尺的海岸线、河流、测深、政治边界、公路和铁路。数据一直是科学研究绕不开的话题，为了方便快捷的获取各种地理数据，许多R包被开发出来，今天介绍一些方便快捷的数据R包。ggOceanMaps，顾名思义，专注于海洋地图，包括海岸线、测深和冰川。ggOceanMaps 包包含一些数据。就这些吧，还有其他的数据，比如OSM，下次再写。
2024-11-29 00:04:10
829
原创 探索天空中的“名字”——用Landsat影像记录你的名字形状！
这些影像不仅能帮助我们分析地球的变化，还能让我们像“天文学家”一样，用卫星图像探索地球的奇妙景观。今天就来试试这个有趣的工具，看看卫星影像中的地球会为你的名字呈现怎样的“形状”。今天我发现了一个特别有趣的工具——NASA官网上有一个功能，允许你输入自己的名字，然后它会根据Landsat卫星影像显示出与你名字形状相符的地形图。神奇的事就发生了——你会看到一个带有你名字形状的Landsat影像！这些影像代表了你的名字对应的地理区域，也许是自然的山脉、河流、森林，甚至是城市的轮廓。wtgis(梧桐GIS)
2024-11-29 00:00:35
1045
原创 使用R语言绘制简单地图的教程
这个组图的图例出现了一些问题，是由于图例太窄放不下文字造成的，另外我们还需要添加一些地图要素，如比例尺，指北针等等。通过上面的代码，事实上我们分别可以得到三个不同版本的地图，分别为“只显示填充”（左），“只显示边界”（中），“边界和填充都显示”（右）tmap可以提供“tmap”类的对象，我们可以把绘制好的地图暂存在tmap对象中，通过制作多个tmap对象来达到绘制组图的目的。今天主要讲的部分是绘制静态地图，使用的R语言绘图包是tmap，关于介绍就不多讲，下面开始代码的讲解，小白也可以放心食用。
2024-11-26 11:14:21
699
原创 GIS空间分析案例---城市公共设施配置与服务评价
然后，为道路要素添加一个表达实际通行时间成本的Time字段（分钟），计算公式为Time= Length÷Speed。考虑某一种公共服务设施（如购物中心、药店、银行网点、派出所、加油站、酒店等）对住宅区人口的服务性质，通过OD成本矩阵分析与最近设施分析功能来评价设施点服务范围的可达性问题，具体内容包括：1) 不同设施点至住宅区的平均成本的分布与比较分析；采用与OD成本矩阵一样的加载位置的方法，将星级酒店和筛选出来的居民点设置为设施点和事故点。使用完全相同的方法，求出植物园和住宅区的OD成本矩阵。
2024-11-14 11:36:51
829
1
原创 R语言实战——一些批量对地理数据进行操作的方法
通过上述代码，我们可以实现对数据的格式转换，中间我们用if函数检查了一下我们的文件夹里面是否存在我们需要转换的数据格式文件，如果没有会报出，这时我们需要对文件夹重新整理。小编在处理的时候，一开始有info文件夹在里头，R语言读不了，后面吧info文件删除之后，就能够开始读取了。在进行数据处理时，tif格式的数据是常用的数据格式，我们从一些网站下面下载的数据可能是GRID格式或者其他格式，小编在处理多年风速数据的时候就遇到了这个问题，如果你需要将多张其他格式的数据转换成tif格式，应该怎么做？
2024-11-05 11:01:42
827
原创 GIS空间分析案例——城市用地适宜性评价分析
"省道距离分析" * 0.14 + "县道距离分析" * 0.14 + "河流距离分析" * 0.0485 + "溪流距离分析" * 0.0485 + "工业区距离分析" * 0.058 + "林区距离分析" * 0.076 + "城镇距离分析" * 0.09 + "村庄距离分析" * 0.09 + "DEM" * 0.155 + "坡度" * 0.155。本案例的思路是：①确定城市用地适应性的评价因子，本案例有交通便捷性、滨水环境、森林环境、城市氛围、地形等五个方面七个指标；
2024-10-18 14:10:13
1512
1
原创 GIS小知识——地理加权回归和地理探测器
由于该方法研究的是空间分异性，因此自变量的要求是类型值（如1，2，3，4，5 五种类别），而不是数量值（如0.1，0.2，0.3，0.4，0.5，...）。但是，GD依赖于空间分异性的假设，而在实际问题中，地理现象的空间分异性往往不是那么明显，这就导致地理探测器得到的结论有的时候并不完全可靠。对于GWR,不管空间分异性是否明显，其都能得到很好的应用。所以在实际研究当中，我们可以将两者结合起来，先利用GD对可能的影响因子进行分异及因子探测，筛选出显著性较高的因子，再利用GWR进一步分析各因子的解释力。
2024-10-17 08:56:58
1907
原创 探索Google Earth Engine：白莲种植范围提取
近年来，随着乡村振兴战略的提出，我国的农业种植模式呈现出多元化的趋势。荷花盛开的7-8月份，荷田、森林、早稻的光谱特征相似且近似平行，但荷田的反射率高于其他两者。步骤四：将F1和F2分别二值化后得到F1’、F2’，与NDVI相乘，得到F3，F3值越大，表示为白莲的概率越大。不同类型的地貌和植被覆盖率不同的地貌往往表现出不同的NDVI值，因此经常使用NDVI提取植被。以上就是今日的分享，欢迎关注梧桐GIS，让我们一起学习，一起进步！计算F2，以提高荷花田和其他植被类型如森林或早稻之间的差异。
2024-10-14 19:13:12
309
原创 jupyter的安装
当我们的Anaconda安装好之后，我们进入软件，开始看看jupyter怎么安装。初次进入软件的朋友需要注册一个账号，使用邮箱即可，该软件会向你的软件发送一个验证码链接，输入验证码，验证完毕之后即可登录成功。为了顺利使用pythonAPI进行编程，听说，我们必须要安装一款强大的工具，它的大名是jupyter。更改文件夹目录，此处大家可以根据自己的兴趣和需求该便软件的安装目录，小编这里是改在F盘新建的Anaconda文件夹内部，主要，安装的文件夹需要保证是空的文件夹，不然会报错，完成后点击“Next”。
2024-10-12 17:01:08
1042
原创 ArcGIS教程——一文制作水墨画地图
吟诗作赋，泼墨山水对现在的小编来说，难度实在有点大，特别是绘画，目前只停留在画火柴人的阶段，虽然前几个月确实有所练习，但无奈毫无艺术细胞，只能作罢。1、将桂林市的DEM数据（共九幅栅格）和桂林市行政边界（shp文件）导入Arcgis-Pro，并加载到地图，可以看见，DEM数据完全覆盖了桂林市，但是分幅之间有明显的的界线，一看就是个九宫格，这是我们需要对数据进行镶嵌。雾凇沆砀，天与云与山与水，上下一白。但只要想做，办法总是有的，虽然没有承诺过什么，但我还是决定，用我的方式，“画”一幅桂林山水。
2024-10-11 10:17:19
491
原创 【B站教程】区位图数据分享
这里给大家分享一下区位图的数据链接，大家多多转发点赞收藏支持一下，数据链接放在这里了。链接: https://pan.baidu.com/s/1y4jq2k24jY4VlSMsw59oOw 提取码: ya5b 复制这段内容后打开百度网盘手机App，操作更方便哦.
2024-10-11 09:47:38
180
原创 探索 Google Earth Engine：连续干旱天数的计算与空间可视化
本文教会你在GEE上实现连续干旱天数的空间可视化，一篇很好的案例！
2024-10-10 23:57:04
1696
1
原创 GIS小技巧——一文教会你安装BlenderGIS插件
不过，原始的Blender是没有地图底图的，要结合三维影像制作三维地形图，需要结合其他的软件，比如ArcGIS Pro，小编之前也曾经介绍过如何结合ArcGIS和Blender开展3D地图绘制，那么是否能够不借助其他的软件，直接在Blender上使用地图呢？关闭文件夹和Blender软件，并再次打开，让我们按照此前的步骤，点击GIS内部的底图，成功加载，让我们将底图换成Esri的（谷歌图源用不了呜呜呜~），点击确定打开。可是啊，我们发现，它报错了，显示的是没有影像，说明我们并没有正确安装影像的ID链接。
2024-10-09 19:09:47
4440
12
原创 ArcGIS实战——一文教会你调整适合中国宝宝体质的标准地图投影参数
各位要出图的小伙伴有没有遇到一个问题，当我们的研究区域是中国时，使用常见的WGS1984投影，会让我们的中国地图看起来不够挺拔，投影使得左右的拉伸比较长，进而让我们的中国地图变得不够美观。一切都完成之后，让我们动动小手，点击确定应用，看看投影之后的效果。
2024-10-07 15:20:07
1356
原创 一文讲清投影，定义投影，投影栅格……
这也就意味着你可以将数据从地理坐标系转换成投影坐标系，也可以从投影坐标系转换成地理坐标系，当然也可以进行相同类型坐标系的转换。这样子显然是不对的。最后来归纳一下定义投影工具和投影工具的区别吧：Ⅰ.当数据只有地理坐标系时，如果我们想把它转换成投影坐标系，可直接使用投影工具，或者可以先用定义投影工具转换成其他的地理坐标系，再用投影工具转换成对应的投影坐标系；Ⅱ.当数据只有投影坐标系时，如果我们要把它转换成其他的投影坐标系，定义投影工具和投影工具都可以实现，但是如果我们要转换成地理坐标系，就只能用投影工具了。
2024-09-10 11:02:17
2192
原创 矢量数据与栅格数据详解-01
这就是今天的大部分内容，事实上一些细节上的事情并没有详细说明，大家可以详细看看栅格数据的压缩编码的方式。今天就到这里结束了。矢量数据和栅格数据都是用来描述实际物理世界的方式，描述世界的方式和我们认知世界的方式是一致的。矢量数据和栅格数据是我们地理人最常见的数据类型之一了，也是地信原理最早介绍的内容之一，就让我们从他们开始。3、注意：通常拓扑关系越多，数据编辑和维护的难度就越大，越复杂，但要进行复杂的空间分析就方便。矢量数据是一种常见的数据结构，常被用于表示空间实体的空间位置以及附加的属性信息和几何关系。
2024-09-10 10:58:23
1348
原创 利用高德+ArcGIS优雅获取任何感兴趣的矢量边界
打开GIS免费神器QGIS，为了实现坐标的转换，我们需要安装一个小插件，在窗口里找到工具（Plugins），下拉找到管理并安装工具（Manang and install plugins），搜索“GeoHey”，点击安装插件即可。今天的分享到这里就结束了，本篇教程，禾穗联合加拿大一枝黄花，倾力奉上，有了它，大家想要哪个地方的shp，就可以轻轻松松解决了，还不快动起手来试试吧。我们的数据源来自高德地图，高德采用的坐标系是火星坐标系，在实际的使用过程中，我们还是要转成WCS1984坐标，才方便后续的使用。
2024-09-06 09:04:12
1396
3
原创 GEE必须会教程——时间序列影像填补方法
第5行，替换无效像素。从2024年11月13日起，如果在GEE上没有云端项目，将限制使用GEE，这个GEE的新闻头条部分读者已经在其他的文章或者是GEE开发官网得到消息了，还没有云端项目的伙伴，赶紧找到网址创建好属于自己的新项目，从项目进入即可！“七”待已久，“夕”望是你，我是禾穗，祝正在学代码的朋友们早日遇见幸福噢~，不过单身的朋友也没关系啦，看看日历，今天没有早八，没有工作，是快乐周六呀！第3-4行，计算当前影像前后两个月的日期，过滤期间所有的影像，之后，计算其平均值，利用mean方法生成平均影像。
2024-08-10 11:35:27
1183
原创 GEE必须会教程——基于Landsat影像构建NDVI时间序列
很久很久以前，小编写了一篇基于MODIS影像构建归一化植被指数的文章，不知道大家还有没有印象，有一段时间没有更新时间序列分析相关的文章了。今天，我们来看看基于Lansat影像，我们来看看在GEE上如何构建NDVI的时间序列。今天写起来有点手痛，小编就只能水一篇代码了。
2024-08-03 12:57:40
542
原创 GEE必须会教程——Savitzky-Golay（S-G）滤波器在MODIS影像时间序列数据的应用
今天比较懒，不想多写文字啦，还记得昨天分享的S-G滤波器，大家有没有去试试。今天布置个任务，同样以昨天的研究区域为例，试试看MODIS时间序列数据如何实现平滑！好了，今天的分享到这里就结束了，更多内容欢迎持续关注小编的公众号“梧桐GIS”，我是禾穗，咱们下期再会啦！结果还比较赏心悦目，好啦好啦，开启元气满满新一天！话不多说，今天直接上代码吧。
2024-07-22 23:16:33
570
原创 GEE必须会教程——Savitzky-Golay（S-G）滤波器在Sentinel影像时间序列数据的应用
此处创建一个包含n天（这里设置为一周时间）的时间序列数据，源于在处理哨兵-2卫星的遥感数据时，由于其两个传感器（Sentinel-2A 和 Sentinel-2B）每天可能会产生多景影像，所以我们需要对一定时间范围内的NDVI计算结果取平均。好了，经过上面的一通操作，我们得到了NDVI数据的S-G滤波结果，基本保留了原时间序列数据的形状、趋势，将高低频噪声进行有效平滑后，我们更能看出NDVI时间序列的变化趋势和大致的周期，是不是很有用呢。，因此S-G滤波器同样也适用于遥感影像的时间序列数据分析。
2024-07-22 13:01:55
1665
2
Sen+MK趋势分析代码
2024-11-28
matlab实现CV变异系数代码
2024-11-28
生态环境评价：基于MODIS数据和R语言的2000-2021年RSEI时空趋势分析实践
2024-11-28
利用Earth Engine API进行西南中国地区植被覆盖率计算与分类
2024-09-20
树轮重建的长江中下游地区夏季降水数据
2024-09-19
气象绘图软件-Panoply软件，快看过来
2024-09-19
“零基础轻松入门：小编精心整理的GEE学习笔记”
2024-02-19
武夷山国家森林公园矢量边界
2024-02-17
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人