notebook.community

Edit and run 





In [2]:



    


%matplotlib inline

import shapely.geometry, shapely.wkt
from shapely.geometry.point import Point
from shapely.geometry.polygon import LinearRing
from shapely.geometry.linestring import LineString
from shapely.geometry.polygon import Polygon
import shapely as sl
import fiona
import numpy as np

import matplotlib
import matplotlib.pyplot as plt
import pylab

from utils.shapely_plot import draw

pylab.rcParams['figure.figsize'] = (17.0, 15.0)



    











In [3]:



    


shp_path = r'../data/OpenStreetMap/australia-oceania-latest-5030073410-water-lines.shp'



    











In [4]:



    


lengths = []
with fiona.collection(shp_path, "r") as input:
    for f in input:
        geom = sl.geometry.shape(f['geometry'])
        lengths.append(geom.length)



    











In [5]:



    


fig = plt.figure(figsize=(17, 17))
ax = plt.axes()
ax.set_aspect('equal', 'datalim')

with fiona.collection(shp_path, "r") as input:
    for f in input:
        geom = sl.geometry.shape(f['geometry'])
        draw(geom, outline='#aaaaff', alpha=0.5)
            
plt.show()



    


















    






---------------------------------------------------------------------------
KeyboardInterrupt                         Traceback (most recent call last)
<ipython-input-5-fe2820fda707> in <module>()
      6     for f in input:
      7         geom = sl.geometry.shape(f['geometry'])
----> 8         draw(geom, outline='#aaaaff', alpha=0.5)
      9 
     10 plt.show()

D:\src\GitHub\paper-osm-2016\paper-osm-2015\notebooks\utils\shapely_plot.pyc in draw(gs, fill, outline, lw, alpha)
     32             draw(g.geoms)
     33         else:
---> 34             draw_(g, fill, outline, lw, alpha)
     35 
     36 def plot(geoms):

D:\src\GitHub\paper-osm-2016\paper-osm-2015\notebooks\utils\shapely_plot.pyc in draw_(g, fill, outline, lw, alpha)
      9         plt.plot(g.x, g.y, 'k,')
     10     elif gType == 'LineString':
---> 11         x, y = g.xy
     12         plt.plot(x, y, 'b-', color=outline, lw=lw)
     13     elif gType == 'Polygon':

