是否可以根据此代码创建热图

问题描述

是否可以将其与热图叠加？或者有没有办法用 cartopy 或 geopandas 创建热图？

可以在here找到数据。

data=pd.read_csv('Data gathered1.csv')
data
x=data['LONGITUDE']
y= data["LATITUDE"]

import matplotlib.pyplot as plt 
import numpy as np
import cartopy.crs as ccrs
import cartopy.io.img_tiles as cimgt
import io
from urllib.request import urlopen,Request
from PIL import Image

def image_spoof(self,tile): # this function pretends not to be a Python script
      url = self._image_url(tile) # get the url of the street map API
      req = Request(url) # start request
      req.add_header('User-agent','Anaconda 3') # add user agent to request
      fh = urlopen(req) 
     im_data = io.BytesIO(fh.read()) # get image
     fh.close() # close url
     img = Image.open(im_data) # open image with PIL
     img = img.convert(self.desired_tile_form) # set image format
     return img,self.tileextent(tile),'lower' # reformat for cartopy

cimgt.OSM.get_image = image_spoof # reformat web request for street map spoofing
osm_img = cimgt.OSM() # spoofed,downloaded street map

fig = plt.figure(figsize=(12,9)) # open matplotlib figure
ax1 = plt.axes(projection=osm_img.crs) # project using coordinate reference system (CRS) of street 
mapcenter_pt = [26.2271,-98.2087] # lat/lon hidalgo
zoom = 0.5 # for zooming out of center point
extent = [center_pt[1]-(zoom*2.0),center_pt[1]+(zoom*2.0),center_pt[0]-zoom,center_pt[0]+zoom] # 
adjust 
to zoom
ax1.set_extent(extent) # set extents
ax1.scatter(x,y,transform=ccrs.PlateCarree())
scale = np.ceil(-np.sqrt(2)*np.log(np.divide(zoom,350.0))) # empirical solve for scale based on zoom
scale = (scale<20) and scale or 19 # scale cannot be larger than 19
ax1.add_image(osm_img,int(scale)) # add OSM with zoom specification
plt.show()

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）