我有描述旋转椭圆的数据(椭圆的中心在纬度经度坐标中,长轴和短轴的长度以公里为单位,以及椭圆的方向角)。我不知道焦点的位置,但假设有一种方法可以找到它们。我想确定一个特定的纬度经度点是否在这个椭圆内。我已经找到了一个很好的方法来确定一个点是否在笛卡尔网格上的椭圆内,但不知道如何处理纬度经度点。任何帮助都将不胜感激。
-Cody O.
我想要将经度/纬度转换为x/y地图投影坐标(以米为单位),但使用cartopy + pyplot而不是底图。 假设这是底图代码: >>> from mpl_toolkits.basemap import Basemap
>>> import numpy as np
>>> import matplotlib.pyplot as plt
>>> # read in topo data (on a regular lat/lon grid)
>>> etopo = np.loadtxt('etopo
我有一个网格化的温度数据和全国各地的气象站,以及它们的纬度和经度。我想找到离气象站最近的网格点。我的网格数据有坐标x,y,其中的纬度和经度是一个函数。
我发现,找到最近的网格点最简单的方法是首先将气象站的纬度和经度(Lat,Lon)转换为x和y值,然后找到最近的网格点。我为一个电台(lat=,lon= )做了以下操作:
import matplotlib.pyplot as plt
from netCDF4 import Dataset as netcdf_dataset
import numpy as np
from cartopy import config
import carto
我已经成功地从天然地球站点创建了一张地图(国家边界和海岸线),但我发现很难绘制天气的经纬度坐标,有些站点很难绘制到地图上。经度和纬度坐标是附加的CSV文件。
下面是目前为止编译的is代码,并生成了地图:
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import cartopy.feature as feature
import cartopy.io.shapereader as shapereader
在这里输入图像描述
countries = shapereader.natural_earth(resolution=&
我有一个有点的GeoAxesSubplot,还有潜在的线或多边形。这个投影可以是任何由拼图支持的投影,包括正字法。
我可以使用不同的转换来绘制图,如所解释的
from matplotlib import pyplot as plt
import cartopy.crs as ccrs
# Projection could be any, e.g. globe or Arctic Stereopolar...
ax = plt.axes(projection=ccrs.Mollweide())
ax.coastlines()
# Plot using the coordinate syste
我画了两组重叠的轴,一组是另一组的放大版。我想在放大的轴的角和它在大轴上所代表的矩形的角之间画线。然而,我画的线条稍微偏离了位置。我试着把它压缩成一个简单的例子:
import cartopy.crs as ccrs
import matplotlib.pyplot as plt
# Create a large figure:
fig = plt.figure(figsize=(10, 10))
# Add an axes set and draw coastlines:
ax1 = plt.axes([0.01, 0.49, 0.8, 0.5], projection=ccrs.Pl
是否有任何支持,允许轴标签极地,Cartopy?我注意到只支持使用ax.gridlines(draw_labels=True)的PlateCarree和Mercator。
此外,ax.set_xticks只允许从一个直角坐标系统转换到另一个直角坐标系统。
下面是一些比较Basemap和Cartopy的代码。我想知道是否有人发明了一种允许标签的解决办法。
import matplotlib.pyplot as plt
import cartopy.crs as ccrs
from mpl_toolkits.basemap import Basemap
import numpy as np
ax