我正在绘制一幅insar图像,并想在其中添加一些地质信息。也就是说,我希望将我的图像(我所做的,一幅由matplotlib在纬度和经度坐标下制作的图片)与That 结合起来。
首先,我想使用 matplotlib工具包中的 basemap 来绘制地形基地图,因为我使用matplotlib绘制卫星图像。我在basemap模块中找到了一些与我想要的类似的东西:
画一个etopo浮雕图片:[http://matplotlib.org/basemap/users/geography.html][2] (对不起,但我没有名声发布作为初学者的图片)。
我想我只需要转移中心。因此,我在给定的示例中模拟了脚本,并对坐标做了一些修改:
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
# setup Lambert Conformal basemap.
# set resolution=None to skip processing of boundary datasets.
m = Basemap(projection='cyl',llcrnrlon=-118.5,llcrnrlat=33.15,urcrnrlon=-117.15,urcrnrlat=34.5,\
resolution='f')
m.shadedrelief()
plt.show()然后我看到了一张丑陋模糊的照片。
后来我在网上搜索,找出了原因。默认的图像大小是5400x2700,这使得它的绘制速度慢,在小规模地图中模糊。所以这似乎是一种不正确的方法。
这是否是为了使用matplotlib获得高分辨率的地形图?我必须开始学习GMT吗?请帮帮我..。
发布于 2018-05-08 21:48:56
尝试像这样的东西,但是从你的区域改变epsg。这个代码在http://spatialreference.org/ref/epsg/
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
map = Basemap(llcrnrlon=-76.5,llcrnrlat=7., urcrnrlon=-75.9,urcrnrlat=7.5, resolution = 'h', epsg=3115)
map.arcgisimage(service='Specialty/DeLorme_World_Base_Map', xpixels = 1500, verbose= True)
#map.arcgisimage(service='NatGeo_World_Map', xpixels = 2000, verbose= False)
plt.show()
plt.close()https://stackoverflow.com/questions/31713415
复制相似问题