Le portail francophone de la géomatique

# Import des packages
import sys, os, subprocess, time
try:
    #pour les utilisateurs ArcGIS 9.2/9.3
    import arcgisscripting
    gp = arcgisscripting.create()
    print "\n" + "Utilisation de ArcGIS 9.2/9.3 avec arcgisscripting..." + "\n"
    gp.AddMessage("\n" + "Utilisation de ArcGIS avec arcgisscripting..." + "\n")

except:
    #pour les utilisateurs ArcGIS 9.0/9.1
    import win32com.client
    gp = win32com.client.Dispatch("esriGeoprocessing.GpDispatch.1")
    print "\n" + "Utilisation de ArcGIS 9.0/9.1 avec win32com..." + "\n"
    gp.AddMessage("\n" + "Utilisation de ArcGIS 9.0/9.1avec win32com..." + "\n")

# Lecture arguments paths
inputfilename = sys.argv[1]
wsize = sys.argv[2]
tfactor = sys.argv[3]
outputworkspace = sys.argv[4]
algo = sys.argv[5]

try:
    (inpath,infile)=os.path.split(inputfilename)
    infile=infile[:len(infile)-4]+'_ALGO_'+str(wsize)+'_'+str(tfactor)+'.tif'
    outputfilename=os.path.join(outputworkspace, infile)
    print outputfilename
    
    paramsfile = open('C:\\IDL\\parameters.txt', "w")
    paramsfile.write(inputfilename)
    paramsfile.write('\n')
    paramsfile.write(wsize)
    paramsfile.write('\n')
    paramsfile.write(tfactor)
    paramsfile.write('\n')
    paramsfile.write(outputfilename)
    paramsfile.close()

        
except:
    gp.AddMessage(gp.GetMessages(2))
    print gp.GetMessages(2)
    
elif algo =="GDAL":
    print "GDAL"
        subprocess.call('C:/gdalalgo.exe', shell=True)
    
elif algo =="IDL":
    print "IDL"
    subprocess.call('C:/gdalalgo.sav', shell=True)
else :
    print "ERROR"
    
time.sleep(2)