LuneraysFlu_step3 - gama-platform/gama GitHub Wiki
This third step illustrates how to load GIS data and to agentify them.
- Define 2 new species that will just be displayed:
road
andbuilding
. - Define new global attributes to load GIS data (shapefile).
- Use the GIS data to create the
road
andbuilding
agents. - Add the
road
andbuilding
agents to the display.
For this step, you will need to add the shapefiles of the roads and buildings inside the includes
folder of the project. The shapefiles (and all the other files) can be found here.
In this model, we have to define two species of agents: the road
agents and the building
ones. These agents will not have a particular behavior, they will just be displayed.
For each of these species, we define an aspect called geom
. As we want to represent the geometry of the agent, we then use the keyword draw
that allows to draw a given geometry. In order to draw the geometry of the agent, we use the attribute shape
(which is a built-in attribute of all agents). The road will be displayed in black and the building in gray.
species road {
aspect geom {
draw shape color: #black;
}
}
species building {
aspect geom {
draw shape color: #gray;
}
}
GAMA allows to automatically read GIS data that are formatted as shapefiles (or as OSM file). In our model, we define 2 shapefiles: one corresponding to the roads and the other ones to the buildings. Note that GAMA is able to manage the projection of the GIS data. In order to set the right size (and position) of the world geometry, we define its value as the envelope of the road shapefile (and no more a square of 1500 meters).
global {
//... other attributes
file roads_shapefile <- file("../includes/roads.shp");
file buildings_shapefile <- file("../includes/buildings.shp");
geometry shape <- envelope(roads_shapefile);
//... init
}
In GAMA, the agentification of GIS data is very straightforward: it only requires to use the create
command with the from
facet to pass the shapefile. Each object of the shapefile will be directly used to instantiate an agent of the specified species. The reading of an attribute in a shapefile is also very simple. It only requires to use the with
facet: the argument of this facet is a dictionary of which the keys are the names of the agent attributes and the value the read
command followed by the name of the shapefile attribute.
In our model, we modify the init section in order to first create the road
agents from the road shapefile, and the building
agents from the building shapefile. Then, when creating people agents, we choose for them a random location inside a random building.
Note that it is possible to execute a sequence of statements at the creation of agents by using a block ({...}) rather than a simple line (;) when using the create
statement.
global {
// world variable definition
init{
create road from: roads_shapefile;
create building from: buildings_shapefile;
create people number:nb_people {
location <- any_location_in(one_of(building));
}
ask nb_infected_init among people {
is_infected <- true;
}
}
}
We used here the one_of operator that returns a random element from a list and the any_location_in operator that returns a random location inside a geometry.
In the map
display, we add the road
and building
species with their geom
aspect just before the people
species (in order to draw the people agents on the top of the roads and buildings).
experiment main_experiment type: gui {
... //parameter definition
output {
... //monitor definition
display map type: opengl{
species road aspect: geom;
species building aspect: geom;
species people aspect: circle;
}
... //chart display definition
}
}
https://github.com/gama-platform/gama.old/blob/GAMA_1.9.2/msi.gama.models/models/Tutorials/Luneray%20flu/models/model3.gaml