staHuri
https://github.com/huifer/planar_algorithm
展开
-
GIT版本管理
基于GIT版本管理Author: HuiFerGit是一个 “分布式版本管理工具”,简单的理解版本管理工具:大家在写东西的时候都用过 “回撤” 这个功能,但是回撤只能回撤几步,假如想要找回我三天之前的修改,光用 “回撤” 是找不回来的。而 “版本管理工具” 能记录每次的修改,只要提交到版本仓库,你就可以找到之前任何时刻的状态(文本状态)。开卷必读如果你还没有使用过GIT或者不知道...原创 2019-10-28 20:38:40 · 314 阅读 · 0 评论 -
【python】图片插入文字
问题如何在图片中插入大量文字并且自动换行效果原始图效果图实现方式from PIL import Image, ImageDraw, ImageFontclass ImgText: font = ImageFont.truetype("micross.ttf", 24) def __init__(self, text): # 预设宽度 可以修改...原创 2018-11-07 13:44:43 · 1280 阅读 · 0 评论 -
高斯消除矩阵
#! /usr/bin/env python# -*- coding: utf-8 -*-#def pprint(A): for i in A: print(i) print("")data = [ [1, 2, 1, 2], [3, 8, 1, 12], [0, 4, 1, 2]]n = len(data)print...原创 2018-11-09 15:59:28 · 870 阅读 · 0 评论 -
【测绘程序设计试题集】 试题06 轨迹数据压缩算法
数据P0,107.605,137.329P1,122.274,169.126P2,132.559,179.311P3,153.324,184.276P4,171.884,174.654P5,186.408,168.634P6,196.566,145.204P7,200.549,127.877P8,211.391,118.179P9,216.318,116.547P10,225...原创 2018-10-19 09:54:10 · 2152 阅读 · 0 评论 -
【测绘程序设计试题集】 试题03 利用线性回归模型进行卫星轨道的预报
数据300,21182.88,-7044.56,14639.48600,21707.87,-6930.28,13906.68900,22207.04,-6828.65,13147.661200,22679.16,-6738.66,12363.841500,23123.06,-6659.23,11556.711800,23537.69,-6589.21,10727.782100,239...原创 2018-10-18 15:12:05 · 1698 阅读 · 0 评论 -
【测绘程序设计试题集】 试题02 矩阵卷积计算
数据N 矩阵10.00 13.50 14.00 13.80 13.90 15.60 13.30 14.50 13.70 14.4013.50 13.30 15.10 16.40 15.40 14.90 11.30 13.50 17.70 13.3015.70 14.00 16.30 18.60 16.80 16.60 12.50 15.50 16.70 14.8016.50 15.90 ...原创 2018-10-18 10:00:37 · 2015 阅读 · 3 评论 -
【测绘程序设计试题集】 试题01 出租车轨迹数据计算
数据问题解#!/usr/bin/env python# -*- coding: utf-8 -*-# @File : Taxi.py# @Author : huifer# @Time : 2018/10/17 20:29import mathimport datetimeclass Taxi(object): def __init__(self, ...原创 2018-10-18 08:19:27 · 5795 阅读 · 9 评论 -
【python】字符串问题
python 字符串问题在arcpy中版本为 python2.x在QGIS中版本为 python2.x 或者 python3.xpython2 和python3 之间的str处理方式经常会导致乱码,故出此文python3版本# 将str或字节并始终返回strdef to_str(bytes_or_str): if isinstance(bytes_or_str, ...原创 2018-09-05 11:48:15 · 267 阅读 · 0 评论 -
【python】多图片合并PDF
python 多图片合并pdf起因一个做美工的朋友需要将多个图片jpg 、png 合并起来,PS操作太慢了所以用了python进行完成这个任务代码#!/usr/bin/env python# -*- coding: utf-8 -*-# @File : 2.py# @Author: huifer# @Date : 2018/8/20from PIL impor...原创 2018-08-20 19:17:37 · 6018 阅读 · 9 评论 -
最小二乘法矩阵
#! /usr/bin/env python# -*- coding: utf-8 -*-import numpy as npdef calc_left_k_mat(k): """ 获得左侧k矩阵 :param k: :return: """ k_mat = [] for i in range(k + 1): no...原创 2018-11-28 09:38:04 · 880 阅读 · 0 评论 -
多边形等分
多边形等分思路前提条件封闭面,不可以有空洞解题思路封闭面中随机构造点利用 ** k-means ** 分组(在此方法中设置分组数量 簇数)计算每一个簇的质心利用质心绘制 voronoi 泰森多边形利用封闭面切割泰森多边形注本文中判断点是否在面中利用了 这篇文章实现#! /usr/bin/env python# -*- coding: utf-8 -*-impor...原创 2019-01-14 19:54:55 · 4013 阅读 · 1 评论 -
GIS就业分析
简介本文针对 北京、上海、深圳、广州、杭州、武汉 这几个城市进行GIS就业分析,分析内容:gis行业需求量,gis开发人员需求占比,gis相关技能。本文数据来源 智联数据下载爬虫构思方案相对简单,直接访问一个接口即可#! /usr/bin/env python# -*- coding: utf-8 -*-# __file__: download_zhilianimport r...原创 2019-05-13 12:15:16 · 1279 阅读 · 2 评论 -
[python] 路径相似度
问题以及解/path/1/b1/path/{a}/b1/path/{a}/b1/c1.输入一个路径的值/path/{正则表达式}/{正则表达式} 设这个字符串为ST12.输入一个被比较的路径值/path/1/b1 :ST2将第一个字符串ST1 进行split("/")将第二个字符串ST2 进行split("/")【“path”,正则1,正则2】列表1【path,1,b...原创 2019-05-09 16:40:39 · 2276 阅读 · 0 评论 -
tensorflow 按区域提取实例
tensorflow 按区域提取实例需求输入一张图片,按照规定正方形尺寸(如30*30 像素) 对输入图片进行提取,类似于口 ,田两个文字 一张图变成四张,四张图能够变成一张图。数据与思路数据资源输入图为534*300提取尺寸为30*30思路根据输入情况来看 534 不能够被30 整除,需要对534这个宽度进行调整 。300则被30整除不需要进行特殊处理。...原创 2019-05-14 11:24:43 · 1012 阅读 · 0 评论 -
[python] 天地图切片下载优化
天地图切片下载优化前提下载天地图切片地图,这篇文章中提供了一种天地图切片下载方案,其方案需要去寻找一个范围(切片的xy索引区间),在寻找切片范围的时候需要在大量的图片中寻找到结果比较麻烦,因此提出如下需求。需求输入经纬度返回切片索引制作本文参考内容上述代码不能正常运行所以对这个代码进行重构#! /usr/bin/env python# -*- coding: utf-8...原创 2019-05-06 10:49:16 · 2263 阅读 · 0 评论 -
[python] 下载天地图切片地图
下载xyz地图资源下列为常用xyz路由地址为了避免图片中出现文字标注(道路名称,建筑物名称等)本文选择天地图tian-vec 作为获取资源对象 var mapUrl = { /**** * 高德地图 * lang可以通过zh_cn设置中文,en设置英文,size基本无作用,scl设置标注还是底图,scl=1代表注记, ...原创 2019-04-24 10:07:41 · 7786 阅读 · 4 评论 -
[python] 旅游
旅游爬虫篇目标网页将这部分内容获取到,主要内容有城市名称,城市id,城市网页进一步增加json数据量,把这个目的地下的 行程路线,景点的路由获取 , 以杭州为例,其他城市替换对应的城市id即可,本文就直接进行字符串替换了路线:http://www.mafengwo.cn/mdd/route/10156.html景点:http://www.mafengwo.cn/jd/10...原创 2019-04-12 15:48:59 · 872 阅读 · 0 评论 -
球面拟合-基于最小二乘法
球面拟合-基于最小二乘法球面公式R2=x2+y2+z2R^2 = x^2 + y^2+z^2R2=x2+y2+z2流程思路最终结果是一个质地均匀的球体,并非椭圆,球面上的每一个点到圆心的距离相等设球体中心为(x0,y0,z0),球体表面点和中心点的关系如下R2=(X−x0)2+(Y−y0)2+(Z−z0)2向量化简x2+y2+z2=2Xx0+2Yy0+2Zz0+r2−...原创 2019-04-11 20:20:25 · 4949 阅读 · 0 评论 -
[python] 层次聚类 Agglomerative Clustering
层次聚类下图中动物 就是一种层次分类主要函数cls = AgglomerativeClustering(n_clusters=group_size,linkage='ward')linkage 参数说明:ward (默认值):每一个类簇的方差最小化average:每一个类簇之间的距离的平均值最小complete:每一个类簇之间的距离最大single:每一个类簇之间的距离最小...原创 2019-04-08 16:01:21 · 14925 阅读 · 1 评论 -
【python】pyecharts使用
pyecharts 使用直方图#!/usr/bin/env python# -*- coding: utf-8 -*-# @File : my_bar.py# @Author : huifer# @Date : 2018/8/19from pyecharts import *import randomclass MyBar(Bar): """直方图"""...原创 2018-08-20 10:27:59 · 639 阅读 · 0 评论 -
【python】pyinstaller 打包
起因明日七夕好友表白在网上找到一个表白工具,地址如下 python表白工具由于需要使用python才可以运行,故而使用pyinstaller 进行打包代码说明本文表白代码来自代码文件路径 import sysimport randomimport pygamefrom pygame.locals import *WIDTH, HEIGHT = 640, 4...原创 2018-08-16 19:40:12 · 532 阅读 · 0 评论 -
【python】gdal 获取相关属性
from osgeo import ogr, osrimport redriver = ogr.GetDriverByName('ESRI Shapefile')daShapefile = r'E:\gitck\lineChick\Shapfile\point_test.shp'dataSource = ogr.Open(daShapefile)daLayer = dataSource...原创 2018-04-26 16:49:23 · 1895 阅读 · 0 评论 -
【python】gdal创建普通点线面
import ostry: from osgeo import ogr, gdal from osgeo import osrexcept: import ogr# 属性reader = [ { 'Name': "a", 'Region': "我是", 'Latitude': "86478.950188794...原创 2018-04-26 14:31:38 · 3098 阅读 · 0 评论 -
selenium 滑块验证初级
滑块验证滑块验证网页 <!DOCTYPE html><html><head> <meta charset="utf-8"> <meta http-equiv="Cache-Control" content="no-cache, no-store, must-revalidate">原创 2018-02-28 15:14:38 · 3564 阅读 · 0 评论 -
python 编码测试
import chardettext = "测试编码"""" 测试文本编码"""# print(text.__str__())# print(type(text))encode_s = chardet.detect(text.encode('utf-8'))['encoding']print(text.encode(encode_s).decode(encode_s))原创 2018-01-30 16:47:52 · 219 阅读 · 0 评论 -
关于数组查询的问题
# 假定AA BB是两组线 以起点,终点线段的形式存放在数组中 寻找 AA BB的两组连线所用到的点 此处所有点都使用gx = [['AAp4', 'AAp3'], ['BBp14', 'BBp13'], ['AAp15y1', 'AAp15'], ['BBp7', 'BBp6'], ['AAp2', 'AAp1'], ['AAp3', 'AAp2'], ['AAp4y2', 'AA原创 2018-01-30 14:20:53 · 392 阅读 · 0 评论 -
python 高德地图面坐标获取
id 和 city获取方式可从高德直接搜索获取到"""URL 说明 前面的一堆内容不改变修改 id=B02430IPMB&city=330702"""lis = [ 'https://gaode.com/service/poiInfo?query_type=IDQ&pagesize=20&pagenum=1&qii=true&cluster_state=5&need_utd=原创 2018-01-19 08:43:53 · 2665 阅读 · 0 评论 -
协程url处理模型
from collections import dequedef urlgenerator(url_list): # url 生成器 for url in url_list: yield (url)class UrlProcessor(object): def __init__(self, url_generator): sel原创 2018-01-19 08:40:23 · 340 阅读 · 0 评论 -
python 多线程爬虫模型
import codecsimport queueimport threadingimport timefrom lxml import etreeimport requestsexitFlag = 0class MyThread(threading.Thread): def __init__(self, threadID, name, q): th原创 2018-01-19 08:17:52 · 387 阅读 · 0 评论 -
【python】根据数据模板创建shapefile
#!/usr/bin/env python# -*- coding: utf-8 -*-# @File : copyShapefile.py# @Author: huifer# @Date : 2018-4-28from os.path import existsimport gdalfrom osgeo import ogrfrom os import removeg...原创 2018-04-28 08:57:39 · 1440 阅读 · 0 评论 -
【python】读取JSON 根据模板创建数据
#!/usr/bin/env python# -*- coding: utf-8 -*-# @File : result_point_data.py# @Author: huifer# @Date : 2018-4-28import jsonfrom os import removefrom os.path import existsfrom osgeo import og...原创 2018-04-28 14:59:48 · 1935 阅读 · 0 评论 -
【python】使用folium excel 绘制point
使用folium excel 绘制point制作内容根据气象台资料获得的点进行绘制对一个特殊的点做特别的标注数据来源#!/usr/bin/env python# -*- coding: utf-8 -*-# @File : map03.py# @Author: huifer# @Date : 2018/6/28import pandas as pdimport ...原创 2018-06-28 14:24:34 · 2479 阅读 · 2 评论 -
【python】关于excel和shp的使用在matplotlib
关于excel和shp的使用在matplotlib使用pandas 对excel进行简单操作使用cartopy 读取shpfile 展示到matplotlib中利用shpfile文件中的一些字段进行一些着色处理#!/usr/bin/env python# -*- coding: utf-8 -*-# @File : map02.py# @Author: huifer# @...原创 2018-06-28 11:13:03 · 1575 阅读 · 0 评论 -
【python】shapefile转换geojson
import shapefileimport codecsfrom json import dumps# read the shapefiledef shp2geo(file="line出产.shp"): reader = shapefile.Reader(file) fields = reader.fields[1:] field_names = [field...原创 2018-06-28 09:28:47 · 3072 阅读 · 0 评论 -
【python】folium 库绘制地图点击框
使用folium 库 生成 地图网页folium 官网import foliumimport pandas as pddef mark_map(data): """ 带有标注的地图 :param data: :return: """ # 地图制作 myMap = folium.Map(location=[20, 0], t...原创 2018-06-28 09:25:48 · 5012 阅读 · 6 评论 -
【python】获取航线信息并且制作成图
获取航线信息并且制作成图航线信息航线信息查询网站 本次实例使用的航班号为 CES5496 查询后在network中可以寻找到如下内容https://zh.flightaware.com/ajax/ignoreall/trackpoll.rvt?token=c35ca45ecbca57cd1ea443d1c65c36426ea06630de026ffd737977e4a40a26ead...原创 2018-07-03 14:26:03 · 6364 阅读 · 2 评论 -
【python】GeoJson和bokeh-1
GeoJson 文档{ "type": "FeatureCollection", "features": [ { "geometry": { "type": "Polygon", "coordinates": [ [ [ 3,原创 2018-07-16 15:04:59 · 422 阅读 · 0 评论 -
【python】geoJson转换Shapefile
#!/usr/bin/env python# -*- coding: utf-8 -*-# @File : C.py# @Author: huifer# @Date : 2018/5/22 0022import gdalimport osrimport ogra = { "type": "FeatureCollection", "features": [ ...原创 2018-05-24 11:29:39 · 6514 阅读 · 7 评论 -
【python】文本文档转换geoJson
数据文件如下 格式版本号=数据产生单位=国土资源部数据产生日期=2013-2-18坐标系=80国家大地坐标系几度分带=3投影类型=高斯克吕格计量单位=米带号=39精度=0.01转换参数=0,0,0,0,0,0,0[地块坐标]23,15824.35,2003-11,双桥乡地块3,面,I-50-77-(25),公共基础设施,@12,1,3753689.660000,3943...原创 2018-05-24 10:40:39 · 2397 阅读 · 0 评论 -
python网页数据转码方式
python 网页数据转码方式import requestsurl = 'http://stockapp.finance.qq.com/mstats/menu_childs.php?id=bd021311'print requests.get(url).text.decode("unicode-escape")原创 2017-12-16 16:45:52 · 1696 阅读 · 0 评论