Python GlobalRasterStats示例

示例#1

def main():
    record = SimulationRecord(100, False)
    record.loadHDF5('../data/anazasi.h5', True, True)

    agentResults = GlobalAgentStats(';')
    agentResults.addAnalysis(AgentNum())
    agentResults.addAnalysis(AgentMean('estimates'))
    agentResults.applyTo(record, 'households.csv', 'Household')

    rasterResults = GlobalRasterStats(';')
    rasterResults.addAnalysis(RasterMean())
    rasterResults.addAnalysis(RasterSum())
    rasterResults.applyTo(record, 'water.csv', 'yield x100')

示例#2

文件： analyzeData.py 项目： digideskio/pandora

def main():
    record = SimulationRecord(100, False)
    record.loadHDF5('../data/anazasi.h5', True, True)

    agentResults = GlobalAgentStats(';')
    agentResults.addAnalysis(AgentNum())
    agentResults.addAnalysis(AgentMean('estimates'))
    agentResults.applyTo(record,'households.csv', 'Household')
    
    rasterResults = GlobalRasterStats(';')
    rasterResults.addAnalysis(RasterMean())
    rasterResults.addAnalysis(RasterSum())
    rasterResults.applyTo(record, 'water.csv', 'yield x100')

示例#3

def main():
    """ this file provides a template that can be used to analyze any Pandora model """
    record = SimulationRecord(1, False)
    record.loadHDF5('data/results.h5', True, True)

    agentResults = GlobalAgentStats(';')
    agentResults.addAnalysis(AgentNum())
    agentResults.addAnalysis(AgentMean('x'))
    agentResults.addAnalysis(AgentMean('y'))
    agentResults.addAnalysis(AgentMean('value'))
    agentResults.addAnalysis(AgentSum('value'))
    agentResults.applyTo(record, 'agents.csv', 'HelloAgent')

    rasterResults = GlobalRasterStats(';')
    rasterResults.addAnalysis(RasterMean())
    rasterResults.addAnalysis(RasterSum())
    rasterResults.applyTo(record, 'values.csv', 'values')

示例#4

文件： pyPandoraAnalysisTemplate.py 项目： digideskio/pandora

def main():
    """ this file provides a template that can be used to analyze any Pandora model """
    record = SimulationRecord(1, False)
    record.loadHDF5('data/results.h5', True, True)

    agentResults = GlobalAgentStats(';')
    agentResults.addAnalysis(AgentNum())
    agentResults.addAnalysis(AgentMean('x'))
    agentResults.addAnalysis(AgentMean('y'))
    agentResults.addAnalysis(AgentMean('value'))
    agentResults.addAnalysis(AgentSum('value'))
    agentResults.applyTo(record, 'agents.csv', 'HelloAgent')

    rasterResults = GlobalRasterStats(';')
    rasterResults.addAnalysis(RasterMean())
    rasterResults.addAnalysis(RasterSum())
    rasterResults.applyTo(record, 'values.csv', 'values')

示例#5

def main():
    record = SimulationRecord(1, False)
    record.loadHDF5('../data/results.h5', True, True)

    agentResults = GlobalAgentStats(';')

    agentResults.addAnalysis(AgentNum())
    agentResults.addAnalysis(AgentMean('strength x100'))
    agentResults.addAnalysis(AgentMean('aggressiveness x100'))
    agentResults.addAnalysis(AgentStdDev('strength x100'))
    agentResults.addAnalysis(AgentStdDev('aggressiveness x100'))
    agentResults.addAnalysis(AgentSum('strength x100'))
    agentResults.addAnalysis(AgentSum('aggressiveness x100'))

    agentResults.applyTo(record, 'farmers.csv', 'Farmer')
    agentResults.applyTo(record, 'herders.csv', 'Herder')

    rasterResults = GlobalRasterStats(';')
    rasterResults.addAnalysis(RasterMean())
    rasterResults.addAnalysis(RasterSum())

    rasterResults.applyTo(record, 'conflicts.csv', 'conflicts')

示例#6

文件： testAnalysis.py 项目： digideskio/pandora

