我已经建立了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 ...