def plot_fence(fence, fence_points, subplot, extent=[13.0, 13.8, 52.3, 52.7]):
proj = cimgt.MapQuestOSM()
ax = plt.subplot(subplot[0], subplot[1], subplot[2], projection=proj.crs)
ax.add_feature(cartopy.feature.BORDERS, linestyle=':')
ax.add_feature(cartopy.feature.LAND)
ax.add_feature(cartopy.feature.COASTLINE)
ax.add_feature(cartopy.feature.OCEAN)
ax.add_feature(cartopy.feature.LAKES)
ax.set_extent(extent)
geoms = []
geo = Polygon(fence)
geoms.append(geo)
geoms2 = []
for cell in fence_points:
g = Polygon(cell)
geoms2.append(g)
ax.title.set_text("Number of hashes: " + str(len(geoms2)))
ax.add_geometries(geoms,
ccrs.PlateCarree(),
facecolor='red',
edgecolor='black',
alpha=1)
ax.add_geometries(geoms2,
ccrs.PlateCarree(),
facecolor='orange',
edgecolor='black',
alpha=0.2)
def test_web_tiles():
extent = [-15, 0.1, 50, 60]
target_domain = sgeom.Polygon([[extent[0], extent[1]],
[extent[2], extent[1]],
[extent[2], extent[3]],
[extent[0], extent[3]],
[extent[0], extent[1]]])
map_prj = cimgt.GoogleTiles().crs
ax = plt.subplot(3, 2, 1, projection=map_prj)
gt = cimgt.GoogleTiles()
gt._image_url = types.MethodType(ctest_tiles.GOOGLE_IMAGE_URL_REPLACEMENT,
gt)
img, extent, origin = gt.image_for_domain(target_domain, 1)
ax.imshow(np.array(img),
extent=extent,
transform=gt.crs,
interpolation='bilinear',
origin=origin)
ax.coastlines(color='white')
ax = plt.subplot(3, 2, 2, projection=map_prj)
qt = cimgt.QuadtreeTiles()
img, extent, origin = qt.image_for_domain(target_domain, 1)
ax.imshow(np.array(img),
extent=extent,
transform=qt.crs,
interpolation='bilinear',
origin=origin)
ax.coastlines(color='white')
ax = plt.subplot(3, 2, 3, projection=map_prj)
mq_osm = cimgt.MapQuestOSM()
img, extent, origin = mq_osm.image_for_domain(target_domain, 1)
ax.imshow(np.array(img),
extent=extent,
transform=mq_osm.crs,
interpolation='bilinear',
origin=origin)
ax.coastlines()
ax = plt.subplot(3, 2, 4, projection=map_prj)
mq_oa = cimgt.MapQuestOpenAerial()
img, extent, origin = mq_oa.image_for_domain(target_domain, 1)
ax.imshow(np.array(img),
extent=extent,
transform=mq_oa.crs,
interpolation='bilinear',
origin=origin)
ax.coastlines()
ax = plt.subplot(3, 2, 5, projection=map_prj)
osm = cimgt.OSM()
img, extent, origin = osm.image_for_domain(target_domain, 1)
ax.imshow(np.array(img),
extent=extent,
transform=osm.crs,
interpolation='bilinear',
origin=origin)
ax.coastlines()
def test_web_tiles():
extent = [-15, 0.1, 50, 60]
target_domain = shapely.geometry.Polygon([[extent[0], extent[1]],
[extent[2], extent[1]],
[extent[2], extent[3]],
[extent[0], extent[3]],
[extent[0], extent[1]]])
ax = plt.subplot(3, 2, 1, projection=ccrs.Mercator())
gt = cimgt.GoogleTiles()
img, extent, origin = gt.image_for_domain(target_domain, 1)
ax.imshow(np.array(img),
extent=extent,
transform=ccrs.Mercator(),
interpolation='bilinear',
origin=origin)
ax.coastlines(color='white')
ax = plt.subplot(3, 2, 2, projection=ccrs.Mercator())
qt = cimgt.QuadtreeTiles()
img, extent, origin = qt.image_for_domain(target_domain, 1)
ax.imshow(np.array(img),
extent=extent,
transform=ccrs.Mercator(),
interpolation='bilinear',
origin=origin)
ax.coastlines(color='white')
ax = plt.subplot(3, 2, 3, projection=ccrs.Mercator())
mq_osm = cimgt.MapQuestOSM()
img, extent, origin = mq_osm.image_for_domain(target_domain, 1)
ax.imshow(np.array(img),
extent=extent,
transform=ccrs.Mercator(),
interpolation='bilinear',
origin=origin)
ax.coastlines()
ax = plt.subplot(3, 2, 4, projection=ccrs.Mercator())
mq_oa = cimgt.MapQuestOpenAerial()
img, extent, origin = mq_oa.image_for_domain(target_domain, 1)
ax.imshow(np.array(img),
extent=extent,
transform=ccrs.Mercator(),
interpolation='bilinear',
origin=origin)
ax.coastlines()
ax = plt.subplot(3, 2, 5, projection=ccrs.Mercator())
osm = cimgt.OSM()
img, extent, origin = osm.image_for_domain(target_domain, 1)
ax.imshow(np.array(img),
extent=extent,
transform=ccrs.Mercator(),
interpolation='bilinear',
origin=origin)
ax.coastlines()
import matplotlib.pyplot as plt
from s2 import *
from shapely.geometry import Polygon
import cartopy.crs as ccrs
import cartopy.io.img_tiles as cimgt
proj = cimgt.MapQuestOSM()
plt.figure(figsize=(20, 20), dpi=200)
ax = plt.axes(projection=proj.crs)
ax.add_image(proj, 12)
ax.set_extent([-51.411886, -50.922470, -30.301314, -29.94364])
region_rect = S2LatLngRect(S2LatLng.FromDegrees(-51.264871, -30.241701),
S2LatLng.FromDegrees(-51.04618, -30.000003))
coverer = S2RegionCoverer()
coverer.set_min_level(8)
coverer.set_max_level(15)
coverer.set_max_cells(500)
covering = coverer.GetCovering(region_rect)
geoms = []
for cellid in covering:
new_cell = S2Cell(cellid)
vertices = []
for i in xrange(0, 4):
vertex = new_cell.GetVertex(i)
latlng = S2LatLng(vertex)
vertices.append((latlng.lat().degrees(), latlng.lng().degrees()))
geo = Polygon(vertices)
geoms.append(geo)
print "Total Geometries: {}".format(len(geoms))
ax.add_geometries(geoms,
ccrs.PlateCarree(),
