Pular para o conteúdo

Recuperar informações de APA do Aimsun Next

O plugin de kernel do Aimsun Next oferece classes para iterar por todos os caminhos armazenados em um arquivo Path Assignment (.apa). A ideia é criar uma classe que derive da classe ForestCalculator e implementa o doPath function. É possível filtrar os caminhos antes de iterar por eles. Este filtro pode ser feito usando a classe AnalysisFilter e configurando os filtros correspondentes. Por exemplo, você pode filtrar caminhos que passam por um conjunto de seções, ou caminhos que entram por um conjunto de seções de entrada, ou apenas caminhos relacionados a um intervalo de tempo especificado.

O doPath função recebe todos os caminhos para uma determinada pathKey: - mInterval -> índice do intervalo de tempo - mVehicleIndex -> índice do tipo de veículo - mOriginIndex -> índice do centroide de origem - mOriginCon -> índice da conexão de origem - mDestinationIndex -> índice do centroide de destino - mDestinationCon -> índice da conexão de destino

E as informações do caminho dentro do parâmetro NSPTreeData. Esta classe contém a porcentagem de demanda atribuída ao caminho, bem como as informações do caminho e outras estatísticas de caminhos, como os veículos atribuídos e os componentes da função.

Uma vez que a classe derivada de ForestCalculator se estiver implementado, este é o código para iterar por todos os caminhos:

#Filter
customFilter=AnalysisFilter()
#PathReader
#baseForest is the apa file
pr=PathReader(baseForest, model, repIdName)
pr.init()
pr.calculate(customFilter)
pr.finishPaths()

Exemplo: Recuperar informações de atribuição de trajetórias

import os.path

class PathReader(ForestCalculator):
    def __init__(self, forest, model, name):
        self.forest = forest
        self.model = model
        self.name = name
        self.geo = forest.getGeometry()
        self.nbLines = 1;
        super(ForestCalculator, self).__init__(forest, model, ForestCalculator.eOnlyAvailablePaths)

    def init(self):
        info= QFileInfo(self.model.getDocumentFileName())
        filename = os.path.join(info.absolutePath(),self.name+'_paths.csv')
        self.file = open(filename, 'w')
        if self.file!=None:
            self.file.write("interval;vehicle;origin;destination;type;path;probability;veh generated;veh arrived;cost;distance;travel time;speed;function component\n")
            self.file.flush()
        else:
            print("Could not open file %s"%filename)

    def doPath(self,pathKey, treeData, connections,aux):
        if self.file!=None:
            #print('entro')
            originId = self.geo.getOriginCentroid(pathKey.mOriginIndex)
            destId = self.geo.getDestinationCentroid(pathKey.mDestinationIndex)
            self.file.write("%d;%d;%d;%d;"%(pathKey.mInterval, pathKey.mVehicleIndex, originId, destId))
            typepath = treeData.pathType()
            typelist = ['None', 'RC', 'PAR', 'ODR', 'TRJ','P2P']
            self.file.write("%s;"%typelist[typepath])
            p = treeData.getPath(self.forest)
            sections = self.forest.getPathSections(self.geo.getSection(treeData.entrancePosition), p, pathKey.mVehicleIndex )
            s='-'.join(map(str, sections))
            self.file.write("%s"%s)
            comp = treeData.getComponents()
            self.file.write(";%f"%comp[1])
            self.file.write(";%f"%comp[0])
            self.file.write(";%f"%comp[2])

            costs =  self.geo.getPathCost( pathKey.mOriginIndex, pathKey.mDestinationIndex, connections, pathKey.mVehicleIndex, pathKey.mInterval )
            cost = costs.mCost
            self.file.write(";%f"%cost)

            distance = costs.mDistance
            self.file.write(";%f"%distance)

            if comp[2]==0.0:
                travel_time = 0.0
            else:
                travel_time = comp[3]/comp[2]

            self.file.write(";%f"%travel_time)

            if travel_time==0.0:
                speed = 0.0*3.6
            else:
                speed = distance*3.6/travel_time

            self.file.write(";%f"%speed)

            components = costs.mCostComponents
            if len(components)>0:
                for component in components:
                    self.file.write(";%f"%component)
            #c = ';'.join(map(str, comp))

            #self.file.write("%s"%c)
            self.file.write("\n")
            self.nbLines =self.nbLines+1
            if self.nbLines%1000==0:
                self.file.flush()

    def finishPaths(self):
        if self.file!=None:
            self.file.close()


