爱吃柚子的梨-优快云博客

原创【超好玩】Aerialod：开启 3D 艺术质感地理地图新时代？

地图不仅仅是一张纸，它是我们探索世界的窗口。’ 在科技飞速发展的今天，有什么好的地图工具可以制作出令人惊叹的地理地图呢？今天我将介绍如何利用 Aerialod 制作带有 3D 艺术质感的地理地图。

2025-02-03 22:24:58 465

原创【巨实用】Git客户端基本操作

本文主要分享Git的一些基本常规操作，手把手教你如何配置~

2025-01-09 22:48:02 561

原创【实用分享】发表在Science 上的一种聚类算法的实现-基于快速搜索和发现密度峰值的聚类方法

本文基于2021年发表在Science上的聚类算法《Clustering by fast search and find of density peaks》的实现，并采用Python语言实现了该聚类算法。采用的Python语言版本为3.8，除去读取GIS数据所用的ogr库和绘图所用的matplotlib、pandas、seaborn库，真正运行算法所必要的库仅为numpy。

2025-01-09 22:20:00 687

翻译【分享】发表在Science上的一种实用聚类算法-基于快速搜索和发现密度峰值的聚类方法

聚类分析的目的是根据元素的相似性将其分类。它的应用范围从天文学到生物信息学、文献计量学和模式识别。我们提出了基于这样一种思想的方法，即簇中心的密度比它们的邻居高，并且与高密度点的距离相对较大。这一思想构成了聚类过程的基础，在这个过程中，聚类的数量会直观地产生，离群值会自动被发现并排除在分析之外，而不管聚类的形状和嵌入的空间的维数如何，聚类都会被识别出来。我们在几个测试用例上展示了该算法的强大功能。

2025-01-09 22:17:41 73

原创【超实用】QGIS矢量数据字段计算器的使用（QGIS Field Calculator）

打开属性表（Attribute Table）,点击红框所指的按钮，便是字段计算器（Field Caculator）：字段计算器界面可以分为三个部分，图中第一部分为字段的选择。如果你要创建一个新的字段，就勾上“Create a new field“，并更新”Output field“相关信息。或者你也可以选择更新现有字段的值，勾上“Update existing field”。第二部分为计算表达式，用于计算字段值。第三部分为函数的选择，用于辅助你利用表达式计算复杂的字段值。

2025-01-04 17:08:57 850

原创【超简单】官方下载全国标准行政规划矢量地图（精确到县）（带审图号）GS(2024) 0650

如果你实在不想动手的话，我也替大家把数据直接下载好了，并提供转换好的shp的版本给到大家，关注我的微信公众号“爱吃柚子的梨”，回复关键字“XZGH”即可~

2024-12-29 23:49:22 1022

原创关于WGS84地理坐标转CGCS2000投影坐标的简单处理办法

首先我是基于WGS84和CGCS2000两个地理坐标系相差极小的事实，对此我默认二者的坐标就是一致的（误差可以忽略）。那么问题就变得十分简单了，只需要现在ArcGIS中定义投影为China Geodetic Coordinate System 2000，然后再投影至相对应的投影坐标系即可！......

2022-06-19 22:41:37 6356

原创关于按照矢量数据的区域划分对栅格进行总和统计的方式

首先，使用ArcGIS中的按表分区统计（Zonal Statistics as Table），矢量数据中的分区属性最好是Unique，生成表格。接着对矢量数据按照分区属性的字段进行连接（joining），这样刚刚生成的统计表格就很好的和矢量数据连接起来，可以在一张属性表中显示和操作了。...

2022-06-18 16:49:53 1857 1

原创【ArcMap学习系列】符号编辑与制作（国界、省界、区县界等线状符号制作方法）

ArcMap中符号制作的入口ArcMap中线状符号制作方法（未完，待续…）如果你认为对你有用，关注我的微信公众号支持我一下吧！~

2022-02-12 11:38:02 7530

原创多年份栅格数据裁剪后范围不一致解决方法

首先使用ArcGIS的clip raster工具重新裁剪得到一份新的栅格数据（用矢量数据进行裁剪），然后在ArcGIS中查看该数据，右击选择导出，导出的时候范围选择基准栅格数据的范围即可。

2021-11-16 16:04:07 3530

原创【FME小计】关于用户参数的设置

工作的时候遇到用户参数相关的需求：用户选择特定年份的时候，就运行特定年份数据。该项需求的实现有两种方案：使用Choice with Alias（带别名的选择），然后将该参数关联到数据表上（Feature to Read）使用字符拼接的方式，先创建年份的用户参数，然后创建一个私有参数，类型为文本，然后选择文本编辑，在里面进行字符的拼接。（注意必须为私有参数）...

