/**
   * Writes the population and their attributes to file.
   *
   * @param outputfolder
   */
  public void writePopulation(String outputfolder) {
    if (this.sc.getPopulation().getPersons().size() == 0 || this.personAttributes == null) {
      throw new RuntimeException("Either no persons or person attributes to write.");
    } else {
      LOG.info("Writing population to file...");
      PopulationWriter pw = new PopulationWriter(this.sc.getPopulation(), this.sc.getNetwork());
      pw.writeV5(outputfolder + "Population.xml");

      LOG.info("Writing person attributes to file...");
      ObjectAttributesXmlWriter oaw = new ObjectAttributesXmlWriter(this.personAttributes);
      oaw.setPrettyPrint(true);
      oaw.writeFile(outputfolder + "PersonAttributes.xml");
    }
  }
  private static void removeActs(String networkFile, String popFile) {
    String popOutFile = popFile + "_removedTransitActs.xml.gz";
    TransitActsRemover remover = new TransitActsRemover();

    Scenario sc = ScenarioUtils.createScenario(ConfigUtils.createConfig());
    new MatsimNetworkReader(sc.getNetwork()).readFile(networkFile);
    new PopulationReader(sc).readFile(popFile);

    for (Person person : sc.getPopulation().getPersons().values()) {
      for (Plan plan : person.getPlans()) {
        remover.run(plan);
      }
    }

    PopulationWriter popWriter = new PopulationWriter(sc.getPopulation(), sc.getNetwork());
    popWriter.write(popOutFile);
  }
Esempio n. 3
0
  public static void main(String[] args) throws Exception {

    String baseDir = "../../../shared-svn/projects/cottbus/data/scenarios/cottbus_scenario/";
    String networkFile = baseDir + "network_wgs84_utm33n.xml.gz";
    boolean useLanduse = true;
    String populationOutputDirectory = baseDir + "cb_spn_gemeinde_nachfrage_landuse_ohneTagebau/";

    //		String shapesOutputDirectory = populationOutputDirectory + "shapes_all_modes/";
    String shapesOutputDirectory = populationOutputDirectory + "shapes/";

    //		String populationOutputFile = populationOutputDirectory +
    // "commuter_population_wgs84_utm33n.xml.gz";
    String populationOutputFile =
        populationOutputDirectory + "commuter_population_wgs84_utm33n_car_only.xml.gz";
    OutputDirectoryLogging.initLoggingWithOutputDirectory(populationOutputDirectory);

    File shapes = new File(shapesOutputDirectory);
    if (!shapes.exists()) {
      shapes.mkdir();
    }

    Config config1 = ConfigUtils.createConfig();
    config1.network().setInputFile(networkFile);
    Scenario sc = ScenarioUtils.loadScenario(config1);

    CommuterDataReader cdr = new CommuterDataReader();
    cdr.addFilterRange(12071000);
    cdr.addFilter("12052000"); // 12052000 == cottbus stadt
    cdr.readFile(
        "../../../shared-svn/studies/countries/de/pendler_nach_gemeinden/brandenburg_einpendler.csv");
    //		cdr.getCommuterRelations().add(new CommuterDataElement("12052000", "12052000", 1000));

    String gemeindenBrandenburgShapeFile =
        "../../../shared-svn/studies/countries/de/brandenburg_gemeinde_kreisgrenzen/gemeinden/dlm_gemeinden.shp";
    ShapeFileReader gemeindenReader = new ShapeFileReader();
    Collection<SimpleFeature> gemeindenFeatures =
        gemeindenReader.readFileAndInitialize(gemeindenBrandenburgShapeFile);

    CommuterDemandWriter cdw =
        new CommuterDemandWriter(
            gemeindenFeatures,
            gemeindenReader.getCoordinateSystem(),
            cdr.getCommuterRelations(),
            MGC.getCRS(TransformationFactory.WGS84_UTM33N));
    // landuse
    if (useLanduse) {
      DgLanduseReader landuseReader = new DgLanduseReader();
      Tuple<Collection<SimpleFeature>, CoordinateReferenceSystem> homeLanduse =
          landuseReader.readLanduseDataHome();
      Tuple<Collection<SimpleFeature>, CoordinateReferenceSystem> workLanduse =
          landuseReader.readLanduseDataWork();
      cdw.addLanduse("home", homeLanduse);
      cdw.addLanduse("work", workLanduse);
    }
    //
    //		cdw.setScalefactor(1.0); // all modes
    cdw.setScalefactor(0.55); // car mode share
    //		cdw.setScalefactor(0.1); //testing

    cdw.computeDemand(sc);
    PopulationWriter populationWriter = new PopulationWriter(sc.getPopulation(), sc.getNetwork());
    populationWriter.write(populationOutputFile);
    log.info("population written to " + populationOutputFile);

    // write some test output
    Config config = ConfigUtils.createConfig();
    config.network().setInputFile(networkFile);
    config.plans().setInputFile(populationOutputFile);
    Scenario baseScenario = ScenarioUtils.loadScenario(config);

    String shapeFilename = shapesOutputDirectory + "commuter_population_home.shp";
    new DgPopulation2ShapeWriter(
            baseScenario.getPopulation(), MGC.getCRS(TransformationFactory.WGS84_UTM33N))
        .write("home", shapeFilename, MGC.getCRS(TransformationFactory.WGS84_UTM33N));
    shapeFilename = shapesOutputDirectory + "commuter_population_work.shp";
    new DgPopulation2ShapeWriter(
            baseScenario.getPopulation(), MGC.getCRS(TransformationFactory.WGS84_UTM33N))
        .write("work", shapeFilename, MGC.getCRS(TransformationFactory.WGS84_UTM33N));

    log.info("done!");
    OutputDirectoryLogging.closeOutputDirLogging();
  }
