- 博客(29)
- 资源 (1)
- 收藏
- 关注
RSS订阅原创 提取栅格影像与点shp相交部分的点，导出为一个新shp
如下图所示，栅格影像中有蓝色和白色两种不同的像素值（比如分别对应0和1），样点则散落在图中不同颜色上，现在我只想提取蓝色区域的样点并导出为一个新的点shp，怎么实现呢？导出一个和原始shp数量一致的shp文件，但其属性表中则增加了各个点对应栅格的数值。然后在属性表中按数值筛选即可，比如要筛选蓝色区域的点就仅提取数值==1的点即可。最后将筛选出的点export成为新数据即可。在 “Analysis” 工具箱中，找到。
2025-03-28 13:24:14
186
原创 关于解决User memory limit exceeded.的几点心得
比如你要下载大区域的多波段图像，每副图像内存以G计算，那就导出至Google drive再下载，这是图个稳定，不至于你下载到一半，报错User memory limit exceeded.一、引用吴秋生老师的话：下载/导出数据量大时，导出至Google drive再下载；下载/导出小数据时，直接导出到本地。一、引用吴秋生老师的话：下载/导出数据量大时，导出至Google drive再下载；但你就下个DEM，每幅影像就几M，那就用python api直接导出到本地，为的是速度，比导到drive快得多！
2025-03-20 16:41:29
358
原创 GEE批量打开asset权限（anyone can read）
我用不同的账号登陆GEE，需要实现asset共享，因此需要随时保持他们是anyone can read的权限。然而，我的asset数量很大，一个个点不现实，而且这个权限还是有效期的，过一阵还得重新点。将asset_id改为自己的id，在python api中运行下述代码即可。通常当asset个数少的时候，可以手动设置权限。所以我得找个办法一键解决。
2025-02-12 20:22:53
177
原创 解决geemap.extract_values_to_points 提取不同点的值时，结果相同
最近卡了一个问题好久好久，我用geemap.extract_values_to_points 函数提取样点对应数值（即feature values）训练模型，结果训练完才发现不同点的feature values大量相同。所以呐，对于新函数一定要认真读参数和对应意义！github的相关讨论区也要关注，今天这个问题就是我在阅读讨论区的时候发现作者的答复才解决的！具体链接如下：common module - geemap。
2025-02-12 15:53:48
183
原创 Google Colab 运行GEE Python 出现“PERMISSION_DENIED“问题
昨天听个讲座讲可以在Google Colab中执行GEE Python代码，并且直接调用Google Drive中储存的GEE结果。这不是免费的翻墙服务器嘛？原来我们还是需要初始化的，可能国内的网络环境还是和国外的不同。但好处是，不需要像geemap那样配置端口，更简洁。
2025-01-12 21:30:30
202
原创 GEE OEEL kriging 包Python API版本（geemap）-克里金插值
众所周知OEEL包是一个好用的GEE libaray, 根据吴秋生老师和OEEL官网的教程，均表示使用下列代码，可以直接在python api中调用JavaScript写的OEEL函数。所以我就把用到的OEEL函数克里金补缺失值的函数源码转为了Python版本。然而我试了好久好久，一直都是在加载的状态....
2024-12-06 15:42:54
306
2
原创 GEE ui界面实现：用户自画多边形， 按面积比例在多边形中自动生成样点，导出多边形和样点shp，以及删除上一组多边形和样点（有视频效果展示）
用户自画多边形， 按面积比例在多边形中自动生成样点，导出多边形。
2024-11-11 14:39:13
458
原创 conda 虚拟环境安装GDAL
后面看了这篇教程解决了问题（https://mp.weixin.qq.com/s/FMabnSCD_NgfeWfsD4bU_Q），尤其是在conda的虚拟环境上安装GDAL，即用conda安装geopandas和cartopy,GDAL即可自动被安装，同时你还将自动安装了许多地理信息处理的依赖包。换了新电脑，要重新安装GDAL。从前是下了GDAL的.wheel文件用pip安装，但平时下轮子的网站现在都打不开，比如。
2024-09-24 09:44:59
1029
2
原创 GEE:多面板同步缩放查看多源数据，并实现交互选点构建NDVI曲线
一. 目标①构建三个面板，分别显示不同来源数据；②面板1显示哨兵数据+面版2显示谷歌高清数据+面板3实现用户任意交互选点，并以该点为中心构建正方形，随后生成该正方形的区域NDVI平均值长时序曲线；③保证前两个面板可以同步缩放查看数据。二.代码三. 结果使用最左侧面板的选点工具在图中点击后，即可得到下图，也就是最右侧的NDVI时间序列。四. 参考GEE中多视窗链接查看地图要素这是本文实现的基础，十分感谢这位大佬~完结撒花~
2024-07-27 17:55:32
445
原创 Landsat 5、7、8和9影像的时间段对比列表
1999年4月15日 - 2003年5月30日。2003年5月31日 - 至今（有SLC故障）水一篇，最近常用，记录在此好查看🤭。
2024-07-26 15:09:25
645
1
原创 GEE：设置ui.Map.Layer上交互矢量边界填充颜色为空，只显示边界
最近在GEE的交互功能鼓捣一些事情，在利用buffer功能实现了通过选点建立一个矩形后，需要将该矩形填充颜色设为空，只留边界。
2024-07-23 18:59:14
291
原创 怎样选博士论文课题
博士论文课题决定了我们未来博士的研究内容，甚至影响着未来的职业选择或学术道路。因此，目前博0的我筛选并学习了b站相关视频，并提炼出有效信息与大家分享。
2024-07-08 14:35:49
1217
原创 逐栅格并行计算偏相关分析
代码分为 ‘process(函数)’ 和 ‘实际运行’ 两个.py文件在pycharm中运行（这样避免了‘溢出’报错），并结合自己情况，一口气输出上述12个组合的结果。自变量（4个）：tem、pre、iet、iep。因变量（3个）：SOS、EOS、LOS。加入偏相关系数显著性（p值）的计算。出处声明：本代码是在借鉴这篇。完结撒花~欢迎交流~
2024-01-31 11:49:36
2346
4
原创 硕士论文结论与展望简要框架
1）不足1+因此，后续研究将，+展望1。2）不足2+因此，后续研究将，+展望2。为探究xxx对xxx的影响，本研究。2）方法2+研究结果表明，+结论2。不足3+因此，未来研究将，+展望3。）方法1+研究结果表明，+结论1。方法3+研究结果表明，+结论3。【简要背景（1句话）
2024-01-25 09:48:48
1113
原创 硕士论文背景与国内外研究进展简要框架
每段结构：目前别人所使用的方法、得到的结果、存在的不足-》突出我的创新点。多少个创新点，就多少个上述结构。每个创新点也是每节的核心。eg:【人类活动-气候变化-植被响应】②为什么【研究的东西是否存在潜在逻辑】（只是为了让自己写起来心中有数）①是什么【研究的东西的基本概念】③重要性【做这个课题的意义】
2024-01-22 11:19:19
561
原创 论文发表必备技能：空间数据（tif图像、shp边界）多子图可视化-基于cartopy与matplotlib的python包
具体地，我有三个因变量分别是SOS, EOS, LOS，四个自变量分别是tem, pre, IET, IEP，因此，我需要画三张大图，每张图内包含四张子图，即三个因变量与四个自变量分别的偏相关系数分布图。其内容涉及到，可视化tif图像、shp边界、绘制colorbar、在图框上显示经纬度等优雅操作。这部分我将colorbar的title写为中文的（可能是因为要毕业了吧）。现新鲜出炉，分享给大家，热忱期盼交流指正~众所周知，我正在写毕业论文，刚画完一张。我的代码总共分为三部分，包括。
2024-01-06 21:52:14
995
1
原创 ax.add_geometries无法加载shp文件
正常运行ax.add_geometries代码，但是shp文件在结果图中一直不显示。一开始还以为是投影有问题，代码中设置的是。最近在写毕业论文，上手练习一下cartopy包。，但后面也尝试了其他投影，依然显示不了。
2024-01-05 10:26:19
808
1
原创 输出两幅相似影像未重叠像素的位置
在gee上下载的物候数据中，EOS波段比SOS波段多了两个像素，这应该就属于系统误差之类的问题了。因此，本文将识别出额外的这两个像素的位置，而后将SOS波段对应位置的数值去除，保证数据的一致性。
2023-10-08 20:23:41
129
1
原创 VNP22Q2 物候数据集简介及波段value转DOY代码
四个物候指标为例。其中LOS为SOS和EOS之差，故而不需要offset计算。因此重点是前三个指标。我的情况：从gee上将2013-2019的数据堆叠为一张图下载，总共38*7=266个波段，目标提取。step1 首先在arcgis将图像裁剪。本文以EOS（对应下表。波段数：每年38个波段。、LOS四个物候指标。空间分辨率：500m。
2023-10-07 14:54:07
526
6
原创 Arcpy批量导出shp文件属性表——使用arcpy.da.SearchCursor函数
Arcpy批量导出shp文件属性表——使用arcpy.da.SearchCursor函数
2022-11-02 11:13:57
2168
原创 Anuspline插值结果问题汇总：在Arcgis中显示错误---有条带；数值范围0-255；插值文件大小0k
Anuspline插值结果问题汇总：在Arcgis中显示错误---有条带；数值范围0-255；插值文件大小0k
2022-09-26 15:28:51
1182
2
原创 Python 自动生成日期（斜杠格式），不自动补零，解决 ‘#‘ is a bad directive in format ‘%Y/%#m/%#d‘报错
Python 自动生成日期（斜杠格式），不自动补零，解决 '#' is a bad directive in format '%Y/%#m/%#d'报错
2022-08-17 16:57:51
2313
原创 循环提取动态命名的dataframe中的某一列：解决SyntaxError: cannot assign to operator 和 如何运用list集成dataframe并提取行列
循环提取动态命名的dataframe中的某一列：解决SyntaxError: cannot assign to operator 和 如何运用list集成dataframe并提取行列
2022-08-08 14:49:00
1517
原创 Citespace安装【Version Check】Your version‘s status cannot be verified due to a network connection issue
解决Citespace安装【Version Check】Your version's status cannot be verified due to a network connection issue问题和citespace（563）问题。
2022-06-09 10:52:55
6834
13
原创 GEE Python API环境配置——解决Jupyter Notebook：KeyError: ‘client_id‘、‘gcloud‘ （超详细内附：如何快速找到ip地址和端口号）
GEE Python API环境配置——解决Jupyter Notebook：KeyError: 'client_id'、'gcloud' （超详细内附：如何快速找到ip地址和端口号）
2022-06-01 20:56:40
4649
13
2024级中国科学院大学（国科大）毕业论文参考文献style文件-Endnote版
2024-03-01
中国气象站点分布点数据图-SHP文件
2022-07-15
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人