After I remove the water from openstreetmap, I lost the administrative boundary and properties. How to save them? - openstreetmap

now I am handling the indonesia's map data from openstreetmap. I get the region datas whose admin_level equals "4". Then I remove the water part through this way,https://stackoverflow.com/questions/25297811/how-can-i-remove-water-from-openstreetmap-ways.But I lose the administrative boundary,either.Who can give me some suggestion to resave them.
My target is just to remove the water part from the original data.
Firstly, I get the pbf data of Indonesia from openstreetmap.
secondly, I get the multipolygons from the pbf through a python script:
thirdly, I remove the water part by ogr2ogr.
ogr2ogr -f GeoJSON ID_filter_level_44.geo.json -clipsrc ID_level_4.geo.json land_polygons.id.shp
And now, I get the file ID_filter_level_44.geo.json, but it's geometry is all "polygon" type. and some diffrent administration are merged.

Related

Importing data from postgres to cytoscape

I have been trying to load some gis data from a postgis database into Cytoscape 3.6. I am trying to get some inDegree and outDegree values I have used the sif file format.
As long as the data is written out in the follow format
source_point\tinteracts with\ttarget_point
Cytoscape is happy to read it.
I am just wondering if there is anyway of including my own metric for the cost of getting between source_point and target_point
Sure! There are several ways to read in text files into Cytoscape -- SIF is just one of them. I would create a file that looks like SIF, but is actually a more complete text file:
Source\tTarget\tScore
source_point\ttarget_point\t1.1
...
And then use the "File->Import Network->File", choose your source and target and leave score as an edge attribute. You can have as many attributes on each line as you want, and can even mix edge attributes, source node attributes, and target node attributes.
-- scooter

Change osm bonus pack input file

I am using osmbonuspack library in my project to find the route between a departure and a destination node and I want to save information related to the traversal time of the edges. So I keep a file which contains the nodes and edges of a specific city (OSM data). My problem is that for a specific route, osmbonuspack returns a road that doesn't agree with my file. I mean that when I search for two nodes of the route in my file, they exist, but there isn't an edge between them, while osmbonuspack has found an edge between them.
I know that this may is happening because I keep information only for a city, so this edge in my file is probably split into more nodes and edges.
But is there any way that I can change an input file or url in osmbonuspack so that it searches for a route in my .osm file?
No way: OSMBonusPack only provides access to online route services.
For offline routing using local map files, you can try GraphHopper Android.

Creating POI on the routes in SUMO

Is it possible to make POIs on the routes of SUMO automatically?
I want to have the locations of Base stations along the routes that cars or pedestrians move. I have manhattan.net.xml and manhattan.rou.xml as network file and route file in SUMO. I want to create the POI file that has POIs along side the routes with some distance from street.
Is it possible to create this POI file automatically?
Thanks.
No, it is not possible to generate this directly from sumo or the sumo-gui.
If you are willing to write a small python script however, you get a lot of support functions in the tools directory of your sumo installation. There is even a script which does almost what you want: tools/route/tracegenerator.py. This one generates points along the route. You will have to add the offset and adapt the output format to generate POIs. Maybe write to the mailing list for more help: http://sumo.dlr.de/wiki/Contact

After osm2pgrouting disappeared street segments

I'm download osm-file with some city. In QGIS this map looks good, all streets is connected. Next, i execute osm2pgrouting and in QGIS (when i connect to postgres db) i watch map without some segments in street. Therefore, find path not worked (i think so). Because of what the problem is and how to solve it?
And more, after i executed osm2pgrouting in QGIS i watch many lines from what that point (like sun, extensively). I'm define coordinates of point and delete all ways with this coordinates.
UPDATE. I'm download two maps. First map is contain town (Russia, Krasnoyarsk) and i got it from openstreetmap.org, chosen area and export it. Second map contain region. I got it from some gis-site, but i forgot the address. Unfortunately I can not post a picture for an example. In map in postgress missing half segments of all. Some hightways segments is missing, but some is not.
UPDATE2. I have osm2pgrouting log:
Trying to load data
Trying to parse data
Reference nd=1162932299 has no corresponding Node Entry (Maybe Node entry after Reference?)
Reference nd=1162931037 has no corresponding Node Entry (Maybe Node entry after Reference?)
...
Split ways
Dropping tables...
Creating tables...
Nodes table created
2create ways failed:
Types table created
Way_tag table created
Relations table created
Relation_ways table created
Classes table created
Adding tag types and classes to database...
Adding relations to database...
Adding nodes to database...
Adding ways to database...
Creating topology...
size of streets: 62931
size of splitted ways : 97921
finished
And i use with osm2pgrouting standart mapconfig.xml.
I noticed that line from the some point and the other end is at the beginning of the missing segment. I think, means that this line is missing segments. But osm file open in QGIS and is correct.

Exporting NetLogo data to graph with nodes and edges

I have created some links between agents (turtles) in NetLogo. This links will change at each time step. My aim is to export this data (i.e., turtles and links b/w them) to graph with vertices (turtles) edges (links), which can be given as input to Gephi. Is it possible to see the changes which occurs in netlogo in the graph when it is linked with Gephi. Can someone help me out. Thanks.
To export your network data in a format usable by Gephi, I would suggest using the nw:save-graphml primitive from NetLogo's NW Extension. This will give produce a file in the GraphML file format, which Gephi can read.
I guess you could re-save your network at each time step and overwrite your file, but I'm not sure if Gephi can display your changes dynamically. And depending on the size of your network, it might be slow.
Are you trying to use Gephi to see how the network changes over time, in a changing network that is generated by NetLogo? That's what #NicolasPayette's answer suggests, so I'll make the same assumption.
Gephi can display "dynamic graphs", i.e. networks that change over time. My understanding is that are two file formats that allow Gephi to import dynamic graphs: GEXF, and a special CSV (comma-separated) format that Gephi calls "Spreadsheet". Nicolas mentioned GraphML, which is a very nice network data format, but it doesn't handle dynamic graphs. And as far as I know, NetLogo doesn't generate GEXF or Gephi's "Spreadsheet" format.
However, the Gephi Spreadsheet format is very simple, and it would not be difficult to write a NetLogo procedure that would write a file in that format. This procedure would write new rows to the "Spreadsheet" CSV file on each NetLogo tick. Then Gephi could read in the file, and you'd be able to move back and forth in time, seeing how the graph changes. (You might need to use a bit of trial and error to figure out how to write Spreadsheet files based on the description on the Gephi site.)
Another option would be to display the evolving graph online using the graphstream protocol. Plugins for NetLogo as well as for gephi provide support for this.