2021-10-21 23:04:14 765 1

原创 IIS网站跨域设置

首先确保IIS相关功能打开了添加跨域响应标头

2021-09-06 10:23:22 1300

原创【IntelliJ IDEA】使用Spring Boot连接oracle数据库出现[08006][17002] IO 错误: Got minus one from a read call

是驱动版本不一致导致的问题，需要更换至所使用的Java版本OJDBC。如下图所示：

2021-09-06 10:15:53 1986

原创【wxPython】创建RibbonGUI

如果想要添加一个RibbonButton，那么顺序必须为Sizer->RibbonPage->RibbonPanel->RibbonButtonBar->RibbonButton对于RibbonButton的事件绑定采用如下代码：elf.Bind(rb.EVT_RIBBONBUTTONBAR_CLICKED, self.SayHello, id=wx.ID_ANY)...

2021-08-19 11:17:18 768

原创【Python报错】UnicodeDecodeError: ‘gbk‘ codec can‘t decode byte 0xac in position 22: illegal multibyte..

【解决办法】在用Open打开文件的时候统一为UTF-8编码即可解决。open(read_file_path, encoding='UTF-8')

2021-08-17 11:15:53 691

原创【Windows 10】我们无法设置移动热点

解决措施打开设备管理器点击网络适配器右击Intel® Dual Band Wireless _AC 8265，点击更新驱动程序点击第二个浏览电脑点击第二个从计算机可用程序中选择随意选一个驱动程序（默认即可），点击下一步便大功告成参考资料https://www.zhihu.com/question/53106181...

2021-08-03 10:00:14 3715 1

原创【聚类算法学习】Clustering by fast search and find of density peaks

该算法的优势与劣势优势能识别非球形的簇（相比较于K-means和K-medoids）参数配置简单（相比较于DBCSAN，该算法仅需关注阶段距离d_c的确定，且该距离的确定作者给出了经验值范围，所以算法的运行几乎依靠默认设定的参数即可）可自动确定聚类中心的个数（该论文被吐槽的一点，目前还是需要结合决策图与降序排序的γ值图人工判定聚类中心个数）...

2021-08-01 11:21:07 365

原创【Python 2.7】Python Exception ＜UnicodeEncodeError＞: ‘ascii‘ codec can‘t encode characters

这是由于对ascii编码错误造成的，可以尝试：在变量后面加.encode(‘gbk’)在变量后面加.encode(‘utf-8’)注意：如果遇到中文输入变量和匹配值不匹配的情况，如x=‘你好’，但是判断x=='你好’的时候为FALSE，这可能是由于半角字符、全角字符或者空格等等原因造成的，建议直接复制或者导入匹配值列表，确保二者信息完全一致。（其实更建议用英文）...

2021-07-14 12:17:09 311

原创【FME报错】Excel Reader: Found named range ‘Print_Area‘ on sheet ‘‘. This is a Microsoft Excel ...

如果出现用FME处理与Excel表格相关的数据时，出现“Excel Reader: Found named range ‘Print_Area’ on sheet ‘’. This is a Microsoft Excel internal name so FME will not show it as a feature type. If you wish to read this named range, please rename it.”的报错信息，目前可以尝试的方法有以下几个：关闭所有Exce

2021-07-13 16:12:06 477

原创【FME & Python】FME中PythonCaller使用小记

打印输出变量使用print无法输出到日志中，使得难以调试，对此可以使用FMELoggler来输出Python中的变量，在PythonCaller文本中插入如下代码：logger = fmeobjects.FMELogFile() logger.logMessageString("Hello world!")参考资料Logging with Python scripts. https://community.safe.com/s/article/logging-with-python-scrip

2021-07-08 15:04:41 730

原创【FME转换器】常用FME模型转换器总结

计算类：StatisticsCalculator：可统计属性最小值、最大值、中位数、总和、计数、平均值和标准差等。条件类：Tester：属性筛选几何操作类：GeometryRemover：将空间数据转为非空间数据，移除其中的空间属性PointOnAreaOverlayer：点LineOnAreaOverlayer：几何判断类：SpatialFilter：根据空间关系如点面之间的相交、包含等筛选要素GeometryValidator：拓扑验证，用于判断数据是否正确合理Geo

2021-07-07 09:18:49 2246

原创【FME & EXCEL】FME玩转EXCEL小记——不定期更新

