年少的我-CSDN博客

原创 Windows10安装ArcGIS10.2出现ArcGIS 10.2 for Desk Setup问题解决方法

安装ArcGIS10.2出现如下问题：可采用如下方法解决：1、右键Windows图标，找到Windows PowerShell，并以管理员权限运行Windows PowerShell。2、在命令行输入该命令：DISM.EXE /Online /Add-Capability /CapabilityName:NetFx3~~~~ /Source:C:\Windows10\sources\sxs现在再安装ArcGIS10.2就不会出现该问题。推荐一篇详细教程的博客，每一步都很详细：h.

2020-12-29 09:36:41 3677 1

原创 NumPy常用方法总结

NumPy是Python数值计算中最为重要的基础包，主要支持高维度数组与矩阵运算，也针对数组运算提供大量的数学函数库numpy简介； numpy常用方法；

2020-12-25 15:39:36 2325 1

原创 ArcGIS教程——以点云数据为输入计算太阳辐照度

Arcgis中的太阳辐射分析工具可以针对特定时间段太阳对某地理位置区域的影响进行分析。如果需要知道特定时刻或特定时间段内，太阳对单点的辐照度，可以选择太阳辐射点分析工具。以下以建筑物屋顶每个位置点进行计算。进行太阳辐射点分析，输入数据为.txt格式的点云坐标，首先在Arcgis中添加数据：接下来，要把txt数据转为shp格式的数据(需要插值为栅格）：导出dbase表，然后指定字...

2020-04-06 10:20:42 2817

原创 Numpy拼接两个数组

连接两个数组有很多方法，当处理数据量较大的时候就需要考虑哪种方式最快了。import numpy as npfrom timeit import repeat# 方法1def f1(): a = np.array([[1, 2, 3], [1, 2, 3]]) b = np.array([[1, 2, 3], [1, 2, 3]]) c = np.row_st...

2020-03-17 17:55:23 4008

原创 Python比较SciPy与sympy库计算定积分的快慢

Python计算定积分有两个库可以可以选择（就我目前所知道的），两个库计算定积分的精度和速度具有显著差异，以下我做了一个简单的对比，以计算函数f(x) = sin(x)*cos(x)在区间[0,pi/2]上的值为例：import mathfrom scipy import integrateimport sympyfrom sympy.abc import xfrom timeit ...

2020-03-16 20:23:12 4217 1

原创 python计算定积分

python有对应的包sympy，用来计算积分。在计算定积分的时候计算函数与积分上下限都需要用sympy库中的函数，否则会出错。下面代码计算的是值：import sympyfrom sympy.abc import x # 使用符号变量的时候，需要先导入符号value = sympy.integrate(sympy.sin(x)*sympy.cos(x), (x, 0, sympy....

2020-03-16 13:13:53 9091 3

原创 Python中使用 cProfile 和flameprof调优程序性能

有时候发现自己写的Python代码运行的很慢，但是又找不到到底是哪里导致程序运行速度缓慢，这个时候需要对程序进行性能分析。Python中的cProfile模块可以生成程序运行的火焰图，检测哪里效率低下，方便更好的修改程序代码。在使用之前需要安装flameprof库和cProfile库，安装的路径自己需要知道。pip install flameprofpip install cProfi...

2020-03-16 12:54:11 2653

原创一个文件夹下所有txt中的数据进行操作

import numpy as npimport ospath = r"F:\CSU_lib\Commercial\A\G" # 文件夹的路径result = [] filelist = os.listdir(path) # 文件夹中的txt文件名for file in filelist: file_name = file.split(".") # 对每一个文件进行处理 ...

2019-09-26 15:04:19 438

原创判断一点是否处于多个三角形内部（python）

二维平面上判断一个点是否在三角形内部有很多种方法。在这我采用的是向量的叉乘，简单来说就是：待计算点P与另外三个顶点A，B，C构成的向量PA，PB，PC两两之间的叉乘方向要一致。详细过程可以参考这篇博客https://www.cnblogs.com/TenosDoIt/p/4024413.html 直接贴代码，我写的这个函数可以一次判断P点是否位于多个三角...

2019-07-28 11:07:10 1077 2

原创二维空间中判断一点是否在三角形内部

平面中判断一个点是否在三角形内部有很多方法，其中一个比较快的方法是：利用叉乘计算PA，PB，PC，是否同向。代码如下：def Triangle(PointA,PointB,PointC,P): # PointA，PointB，PointC为三角形的三个顶点，P为待计算点 PExtend = np.tile(P,(PointA.shape[0],3)) # 扩展P为与A同等维度 ...

2019-07-16 11:05:41 756 1

原创每隔15min生成一个时刻

按时间间隔生成一个时间序列Time_list = [] Hour,Min,Value = 10.5,0,0 #设置起始时刻为10:30while Value < 30: #每隔30分钟生成一个值 if Min >= 60: Min -= 60 Hour += 1 Min += 15 V...

2019-01-15 10:38:01 804

原创三维线面交点

三维空间中线面交点有两种方法可以求出，我所用的是第二种方法：向量法。具体代码如下：def Line_intersection(P1, P2, coeffi): #P1、P2为存贮平面外两点xyz坐标的列表，coeffi为存贮平面法向量a、b、c、d（常数项）的列表ax+by+cZ+d=0 P12 = [P2[0]-P1[0],P2[1]-P1[1],P2[2]-P1[...

2019-01-02 19:29:40 3628

DDremote的博客