- 博客(83)
- 资源 (2)
- 收藏
- 关注
RSS订阅
原创 决策树经典算法ID3——我的第一篇博客
简介 决策树是机器学习中一种常见的分类算法，属于有监督学习算法(至于什么是有监督学习，什么是无监督学习读者可以自行百度)。决策树算法有多种，ID3算法是其中一种经典的决策树算法，这种算法的核心是信息熵（至于什么是信息熵，后面会进行详细介绍）。现在已经商用的决策树算法C45，C50等都是在ID3的基础上进行改进优化而来。一些术语1、分类特征比如一个数据集有A、B、C、D、E五个属性（或者说五个字段），这五个属性不是相互独立的，而是存在某种联系，基于这种联系可以通过其中某些属性的取值，确定另外一个或
2021-02-06 15:53:02
6231
5
原创 geopandas：数据读取、坐标设置与转换、空间查询连接、数据输出、距离计算
import geopandas as gpimport pandas as pdpolygonShp=gp.read_file("./shp/globeland30-10.shp")polygonShp OBJECTID Id gridcode Shape_Leng Shape_Area geometry 0 1 1
2022-05-16 20:34:56
4211
原创 点集凸包算法python实现（二）
算法逻辑在点集凸包算法python实现这篇博客中介绍了一种凸包算法，这种算法中凸包点搜索的过程较为麻烦，主要是因为计算点集连线与X轴的夹角需要考虑到四个不同象限，在这里通过计算向量夹角的方式，对凸包点集的搜索过程进行了优化算法逻辑如下图所示：python代码实现from osgeo import ogr, gdalconst, osrimport numpy as npdef getAlfa(startPoint, endPoint): dx = endPoint[0] - sta
2022-04-22 22:16:09
1299
原创 点集凸包算法python实现
什么是凸包？凸包定义点集p的凸包是指一个最小凸多边形（内角均小于180°），满足p中的点或者在多边形边上或者在其内下图中的红色线段表示的多边形就是点集p={p0,p1,p2,p3,…………,p12}的凸包通俗理解一组平面上的点，求一个包含所有点的最小的凸多边形这可以形象地想成这样：在地上放置一些不可移动的木桩(代表点集中的点)，用一根绳子把他们尽量紧地圈起来，这就是凸包凸包有什么特点？整个凸包都在任意一条边的一侧凸包任意两点的中点都在凸包内凸包内的任意点集的加权平均（凸组合）都
2022-04-22 15:58:41
4129
原创 python gdal 栅格影像拼接工具
先上代码from osgeo import gdal, gdalconst, osrimport argparseimport osif __name__ == "__main__": parser = argparse.ArgumentParser( prog='raster files mosaic', description='对输入目录中的栅格影像进行拼接，要求这些影像具有相同的空间参考') parser.add_argument('inputDir',
2022-04-18 23:43:32
2957
原创 Arcpy自动批量化制图：修改图层数据源 修改布局元素 导出布局地图
先上代码# coding:utf-8_import arcpyimport osif __name__ == "__main__": tifdir = r'I:\DInSAR\PJQZ\result' shpdir = r'I:\DInSAR\reference' tifs = [] with open("./filter_list.txt", mode='r') as f: tifs = f.readlines() tifs = [tif.
2022-04-14 12:09:48
1431
原创 python gdal多值提取至点工具
提取shp点文件栅格目录中的栅格数值序列import osfrom osgeo import gdal, gdalconst, ogrimport pandas as pdclass ExtractByPoint(object): def __init__(self, x, y, tifFile) -> None: self.x = x self.y = y self.tifFile = tifFile# 地理坐标转像素坐标
2022-04-12 20:36:04
3902
原创 C++ GDAL 遥感影像K-means聚类
#include <gdal_priv.h>#include <string>#include <iostream>#include <ctime>#include <vector>#include <map>#include <algorithm>using namespace std;class Kmeans{private: GDALDataset* m_ds;//待分类的影像 int m_x
2022-04-12 11:05:36
609
原创 Arcpy计算字段
# coding=utf-8# 字段分类import arcpyimport os# 获取当前shp的所有字段def getFieldNames(fc): fieldList = [] for f in arcpy.ListFields(fc): fieldList.append(str(f.name)) return fieldListRecalss_area = """def Reclass(area1): if (area1 <
2022-04-10 10:16:32
1192
原创 Arcpy读取SHP
#coding:utf-8_import arcpyimport osdef sim(shpFile1, shpFile2): ds1 = arcpy.CopyFeatures_management(shpFile1, arcpy.Geometry()) ds2 = arcpy.CopyFeatures_management(shpFile2, arcpy.Geometry()) e1 = ds1[0] e2 = ds2[0] # 求两个椭圆的相交部分
2022-04-08 14:04:23
1298
原创 QGIS插件开发详细教程
为了提高生产力，接触到了QGIS插件开发。从相关的知识准备，到实际插件开发，断断续续大约花了一周的时间，最后也算是实现了自己的目的，感觉有必要记录一下其中遇到的坑。1 知识准备进行QGIS插件开发，总体上需要两个方面的知识储备熟悉QGIS提供的各种类库和函数接口这方面最好的资料是官方文档,内容详细而且比较有条理。如果以后经常接触到QGIS二次开发的话，建议吧整个文档好好看一遍，如果时间有限，也可以根据自己的需要有针对性的阅读了解PyQt的基本使用QT是C++编写的跨平台GUI库，PyQt
2022-03-23 10:47:57
10256
6
原创 PyQt QGraphicsView组件视图导出为图片
PyQt QGraphicsView组件视图导出为图片可以通过调用QGraphicsView.render()接口将组件视图渲染到QPixmap或QImage上，再调用QPixmap或QImage的**save()**方法保存为图片rect = QGraphicsView.viewport(self.dlg.gvPointRecords).rect()pixmap = QPixmap(rect.size())painter = QPainter(pixmap)painter.begin(pixma
2022-03-22 12:02:49
1463
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人