目录前言本文记述FME使用过程中，对EXCEL操控心得，皆为平时的有感而发，故会不定时的更新。模版的定义与使用有如下一份成绩表，先存在两个问题，一个是所有同学的成绩少加了5分，另一个是登记表的位置不太对，应该是在A-3的位置。我们通过FME来尝试解决以上两个问题。1. 通过输入数据确定属性列表如下图所示，FME中读取输入数据，将第二行数据作为属性列表。点击Parameters可以修改输入时的参数，如下图红框所示，通过修改行列号可以确定属性列表的读取范围以及属性的类型等：2. 操控数据

2021-07-05 21:51:35 3686 2

原创 ArcGIS报错the database client software failed to load不能连接到oracle数据库,因为数据库客户端软件加载失败

原因很简单，就是没有按照oracle数据库导致的，去官网下载一个数据库即可，我是ArcGIS 10.2版本的下载了Oracle 11g以后重启电脑成功解决问题！

2021-06-26 15:14:20 2105 1

原创【ArcGIS学习系列】高级编辑功能介绍——持续不定期更新

写在前面随着研究生阶段学习的深入，逐渐地用到了ArcGIS中一些以前从未想过会使用的功能，以下是对我用到的关于编辑器功能的介绍。合并（Merge）以线要素为例，该功能可以在编辑状态下将多条线段合并为一个要素。使用也很简单，选中多条线段，在Editor下拉列表中选择Merge。打断相交线（Planarize Lines）该功能可以将多条在空间上相交的线段在其相交点处一次性全部打断，这样你可以不用重复一天天的手动分割线段。使用上，先选择多条线段，再在Advanced Editing中

2021-06-10 15:50:11 2937

原创 shapely.errors.TopologicalError:...Likely cause is invalidity of the geometry

这个报错的原因是由于输入的矢量文件中包含了无效的多边形导致的，解决办法是使用ArcGIS Pro的Check geometry和Repair geometry工具，其中Validation Method选择OGC。

2021-05-11 16:37:29 875

原创【Python实现Select by Location】使用geopandas对多边形进行buffer筛选

import geopandas as gpddist_buffer = 500jy_12_poly_gdf = gpd.read_file(jy_12_split_file)for i in range(100): nei_poly = jy_12_poly_gdf[jy_12_poly_gdf.within(jy_12_poly_gdf['geometry'][i].buffer(dist_buffer))]

2021-05-05 23:30:33 878

原创【numpy计算欧几何距离】多点至单点的距离

方式一下面数组中，第一列是表示点的ID号，第二、三列分别表示x和y坐标。import numpy as npa = np.array([[0, 1, 2], [1, 2, 3], [2, 3, 4]])b = np.array([[1, 2, 3]])d_arr = np.sqrt(np.power(a[:, 1] - b[:, 1], 2) + np.power(a[:, 2] - b[:, 2], 2))print(d_arr)方式二numpy.linalg.norm...

2021-05-05 21:31:50 879

原创【PostgreSQL】ERROR: 错误: 关系 “weather“ 不存在

错误的SQL查询语句：SELECT * FROM Weather引发报错ERROR: 错误: 关系 “weather” 不存在解决办法：给Weather添加双引号，如果你创建表的时候使用了双引号，那么在引用的时候也要添加。SELECT * FROM "Weather"...

2021-04-21 11:22:37 2169

原创【Python】ogr使用小记

图层的属性筛选 # 通过属性进行要素筛选 lyr.SetAttributeFilter(field_name + "=" + str(val)) # 还原图层值筛选前的状态 lyr.SetAttributeFilter(None)ogr中的geometry与shapely互相转换test_poly_file = r"F:\Temp\temp_gis_data\testSplittingPolygon.shp"ds = shp_driver.Open(test.

2021-04-21 09:54:50 620

原创【pyinstaller】modulenotfounderror no module named ‘rasterio._shim‘、‘rasterio.control‘等

对于出现no module not found error的情况一般的解决办法就是使用hidden imports参数当前对于这种问题的解决方法就是将所有缺少的模块用hidden imports给加进去，当然手动加有些麻烦，故参考博文https://stackoverflow.com/questions/53149750/something-wrong-with-how-im-bundling-rasterio-into-an-executable给出的解决方法：# -*- mode: pyt

2021-04-18 22:30:29 2443 4

原创 conda安装rasterio失败原因分析

conda提示是使用的Python版本不对，当我从3.8改成了3.7的时候conda安装rasterio倒是成功了。

2021-04-18 19:34:46 1173 1

原创【Python】geos_c.dll缺失解决方法

