/edl-scripts/bil2tif.py

http://alageospatialportal.googlecode.com/ · Python · 16 lines · 9 code · 3 blank · 4 comment · 3 complexity · 2fe298af5b46335aab46bd33a15e72a4 MD5 · raw file 

    

  1. #!/usr/bin/python
    2. 

  2. import os
    3. 

  3. import shutil
    4. 

  4. import edlconfig
    5. 

  5. 

  6. for root, dirs, files in os.walk(edlconfig.dataset):
    7. 

  7. 	for name in files:
    8. 

  8. 		if ".BIL" in name:
    9. 

  9. 			layername = name.replace(".BIL","")
    10. 

  10. 			print(layername)		
    11. 

  11. 			os.system(edlconfig.gdalapps+"/gdalwarp -s_srs '" + edlconfig.s_srs + "' -t_srs '" + edlconfig.t_srs + "' -srcnodata \"0\" -dstnodata \"-9999\" -of GTiff " + os.path.join(root,name) + " "  + os.path.join(root,layername) + ".tif")
    12. 

  12. #			os.system(edlconfig.gdalapps+"/gdal_translate -of GTiff " + os.path.join(root,name) + " "  + os.path.join(root,layername) + ".tif")
    13. 

  13. 

  14. 	#		os.system(edlconfig.gdalapps+"/gdalwarp -t_srs '" + edlconfig.t_srs + "' " +  os.path.join(root,layername) + "GDA.tif " + os.path.join(root,layername) + ".tif")
    15. 

  15. 

  16. 	#		os.system("rm " + os.path.join(root,layername) + "GDA.tif")
    17.