from PyRaster import *
import fileInput
from string import replace,joinfields
from os import path
#----------------------------------------------------
concFile="//Winfs/data/proj/Luftkvalitet/Airviro_projects/Chile_Santiago/Transsantiago/backgroundInterpolation/timeSeries/O3_timeseries_7stations.txt"
stationFile="//Winfs/data/proj/Luftkvalitet/Airviro_projects/Chile_Santiago/Transsantiago/backgroundInterpolation/stations/estacionesO3.txt"
substanceName="O3"
outputDirectory="D:/Projekt/Transsantiago/backgroundInterpolation/interpolatedAirviroResults/Airviro"
#interpolation parameters
power=1
smoothFactor=0
IDWweight=1
nNeighbours='all'
#-----------------------------------------------------
stations=fileInput.file2list(stationFile)
concData=TimeSeries()
concData.read(concFile)
concData.clean()
averageRaster=raster(329800,6278200,162,154,200,-9999)
tsOutFileName=path.join(outputDirectory,substanceName+'_timeSeries.txt')
ts_fid=open(tsOutFileName,'w')
header=""
stationX=[]
stationY=[]
stationNames=[]
for row in range(1,len(stations)):  
    stationX.append(float(stations[row][1]))
    stationY.append(float(stations[row][2]))
    stationNames.append(stations[row][0]+"_med")
    stationNames.append(stations[row][0]+"_int")
Esempio n. 2
0
import fileInput
from string import replace
from os import path
#----------------------------------------------------
concFile="Z:/Airviro_projects/Chile_Santiago/Transsantiago/backgroundInterpolation/PM10_Red_Macam.txt"
stationFile="Z:/Airviro_projects/Chile_Santiago/Transsantiago/backgroundInterpolation/estaciones.txt"
dummyStationFile="Z:/Airviro_projects/Chile_Santiago/Transsantiago/backgroundInterpolation/dummyStationsDistant.txt"
substanceName="PM10"
emisRasterFile="Z:/Airviro_projects/Chile_Santiago/Transsantiago/backgroundInterpolation/santiagopm10emissionsraster.txt"
outputDirectory="D:/Projekt/Santiago/testResults"
#interpolation parameters
power=1
smoothFactor=0
IDWweight=0.5
#-----------------------------------------------------
stations=fileInput.file2list(stationFile)
dstations=fileInput.file2list(dummyStationFile)
concData=TimeSeries()
concData.read(concFile)
emisRast=raster()

emisRast.readAscii(emisRasterFile)

stationX=[]
stationY=[]
for row in range(1,len(stations)):  
    stationX.append(float(stations[row][1]))
    stationY.append(float(stations[row][2]))

dstationX=[]
dstationY=[]