C:\Anaconda\lib\site-packages\shapely\geometry\linestring.pyc in xy(self)
    108           [0.0, 1.0]
    109         """
--> 110         return self.coords.xy
    111 
    112     def parallel_offset(

C:\Anaconda\lib\site-packages\shapely\coords.pyc in xy(self)
    173             lgeos.GEOSCoordSeq_getX(self._cseq, i, byref(temp))
    174             x.append(temp.value)
--> 175             lgeos.GEOSCoordSeq_getY(self._cseq, i, byref(temp))
    176             y.append(temp.value)
    177         return x, y

KeyboardInterrupt: 









    
























In [6]:



    


fig = plt.figure(figsize=(17, 2))
ax = plt.axes()

ax.set_xscale('log')
bins = 10 ** np.linspace(np.log10(0.0001), np.log10(0.5), 50)

#bins = 100


_ = plt.hist(lengths, bins)



    


















    
























In [7]:



    


def pair(list):
    '''Iterate over pairs in a list -> pair of points '''
    for i in range(1, len(list)):
        yield list[i-1], list[i]

for seg_start, seg_end in pair(geom.coords):
    line_start = Point(seg_start)
    line_end = Point(seg_end)
    segment = LineString([line_start.coords[0],line_end.coords[0]])

    print segment



    


















    






LINESTRING (150.6029044 -29.5390386, 150.6023251 -29.5395613)
LINESTRING (150.6023251 -29.5395613, 150.6014893 -29.5397144)
LINESTRING (150.6014893 -29.5397144, 150.5999433 -29.5398787)
LINESTRING (150.5999433 -29.5398787, 150.59857 -29.5398413)
LINESTRING (150.59857 -29.5398413, 150.5974327 -29.5394493)
LINESTRING (150.5974327 -29.5394493, 150.5964682 -29.5399011)
LINESTRING (150.5964682 -29.5399011, 150.5938718 -29.5399197)
LINESTRING (150.5938718 -29.5399197, 150.5924116 -29.54014)
LINESTRING (150.5924116 -29.54014, 150.5906961 -29.5411145)
LINESTRING (150.5906961 -29.5411145, 150.5895803 -29.5415252)
LINESTRING (150.5895803 -29.5415252, 150.5881845 -29.5415775)
LINESTRING (150.5881845 -29.5415775, 150.5875204 -29.5419172)
LINESTRING (150.5875204 -29.5419172, 150.5859539 -29.5441014)
LINESTRING (150.5859539 -29.5441014, 150.5854808 -29.5443777)
LINESTRING (150.5854808 -29.5443777, 150.5848156 -29.544471)
LINESTRING (150.5848156 -29.544471, 150.58413 -29.544176)
LINESTRING (150.58413 -29.544176, 150.5830561 -29.5432576)
LINESTRING (150.5830561 -29.5432576, 150.5823241 -29.543365)
LINESTRING (150.5823241 -29.543365, 150.5820931 -29.5434259)
LINESTRING (150.5820931 -29.5434259, 150.5815363 -29.5434861)
LINESTRING (150.5815363 -29.5434861, 150.5800603 -29.543589)
LINESTRING (150.5800603 -29.543589, 150.5794739 -29.5434819)
LINESTRING (150.5794739 -29.5434819, 150.5790389 -29.5431046)
LINESTRING (150.5790389 -29.5431046, 150.5788485 -29.54265)
LINESTRING (150.5788485 -29.54265, 150.5787819 -29.5415684)
LINESTRING (150.5787819 -29.5415684, 150.5786502 -29.5411666)
LINESTRING (150.5786502 -29.5411666, 150.5780438 -29.5404075)
LINESTRING (150.5780438 -29.5404075, 150.5775119 -29.540075)
LINESTRING (150.5775119 -29.540075, 150.5764068 -29.5397144)
LINESTRING (150.5764068 -29.5397144, 150.5755804 -29.5396274)
LINESTRING (150.5755804 -29.5396274, 150.5741634 -29.5396444)
LINESTRING (150.5741634 -29.5396444, 150.5725238 -29.539636)
LINESTRING (150.5725238 -29.539636, 150.568814 -29.5392372)
LINESTRING (150.568814 -29.5392372, 150.5655498 -29.5399611)
LINESTRING (150.5655498 -29.5399611, 150.5644609 -29.5405578)
LINESTRING (150.5644609 -29.5405578, 150.5639461 -29.5407207)
LINESTRING (150.5639461 -29.5407207, 150.5636207 -29.5408688)
LINESTRING (150.5636207 -29.5408688, 150.5634239 -29.5411633)
LINESTRING (150.5634239 -29.5411633, 150.5625248 -29.5425574)
LINESTRING (150.5625248 -29.5425574, 150.5613983 -29.5432854)
LINESTRING (150.5613983 -29.5432854, 150.5603839 -29.5433461)
LINESTRING (150.5603839 -29.5433461, 150.5582582 -29.5435279)
LINESTRING (150.5582582 -29.5435279, 150.5565751 -29.543491)
LINESTRING (150.5565751 -29.543491, 150.5546662 -29.5437659)
LINESTRING (150.5546662 -29.5437659, 150.5519558 -29.5448864)
LINESTRING (150.5519558 -29.5448864, 150.5497792 -29.5456793)
LINESTRING (150.5497792 -29.5456793, 150.5492736 -29.5460531)
LINESTRING (150.5492736 -29.5460531, 150.5480894 -29.5472746)
LINESTRING (150.5480894 -29.5472746, 150.5473742 -29.5477487)
LINESTRING (150.5473742 -29.5477487, 150.5464372 -29.5481247)
LINESTRING (150.5464372 -29.5481247, 150.5451575 -29.548352)
LINESTRING (150.5451575 -29.548352, 150.5446927 -29.5485778)
LINESTRING (150.5446927 -29.5485778, 150.5435908 -29.5499039)
LINESTRING (150.5435908 -29.5499039, 150.5431748 -29.5510257)
LINESTRING (150.5431748 -29.5510257, 150.5424136 -29.5517789)
LINESTRING (150.5424136 -29.5517789, 150.5412951 -29.552733)
LINESTRING (150.5412951 -29.552733, 150.538905 -29.5543691)
LINESTRING (150.538905 -29.5543691, 150.5381473 -29.5548877)
LINESTRING (150.5381473 -29.5548877, 150.5368614 -29.556216)
LINESTRING (150.5368614 -29.556216, 150.5364164 -29.5570079)
LINESTRING (150.5364164 -29.5570079, 150.5362059 -29.5580884)
LINESTRING (150.5362059 -29.5580884, 150.5356324 -29.5594509)
LINESTRING (150.5356324 -29.5594509, 150.5352963 -29.5605505)
LINESTRING (150.5352963 -29.5605505, 150.5353428 -29.5615087)
LINESTRING (150.5353428 -29.5615087, 150.535482 -29.5623339)
LINESTRING (150.535482 -29.5623339, 150.5355675 -29.5634939)
LINESTRING (150.5355675 -29.5634939, 150.5355059 -29.5643633)
LINESTRING (150.5355059 -29.5643633, 150.5345646 -29.5665393)
LINESTRING (150.5345646 -29.5665393, 150.5343026 -29.5669268)
LINESTRING (150.5343026 -29.5669268, 150.5336843 -29.5675378)
LINESTRING (150.5336843 -29.5675378, 150.5330551 -29.5678271)
LINESTRING (150.5330551 -29.5678271, 150.5320174 -29.5688942)
LINESTRING (150.5320174 -29.5688942, 150.530569 -29.5723284)
LINESTRING (150.530569 -29.5723284, 150.530254 -29.5748245)
LINESTRING (150.530254 -29.5748245, 150.5292811 -29.5757688)
LINESTRING (150.5292811 -29.5757688, 150.5257756 -29.5763314)
LINESTRING (150.5257756 -29.5763314, 150.5248355 -29.5768355)
LINESTRING (150.5248355 -29.5768355, 150.524645 -29.5774527)
LINESTRING (150.524645 -29.5774527, 150.5248695 -29.5791129)
LINESTRING (150.5248695 -29.5791129, 150.5254256 -29.5797646)
LINESTRING (150.5254256 -29.5797646, 150.5262241 -29.5821161)
LINESTRING (150.5262241 -29.5821161, 150.5268304 -29.5839193)
LINESTRING (150.5268304 -29.5839193, 150.5270828 -29.5848381)
LINESTRING (150.5270828 -29.5848381, 150.5273189 -29.5854374)
LINESTRING (150.5273189 -29.5854374, 150.5275871 -29.5862071)
LINESTRING (150.5275871 -29.5862071, 150.5279519 -29.5871261)
LINESTRING (150.5279519 -29.5871261, 150.5280807 -29.587476)
LINESTRING (150.5280807 -29.587476, 150.5281933 -29.5878212)
LINESTRING (150.5281933 -29.5878212, 150.5283435 -29.588325)
LINESTRING (150.5283435 -29.588325, 150.5283221 -29.588451)
LINESTRING (150.5283221 -29.588451, 150.5282148 -29.5887682)
LINESTRING (150.5282148 -29.5887682, 150.5274852 -29.5892813)
LINESTRING (150.5274852 -29.5892813, 150.5264499 -29.5897991)
LINESTRING (150.5264499 -29.5897991, 150.525495 -29.5902703)
LINESTRING (150.525495 -29.5902703, 150.5244275 -29.590816)
LINESTRING (150.5244275 -29.590816, 150.5241003 -29.590858)
LINESTRING (150.5241003 -29.590858, 150.5237784 -29.590858)
LINESTRING (150.5237784 -29.590858, 150.5229523 -29.5907694)
LINESTRING (150.5229523 -29.5907694, 150.5219384 -29.5906294)
LINESTRING (150.5219384 -29.5906294, 150.5212893 -29.5904475)
LINESTRING (150.5212893 -29.5904475, 150.5204417 -29.5901303)
LINESTRING (150.5204417 -29.5901303, 150.5201252 -29.5899437)
LINESTRING (150.5201252 -29.5899437, 150.5197444 -29.5895006)
LINESTRING (150.5197444 -29.5895006, 150.5194547 -29.5891787)
LINESTRING (150.5194547 -29.5891787, 150.518768 -29.5883903)
LINESTRING (150.518768 -29.5883903, 150.5182584 -29.5877792)
LINESTRING (150.5182584 -29.5877792, 150.5176469 -29.5872568)
LINESTRING (150.5176469 -29.5872568, 150.5170783 -29.5866923)
LINESTRING (150.5170783 -29.5866923, 150.5167081 -29.5863331)
LINESTRING (150.5167081 -29.5863331, 150.516118 -29.5860019)
LINESTRING (150.516118 -29.5860019, 150.515426 -29.585694)
LINESTRING (150.515426 -29.585694, 150.5145087 -29.5852741)
LINESTRING (150.5145087 -29.5852741, 150.5140259 -29.5852088)
LINESTRING (150.5140259 -29.5852088, 150.5137148 -29.5852648)
LINESTRING (150.5137148 -29.5852648, 150.5132695 -29.5854747)
LINESTRING (150.5132695 -29.5854747, 150.5129316 -29.5855307)
LINESTRING (150.5129316 -29.5855307, 150.5125185 -29.585568)
LINESTRING (150.5125185 -29.585568, 150.5123468 -29.5856333)
LINESTRING (150.5123468 -29.5856333, 150.5121644 -29.5857966)
LINESTRING (150.5121644 -29.5857966, 150.511923 -29.5860299)
LINESTRING (150.511923 -29.5860299, 150.511805 -29.5862864)
LINESTRING (150.511805 -29.5862864, 150.511746 -29.5866783)
LINESTRING (150.511746 -29.5866783, 150.5117782 -29.5869302)
LINESTRING (150.5117782 -29.5869302, 150.5119767 -29.587476)
LINESTRING (150.5119767 -29.587476, 150.512497 -29.5886142)
LINESTRING (150.512497 -29.5886142, 150.5125829 -29.5889688)
LINESTRING (150.5125829 -29.5889688, 150.5125614 -29.589202)
LINESTRING (150.5125614 -29.589202, 150.512497 -29.5893933)
LINESTRING (150.512497 -29.5893933, 150.5123468 -29.5895799)
LINESTRING (150.5123468 -29.5895799, 150.5118801 -29.589995)
LINESTRING (150.5118801 -29.589995, 150.5116709 -29.5902283)
LINESTRING (150.5116709 -29.5902283, 150.5115314 -29.5904335)
LINESTRING (150.5115314 -29.5904335, 150.51151 -29.5905921)
LINESTRING (150.51151 -29.5905921, 150.5115314 -29.590816)
LINESTRING (150.5115314 -29.590816, 150.5115851 -29.5910026)
LINESTRING (150.5115851 -29.5910026, 150.5116924 -29.5911845)
LINESTRING (150.5116924 -29.5911845, 150.5119391 -29.5915484)
LINESTRING (150.5119391 -29.5915484, 150.5128833 -29.59243)
LINESTRING (150.5128833 -29.59243, 150.5140634 -29.5932323)
LINESTRING (150.5140634 -29.5932323, 150.5143102 -29.5935029)
LINESTRING (150.5143102 -29.5935029, 150.5145302 -29.5938901)
LINESTRING (150.5145302 -29.5938901, 150.5145677 -29.5941233)
LINESTRING (150.5145677 -29.5941233, 150.5145355 -29.5943892)
LINESTRING (150.5145355 -29.5943892, 150.5143102 -29.5949256)
LINESTRING (150.5143102 -29.5949256, 150.5139508 -29.5956066)
LINESTRING (150.5139508 -29.5956066, 150.5133393 -29.5963343)
LINESTRING (150.5133393 -29.5963343, 150.5128672 -29.5970246)
LINESTRING (150.5128672 -29.5970246, 150.512658 -29.5973465)
LINESTRING (150.512658 -29.5973465, 150.51254 -29.597645)
LINESTRING (150.51254 -29.597645, 150.5122986 -29.5984846)
LINESTRING (150.5122986 -29.5984846, 150.5121215 -29.5990956)
LINESTRING (150.5121215 -29.5990956, 150.5119338 -29.5999539)
LINESTRING (150.5119338 -29.5999539, 150.5118426 -29.6005229)
LINESTRING (150.5118426 -29.6005229, 150.5115958 -29.6017543)
LINESTRING (150.5115958 -29.6017543, 150.5114188 -29.6027058)
LINESTRING (150.5114188 -29.6027058, 150.5113383 -29.6032142)
LINESTRING (150.5113383 -29.6032142, 150.511274 -29.6044175)
LINESTRING (150.511274 -29.6044175, 150.511231 -29.6046881)
LINESTRING (150.511231 -29.6046881, 150.5110916 -29.604884)
LINESTRING (150.5110916 -29.604884, 150.5107804 -29.6051918)
LINESTRING (150.5107804 -29.6051918, 150.5104693 -29.605341)
LINESTRING (150.5104693 -29.605341, 150.5102333 -29.6054343)
LINESTRING (150.5102333 -29.6054343, 150.5096163 -29.6055229)
LINESTRING (150.5096163 -29.6055229, 150.509155 -29.6054483)
LINESTRING (150.509155 -29.6054483, 150.508286 -29.6051358)
LINESTRING (150.508286 -29.6051358, 150.5073579 -29.6048047)
LINESTRING (150.5073579 -29.6048047, 150.5064192 -29.6046927)
LINESTRING (150.5064192 -29.6046927, 150.5056252 -29.6046414)
LINESTRING (150.5056252 -29.6046414, 150.5047186 -29.6046554)
LINESTRING (150.5047186 -29.6046554, 150.5038711 -29.6046741)
LINESTRING (150.5038711 -29.6046741, 150.5032756 -29.6046787)
LINESTRING (150.5032756 -29.6046787, 150.5027231 -29.604814)
LINESTRING (150.5027231 -29.604814, 150.5022886 -29.6049772)
LINESTRING (150.5022886 -29.6049772, 150.502133 -29.6051172)
LINESTRING (150.502133 -29.6051172, 150.5020096 -29.6052804)
LINESTRING (150.5020096 -29.6052804, 150.5019881 -29.6055649)
LINESTRING (150.5019881 -29.6055649, 150.502015 -29.6056908)
LINESTRING (150.502015 -29.6056908, 150.5020686 -29.6058634)
LINESTRING (150.5020686 -29.6058634, 150.502192 -29.60605)
LINESTRING (150.502192 -29.60605, 150.5023797 -29.6061573)
LINESTRING (150.5023797 -29.6061573, 150.502589 -29.6062319)
LINESTRING (150.502589 -29.6062319, 150.5027821 -29.6062785)
LINESTRING (150.5027821 -29.6062785, 150.5032327 -29.6063811)
LINESTRING (150.5032327 -29.6063811, 150.5034955 -29.606535)
LINESTRING (150.5034955 -29.606535, 150.5036887 -29.606675)
LINESTRING (150.5036887 -29.606675, 150.5038335 -29.6068662)
LINESTRING (150.5038335 -29.6068662, 150.5039837 -29.6070481)
LINESTRING (150.5039837 -29.6070481, 150.5040588 -29.607202)
LINESTRING (150.5040588 -29.607202, 150.5041607 -29.6074725)
LINESTRING (150.5041607 -29.6074725, 150.5042305 -29.6077803)
LINESTRING (150.5042305 -29.6077803, 150.5044236 -29.6087644)
LINESTRING (150.5044236 -29.6087644, 150.5045792 -29.609562)
LINESTRING (150.5045792 -29.609562, 150.5046757 -29.6101263)
LINESTRING (150.5046757 -29.6101263, 150.5047025 -29.6108912)
LINESTRING (150.5047025 -29.6108912, 150.5046427 -29.611674)
LINESTRING (150.5046427 -29.611674, 150.5044987 -29.6124255)
LINESTRING (150.5044987 -29.6124255, 150.5043002 -29.6129759)
LINESTRING (150.5043002 -29.6129759, 150.5038925 -29.613862)
LINESTRING (150.5038925 -29.613862, 150.5031254 -29.6148507)
LINESTRING (150.5031254 -29.6148507, 150.5026372 -29.6153497)
LINESTRING (150.5026372 -29.6153497, 150.5021813 -29.615886)
LINESTRING (150.5021813 -29.615886, 150.5016448 -29.6164876)
LINESTRING (150.5016448 -29.6164876, 150.5012908 -29.6167674)
LINESTRING (150.5012908 -29.6167674, 150.5009099 -29.6169446)
LINESTRING (150.5009099 -29.6169446, 150.5004915 -29.6172338)
LINESTRING (150.5004915 -29.6172338, 150.5001535 -29.6175835)
LINESTRING (150.5001535 -29.6175835, 150.4999282 -29.6180126)
LINESTRING (150.4999282 -29.6180126, 150.4994776 -29.6187354)
LINESTRING (150.4994776 -29.6187354, 150.4991182 -29.6195795)
LINESTRING (150.4991182 -29.6195795, 150.4987158 -29.6200971)
LINESTRING (150.4987158 -29.6200971, 150.4980989 -29.6207593)
LINESTRING (150.4980989 -29.6207593, 150.4974284 -29.6213143)
LINESTRING (150.4974284 -29.6213143, 150.4965379 -29.6218972)
LINESTRING (150.4965379 -29.6218972, 150.4950734 -29.62269)
LINESTRING (150.4950734 -29.62269, 150.4944082 -29.6230164)
LINESTRING (150.4944082 -29.6230164, 150.4927613 -29.62359)
LINESTRING (150.4927613 -29.62359, 150.4917743 -29.6238091)
LINESTRING (150.4917743 -29.6238091, 150.4913934 -29.6237858)
LINESTRING (150.4913934 -29.6237858, 150.4910501 -29.6236926)
LINESTRING (150.4910501 -29.6236926, 150.4909321 -29.6235853)
LINESTRING (150.4909321 -29.6235853, 150.4908194 -29.6233475)
LINESTRING (150.4908194 -29.6233475, 150.4907121 -29.623105)
LINESTRING (150.4907121 -29.623105, 150.490637 -29.6228299)
LINESTRING (150.490637 -29.6228299, 150.4903474 -29.6224801)
LINESTRING (150.4903474 -29.6224801, 150.4897948 -29.6220138)
LINESTRING (150.4897948 -29.6220138, 150.4893335 -29.621636)
LINESTRING (150.4893335 -29.621636, 150.4887058 -29.6213096)
LINESTRING (150.4887058 -29.6213096, 150.4882499 -29.621249)
LINESTRING (150.4882499 -29.621249, 150.4874828 -29.6212257)
LINESTRING (150.4874828 -29.6212257, 150.4867639 -29.6213003)
LINESTRING (150.4867639 -29.6213003, 150.4858466 -29.6213376)
LINESTRING (150.4858466 -29.6213376, 150.4843821 -29.6213656)
LINESTRING (150.4843821 -29.6213656, 150.4835238 -29.6212303)
LINESTRING (150.4835238 -29.6212303, 150.4824778 -29.6209039)
LINESTRING (150.4824778 -29.6209039, 150.4813727 -29.6205122)
LINESTRING (150.4813727 -29.6205122, 150.480155 -29.6200412)
LINESTRING (150.480155 -29.6200412, 150.4794039 -29.6197334)
LINESTRING (150.4794039 -29.6197334, 150.4779073 -29.6192484)
LINESTRING (150.4779073 -29.6192484, 150.4776122 -29.6190992)
LINESTRING (150.4776122 -29.6190992, 150.4773923 -29.6189219)
LINESTRING (150.4773923 -29.6189219, 150.4770919 -29.6184369)
LINESTRING (150.4770919 -29.6184369, 150.47677 -29.6177001)
LINESTRING (150.47677 -29.6177001, 150.4763677 -29.6172198)
LINESTRING (150.4763677 -29.6172198, 150.4761316 -29.6166135)
LINESTRING (150.4761316 -29.6166135, 150.4761638 -29.6162311)
LINESTRING (150.4761638 -29.6162311, 150.4759654 -29.6154243)
LINESTRING (150.4759654 -29.6154243, 150.4755308 -29.6142584)
LINESTRING (150.4755308 -29.6142584, 150.4752412 -29.6133816)
LINESTRING (150.4752412 -29.6133816, 150.4750105 -29.6130598)
LINESTRING (150.4750105 -29.6130598, 150.4746672 -29.6128033)
LINESTRING (150.4746672 -29.6128033, 150.4743024 -29.6126774)
LINESTRING (150.4743024 -29.6126774, 150.4736694 -29.6124209)
LINESTRING (150.4736694 -29.6124209, 150.4730417 -29.6121131)
LINESTRING (150.4730417 -29.6121131, 150.4723068 -29.6115627)
LINESTRING (150.4723068 -29.6115627, 150.4715719 -29.6109844)
LINESTRING (150.4715719 -29.6109844, 150.470837 -29.6104574)
LINESTRING (150.470837 -29.6104574, 150.4701611 -29.6100843)
LINESTRING (150.4701611 -29.6100843, 150.469748 -29.6098791)
LINESTRING (150.469748 -29.6098791, 150.4691579 -29.6097299)
LINESTRING (150.4691579 -29.6097299, 150.4687609 -29.6097112)
LINESTRING (150.4687609 -29.6097112, 150.4680367 -29.6097578)
LINESTRING (150.4680367 -29.6097578, 150.46757 -29.6098231)
LINESTRING (150.46757 -29.6098231, 150.4659178 -29.6098231)
LINESTRING (150.4659178 -29.6098231, 150.4650005 -29.6099257)
LINESTRING (150.4650005 -29.6099257, 150.4643782 -29.610019)
LINESTRING (150.4643782 -29.610019, 150.4639866 -29.6102102)
LINESTRING (150.4639866 -29.6102102, 150.4637345 -29.6104714)
LINESTRING (150.4637345 -29.6104714, 150.4636004 -29.6108772)
LINESTRING (150.4636004 -29.6108772, 150.4637398 -29.6115394)
LINESTRING (150.4637398 -29.6115394, 150.4638525 -29.6120478)
LINESTRING (150.4638525 -29.6120478, 150.4643407 -29.6130038)
LINESTRING (150.4643407 -29.6130038, 150.4648825 -29.6139972)
LINESTRING (150.4648825 -29.6139972, 150.4651239 -29.6147061)
LINESTRING (150.4651239 -29.6147061, 150.4653438 -29.615816)
LINESTRING (150.4653438 -29.615816, 150.4653116 -29.6161332)
LINESTRING (150.4653116 -29.6161332, 150.4651721 -29.6164176)
LINESTRING (150.4651721 -29.6164176, 150.4647913 -29.6168094)
LINESTRING (150.4647913 -29.6168094, 150.4645928 -29.6169913)
LINESTRING (150.4645928 -29.6169913, 150.4638632 -29.6176442)
LINESTRING (150.4638632 -29.6176442, 150.4633643 -29.6179146)
LINESTRING (150.4633643 -29.6179146, 150.462844 -29.6180825)
LINESTRING (150.462844 -29.6180825, 150.4625811 -29.6180685)
LINESTRING (150.4625811 -29.6180685, 150.4623236 -29.6179566)
LINESTRING (150.4623236 -29.6179566, 150.4620018 -29.6176115)
LINESTRING (150.4620018 -29.6176115, 150.4615082 -29.6172244)
LINESTRING (150.4615082 -29.6172244, 150.4607626 -29.6169073)
LINESTRING (150.4607626 -29.6169073, 150.460253 -29.616427)
LINESTRING (150.460253 -29.616427, 150.4596682 -29.6158673)
LINESTRING (150.4596682 -29.6158673, 150.459502 -29.6153543)
LINESTRING (150.459502 -29.6153543, 150.4592659 -29.6148973)
LINESTRING (150.4592659 -29.6148973, 150.4588314 -29.6145942)
LINESTRING (150.4588314 -29.6145942, 150.4586168 -29.6145149)
LINESTRING (150.4586168 -29.6145149, 150.4584559 -29.6145009)
LINESTRING (150.4584559 -29.6145009, 150.4582628 -29.6145102)
LINESTRING (150.4582628 -29.6145102, 150.4575815 -29.6147527)
LINESTRING (150.4575815 -29.6147527, 150.4563799 -29.6151631)
LINESTRING (150.4563799 -29.6151631, 150.4556771 -29.6153777)
LINESTRING (150.4556771 -29.6153777, 150.4534992 -29.6161705)
LINESTRING (150.4534992 -29.6161705, 150.4525711 -29.6165809)
LINESTRING (150.4525711 -29.6165809, 150.450479 -29.6178727)
LINESTRING (150.450479 -29.6178727, 150.4493739 -29.6184136)
LINESTRING (150.4493739 -29.6184136, 150.4485264 -29.6187914)
LINESTRING (150.4485264 -29.6187914, 150.4465952 -29.6197054)
LINESTRING (150.4465952 -29.6197054, 150.4460051 -29.6197707)
LINESTRING (150.4460051 -29.6197707, 150.445533 -29.6198966)
LINESTRING (150.445533 -29.6198966, 150.4451146 -29.6200365)
LINESTRING (150.4451146 -29.6200365, 150.4449268 -29.6201578)
LINESTRING (150.4449268 -29.6201578, 150.4448249 -29.6203117)
LINESTRING (150.4448249 -29.6203117, 150.4447337 -29.6205308)
LINESTRING (150.4447337 -29.6205308, 150.4446908 -29.6207547)
LINESTRING (150.4446908 -29.6207547, 150.444664 -29.6211464)
LINESTRING (150.444664 -29.6211464, 150.4448893 -29.6217107)
LINESTRING (150.4448893 -29.6217107, 150.4451039 -29.6224008)
LINESTRING (150.4451039 -29.6224008, 150.4454625 -29.6236175)
LINESTRING (150.4454625 -29.6236175, 150.4458012 -29.6252594)
LINESTRING (150.4458012 -29.6252594, 150.4461338 -29.6264065)
LINESTRING (150.4461338 -29.6264065, 150.4461016 -29.6275956)
LINESTRING (150.4461016 -29.6275956, 150.4460426 -29.6286774)
LINESTRING (150.4460426 -29.6286774, 150.4458871 -29.6295727)
LINESTRING (150.4458871 -29.6295727, 150.4455437 -29.6305053)
LINESTRING (150.4455437 -29.6305053, 150.4450824 -29.6319741)
LINESTRING (150.4450824 -29.6319741, 150.4447981 -29.632557)
LINESTRING (150.4447981 -29.632557, 150.4446157 -29.6327715)
LINESTRING (150.4446157 -29.6327715, 150.444326 -29.6329347)
LINESTRING (150.444326 -29.6329347, 150.4434677 -29.633331)
LINESTRING (150.4434677 -29.633331, 150.4424324 -29.6335315)
LINESTRING (150.4424324 -29.6335315, 150.4414561 -29.6335502)
LINESTRING (150.4414561 -29.6335502, 150.4407533 -29.6333823)
LINESTRING (150.4407533 -29.6333823, 150.4398575 -29.6331492)
LINESTRING (150.4398575 -29.6331492, 150.4379692 -29.6325197)
LINESTRING (150.4379692 -29.6325197, 150.4358663 -29.6321047)
LINESTRING (150.4358663 -29.6321047, 150.4352655 -29.6319182)
LINESTRING (150.4352655 -29.6319182, 150.4348793 -29.6316897)
LINESTRING (150.4348793 -29.6316897, 150.4346003 -29.6314472)
LINESTRING (150.4346003 -29.6314472, 150.4343053 -29.6308877)
LINESTRING (150.4343053 -29.6308877, 150.4337527 -29.6294468)
LINESTRING (150.4337527 -29.6294468, 150.4333182 -29.6285748)
LINESTRING (150.4333182 -29.6285748, 150.4332163 -29.6279593)
LINESTRING (150.4332163 -29.6279593, 150.4331841 -29.62622)
LINESTRING (150.4331841 -29.62622, 150.4330232 -29.6254879)
LINESTRING (150.4330232 -29.6254879, 150.4329213 -29.625278)
LINESTRING (150.4329213 -29.625278, 150.432712 -29.6250216)
LINESTRING (150.432712 -29.6250216, 150.4314139 -29.6236552)
LINESTRING (150.4314139 -29.6236552, 150.4303517 -29.62269)
LINESTRING (150.4303517 -29.62269, 150.4299333 -29.6224288)
LINESTRING (150.4299333 -29.6224288, 150.4294076 -29.6222516)
LINESTRING (150.4294076 -29.6222516, 150.4290482 -29.6221537)
LINESTRING (150.4290482 -29.6221537, 150.4287531 -29.6222096)
LINESTRING (150.4287531 -29.6222096, 150.4283454 -29.6222842)
LINESTRING (150.4283454 -29.6222842, 150.4280235 -29.6224428)
LINESTRING (150.4280235 -29.6224428, 150.4276105 -29.6227925)
LINESTRING (150.4276105 -29.6227925, 150.4267683 -29.6237159)
LINESTRING (150.4267683 -29.6237159, 150.425851 -29.6245599)
LINESTRING (150.425851 -29.6245599, 150.4243597 -29.6255951)
LINESTRING (150.4243597 -29.6255951, 150.4235282 -29.6261454)
LINESTRING (150.4235282 -29.6261454, 150.4223534 -29.6274138)
LINESTRING (150.4223534 -29.6274138, 150.4214835 -29.6286443)
LINESTRING (150.4214835 -29.6286443, 150.4208781 -29.6290598)
LINESTRING (150.4208781 -29.6290598, 150.4202076 -29.6293676)
LINESTRING (150.4202076 -29.6293676, 150.4198428 -29.6293815)
LINESTRING (150.4198428 -29.6293815, 150.4194673 -29.6292323)
LINESTRING (150.4194673 -29.6292323, 150.4188289 -29.6289759)
LINESTRING (150.4188289 -29.6289759, 150.4181798 -29.6289199)
LINESTRING (150.4181798 -29.6289199, 150.4176058 -29.6289526)
LINESTRING (150.4176058 -29.6289526, 150.4166939 -29.6291344)
LINESTRING (150.4166939 -29.6291344, 150.4160073 -29.6293582)
LINESTRING (150.4160073 -29.6293582, 150.4154815 -29.6295447)
LINESTRING (150.4154815 -29.6295447, 150.4148217 -29.630025)
LINESTRING (150.4148217 -29.630025, 150.4143496 -29.6305753)
LINESTRING (150.4143496 -29.6305753, 150.4140546 -29.6312048)
LINESTRING (150.4140546 -29.6312048, 150.4139259 -29.6319088)
LINESTRING (150.4139259 -29.6319088, 150.4135825 -29.6330279)
LINESTRING (150.4135825 -29.6330279, 150.4133626 -29.6334802)
LINESTRING (150.4133626 -29.6334802, 150.4130783 -29.6337227)
LINESTRING (150.4130783 -29.6337227, 150.4121824 -29.6340771)
LINESTRING (150.4121824 -29.6340771, 150.4117801 -29.6341843)
LINESTRING (150.4117801 -29.6341843, 150.4113134 -29.6341843)
LINESTRING (150.4113134 -29.6341843, 150.4102888 -29.6340957)
LINESTRING (150.4102888 -29.6340957, 150.409232 -29.633732)
LINESTRING (150.409232 -29.633732, 150.408202 -29.6332751)
LINESTRING (150.408202 -29.6332751, 150.4078265 -29.6330512)
LINESTRING (150.4078265 -29.6330512, 150.4075368 -29.6326456)
LINESTRING (150.4075368 -29.6326456, 150.4073598 -29.6318762)
LINESTRING (150.4073598 -29.6318762, 150.4073652 -29.6309063)
LINESTRING (150.4073652 -29.6309063, 150.4074993 -29.6304447)
LINESTRING (150.4074993 -29.6304447, 150.4076656 -29.630039)
LINESTRING (150.4076656 -29.630039, 150.4078158 -29.6298945)
LINESTRING (150.4078158 -29.6298945, 150.408363 -29.6294049)
LINESTRING (150.408363 -29.6294049, 150.4089852 -29.6289339)
LINESTRING (150.4089852 -29.6289339, 150.4094037 -29.6284303)
LINESTRING (150.4094037 -29.6284303, 150.4100581 -29.6277681)
LINESTRING (150.4100581 -29.6277681, 150.4106428 -29.6275397)
LINESTRING (150.4106428 -29.6275397, 150.4115065 -29.6272179)
LINESTRING (150.4115065 -29.6272179, 150.412059 -29.6270547)
LINESTRING (150.412059 -29.6270547, 150.4124775 -29.6268728)
LINESTRING (150.4124775 -29.6268728, 150.4126169 -29.6267143)
LINESTRING (150.4126169 -29.6267143, 150.412751 -29.6265231)
LINESTRING (150.412751 -29.6265231, 150.4128047 -29.6262853)
LINESTRING (150.4128047 -29.6262853, 150.4128101 -29.6260241)
LINESTRING (150.4128101 -29.6260241, 150.4126652 -29.6252874)
LINESTRING (150.4126652 -29.6252874, 150.412423 -29.6247457)
LINESTRING (150.412423 -29.6247457, 150.4121127 -29.6243594)
LINESTRING (150.4121127 -29.6243594, 150.4117908 -29.6240516)
LINESTRING (150.4117908 -29.6240516, 150.4113456 -29.6237765)
LINESTRING (150.4113456 -29.6237765, 150.4107448 -29.623562)
LINESTRING (150.4107448 -29.623562, 150.4101761 -29.6234547)
LINESTRING (150.4101761 -29.6234547, 150.4094144 -29.6234314)
LINESTRING (150.4094144 -29.6234314, 150.4085132 -29.6234221)
LINESTRING (150.4085132 -29.6234221, 150.4077568 -29.6234081)
LINESTRING (150.4077568 -29.6234081, 150.4069682 -29.6233195)
LINESTRING (150.4069682 -29.6233195, 150.4065391 -29.6231143)
LINESTRING (150.4065391 -29.6231143, 150.4062494 -29.6228625)
LINESTRING (150.4062494 -29.6228625, 150.4058041 -29.6224708)
LINESTRING (150.4058041 -29.6224708, 150.4053642 -29.6222609)
LINESTRING (150.4053642 -29.6222609, 150.4041251 -29.6218552)
LINESTRING (150.4041251 -29.6218552, 150.402387 -29.621207)
LINESTRING (150.402387 -29.621207, 150.4019203 -29.6209505)
LINESTRING (150.4019203 -29.6209505, 150.401636 -29.6207034)
LINESTRING (150.401636 -29.6207034, 150.4014429 -29.6202883)
LINESTRING (150.4014429 -29.6202883, 150.4011264 -29.6198826)
LINESTRING (150.4011264 -29.6198826, 150.4005631 -29.6194489)
LINESTRING (150.4005631 -29.6194489, 150.4000588 -29.6192344)
LINESTRING (150.4000588 -29.6192344, 150.3996458 -29.6192857)
LINESTRING (150.3996458 -29.6192857, 150.3993239 -29.6195562)
LINESTRING (150.3993239 -29.6195562, 150.3990396 -29.6199013)
LINESTRING (150.3990396 -29.6199013, 150.3988948 -29.6203583)
LINESTRING (150.3988948 -29.6203583, 150.3989323 -29.6209412)
LINESTRING (150.3989323 -29.6209412, 150.399163 -29.6217713)
LINESTRING (150.399163 -29.6217713, 150.3994473 -29.6226713)
LINESTRING (150.3994473 -29.6226713, 150.3996297 -29.6232962)
LINESTRING (150.3996297 -29.6232962, 150.3995975 -29.6236646)
LINESTRING (150.3995975 -29.6236646, 150.3994205 -29.6242195)
LINESTRING (150.3994205 -29.6242195, 150.3991308 -29.6245459)
LINESTRING (150.3991308 -29.6245459, 150.3986426 -29.624849)
LINESTRING (150.3986426 -29.624849, 150.3982027 -29.6250449)
LINESTRING (150.3982027 -29.6250449, 150.397838 -29.6251242)
LINESTRING (150.397838 -29.6251242, 150.3974517 -29.6251102)
LINESTRING (150.3974517 -29.6251102, 150.3972425 -29.6249516)
LINESTRING (150.3972425 -29.6249516, 150.3968134 -29.6241402)
LINESTRING (150.3968134 -29.6241402, 150.396116 -29.6226107)
LINESTRING (150.396116 -29.6226107, 150.3956505 -29.6215012)

















In [8]:



    


# get first geom

shp_path = r'../data/OpenStreetMap/australia-oceania-latest-5030073410-water-lines.shp'
    
g = None
with fiona.collection(shp_path, "r") as input:
    for f in input:
        g = sl.geometry.shape(f['geometry'])
        break

print(g.length)
g



    


















    






0.106613844241










    
Out[8]:























In [9]:



    


def getDistance(coord1, coord2):
    return Point(coord1).distance(Point(coord2))

def explodeLineString(lineString, segmentLength):
    length = lineString.length
    coords = lineString.coords
  
    # first point
    pt1 = coords[0]
    ipt2 = 1
    
    currentLength = 0.0
    
    segments = []
    segmentId = 0
  
    while currentLength < length and ipt2 < len(coords):
        segmentPoints = []
        segmentPoints.append(pt1)
        
        currentSegmentLength = 0.0
        
        # print('ipt2:{0}, N:{1} Lc:{2}, L:{3}'
        #      .format(ipt2, len(coords), currentLength, length))
        
        while currentSegmentLength < segmentLength:
            if ipt2 == len(coords):
                break
                
            pt2 = coords[ipt2]
            distance = getDistance(pt1, pt2)
      
            # print('d:{0}, Lc: {1} Ls: {2} ipt2: {3}'
            #      .format(distance, currentSegmentLength, segmentLength, ipt2))

            if currentSegmentLength + distance > segmentLength: # split
                # hypotenuse lengths
                h1 = segmentLength - currentSegmentLength # short
                h2 = distance # long
        
                # cathetus2
                dx2 = pt2[0] - pt1[0]
                dy2 = pt2[1] - pt1[1]

                # cathetus1
                ratio = h1 / h2
                x1 = pt1[0] + dx2 * ratio
                y1 = pt1[1] + dy2 * ratio

                currentSegmentLength += h1
        
                pt1 = [x1, y1]
                segmentPoints.append(pt1)
            else: # next
                currentSegmentLength += distance

                segmentPoints.append(pt2)
                pt1 = pt2
                ipt2 += 1
        
        currentLength += currentSegmentLength
        segmentId += 1
    
        segments.append(LineString(segmentPoints))
        
    return segments



    











In [10]:



    


print(getDistance(g.coords[0], g.coords[1]))
print(g.length)



    


















    






0.000502866682127
0.106613844241

















In [11]:



    


segments = explodeLineString(g, 0.02)



    











In [12]:



    


len(segments)



    


















    
Out[12]:








6

















In [13]:



    


shp_path = r'../data/OpenStreetMap/australia-oceania-latest-5030073410-water-lines.shp'
    
lengths_exploded = []
with fiona.collection(shp_path, "r") as input:
    for f in input:
        segments = explodeLineString(sl.geometry.shape(f['geometry']), 0.02)
        for s in segments:
            lengths_exploded.append(s.length)



    











In [14]:



    


print('Total length: ' + str(sum(lengths)))



    


















    






Total length: 576.87027753

















In [15]:



    


print(len(lengths))
print(len(lengths_exploded))



    


















    






9576
34936

















In [16]:



    


f, (ax1, ax2) = plt.subplots(1,2, figsize=(20,5))
ax1.set_xlim(0, 0.5)
_ = ax1.hist(lengths, 100)

ax2.set_xlim(0, 0.01)
_ = ax2.hist(lengths_exploded, 100)
plt.show()



    


















    
























In [18]:



    


f, (ax1, ax2) = plt.subplots(1,2, figsize=(20,5))

bins = 10 ** np.linspace(np.log10(0.0001), np.log10(0.5), 50)
_ = ax1.hist(lengths, bins)
ax1.set_xscale('log')

bins = 10 ** np.linspace(np.log10(0.0001), np.log10(0.5), 50)
_ = ax2.hist(lengths_exploded, bins)
ax2.set_xscale('log')

plt.show()



    


















    
























In [17]:



    


fig = plt.figure(figsize=(17, 17))
ax = plt.axes()
ax.set_aspect('equal', 'datalim')

# original
# draw(g, outline='#000000')

cmap = matplotlib.colors.ListedColormap (np.random.rand ( 256,3))

# segments
for (i, s) in enumerate(segments):
    draw(s, outline=cmap.colors[i])
            
plt.show()



    


















    
























In [19]:



    


from shapely.geometry import mapping

# take polylines and split them into segments of a given length

shp_path = r'../data/OpenStreetMap\australia-oceania-latest-5030073410-water-lines.shp'
shp_segments_path = r'../data/OpenStreetMap\australia-oceania-latest-5030073410-water-lines-segments.shp'

step = 0.02 # degrees

with fiona.collection(shp_path, "r") as input:
    schema = input.schema
    schema['properties']['segment'] = 'int'
    
    with fiona.open(shp_segments_path, 'w', 'ESRI Shapefile', schema) as output:
        for f in input:
            geom = sl.geometry.shape(f['geometry'])
            segments = explodeLineString(geom, step)
            
            for (i, s) in enumerate(segments):
                properties = f['properties']
                properties['segment'] = i
                output.write({'geometry': mapping(s), 'properties': properties})

Content source: gena/paper-osm-2015

Similar notebooks:

notebook.community | gallery | about