def main():
    record = SimulationRecord(1, False)
    record.loadHDF5('../data/results.h5', True, True)

    agentResults = GlobalAgentStats(';')

    agentResults.addAnalysis(AgentNum())
    agentResults.addAnalysis(AgentMean('strength x100'))
    agentResults.addAnalysis(AgentMean('aggressiveness x100'))
    agentResults.addAnalysis(AgentStdDev('strength x100'))
    agentResults.addAnalysis(AgentStdDev('aggressiveness x100'))
    agentResults.addAnalysis(AgentSum('strength x100'))
    agentResults.addAnalysis(AgentSum('aggressiveness x100'))

    agentResults.applyTo(record,'farmers.csv', 'Farmer')
    agentResults.applyTo(record,'herders.csv', 'Herder')

    rasterResults = GlobalRasterStats(';')
    rasterResults.addAnalysis(RasterMean())
    rasterResults.addAnalysis(RasterSum())

    rasterResults.applyTo(record, 'conflicts.csv', 'conflicts')

示例#7

文件： testAnalysis.py 项目： alcaina/pandora

#!/usr/bin/python

import sys, random
sys.path.append('../../../pyPandora/')
sys.path.append('../../../')
from pyPandora import SimulationRecord, GlobalAgentStats, AgentNum, AgentMean, AgentSum, AgentStdDev, GlobalRasterStats, RasterMean, RasterSum

record = SimulationRecord(1, False)
record.loadHDF5('../data/results.h5', True, True)

agentResults = GlobalAgentStats(';')

agentResults.addAnalysis(AgentNum())
agentResults.addAnalysis(AgentMean('strength x100'))
agentResults.addAnalysis(AgentMean('aggressiveness x100'))
agentResults.addAnalysis(AgentStdDev('strength x100'))
agentResults.addAnalysis(AgentStdDev('aggressiveness x100'))
agentResults.addAnalysis(AgentSum('strength x100'))
agentResults.addAnalysis(AgentSum('aggressiveness x100'))

agentResults.applyTo(record,'farmers.csv', 'Farmer')
agentResults.applyTo(record,'herders.csv', 'Herder')

rasterResults = GlobalRasterStats(';')
rasterResults.addAnalysis(RasterMean())
rasterResults.addAnalysis(RasterSum())

rasterResults.applyTo(record, 'conflicts.csv', 'conflicts')

示例#8

        for i in range(0, 10):
            newAgent = MyAgent('id_' + str(i))
            self.addAgent(newAgent)
            newAgent.setRandomPosition()
            newAgent._value = random.randint(0, 10)


myConfig = Config(SizeInt(32, 32), 10)
myWorld = MyWorld(myConfig)
myWorld.initialize()
myWorld.run()

# analysis

record = SimulationRecord(1, False)
record.loadHDF5('data/results.h5', True, True)

agentResults = GlobalAgentStats(';')
agentResults.addAnalysis(AgentNum())
agentResults.addAnalysis(AgentMean('x'))
agentResults.addAnalysis(AgentMean('y'))
agentResults.addAnalysis(AgentMean('value'))
agentResults.addAnalysis(AgentSum('value'))
agentResults.applyTo(record, 'agents.csv', 'id')
#agentResults.addAnalysis(AgentHDFtoSHP('shp/agents.shp', -1))

rasterResults = GlobalRasterStats(';')
rasterResults.addAnalysis(RasterMean())
rasterResults.addAnalysis(RasterSum())
rasterResults.applyTo(record, 'resources.csv', 'test')

示例#9

文件： testAnalysis.py 项目： alcaina/pandora

		for i in range (0, 10):
			newAgent = MyAgent('id_'+str(i))
			self.addAgent(newAgent)
			newAgent.setRandomPosition()
			newAgent._value = random.randint(0,10)

mySimulation = Simulation(SizeInt(32,32), 10)
myWorld = MyWorld(mySimulation)
myWorld.initialize()
myWorld.run()

# analysis

record = SimulationRecord(1, False)
record.loadHDF5('data/results.h5', True, True)

agentResults = GlobalAgentStats(';')
agentResults.addAnalysis(AgentNum())
agentResults.addAnalysis(AgentMean('x'))
agentResults.addAnalysis(AgentMean('y'))
agentResults.addAnalysis(AgentMean('value'))
agentResults.addAnalysis(AgentSum('value'))
agentResults.applyTo(record, 'agents.csv', 'id')
#agentResults.addAnalysis(AgentHDFtoSHP('shp/agents.shp', -1))

rasterResults = GlobalRasterStats(';')
rasterResults.addAnalysis(RasterMean())
rasterResults.addAnalysis(RasterSum())
rasterResults.applyTo(record, 'resources.csv', 'test')