generated from matsim-scenarios/matsim-scenario-template
-
Notifications
You must be signed in to change notification settings - Fork 2
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
Merge pull request #11 from matsim-scenarios/v1.1_emissions
V1.1 emissions
- Loading branch information
Showing
5 changed files
with
206 additions
and
122 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
115 changes: 0 additions & 115 deletions
115
src/test/java/org/matsim/run/EmissionAnalysisOutputTest.java
This file was deleted.
Oops, something went wrong.
186 changes: 186 additions & 0 deletions
186
src/test/java/org/matsim/run/EmissionAndNoiseAnalysisOutputTest.java
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,186 @@ | ||
| package org.matsim.run; | ||
|
|
||
| import org.junit.jupiter.api.Assertions; | ||
| import org.junit.jupiter.api.Disabled; | ||
| import org.junit.jupiter.api.Test; | ||
| import org.junit.jupiter.api.extension.RegisterExtension; | ||
| import org.junit.jupiter.api.io.TempDir; | ||
| import org.matsim.api.core.v01.Coord; | ||
| import org.matsim.api.core.v01.Id; | ||
| import org.matsim.api.core.v01.TransportMode; | ||
| import org.matsim.api.core.v01.population.*; | ||
| import org.matsim.application.MATSimApplication; | ||
| import org.matsim.core.config.Config; | ||
| import org.matsim.core.config.ConfigUtils; | ||
| import org.matsim.core.population.PersonUtils; | ||
| import org.matsim.core.population.PopulationUtils; | ||
| import org.matsim.core.utils.io.IOUtils; | ||
| import org.matsim.dashboards.LausitzSimWrapperRunner; | ||
| import org.matsim.simwrapper.SimWrapperConfigGroup; | ||
| import org.matsim.testcases.MatsimTestUtils; | ||
|
|
||
| import java.io.BufferedReader; | ||
| import java.io.IOException; | ||
| import java.nio.file.Path; | ||
| import java.util.HashMap; | ||
| import java.util.Map; | ||
| import java.util.Set; | ||
|
|
||
| import static org.matsim.application.ApplicationUtils.globFile; | ||
|
|
||
| class EmissionAndNoiseAnalysisOutputTest { | ||
|
|
||
| @RegisterExtension | ||
| public MatsimTestUtils utils = new MatsimTestUtils(); | ||
|
|
||
| @TempDir | ||
| public Path p; | ||
|
|
||
| private final static Id<Person> carPersonId = Id.createPersonId("Hoyerswerda-Cottbus_CAR"); | ||
|
|
||
| private static Config config = ConfigUtils.loadConfig(String.format("input/v%s/lausitz-v%s-10pct.config.xml", LausitzScenario.VERSION, LausitzScenario.VERSION)); | ||
|
|
||
| @Disabled("Test is used to secure functionality of emission analysis. As the analysis needs" + | ||
| "a lot of RAM, it is disabled and only run manually. -sme0924") | ||
| @Test | ||
| void runEmissionAnalysisOutputTest() throws IOException { | ||
| // Here, we do want dashboards in general (emission dashboard), but do not want to waste computation time on the standard dashboards | ||
| ConfigUtils.addOrGetModule(config, SimWrapperConfigGroup.class).exclude = Set.of("TripDashboard", "OverviewDashboard", "StuckAgentDashboard", | ||
| "TrafficCountsDashboard", "TrafficDashboard", "PublicTransitDashboard"); | ||
|
|
||
| Path inputPath = p.resolve("test-population.xml.gz"); | ||
| createTestPopulation(inputPath); | ||
|
|
||
| assert MATSimApplication.execute(LausitzScenario.class, config, | ||
| "--1pct", | ||
| "--iterations", "0", | ||
| "--output", utils.getOutputDirectory(), | ||
| "--config:plans.inputPlansFile", inputPath.toString(), | ||
| "--config:controller.overwriteFiles=deleteDirectoryIfExists") == 0 : "Must return non error code"; | ||
|
|
||
|
|
||
| Path csvPath = globFile(Path.of(utils.getOutputDirectory() + "/analysis/emissions"), "*emissions_per_link.csv*"); | ||
|
|
||
| Map<String, Double[]> nonZeroLinks = new HashMap<>(); | ||
|
|
||
| try { | ||
| BufferedReader reader = IOUtils.getBufferedReader(csvPath.toUri().toURL()); | ||
| String line; | ||
|
|
||
| // skip header | ||
| reader.readLine(); | ||
|
|
||
| while ((line = reader.readLine()) != null) { | ||
| String[] parts = line.split(","); | ||
|
|
||
| // if first value (CO) is zero, all others are too | ||
| if (Double.parseDouble(parts[1]) == 0.) { | ||
| continue; | ||
| } | ||
|
|
||
| Double[] values = new Double[23]; | ||
|
|
||
| for (int i = 1; i < parts.length; i++) { | ||
| values[i - 1] = Double.parseDouble(parts[i]); | ||
| } | ||
|
|
||
| nonZeroLinks.put(parts[0], values); | ||
| } | ||
| } finally { | ||
|
|
||
| } | ||
|
|
||
| Assertions.assertFalse(nonZeroLinks.isEmpty()); | ||
| Assertions.assertTrue(nonZeroLinks.containsKey("28922425#0")); | ||
| Assertions.assertTrue(nonZeroLinks.containsKey("-686055693#1")); | ||
| // CO | ||
| Assertions.assertEquals(830.9171, nonZeroLinks.get("28922425#0")[0], 0.0001); | ||
| Assertions.assertEquals(700.632, nonZeroLinks.get("-686055693#1")[0], 0.001); | ||
| // CO2_TOTAL | ||
| Assertions.assertEquals(8029.8788, nonZeroLinks.get("28922425#0")[1], 0.0001); | ||
| Assertions.assertEquals(6114.849, nonZeroLinks.get("-686055693#1")[1], 0.001); | ||
| } | ||
|
|
||
| @Disabled("Test is used to secure functionality of noise analysis. As the analysis needs" + | ||
| "a lot of RAM, it is disabled and only run manually. -sme0924") | ||
| @Test | ||
| void runNoiseAnalysisOutputTest() throws IOException { | ||
| // Here, we do want dashboards in general (emission dashboard), but do not want to waste computation time on the standard dashboards | ||
| ConfigUtils.addOrGetModule(config, SimWrapperConfigGroup.class).defaultDashboards = SimWrapperConfigGroup.Mode.disabled; | ||
|
|
||
| Path inputPath = p.resolve("test-population.xml.gz"); | ||
| createTestPopulation(inputPath); | ||
|
|
||
| assert MATSimApplication.execute(LausitzScenario.class, config, | ||
| "--1pct", | ||
| "--iterations", "0", | ||
| "--output", utils.getOutputDirectory(), | ||
| "--config:plans.inputPlansFile", inputPath.toString(), | ||
| "--config:controller.overwriteFiles=deleteDirectoryIfExists") == 0 : "Must return non error code"; | ||
|
|
||
| new LausitzSimWrapperRunner().execute(utils.getOutputDirectory(), "--noise", "--shp", "../../../../../../../input/shp/lausitz.shp"); | ||
|
|
||
| Path csvPath = globFile(Path.of(utils.getOutputDirectory() + "/analysis/noise-noise"), "*emission_per_day.csv*"); | ||
|
|
||
| Map<String, Double> nonZeroLinks = new HashMap<>(); | ||
|
|
||
| try { | ||
| BufferedReader reader = IOUtils.getBufferedReader(csvPath.toUri().toURL()); | ||
| String line; | ||
|
|
||
| // skip header | ||
| reader.readLine(); | ||
|
|
||
| while ((line = reader.readLine()) != null) { | ||
| String[] parts = line.split(","); | ||
|
|
||
| if (Double.parseDouble(parts[1]) == 0.) { | ||
| continue; | ||
| } | ||
|
|
||
| double value = Double.parseDouble(parts[1]); | ||
|
|
||
| nonZeroLinks.put(parts[0], value); | ||
| } | ||
| } finally { | ||
|
|
||
| } | ||
|
|
||
| Assertions.assertFalse(nonZeroLinks.isEmpty()); | ||
| Assertions.assertTrue(nonZeroLinks.containsKey("28922425#0")); | ||
| Assertions.assertTrue(nonZeroLinks.containsKey("-686055693#1")); | ||
| Assertions.assertEquals(71.21, Math.round(nonZeroLinks.get("28922425#0") * 100.) / 100. , 0.001); | ||
| Assertions.assertEquals(69.72, Math.round(nonZeroLinks.get("-686055693#1") * 100.) / 100., 0.001); | ||
| } | ||
|
|
||
| private void createTestPopulation(Path inputPath) { | ||
| Population population = PopulationUtils.createPopulation(config); | ||
| PopulationFactory fac = population.getFactory(); | ||
| Person person = fac.createPerson(carPersonId); | ||
| Plan plan = PopulationUtils.createPlan(person); | ||
|
|
||
| // home in hoyerswerda, nearest link 28922425#0 | ||
| Activity home = fac.createActivityFromCoord("home_2400", new Coord(863538.13,5711028.24)); | ||
| home.setEndTime(8 * 3600); | ||
| Activity home2 = fac.createActivityFromCoord("home_2400", new Coord(863538.13,5711028.24)); | ||
| home2.setEndTime(19 * 3600); | ||
| // work in hoyerswerda, nearest link(s): 686055693#1, -686055693#1 | ||
| Activity work = fac.createActivityFromCoord("work_2400", new Coord(863866.47,5710961.86)); | ||
| work.setEndTime(17 * 3600 + 25 * 60); | ||
|
|
||
| Leg leg = fac.createLeg(TransportMode.car); | ||
|
|
||
| plan.addActivity(home); | ||
| plan.addLeg(leg); | ||
| plan.addActivity(work); | ||
| plan.addLeg(leg); | ||
| plan.addActivity(home2); | ||
|
|
||
| person.addPlan(plan); | ||
| PersonUtils.setIncome(person, 1000.); | ||
| person.getAttributes().putAttribute("subpopulation", "person"); | ||
| population.addPerson(person); | ||
|
|
||
| new PopulationWriter(population).write(inputPath.toString()); | ||
| } | ||
| } |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters