drewfustin / isocronut Goto Github PK
View Code? Open in Web Editor NEWFor a given geospatial location, calculate an isochrone (same time) contour around it.
isocronut's People
Contributors
Stargazers
Watchers
Forkers
akira1996 mwoytowich kaynewest ddd1600 ddclark quevedin gratianus wolframteetz kdmon ashenfad chrisschwer thetravelingsalesman dindinet moun spencer-scott jd-carroll patihe minsu-daniel-kim fototo oneslowpony fcpenha abk-chicago pilgrim2go dlmccaffrey profshen ankush23288 recerca zombiz pippo-sci blumonkey altophoto realmassive dairejpwalsh kostalihas bnomb3 mathematiguy xinyuez alptuzel gisdevelope allorimd alex4impact mayank06 millacy isallisocronut's Issues
Doubt on isochrones for 2 hours
Something's weird about the isochrones... the 90 and 120 minute isochrones are almost identical to the west although there is a highway in that direction. I strongly doubt that these are correct. Need to investigate the cause.
(html file with 60,90,120min isochrones and 48 directional vectors)
Unfortunately, my daily limit's gone again --
<!DOCTYPE html>
<html>
<head>
<meta name="viewport" content="initial-scale=1.0, user-scalable=no">
<meta charset="utf-8">
<title>Isochrone</title>
<style>
html, body, #map-canvas {
height: 100%;
margin: 0px;
padding: 0px
}
</style>
<script src="https://maps.googleapis.com/maps/api/js?v=3.exp&signed_in=true"></script>
<script>
function initialize() {
var mapOptions = {
zoom: 8,
center: new google.maps.LatLng(48.1089197,11.451958),
mapTypeId: google.maps.MapTypeId.TERRAIN
};
var map = new google.maps.Map(document.getElementById('map-canvas'), mapOptions);
var marker = new google.maps.Marker({
position: new google.maps.LatLng(48.1089197,11.451958),
map: map
});
var isochrone;
var isochroneCoords = [
new google.maps.LatLng(48.69084,11.45411),
new google.maps.LatLng(48.6821922,11.611889),
new google.maps.LatLng(48.6821922,11.611889),
new google.maps.LatLng(48.64269,11.78697),
new google.maps.LatLng(48.4982456,11.7956948),
new google.maps.LatLng(48.564644,11.9831995),
new google.maps.LatLng(48.5040218,12.0576007),
new google.maps.LatLng(48.4573868,12.1469998),
new google.maps.LatLng(48.28024,12.07553),
new google.maps.LatLng(48.3974082,12.7495727),
new google.maps.LatLng(48.2167837,12.1330527),
new google.maps.LatLng(48.1665,12.12216),
new google.maps.LatLng(48.1101058,12.1113211),
new google.maps.LatLng(48.0780006,11.9669667),
new google.maps.LatLng(48.0053526,12.0136777),
new google.maps.LatLng(47.9346692,12.014565),
new google.maps.LatLng(47.81522,12.19909),
new google.maps.LatLng(47.8640182,12.0093171),
new google.maps.LatLng(47.7172118,11.9563338),
new google.maps.LatLng(47.7505202,11.85938),
new google.maps.LatLng(47.74573,11.76186),
new google.maps.LatLng(47.7168545,11.719957),
new google.maps.LatLng(47.75063,11.5962),
new google.maps.LatLng(47.72323,11.52728),
new google.maps.LatLng(47.7058469,11.409518),
new google.maps.LatLng(47.62628,11.35382),
new google.maps.LatLng(47.6341305,11.233276),
new google.maps.LatLng(47.5975289,11.1884435),
new google.maps.LatLng(47.6645158,11.0703333),
new google.maps.LatLng(47.7863276,11.0604693),
new google.maps.LatLng(47.8739757,11.0244206),
new google.maps.LatLng(47.9118491,11.0575774),
new google.maps.LatLng(47.9139271,10.8490461),
new google.maps.LatLng(47.8855514,10.6467416),
new google.maps.LatLng(48.0051291,10.5949358),
new google.maps.LatLng(48.1005133,10.561436),
new google.maps.LatLng(48.17436,10.68676),
new google.maps.LatLng(48.2272487,10.7888558),
new google.maps.LatLng(48.3265549,10.6511387),
new google.maps.LatLng(48.4659587,10.5135426),
new google.maps.LatLng(48.5297439,10.6963314),
new google.maps.LatLng(48.5617904,10.6791623),
new google.maps.LatLng(48.56507,10.91982),
new google.maps.LatLng(48.5628425,11.0627245),
new google.maps.LatLng(48.57973,11.14955),
new google.maps.LatLng(48.5626401,11.2630624),
new google.maps.LatLng(48.5561561,11.3652545),
];
isochrone = new google.maps.Polygon({
paths: isochroneCoords,
strokeColor: '#000',
strokeOpacity: 0.5,
strokeWeight: 1,
fillColor: '#0000FF',
fillOpacity: 0.1
});
isochrone.setMap(map);
var isochrone2;
var isochroneCoords2 = [
new google.maps.LatLng(49.01559,11.63351),
new google.maps.LatLng(48.9352388,11.7886699),
new google.maps.LatLng(49.0347825,12.0631613),
new google.maps.LatLng(48.9854354,12.1959296),
new google.maps.LatLng(48.77459,12.23209),
new google.maps.LatLng(48.7879195,12.4686735),
new google.maps.LatLng(48.80049,12.83738),
new google.maps.LatLng(48.5510303,12.556838),
new google.maps.LatLng(48.3889649,12.4925531),
new google.maps.LatLng(48.31527,12.64295),
new google.maps.LatLng(48.1847363,12.4313311),
new google.maps.LatLng(48.1005422,12.4657629),
new google.maps.LatLng(48.0077528,12.4679199),
new google.maps.LatLng(47.87637,12.69482),
new google.maps.LatLng(47.7629182,12.6453054),
new google.maps.LatLng(47.7070091,12.4768688),
new google.maps.LatLng(47.6197486,12.3883036),
new google.maps.LatLng(47.5761214,12.2262449),
new google.maps.LatLng(47.5683328,12.062927),
new google.maps.LatLng(47.6892024,11.7704106),
new google.maps.LatLng(47.6446218,11.7433702),
new google.maps.LatLng(47.3445665,11.6890377),
new google.maps.LatLng(47.3445665,11.6890377),
new google.maps.LatLng(47.6825636,11.5765039),
new google.maps.LatLng(47.6825636,11.5765039),
new google.maps.LatLng(47.4433478,11.1744672),
new google.maps.LatLng(47.4617723,11.0053989),
new google.maps.LatLng(47.4035774,10.8806966),
new google.maps.LatLng(47.6316531,10.8201512),
new google.maps.LatLng(47.61059,10.71506),
new google.maps.LatLng(47.6366,10.54275),
new google.maps.LatLng(47.6644577,10.3069488),
new google.maps.LatLng(47.7724165,10.2581077),
new google.maps.LatLng(47.789715,9.8857952),
new google.maps.LatLng(47.9675891,9.948535),
new google.maps.LatLng(48.069478,9.9515912),
new google.maps.LatLng(48.223185,10.1784321),
new google.maps.LatLng(48.3497101,10.0558781),
new google.maps.LatLng(48.5405915,9.9226157),
new google.maps.LatLng(48.6223985,10.2391536),
new google.maps.LatLng(48.6641069,10.3632798),
new google.maps.LatLng(48.7518434,10.4708748),
new google.maps.LatLng(48.8600826,10.5008528),
new google.maps.LatLng(48.93889,10.71853),
new google.maps.LatLng(48.9532292,10.9188556),
new google.maps.LatLng(48.8405607,11.218338),
new google.maps.LatLng(49.2316261,11.2074204),
new google.maps.LatLng(49.09484,11.45046),
];
isochrone2 = new google.maps.Polygon({
paths: isochroneCoords2,
strokeColor: '#000',
strokeOpacity: 0.5,
strokeWeight: 1,
fillColor: '#000',
fillOpacity: 0.25
});
isochrone2.setMap(map);
var isochrone3;
var isochroneCoords3 = [
new google.maps.LatLng(49.26565,11.45254),
new google.maps.LatLng(49.25586,11.68265),
new google.maps.LatLng(49.216316,11.8367523),
new google.maps.LatLng(49.186964,12.1129601),
new google.maps.LatLng(49.10618,12.33497),
new google.maps.LatLng(49.02192,12.52515),
new google.maps.LatLng(48.9199723,12.6965073),
new google.maps.LatLng(48.8169289,12.8017088),
new google.maps.LatLng(48.6804397,12.9679205),
new google.maps.LatLng(48.5328439,13.0575411),
new google.maps.LatLng(48.3720899,13.1193429),
new google.maps.LatLng(48.2452541,13.1805614),
new google.maps.LatLng(48.097362,13.1838851),
new google.maps.LatLng(47.9448083,13.1724472),
new google.maps.LatLng(47.7967697,13.1164647),
new google.maps.LatLng(47.65502,13.03983),
new google.maps.LatLng(47.5251713,12.9483501),
new google.maps.LatLng(47.4035357,12.8251358),
new google.maps.LatLng(47.2857166,12.6596103),
new google.maps.LatLng(47.1893922,12.4765612),
new google.maps.LatLng(47.1427419,12.2762211),
new google.maps.LatLng(47.037923,12.100363),
new google.maps.LatLng(46.990808,11.9637515),
new google.maps.LatLng(46.9702609,11.594396),
new google.maps.LatLng(47.0028056,11.5042219),
new google.maps.LatLng(46.8684818,11.3173944),
new google.maps.LatLng(46.9917225,11.0124301),
new google.maps.LatLng(47.0211632,10.7792317),
new google.maps.LatLng(47.1202927,10.6461183),
new google.maps.LatLng(47.1202927,10.6461183),
new google.maps.LatLng(47.4098581,10.2797448),
new google.maps.LatLng(47.4656903,10.1044633),
new google.maps.LatLng(47.5211163,9.9694213),
new google.maps.LatLng(47.6475674,9.8293321),
new google.maps.LatLng(47.776013,9.7674612),
new google.maps.LatLng(47.9226998,9.7532823),
new google.maps.LatLng(48.0951473,9.7901525),
new google.maps.LatLng(48.2494496,9.7289785),
new google.maps.LatLng(48.39828,9.76904),
new google.maps.LatLng(48.5430757,9.7913518),
new google.maps.LatLng(48.6799436,9.9370948),
new google.maps.LatLng(48.825928,10.1004936),
new google.maps.LatLng(48.9176446,10.2259365),
new google.maps.LatLng(49.0195622,10.3680885),
new google.maps.LatLng(49.0872864,10.5405129),
new google.maps.LatLng(49.1751,10.77501),
new google.maps.LatLng(49.2322124,10.9798272),
new google.maps.LatLng(49.2533173,11.2388746),
];
isochrone3 = new google.maps.Polygon({
paths: isochroneCoords3,
strokeColor: '#000',
strokeOpacity: 0.5,
strokeWeight: 1,
fillColor: '#000',
fillOpacity: 0.25
});
isochrone3.setMap(map);
}
google.maps.event.addDomListener(window, 'load', initialize);
</script>
</head>
<body>
<div id="map-canvas"></div>
</body>
</html>
Can't start from local on my PC. The project might update "map"
Can't start from local on my PC.
I have config api key as well in cfg file, however, I still see above.
I think the map should be updated.
Javascript implementation issues
If the complexity of the requested isochrones is high, Chrome sais the script crashes and Firefox repeatedly throws unresponsive script errors. Would be great if we could fix this.
Example
var m = [20, 50, 100, 300, 400, 500, 600, 700, 800, 900, 1000];
var s = [20, 50, 100, 100, 100, 100, 100, 100, 100, 100, 100];
(Complete "time-out-problematic" html without token, at least on my mac it times out)
<!DOCTYPE html>
<html>
<head>
<meta name="viewport" content="initial-scale=1.0, user-scalable=no">
<meta charset="utf-8">
<title>Isochrone</title>
<script src="https://api.tiles.mapbox.com/mapbox.js/v2.2.2/mapbox.js"></script>
<script src="https://api.tiles.mapbox.com/mapbox.js/plugins/turf/v2.0.0/turf.min.js" charset="utf-8"></script>
<script src="https://ajax.googleapis.com/ajax/libs/jquery/1.9.1/jquery.min.js"></script>
<link href='https://api.tiles.mapbox.com/mapbox.js/v2.2.2/mapbox.css' rel='stylesheet' />
<style>
body { margin:0; padding:0; }
#map { position:absolute; top:0; bottom:0; width:100%; }
</style>
</head>
<body>
<div id="map"></div>
<script>
var accessToken = ' $$$ TOKEN $$$ ';
L.mapbox.accessToken = accessToken;
var map = L.mapbox.map('map', 'mapbox.streets')
.setView([48.10892,11.451958], 15);
map.on('click', function(e) {
create_isochrones(map, [e['latlng']['lat'], e['latlng']['lng']]);
});
function create_isochrones(map, origin) {
map.setView(origin, 6); // Zoom out
var well = {
type: 'Feature',
geometry: {
type: 'Point',
coordinates: [
origin[0],
origin[1]
]
},
properties: {}
};
L.mapbox.featureLayer(turf.flip(well)).addTo(map);
var m = [20, 50, 100, 300, 400, 500, 600, 700, 800, 900, 1000];
var s = [20, 50, 100, 100, 100, 100, 100, 100, 100, 100, 100];
for (var i in m) {
if (i > 0) {
var circle = turf.buffer(well, m[i], 'miles');
var extent = turf.extent(circle);
var additions = turf.random('points', s[i], {bbox: extent});
points['features'].push.apply(points['features'], additions['features']);
} else {
var circle = turf.buffer(well, m[i], 'miles');
var extent = turf.extent(circle);
var points = turf.random('points', s[i], {bbox: extent});
};
};
var deferreds = [];
var destinations = {
type: 'FeatureCollection',
features: []
};
for (i = 0; i < points.features.length; i++) {
var url = 'https://api.tiles.mapbox.com/v4/directions/mapbox.driving/' + well.geometry.coordinates[1].toString() + ',' + well.geometry.coordinates[0].toString() + ';' + points.features[i].geometry.coordinates[1].toString() + ',' + points.features[i].geometry.coordinates[0].toString() + '.json?access_token=' + L.mapbox.accessToken;
deferreds.push(
$.ajax({
url: url,
contentType: 'text/plain',
type: 'GET',
xhrFields: {
withCredentials: false
},
})
.success(function(data) {
var min_dur = 9999999;
for (r = 0; r < data['routes'].length; r++) {
min_dur = Math.min(min_dur, data['routes'][r]['duration'] / 60);
};
if (data['destination']) {
var d = {
type: 'Feature',
geometry: {
type: 'Point',
coordinates: [
data['destination']['geometry']['coordinates'][1],
data['destination']['geometry']['coordinates'][0]
]
},
properties: {
'z': min_dur
}
};
destinations['features'].push(d);
};
})
.fail(function() {
alert("Ajax failed to fetch data")
})
);
}
$.when.apply($, deferreds).done(
function() {
var breaks = [10, 20, 30, 60, 90, 120, 180, 240, 300, 420, 540, 660];
var colors = {};
colors[breaks[0]] = 'rgb(255, 255, 255)';
colors[breaks[1]] = 'rgb(255, 255, 191)';
colors[breaks[2]] = 'rgb(255, 191, 191)';
colors[breaks[3]] = 'rgb(255, 191, 127)';
colors[breaks[4]] = 'rgb(255, 127, 127)';
colors[breaks[5]] = 'rgb(255, 127, 63)';
colors[breaks[6]] = 'rgb(255, 63, 63';
colors[breaks[7]] = 'rgb(255, 63, 0)';
colors[breaks[8]] = 'rgb(255, 0, 0)';
colors[breaks[9]] = 'rgb(191, 0, 0)';
colors[breaks[10]] = 'rgb(127, 0, 0)';
colors[breaks[11]] = 'rgb(63, 0, 0)';
colors[breaks[12]] = 'rgb(0, 0, 0)';
var isolined = turf.isolines(destinations, 'z', 500, breaks);
for (i = isolined.features.length - 1; i >= 0; i--) {
var polygon_options = {
color: '#000',
opacity: 0.5,
weight: 1,
// fillColor: colors[isolined.features[i].properties.z],
fillColor: 'rgb(0, 0, 0)',
fillOpacity: 0.1
};
L.polygon(isolined.features[i].geometry.coordinates, polygon_options).addTo(map);
}
}
);
};
</script>
</body>
</html>
Cannot run your example - ConfigParser.NoSectionError: No section: 'api'
Import checks
(...)
import hmac
import base64
import urlparse
import ConfigParser
import simplejson
import urllib2
import time
import datetime
import mathimport isocronut
origin = '111 W Washington, Chicago'
duration = 10isochrone = isocronut.get_isochrone(origin, duration)
Traceback (most recent call last):
File "", line 1, in
File "isocronut.py", line 325, in get_isochrone
origin_geocode = geocode_address(origin, access_type, config_path)
File "isocronut.py", line 170, in geocode_address
key = config.get('api', 'api_number')
File "/usr/local/Cellar/python/2.7.10_2/Frameworks/Python.framework/Versions/2.7/lib/python2.7/ConfigParser.py", line 607, in get
raise NoSectionError(section)
ConfigParser.NoSectionError: No section: 'api'
Any idea?
unknown0023dffeacbb:isocronut-master sir$ python --version
Python 2.7.10
Configparser error after 2to3 conversion
Hi Drew,
I was really hoping to get this working, as it looks like a great solution. Alas, after much trial and error and small adjustments to get this running with Python 3.7, I have hit a wall:
Traceback (most recent call last): File "<stdin>", line 1, in <module> File "C:\Program Files (x86)\Python37-32\Scripts\isocronut.py", line 393, in generate_isochrone_map origin_geocode = geocode_address(origin, access_type, config_path) File "C:\Program Files (x86)\Python37-32\Scripts\isocronut.py", line 170, in geocode_address key = config.get('api', 'api_number') File "C:\Program Files (x86)\Python37-32\lib\configparser.py", line 780, in get d = self._unify_values(section, vars) File "C:\Program Files (x86)\Python37-32\lib\configparser.py", line 1146, in _unify_values raise NoSectionError(section) from None configparser.NoSectionError: No section: 'api'
Any ideas what's going on there? I'm stumped.
All the best,
Johnny
Map marker not within isochrone
Hey Drew,
thanks for realeasing this great project! I just have one question: Sometimes isonlines are not calculated based on the marker but close to the coordinates of the marker. Do you have an idea how to tackle this problem?
UnicodeEncodeError when using outside of US
For the error see e.g. http://www.markhneedham.com/blog/2015/05/21/python-unicodeencodeerror-ascii-codec-cant-encode-character-uxfc-in-position-11-ordinal-not-in-range128/
import isocronut
origin = 'Martinsried'
duration = 60
isochrone = isocronut.generate_isochrone_map(origin, duration)
Traceback (most recent call last):
File "", line 1, in
File "isocronut.py", line 399, in generate_isochrone_map
iso = get_isochrone(origin, duration, number_of_angles, tolerance, access_type, config_path)
File "isocronut.py", line 369, in get_isochrone
iso[i] = geocode_address(data[0][i], access_type, config_path)
File "isocronut.py", line 181, in geocode_address
key))
UnicodeEncodeError: 'ascii' codec can't encode character u'\xfc' in position 1: ordinal not in range(128)
Exception: Error. Google Maps API return status: ZERO_RESULTS for geocode_address
Every few hundred location calls I get this error.
import sys
reload(sys)
<module 'sys' (built-in)>
sys.setdefaultencoding('utf8')
origin = 'Martinsried'
duration = 90
number_of_angles=48
tolerance = 5
tolerance = 2
import isocronut
isochrone = isocronut.generate_isochrone_map(origin, duration, number_of_angles, tolerance)
Iteration 0
Iteration 1
Iteration 2
Iteration 3
Iteration 4
Iteration 5
Iteration 6
Traceback (most recent call last):
File "", line 1, in
File "isocronut.py", line 413, in generate_isochrone_map
File "isocronut.py", line 383, in get_isochrone
for i in range(number_of_angles):
File "isocronut.py", line 222, in geocode_address
if not d['status'] == 'OK':
Exception: Error. Google Maps API return status: ZERO_RESULTS
See related issues :
I think one needs to provide a radius iteration to find a place with a name nearby sometimes. See e.g.
http://stackoverflow.com/questions/11754114/why-is-google-place-api-giving-me-zero-results
Google APIs
Which Google APIs do I need to activate to use the service?
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.