我已经建立了lat和lon点的numpy meshgrid,其形状为1750 X 1750及其对应的数据(降雨)形状相同.我需要找到从形状文件读取的多边形内数据点的平均值. 如果我的逻辑是正确的,我必须找出落在多边形内的经纬度的索引,然后使用该索引过滤数据,然后进行平均. 我已经使用Inpolygon函数在Matlab中成功实现了它,但是我想用Python重写代码.我使用过matplotlib Path.contains_points函数,但是在numpy ND数组上不起作用. 有人可以建议一种合适的方...
我正在尝试学习matplotlib和映射函数Basemap,并试图在地图上绘制经纬度的简单列表.但是,将底图坐标转换为简单的x / y坐标会在地图上疯狂地绘制点,我不知道为什么.我的代码如下：locationData = parseFile(locationFile)fig = plt.figure(figsize=(8,8)) m = Basemap(projection='merc', resolution='h',llcrnrlat=49, llcrnrlon=-13.5,urcrnrlat=59.5, urcrnrlon=4) m.drawcoastlines() m.drawmapboundary(fill_color='aqua') m.fi...
我有一个包含纬度和经度的地方的数据框.想象一下,例如城市.df = pd.DataFrame([{'city':"Berlin", 'lat':52.5243700, 'lng':13.4105300},{'city':"Potsdam", 'lat':52.3988600, 'lng':13.0656600},{'city':"Hamburg", 'lat':53.5753200, 'lng':10.0153400}]);现在我试图让所有城市都在另一个城市的半径范围内.假设距离柏林500公里,距汉堡500公里等所有城市.我会通过复制原始数据帧并使用距离函数连接来完成此操作. 中间结果有点像这...
我在两个独立的数组中有纬度和经度坐标：a = np.array([71,75]) b = np.array([43,42])如何轻松找到由这些坐标组成的所有可能点？ 我一直在搞乱itertools.combinations：In [43]:list(itertools.combinations(np.concatenate([a,b]), r=2))Out[43]:[(71, 75), (71, 43), (71, 42), (75, 43), (75, 42), (43, 42)]但这对我不起作用,因为点(71,75)和(43,42)是纬度/纬度和经度/经度对. 我想拥有的是：Out[43]:[(71, 43), (71, 42), (7...
我正在绘制龙卷风的起点/终点位置. csv文件具有以下数据：TouchDownLat TouchDownLong LiftoffLat LiftoffLong 31.53 -97.15 31.74 -96.88 46.45 -100.67 46.67 -100.47 43.1 -83.85 43.17 -83.42等等… 我所做的是每个纬度和经度,并将其分成一个numpy数组,如下：import matplotlib.pyplot as plt import csv import numpy as np from mpl_toolkits.basemap import Base...
我有一个超过500点的列表,在纬度和经度中给出.这些点代表陨石坑,我想绘制这些陨石坑的热图.例如,我想要一个有很多陨石坑的区域被认为是“热”,而更少的陨石坑被认为是“冷”.我使用SciPy查看了KDE,并尝试在Mathematica中使用ListSliceDensityPlot3D,但我无法创建足够的图形. 我将每个点从纬度/经度转换为笛卡尔[x,y,z]坐标,并将它们绘制在球体的表面上,但我不知道如何获取点列表并计算给定的密度区域,然后将其绘制在3D表面上. 我的...
我有一个用起点和终点定义的线段：A: x1 = 10.7196405787775 y1 = 59.9050401935882B: x2 = 10.7109989561813 y2 = 59.9018650448204其中x定义经度,y定义纬度. 我也有一点意见：P: x0 = 10.6542116666667 y0 = 59.429105如何计算线段与点之间的最短距离？我知道如何在笛卡尔坐标中执行此操作,但不能在长/纬度坐标中执行此操作.解决方法:使用有用的Python地理编码库geopy,以及从Chris Veness’s geodesy formulae开始的大圆中点的...
我无法弄清楚python是否是插入我拥有的数据集的最佳选择**lat lon time**x1 y1 3:02(t1)x2 y2 3:05(t2)x3 y3 3:10(t3)x4 y4 3:13(t4) 我留下空间来强调我正在寻找的东西.我希望以1秒的间隔填补缺失的时间.所以我认为我需要做的是t4-t1以获得11个相等的空格,将时间转换为秒以使计算更容易.然后图表看起来像这样：**lat lon time**x1 y1 3:02(t1)a b 3:03c d 3:04x2 y2 ...
我有一个DataFrame,结果如下： 我想将那些坐标列从WGS84转换为Lon& Lat并最后在我的数据框中添加这些新列： 对于转换,我使用以下代码,但我认为应该有一个更好的方法,不将坐标列转换为列表并创建一个新的DataFrame.import pyproj as pp from mpl_toolkits.basemap import Basemap import pandas as pdcx =dfb.COORDENADA_X.tolist() cy =dfb.COORDENADA_Y.tolist()utm15_wgs84 = pp.Proj(init='epsg:32615')for ix, iy in zip(cx, c...
当给定每端点的纬度和经度时,我有一个短函数来计算线的中点.简单来说,当经度大于-90度或小于90度时,它可以正常工作.对于地球的另一半,它提供了一些随机的结果. 该代码是在http://www.movable-type.co.uk/scripts/latlong.html提供的javascript的python转换,并且似乎符合更正的版本here和here.当与两个stackoverflow版本进行比较时,我承认我不用C#或Java编码,但我可以’ t发现我的错误. 代码如下：#!/usr/bin/pythonimport mathdef ...