- 博客(60)
- 问答 (4)
- 收藏
- 关注
RSS订阅原创 PermissionError: [Errno 13] Permission denied: ‘E:\\python\\sta\\dist\\tkinterexercise.exe.notanexec
打包成exe提示PermissionError: [Errno 13] Permission denied: ‘E:\python\sta\dist\tkinterexercise.exe.notanexecutable’要关闭所有保护电脑的防火墙啊什么病毒和威胁防护360啊等等。
2023-07-27 11:23:00
135
原创 vfp常用统计命令汇总（随时增加）
select A.value/B.s as 比率 from A,(select 类别,sum(value) as s from A group by 类别)B where A.类别 =B.类别。select 类别,sum(area)/(select sum(area) from 表名) as 比率 from 表名 group by 类别。select distinct 字段名1,字段名2,sum(面积),count() from 表名 group by 字段名1,字段名2。
2023-06-06 17:43:32
1111
原创 将度分秒换算为度
import pandas as pdimport redata = pd.read_excel(r'F:\2z\guangxi\经纬度.xlsx')y=[]for i in data['Latitud']: a=re.split('°|\′|\″|\'|\"|’',i)#按度分秒分割，用竖线间隔多个分隔符，部分符号需要使用反斜杠 a=[x for x in a if x!=' ' and x!=''] fh=int(a[0])/abs(int(a[0]))
2022-04-19 21:41:30
853
原创 R语言Tips
for (i in A){if (i <0.5 ){print(i)}}#for 循环里嵌套If木有接触过R要是写错了请大佬们指点
2022-04-06 16:56:38
430
原创 有关OSMnx
1安装安装之前要先安装geopandas，之前有些过相关安装方法pip install osmnx（partially initialized module ‘osmnx’ has no attribute ‘graph_from_point’2import osmnx G = ox.graph_from_bbox(24, 21, 107, 103, network_type ='drive') G_projected = ox.project_graph(G) ox.plot_graph(G_pr
2022-01-06 22:09:22
1147
原创 批量自然断点分类法
import pandas as pdimport jenkspy data = pd.read_excel('F:\\data.xls') #读取数据list1=data.columns.values #获取列索引名称c=list1[8:] #获取需要进行自然断点分类的数据的索引名for i in c: breaks = jenkspy.jenks_breaks(data[i], nb_class=5) #计算自然断点分类的数值 data['class_'+i]='' #添
2022-01-03 22:06:33
3703
11
原创 arcpy批量相交矢量数据
import arcpyfrom arcpy import envenv.workspace = "F:\1project\新建文件地理数据库.gdb"lists=["na13","na18","na8","pl18","pl13","pl8"]for i in lists:... arcpy.Intersect_analysis ([i, "village"], i+"xj", "ALL", "", "")
2021-12-20 14:48:48
646
原创 seaborn绘制上三角热力图
import seaborn as snsimport matplotlib.pyplot as pltimport pandas as pdimport numpy as npdf = pd.read_excel('F:/第一个时期.xls')data = df.corr()mask = np.zeros_like(data)mask[np.tril_indices_from(mask)] = True #根据官网的triu_是下三角推测u为upper,相应的掩膜去下三角就得是lower的
2021-10-05 20:10:25
1208
原创 python读取dta文件并转化为dataframe
码住，学习import osimport pandas as pddef Load_dta(filename): import sys reader = pd.read_stata(filename,iterator =True) df= pd.DataFrame() try: chunk = reader.get_chunk(100*1000) while len(chunk)>0: df=df.appen
2021-10-03 14:16:35
1770
原创 arcpy批量按区统计
env.workspace=("F:\\drive\\")>>> rasters=arcpy.ListRasters("*")>>> for factor in rasters:... indata="F:\\drive\\"+factor... out="F:\\drive0821\\"+str(factor)[0:4]+".dbf"... zone=ZonalStatisticsAsTable("nn","OBJECTID",inda
2021-08-21 20:26:25
614
原创 离线安装pyinstaller
先下载依赖的文件future ,pefile,altgraph,pywin32,pyinstaller-hooks-contrib,讲他们都直接pip安装，最后pip install 文件绝对路径
2021-08-21 12:54:22
774
1
原创 安装geopandas
1首先安装GDAL下载网址https://www.lfd.uci.edu/~gohlke/pythonlibs/#gdal，选择与自己Python相匹配包anaconda prompt 安装gdal2下载相应版本的Fiona，pip安装fiona3pip install geopandas 即可
2021-08-15 15:53:03
63
原创 arcpy合并shp文件
arcpy.Merge_management(["native8","plantation8"], "F:\\forest8.shp")
2021-07-20 14:36:22
489
原创 arcpy删除字段
arcpy.DeleteField_management("F:\\paper\\data\\forest\\native.shp",["ZMJ","POLY_ID"])
2021-07-19 22:44:36
470
原创 arcpy按掩膜提取
>>> import arcpy>>> from arcpy.sa import *>>> arcpy.env.workspace = "D:\\Downloads\\Downloads">>> >>> shppath="F:\\1project\\paper\\data\\sheng\\gx1.shp">>> rasterpath="D:\\Downloads\\Downloads\\gx
2021-07-08 18:24:20
633
原创 gis等高线注记掩膜
1右击——属性——标注——标注此图层中的要素；选择标注字段；选择标注位置——确定2右击——将标注转换为注记——转换3工具箱——制图工具——掩膜工具——要素轮廓线掩膜4数据框右击——高级绘制选项——勾选被掩膜图层5将生成的掩膜的图层关闭...
2021-07-06 20:16:35
648
原创 读取表格并写入
读取表格，按条件筛选并且写入表格条件是将一组浮点数据分组，并计数没有研究怎么追加，就直接新生成了一个表格，也许不太规范，但是就能用就成import pandas as pdN = pd.read_excel('花坪2.xlsx')list1=list(N['胸径'])list2=[]for i in list1: if int(i)<i<int(i)+1: list2.append(int(i)) if i== int(i): li
2021-06-15 22:36:50
61
原创 arcpy添加字段
import arcpy... from arcpy import env... env.workspace = "F:\\1project\\beibuwan\\data\\tdlxsg\\un2000\\"... shps= arcpy.ListFeatureClasses()... for shp in shps:... print shp... arcpy.AddField_management(shp, "area", "DOUBLE")#要注意层级关系，用TAb空行
2021-06-10 00:32:00
311
原创 驱动页面出图
1菜单栏右击——数据驱动页面——点击数据驱动页面设置——勾选启用数据驱动页面将标题与显示页面动态一致：点击页面文本——数据驱动页面名称数据框属性——裁剪至当前驱动页面范围，实现按区县范围显示导出地图2固定比例尺制图工具——数据驱动页面——格网索引要素设置数据驱动页面，范围要修改为100%插入数据框，将格网索引放在新数据框内，属性——范围指示器，红色的框即指示当前显示位置导出地图...
2021-05-31 18:08:57
483
原创 地图数据符号化
自定义——样式管理器——样式——创建新样式以新建填充符号为例选择填充符号，在右边空白处右击新建——填充符号，可以自定义可以自定义多个填充符号，在类别那里设置与属性表字段名相同，如性别，名称则命名为属性数据，如=男，女，可以实现自动匹配（属性——符号系统——类别——与样式中的符号匹配——选择新建的符号库——符号匹配）比如标记填充物——标记——编辑符号——选择字符标记（可以选择格网或随机）——点击确定，可以修改符号名称符号选择器里就出现了刚刚制作的样式...
2021-05-31 16:51:18
470
原创 文档跳转到指定内容
1将光标放到需要跳转的位置，插入——书签2在需要点击跳转的地方，右击——超链接——本文档中的位置，选择连接的书签，确定3生成一个蓝色带下划线的1，可以自定义字体
2021-05-28 18:10:38
863
原创 arcgis安装步骤及问题总结
安装了二十多次各种版本，总结一下遇到的问题，大致步骤如下：安装前，首先要确定防火墙、病毒防护中心和各种杀毒卫士已关闭，安装路径设置为英文。1解压ISO文件后双击ESRI.exe2点击Arcgis License Manager一栏的set up安装许可证，安装完毕后，在弹出的License Manager窗口，选择停止服务。3点击Arcgis for Desktop一栏的setup，一直下一步，直到出现administrator放着不管。4打开破解文件夹里的service.txt将第三行主机名改
2021-05-21 13:53:02
1950
原创 gis服务区分析
步行十分钟可达到的范围将交通路网进行打断——数据管理工具——要素——要素转线，即可将相交的线拆成一段一段右键数据——新建网络数据集加载network analyst——新建服务区打开网络分析窗口，设施点右击加载位置，加载locate,确定设置服务区属性，设置500米1000米的范围点击网络求解...
2021-05-16 19:02:17
3734
7
原创 arcmap打开nc数据
打开逐月平均气温数据多维工具——创建NetCDF栅格图层，加载nc文件，确定文件右击属性——值那里可以选择想要的某一月数据
2021-05-13 10:02:04
2991
原创 python绘制密度散点图
头大，外行人做个图咋这么难，趋势线还没有研究出来怎么加上去，哎import matplotlib.pyplot as plt from scipy.stats import gaussian_kdefrom mpl_toolkits.axes_grid1 import make_axes_locatableimport numpy as npimport pandas as pdfrom dbfread import DBFdata=DBF('native9.dbf')#数据量超级超级大，e
2021-05-04 20:03:26
4240
10
空空如也
wps里怎么插html格式文件
2022-06-03
自然断点法分类结果不一样
2022-01-05
pyreadstat读取stata数据类型是元组吗
2021-10-02
pyinstaller打包的exe在自己电脑里可以运行，其它电脑里就不行，而且文件很大
2021-08-23
TA创建的收藏夹 TA关注的收藏夹
TA关注的人