replicacion1626=model.getCatalog().find( 1626 )
if (replicacion1626.getSimulationStatus() == GKGenericExperiment.eLoaded ):

    #Forest
    baseForest=replicacion1626.getPathsForest()

    #Filter
    customFilter=AnalysisFilter()
    #Filter could no applied iif it is passing to the reader without conditions. so this part could be skipped.
    intervals=[0]
    customFilter.setIntervals(intervals)
    userClass=0
    customFilter.setUserClass(userClass)
    originCentroids=[1586]
    customFilter.setOriginCentroids(originCentroids)
    destinationCentroids=[1589]
    customFilter.setDestinationCentroids(destinationCentroids)

    #Name
    repIdName='Replication'+str(replicacion1626.getId())
    #PathReader
    pr=PathReader(baseForest, model, repIdName)
    pr.init()
    pr.calculate(customFilter)
    pr.finishPaths()
else:   
    print( 'Replication %s is not simulated yet.' %repIdName)

print( 'info in file: %s_paths.csv.' %repIdName)

O Aimsun Next recupera informações estáticas APA de transporte público

Analogamente ao kernel do Aimsun Next, o plugin Macro PT oferece classes para iterar sobre todos os caminhos armazenados em um arquivo de Path Assignment gerado a partir de uma Transit Assignment. A ideia é criar uma classe que derive da classe PTForestCalculator e implementa o doPTODPair function. É possível filtrar os caminhos antes de iterar por eles. Este filtro pode ser feito usando a classe PTForestFilter e definindo os filtros correspondentes. Por exemplo, você pode filtrar caminhos que passam por um conjunto de seções ou linhas, ou apenas caminhos relacionados a uma classe de usuário especificada.

O doPTODPair a função recebe uma instância PTODPair contendo todas as informações relacionadas à demanda, aos caminhos e às suas porcentagens.

Uma vez que a classe derivada de ForestCalculator se estiver implementado, este é o código para iterar por todos os caminhos:

#Filter
filter=PTForestFilter()
#PathReader
#forest is the apa file
pr=PathReader(forest, filter, model, repIdName)
pr.visit(forest)

Exemplo: Recuperar informações de custos de atribuição de caminhos

import os.path

class PTPathReader(PTForestCalculator):

    def __init__(self, forest, filter, model, name ):
        self.file = None
        self.model = model
        self.name = name
        self.forest = forest
        super(PTForestCalculator, self).__init__(forest, filter, model, False)

    def init(self):
        info= QFileInfo(self.model.getDocumentFileName())
        filename = os.path.join(info.absolutePath(),self.name+'_paths.csv')
        self.file = open(filename, 'w')
        if self.file!=None:
            self.file.write("user;origin;destination;volume;cost;time;distance;\n")
            self.file.flush()
        else:
            print("Could not open file %s"%filename)

    def doPTODPair(self,ptODpair):

        if self.file!=None:

            userId = self.forest.getUsers()[ ptODpair.user() ]
            originId = self.forest.getOriginCentroidId( ptODpair.origin() )
            destinationId = self.forest.getDestinationCentroidId( ptODpair.destination() )

            for strategy in ptODpair.strategies():

                for path in strategy.paths():

                    self.file.write("%d;"% userId )
                    self.file.write("%d;"% originId )
                    self.file.write("%d;"% destinationId )
                    self.file.write("%f;"%( path.factor()*ptODpair.volume() ) )
                    self.file.write("%f;"%path.cost())
                    self.file.write("%f;"%path.inVehicleTime())
                    self.file.write("%f;"%path.inVehicleDistance())
                    self.file.write("\n")

            self.file.flush()

    def postProcess(self):
        if self.file!=None:
            self.file.close()

EXEMPLO DE USO

experiment=model.getCatalog().findByName("Transit Assignment Experiment") # SET THE REPLICATION ID HERE 
forest=experiment.getForest()

#Filter
filter = PTForestFilter()
filter.setUserClass( model.getCatalog().findByName("Transit Users - Work").getId() )
destinationCentroids=[model.getCatalog().findByName("Mall West").getId()] # SET THE DESTINATION CENTROIDS HERE
filter.setDestinationCentroids(destinationCentroids)

#PTPathReader
pr=PTPathReader(forest, filter, model, experiment.getName())
pr.visit(forest)

print("Done")