打开下面网址，或者自行百度下载geos_c.dllhttps://www.dll-files.com/download/d8b5101f07394b4562ef673869395443/geos_c.dll.html?c=eXJISm5CSFNJREh3aXN4WVllQy9ydz09放到可执行文件旁或者需要对应的目录如果你认为对你有用，或者想要与我有更多交流，可以关注我的微信公众号支持我哦！~...

2021-04-15 23:48:40 2447 1

原创【Pyinstaller/cx_freeza已解决】oserror: could not find or load spatialindex_c-64.dll

出现这个问题的原因是geopandas引用的库rtree找不到spatialindex_c-64.dll导致的，那么解决这个问题的直接办法就是修改源代码，首先根据报错提示（cx_freeza）我找到出错的原文件路径：D:\ProgramsInD\anaxonda3\Lib\site-packages\rtree\core.py报错显示出错的行数是129行，找到该行代码：既然是找到不的原因导致的，那就给它这个文件的路径便是了，直接将源代码改为： if not rt: if o

2021-04-15 15:04:09 3181 6

原创【英语作文】常见错误记录

Right: return to the upper riverWrong: return the upper river返回某地要加toRight: approximate尽量别用缩写approx在作文中Right: enter a lower riverWrong: enter to a lower riverenter是及物动词后面直接跟名词Right: live here for 10 yearsWrong: live in here for 10 yearsHere是副词，l.

2021-04-06 14:48:49 244

原创【Delaunay三角网】scipy.spatial.Delaunay函数解析及其与Geopandas结合绘制Delaunay三角网

常用参数解析points：包含x,y坐标的点，用于生成Delaunay三角网。格式：[[x1, y1], [x2, y2], [x3, y3], …]incremental: 是否允许继续加点，默认False常用属性解析points：即你输入的那些点simplices：Delaunay函数生成的那些三角形，每个三角形用其三个顶点的顶点号表示，每个三角形的顶点号按照逆时针排列neighbors：每个三角形的邻居三角形，每个三角形至多三个邻居，第k个邻居在第k个顶点的对面（就是

2021-04-04 11:29:42 5338 1

原创【VSCode】设置敲写函数自动补全括号

打开设置搜索addBrackets，勾上选项重启VSCode测试是否生效如果你认为对你有用，关注我的微信公众号支持我一下吧！~

2021-03-23 22:55:15 16209 12

原创【ArcGIS系列】道路距离地图制作（Roads to Distance）

简介简单的说就是做出一个城市的道路距离分布图，里面显示了每一个点距离最近的城市道路的距离，并将所有的距离值归化至0-255的区间，效果如下：数据准备某一地区的行政边界数据某一地区的路网数据（可选）将数据进行投影本文选择将地理坐标系转换为投影坐标系备注：UTM如何快速找到自己所属的精度分区？点击下面的网址在地图上直接查看。https://mangomap.com/robertyoung/maps/69585/what-utm-zone-am-i-in-#计算路网距离在T

2021-03-10 21:48:34 15105 9

原创【OpenStreetMap】任意城市道路数据下载（附带数据解释）

准备工作（可选）所需城市的行政边界数据，如江阴市：QGIS 2.x版本，2点几无所谓，本文用的是2.18版本（下载地址：https://qgis.org/downloads/）QGIS中下载OpenStreetMap道路数据打开QGIS 2.18打开QGIS中默认的OpenStreetMap插件数据下载QGIS解压OSM数据目录中点击数据勾选想要查看的数据，如点击第一个另存为矢量shp文件后即可在ArcMap中查看使用路网数据解释在ArcMap中打

2021-03-08 13:37:51 6949 11

原创【读书笔记：算法小抄】动态规划解决编辑距离问题——Java实现

简介可以对一个字符串进行三种操作：插入一个字符，删除一个字符，替换一个字符，给你两个字符串s1和s2，请计算将s1转换成s2最少需要多少次操作。代码import static java.lang.System.out;public class Helloworld { public static int minDistance(String s1, String s2){ int m = s1.length(), n = s2.length(); // 定义

2021-03-08 10:30:17 220

原创使用GEE（Google Earth Engine）下载STRM30m分辨率高程数据

简介本文介绍如何使用GEE（Google Earth Engine）下载STRM30m分辨率高程数据并用ArcMap加载行政边界裁剪栅格图层。准备工作申请Google账号注册Google Drive在GEE上进行注册（https://earthengine.google.com/）搜索STRM 30m，选择红框中的那个：（https://developers.google.com/earth-engine/datasets/catalog/USGS_SRTMGL1_003）编写代码下载到

2021-03-07 16:54:15 7634 9