Esempio n. 4
0
  private void run() {

    // Braunschweig, Stadt - BS
    // Wolfsburg, Stadt - WB
    // Gifhorn - GH
    // Wolfenbüttel - WL
    // Helmstedt - HS
    // Peine - PE
    // Salzgitter, Stadt - SG
    // Börde - BR
    // Altmarkkreis Salzwedel - AS
    // Region Hannover - RH
    // Goslar - GL
    // Harz - HZ
    // Hildesheim - HH
    // Magdeburg, Landeshauptstadt - MB
    // Celle - CL
    // Uelzen - UL
    // Stendal - SD
    // Göttingen - GG
    // Salzlandkreis - SL
    // Jerochower Land - JL
    // Osnabrück, Stadt - OB
    // Osterode am Harz - OR

    // Braunschweig, Stadt - Braunschweig, Stadt | work
    createWorkers("BS", "BS", 62158 * config.getScalefactor(), 0.45, 0.2, 0.2, "03101", "03101");

    // Braunschweig, Stadt - Braunschweig, Stadt | school
    createPupils("BS", "BS", 24907 * config.getScalefactor(), 0.1, 0.13, 0.22, "03101", "03101");

    // Braunschweig, Stadt - Braunschweig, Stadt | university
    createStudents("BS", "BS", 29455 * config.getScalefactor(), 0.1, 0.13, 0.22, "03101", "03101");

    // Wolfsburg, Stadt - Wolfsburg, Stadt | work
    createWorkers(
        "WB", "WB", (41470 - 18313) * config.getScalefactor(), 0.4, 0.2, 0.2, "03103", "03103");
    createVWWorkers("WB", "03103", 11207, 1260 + 5328, 658, 0.4, 0.2, 0.2);
    //		3103	11207	1260	5328	518	18313

    // Wolfsburg, Stadt - Wolfsburg, Stadt | school
    createPupils("WB", "WB", 13867 * config.getScalefactor(), 0.1, 0.13, 0.27, "03103", "03103");

    // Wolfsburg, Stadt - Wolfsburg, Stadt | university
    createStudents("WB", "WB", 8823 * config.getScalefactor(), 0.01, 0.4, 0.25, "03103", "03103");

    // Gifhorn - Wolfsburg, Stadt | work
    createWorkers(
        "GH", "WB", (26484 - 15152) * config.getScalefactor(), 0.8, 0.0, 0.0, "03151", "03103");
    createVWWorkers("GH", "03151", 9156, 1050 + 4758, 518, 0.8, 0.0, 0.0);
    // 3151	9156	1050	4758	188	15152

    // Gifhorn - Gifhorn | work
    createWorkers("GH", "GH", 26414 * config.getScalefactor(), 0.63, 0.15, 0.15, "03151", "03151");

    // Gifhorn - Gifhorn | school
    createPupils("GH", "GH", 26429 * config.getScalefactor(), 0.2, 0.3, 0.2, "03151", "03151");

    // Wolfenbüttel - Braunschweig, Stadt | work
    createWorkers("WL", "BS", 13304 * config.getScalefactor(), 0.8, 0.0, 0.0, "03158", "03101");

    // Helmstedt - Wolfsburg, Stadt | work
    createWorkers(
        "HS", "WB", (12731 - 7908) * config.getScalefactor(), 0.8, 0.0, 0.0, "03154", "03103");
    createVWWorkers("HS", "03154", 4704, 536 + 2668, 188, 0.8, 0.0, 0.0);
    //		3154	4704	536	2668	0	7908

    // Braunschweig, Stadt - Wolfsburg, Stadt | work
    createWorkers(
        "BS", "WB", (10273 - 5006) * config.getScalefactor(), 0.65, 0.01, 0.0, "03101", "03103");
    createVWWorkers("BS", "03101", 4085, 127 + 693, 101, 0.65, 0.01, 0.0);
    // 3101	4085	127	693	101	8107

    // Peine - Braunschweig, Stadt | work
    createWorkers("PE", "BS", 9089 * config.getScalefactor(), 0.8, 0.0, 0.0, "03157", "03101");

    // Gifhorn - Braunschweig, Stadt | work
    createWorkers("GH", "BS", 7586 * config.getScalefactor(), 0.8, 0.0, 0.0, "03151", "03101");

    // Salzgitter, Stadt - Braunschweig, Stadt | work
    createWorkers("SG", "BS", 5572 * config.getScalefactor(), 0.8, 0.0, 0.0, "03102", "03101");

    // Helmstedt - Braunschweig, Stadt | work
    createWorkers("HS", "BS", 4618 * config.getScalefactor(), 0.8, 0.0, 0.0, "03154", "03101");

    // Börde - Wolfsburg, Stadt | work
    createWorkers(
        "BR", "WB", (3685 - 856) * config.getScalefactor(), 0.8, 0.0, 0.0, "15085", "03103");
    createVWWorkers("BR", "15085", 396, 460, 0, 0.8, 0.0, 0.0);
    //		15362	396	0	460	0	856

    // Altmarkkreis Salzwedel - Wolfsburg, Stadt | work
    createWorkers(
        "AS", "WB", (3305 - 740) * config.getScalefactor(), 0.8, 0.0, 0.0, "15081", "03103");
    createVWWorkers("AS", "15081", 290, 450, 0, .8, .0, .0);
    //		15370	290	0	450	0	740

    // Region Hannover - Wolfsburg, Stadt | work
    createWorkers(
        "RH", "WB", (2850 - 1148) * config.getScalefactor(), 0.8, 0.0, 0.0, "03241", "03103");
    createVWWorkers("RH", "03241", 1004, 144, 0, 0.8, 0.0, .0);
    //		3241	1004	0	144	0	1148

    // Region Hannover - Braunschweig, Stadt | work
    createWorkers("RH", "BS", 2833 * config.getScalefactor(), 0.8, 0.0, 0.0, "03241", "03101");

    // Wolfenbüttel - Wolfsburg, Stadt | work
    createWorkers(
        "WL", "WB", (2676 - 1191) * config.getScalefactor(), 0.8, 0.0, 0.0, "03158", "03103");
    createVWWorkers("WL", "03158", 943, 248, 0, .8, .0, .0);
    //		3158	943	0	248	0	1191

    // Braunschweig, Stadt - Gifhorn | work
    createWorkers("BS", "GH", 2635 * config.getScalefactor(), 0.8, 0.0, 0.0, "03101", "03151");

    // Braunschweig, Stadt - Wolfburg, Stadt | work
    createWorkers("BS", "WB", 2310 * config.getScalefactor(), 0.8, 0.0, 0.0, "03101", "03103");

    // Braunschweig, Stadt - Goslar | work
    createWorkers("BS", "GL", 1960 * config.getScalefactor(), 0.8, 0.0, 0.0, "03101", "03153");

    // Wolfsburg, Stadt - Gifhorn | work
    createWorkers("WB", "GH", 1844 * config.getScalefactor(), 0.8, 0.0, 0.0, "03103", "03151");

    // Peine - Wolfburg, Stadt | work
    createWorkers(
        "PE", "WB", (1756 - 782) * config.getScalefactor(), 0.8, 0.0, 0.0, "03157", "03103");
    createVWWorkers("PE", "03157", 551, 231, 0, .8, .0, .0);
    //		3157	551	0	231	0	782

    // Harz - Braunschweig, Stadt | work
    createWorkers("HZ", "BS", 1685 * config.getScalefactor(), 0.8, 0.0, 0.0, "15085", "03101");

    // Altmarkkreis Salzwedel - Gifhorn | work
    createWorkers("AS", "GH", 1649 * config.getScalefactor(), 0.8, 0.0, 0.0, "15081", "03151");

    // Börde - Braunschweig, Stadt | work
    createWorkers("BR", "BS", 1260 * config.getScalefactor(), 0.8, 0.0, 0.0, "15085", "03101");

    // Hildesheim - Braunschweig, Stadt | work
    createWorkers("HH", "BS", 1195 * config.getScalefactor(), 0.8, 0.0, 0.0, "03254", "03101");

    // Salzgitter, Stadt - Wolfburg, Stadt | work
    createWorkers("SG", "WB", 0 * config.getScalefactor(), 0.8, 0.0, 0.0, "03102", "03103");
    createVWWorkers("SG", "03102", 180, 106, 0, .8, .0, .0);
    // 3102	180	0	106	658	944

    // Magdeburg, Landeshauptstadt - Wolfburg, Stadt | work
    createWorkers("MB", "WB", 847 * config.getScalefactor(), 0.8, 0.0, 0.0, "15003", "03103");
    createVWWorkers("MB", "15003", 256, 140, 0, .8, .0, .0);

    // 15003	256	0	140	0	396

    // Peine - Gifhorn | work
    createWorkers("PE", "GH", 834 * config.getScalefactor(), 0.8, 0.0, 0.0, "03157", "03151");

    // Celle - Gifhorn | work
    createWorkers("CL", "GH", 815 * config.getScalefactor(), 0.8, 0.0, 0.0, "03351", "03151");

    // Helmstedt - Gifhorn | work
    createWorkers("HS", "GH", 735 * config.getScalefactor(), 0.8, 0.0, 0.0, "03154", "03151");

    // Region Hannover - Gifhorn | work
    createWorkers("RH", "GH", 680 * config.getScalefactor(), 0.8, 0.0, 0.0, "03241", "03151");

    // Celle - Wolfburg, Stadt | work
    createWorkers(
        "CL", "WB", (607 - 315) * config.getScalefactor(), 0.8, 0.0, 0.0, "03351", "03103");
    createVWWorkers("CL", "03351", 155, 160, 0, .8, .0, .0);
    //		3351	155	0	160	0	315

    // Magdeburg, Landeshauptstadt - Braunschweig, Stadt | work
    createWorkers("MB", "BS", 511 * config.getScalefactor(), 0.8, 0.0, 0.0, "15003", "03101");

    // Wolfenbüttel - Gifhorn | work
    createWorkers("WB", "GH", 496 * config.getScalefactor(), 0.8, 0.0, 0.0, "03158", "03151");

    // Harz - Wolfburg, Stadt | work
    createWorkers("HZ", "WB", 476 * config.getScalefactor(), 0.8, 0.0, 0.0, "15085", "03103");

    // Uelzen - Gifhorn | work
    createWorkers("UL", "GH", 474 * config.getScalefactor(), 0.8, 0.0, 0.0, "03360", "03151");

    // Stendal - Wolfburg, Stadt | work
    createWorkers("SD", "WB", 456 * config.getScalefactor(), 0.8, 0.0, 0.0, "15090", "03103");

    // Celle - Braunschweig, Stadt | work
    createWorkers("CL", "BS", 440 * config.getScalefactor(), 0.8, 0.0, 0.0, "03351", "03101");

    // Börde - Gifhorn | work
    createWorkers("BR", "GH", 400 * config.getScalefactor(), 0.8, 0.0, 0.0, "15083", "03151");

    // Hildesheim - Wolfburg, Stadt | work
    createWorkers(
        "HH", "WB", (391 - 107) * config.getScalefactor(), 0.8, 0.0, 0.0, "03254", "03103");
    createVWWorkers("HH", "03103", 107, 0, 0, .8, .0, .0);
    //		3254	107	0	0	0	107

    // Goslar - Wolfburg, Stadt | work
    createWorkers("GL", "WB", 384 * config.getScalefactor(), 0.8, 0.0, 0.0, "03153", "03103");

    // Göttingen - Braunschweig, Stadt | work
    createWorkers("GG", "BS", 369 * config.getScalefactor(), 0.8, 0.0, 0.0, "03152", "03101");

    // Uelzen - Wolfburg, Stadt | work
    createWorkers("UL", "WB", 338 * config.getScalefactor(), 0.8, 0.0, 0.0, "03360", "03103");

    // Altmarkkreis Salzwedel - Braunschweig, Stadt | work
    createWorkers("AS", "BS", 245 * config.getScalefactor(), 0.8, 0.0, 0.0, "15081", "03101");

    // Salzlandkreis - Wolfburg, Stadt | work
    createWorkers("SL", "WB", 214 * config.getScalefactor(), 0.8, 0.0, 0.0, "15089", "03103");

    // Salzlandkreis - Braunschweig, Stadt | work
    createWorkers("SL", "BS", 185 * config.getScalefactor(), 0.8, 0.0, 0.0, "15089", "03101");

    // Salzgitter, Stadt - Gifhorn | work
    createWorkers("SG", "GH", 174 * config.getScalefactor(), 0.8, 0.0, 0.0, "03102", "03151");

    // Jerichower Land - Wolfburg, Stadt | work
    createWorkers("JL", "WB", 170 * config.getScalefactor(), 0.8, 0.0, 0.0, "15086", "03103");

    // Harz - Gifhorn | work
    createWorkers("HZ", "GH", 148 * config.getScalefactor(), 0.8, 0.0, 0.0, "15085", "03151");

    // Osnabrück, Stadt - Wolfburg, Stadt | work
    createWorkers("OB", "WB", 136 * config.getScalefactor(), 0.8, 0.0, 0.0, "03404", "03103");

    // Osterode am Harz - Braunschweig, Stadt | work
    createWorkers("OR", "BS", 129 * config.getScalefactor(), 0.8, 0.0, 0.0, "03156", "03101");

    // Jerichower Land - Braunschweig, Stadt | work
    createWorkers("JL", "BS", 117 * config.getScalefactor(), 0.8, 0.0, 0.0, "15086", "03101");

    // Magdeburg, Landeshauptstadt - Gifhorn | work
    createWorkers("MB", "GH", 114 * config.getScalefactor(), 0.8, 0.0, 0.0, "15003", "03151");

    // Uelzen - Braunschweig, Stadt | work
    createWorkers("UL", "BS", 110 * config.getScalefactor(), 0.8, 0.0, 0.0, "03360", "03101");

    // Stendal - Braunschweig, Stadt | work
    createWorkers("SD", "BS", 108 * config.getScalefactor(), 0.8, 0.0, 0.0, "15090", "03101");

    // Goslar - Gifhorn | work
    createWorkers("GL", "GH", 101 * config.getScalefactor(), 0.8, 0.0, 0.0, "03153", "03151");

    createA2TransitTruckers(Math.round(15000 * config.getScalefactor()));
    createVWTruckers(Math.round(900 * config.getScalefactor()));
    replaceDoubtfulLegsByOtherMode();

    System.out.println("generated Agents: " + commuterCounter);
    System.out.println("VW Workers: " + vwWorkerCounter);
    System.out.println("Workers: " + workerCounter);
    createAgentGroupNearTransitstrops(scenario, 1500, config.getTransitSchedule());
    replaceSptByPtp();
    PopulationWriter pw = new PopulationWriter(scenario.getPopulation(), scenario.getNetwork());
    pw.write(config.getPlansOutputString());
    new ObjectAttributesXmlWriter(scenario.getPopulation().getPersonAttributes())
        .writeFile(config.getObjectAttributes());
  }
