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")
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=[]