def main():
lat, lon, radius = -66.10165712, 136.51515176, 55.957758952975269
filename = 'Torres-Strait.png'
land = earth.ReadShapes('shapes/ne_10m_land')
klass = crs.AzimuthalEquidistant
proj, land_projected = earth.Project(lat, lon, land, klass)
circle = earth.Circle(lat, lon, radius, proj)
pyplot.figure(figsize=(6, 6))
ax = pyplot.axes(projection=proj)
ax.set_extent([123, 153, -32, 0])
land_feature = feature.ShapelyFeature(land_projected,
proj,
edgecolor='none',
facecolor='black')
circle_feature = feature.ShapelyFeature([circle],
proj,
linewidth=1,
edgecolor='blue',
facecolor='none')
intersection_feature = feature.ShapelyFeature(
[g.intersection(circle) for g in land_projected],
proj,
linewidth=1,
edgecolor='yellow',
facecolor='none')
ax.add_feature(land_feature)
ax.add_feature(circle_feature)
ax.add_feature(intersection_feature)
ax.gridlines(edgecolor='#bbbbbb', linestyle='-', linewidth=0.5)
pyplot.savefig(filename, bbox_inches='tight', pad_inches=0.2)
def main():
land = earth.ReadShapes('shapes/ne_10m_land')
with open('small_circle_10m.txt', 'w') as log:
for lat0, lat1, lon0, lon1, lat_step in [
(-61, -63, 110, 113, 0.1),
(-64.3, -66.7, 135.0, 140.2, 0.1),
(-9, -22, -140, -155, 0.5)]:
earth.EvaluateForMeshPoints(
LargestCircleThroughOcean2,
lat0, lat1, lon0, lon1,
lat_step, land, log)
def main():
land = earth.ReadShapes('shapes/ne_10m_land')
with open('great_circle_10m_a.txt', 'w') as log:
for lat0, lat1, lon0, lon1 in [
(4, 8, -108, -118), (23, 25, 78, 81), (51, 64, 45, 50),
(49, 53, -3, -6), (42, 45, 153, 157), (76, 79, 98, 103),
(0, 7, -151, -161)]:
earth.EvaluateForMeshPoints(
earth.GreatCircleLandFraction,
lat0, lat1, lon0, lon1,
0.5, land, log)
print >>log
def main():
lat, lon, radius = 50, 20, 90
lands = earth.ReadShapes('shapes/ne_110m_land')
proj, lands_projected = earth.Project(lat, lon, lands)
pyplot.figure(figsize=(6, 6))
ax = pyplot.axes(projection=proj)
lands_feature = feature.ShapelyFeature(lands_projected,
proj,
edgecolor='none',
facecolor='black')
ax.add_feature(lands_feature)
ax.gridlines(edgecolor='#bbbbbb', linestyle='-', linewidth=1)
pyplot.savefig('Centre-in-Krakow.png', bbox_inches='tight')
def main():
land = earth.ReadShapes('shapes/ne_10m_land')
with open('great_circle_10m_b.txt', 'w') as log:
for lat0, lat1, lon0, lon1 in [(3.7, 6.2, -112.5, -116.5),
(23.5, 24.5, 78.5, 80.5),
(56.5, 59, 46, 49),
(50.5, 51.5, -3.5, -6.5),
(43, 44.5, 153.5, 155.5),
(76.5, 78.5, 98, 101),
(3.2, 5.2, -153, -156)]:
earth.EvaluateForMeshPoints(earth.GreatCircleLandFraction, lat0,
lat1, lon0, lon1, 0.1, land, log)
print >> log
def main():
land = earth.ReadShapes('shapes/ne_10m_land')
for x0 in [(5.2, -114.45), (24.1, 79.17), (57.5, 46.75), (50.7, -4.76),
(44.1, 154.79), (77.4, 100.5)]:
print optimize.minimize(LandFraction,
x0=x0,
args=(land, ),
method=METHOD)
print optimize.minimize(NegativeLandFraction,
x0=(3.4, -154.14),
args=(land, ),
method=METHOD)
for x0 in [(-65.9, 136.73), (-15.0, -149.64)]:
print optimize.minimize(NegativeCircleRadius,
x0=x0,
args=(land, ),
method=METHOD)
def main():
land = earth.ReadShapes('shapes/ne_110m_land')
with open('small_circle_110m.txt', 'w') as log:
earth.EvaluateForMeshPoints(earth.LargestCircleThroughOcean, -90, 90,
-180, +180, 1., land, log)