爱吃柚子的梨-CSDN博客

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

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

2022-06-19 22:41:37 5777 1

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

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

2022-06-18 16:49:53 1555

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

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

2022-02-12 11:38:02 5894

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

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

2021-11-16 16:04:07 3093

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

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

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

原创 IIS网站跨域设置

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

2021-09-06 10:23:22 1234

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

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

2021-09-06 10:15:53 1755

原创【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 556

原创【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 572

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

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

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

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

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

2021-08-01 11:21:07 298

原创【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 241

原创【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 364

原创【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 567

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

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

2021-07-07 09:18:49 1849

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

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

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

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

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

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

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

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

2021-06-10 15:50:11 2445

原创 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 700

原创【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 755

原创【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 747

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

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

2021-04-21 11:22:37 1977

原创【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 549

原创【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 1863 3

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

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

2021-04-18 19:34:46 957 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 2247 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 3028 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 149

原创【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 4511 1

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

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

2021-03-23 22:55:15 13468 11

原创【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 12495 5

原创【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 6294 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 179

原创使用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 6526 8

原创【读书笔记：算法小抄】动态规划解决最长公共子序列问题——Java实现

描述LCS：最长公共子序列问题描述：求两个字符串的最长公共子序列的长度，例如输入str1 = “abcde”, str2 = “aceb”，算法应该输出3。代码import static java.lang.System.out;public class Helloworld { public static int LCS(String s1, String s2){ char[] s1_arr = s1.toCharArray(); char[] s2

2021-03-05 11:41:31 118 1

原创【读书笔记：算法小抄】斜着遍历二维数组

简介本文参考书P113页代码，实现了右上和右下倾斜遍历二维数组，即：代码import static java.lang.System.out;public class Helloworld { public static void main(String[] args){ int n = 5; // 定义一个5*5的矩阵 out.println("右上倾斜遍历行列号"); for (int l = 1; l <= n; l++

2021-03-04 19:24:27 351 2

原创【读书笔记：算法小抄】动态规划之最大子数组问题——Java实现

简介最大子数组问题输入一个整数数组nums，请你在其中找一个和最大的子数组，返回这个子数组的和。例如输入nums = {-3, 1, 3, -1, 2, -4, -2}，返回，因为最大子数组为{1, 3, -1, 2}代码import static java.lang.System.out;public class Helloworld { public static int maxSubArray(int[] nums){ if (nums.length == 0)

2021-03-03 21:43:10 269 1

原创【读书笔记：算法小抄】二维递增子序列：信封嵌套问题——Java实现

Java常识排序基础对二维数组或者其他对象进行排序，常用Arrays.sort(T[] a, Comparator<? super T> c)。假如传入的是一个二维数组int[][] arr，那么被比较的对象就是arr[0]，arr[1]一类的。实现Comparator接口通常需要定义Compare函数，且看Compare函数的文档解释：int compare(T o1, T o2)Compares its two arguments for order. Returns a neg

2021-02-26 17:02:00 362

原创【读书笔记：算法小抄】最长递增子序列——Java实现

简介最长递增子序列：输入一个无序的整数数组，请你找到其中最长递增子序列的长度。例如，输入nums=[10, 9, 2, 5, 3, 7, 101, 18]，其中最长的递增子序列是[2, 3, 7, 101]，算法输出4。该题有两种主流的解法，一种是使用动态规划法，另一种为二分搜索解法。动态规划解法代码package com.company;import static java.lang.System.out;import java.util.*;public class Main {

2021-02-25 10:55:01 184