Esempio n. 5
0
  private void run() {
    this.shapeMap = readShapeFile(KREISE, "Nr");
    this.buildingsMap = readShapeFile(BUILDINGS, "osm_id");
    this.shops = readFacilityLocations(SHOPS);
    this.kindergartens = readFacilityLocations(KINDERGARTEN);

    // CB-CB (Stadtverkehr: Modal Split mit 50% IV/ÖV)
    double commuters = 22709 * SCALEFACTOR;
    createPersons("12052000", "12052000", commuters, 0.5);

    // CB-LDS (Regionalverkehr: Modal Split mit 80% IV)
    commuters = 399 * SCALEFACTOR;
    createPersons("12052000", "12061000", commuters, 0.8);

    // CB-LOS
    commuters = 139 * SCALEFACTOR;
    createPersons("12052000", "12067000", commuters, 0.8);

    // CB-SPN
    commuters = 4338 * SCALEFACTOR;
    createPersons("12052000", "12071000", commuters, 0.8);

    // LDS-CB
    commuters = 1322 * SCALEFACTOR;
    createPersons("12061000", "12052000", commuters, 0.8);

    // LDS-SPN
    commuters = 522 * SCALEFACTOR;
    createPersons("12061000", "12071000", commuters, 0.8);

    // LOS-CB
    commuters = 382 * SCALEFACTOR;
    createPersons("12067000", "12052000", commuters, 0.8);

    // LOS-SPN
    commuters = 449 * SCALEFACTOR;
    createPersons("12067000", "12071000", commuters, 0.8);

    // SPN-CB
    commuters = 11869 * SCALEFACTOR;
    createPersons("12071000", "12052000", commuters, 0.8);

    // SPN-LDS
    commuters = 408 * SCALEFACTOR;
    createPersons("12071000", "12061000", commuters, 0.8);

    // SPN-LOS
    commuters = 466 * SCALEFACTOR;
    createPersons("12071000", "12067000", commuters, 0.8);

    // SPN-SPN
    commuters = 22524 * SCALEFACTOR;
    createPersons("12071000", "12071000", commuters, 0.8);

    /*
     * Achtung! 3 Nullen extra hier ggü. ISIS Daten
     */

    System.out.println(this.shapeMap.keySet());

    createCarSharingVehicleLocationsFile();

    PopulationWriter pw = new PopulationWriter(scenario.getPopulation(), scenario.getNetwork());
    pw.write(PLANSFILEOUTPUT);
  }