I would like to convert jp2 files (sentinel-2) to geotiff format using python 2.7 on win 10. I saw some solutions using python and GDAL but they did not work for me. The only thing that did work was to open OSGeo4W shell, go to the desired directory using
cd directory (from the shell) and then run (in the shell)
for %i in (*jp2) do gdal_translate -of GTiff -co TFW=YES %i %~ni_tiff.tif. The problem is that I need to perform it on multiple folders and also there are some processes I need to perform on each file before and after the conversion, the conversion from jp2 to geotiff is just a small part of the process. So I would like to know if it is possible to call OSGeo4W shell command from within a python script. My desired results would be a python script that at some part of the script runs a command on the OSGeo shell, each time for a different folder (I can store the folders in a list if necessary) and continues to process the products of the OSGeo4W command shell.
I read solutions regarding command line (cmd) using
subprocees but none of them referred to executing commands in a different platform than win command line.
The script according to my logic would be something like that:
import ... #src is a list of directories that holds directories with jp2 files src = ["c:\jp2_files_1", "c:\jp2_files_2","c:\jp2_files_3","c:\jp2_files_4"] for i in (src): os.chdir(i) SomethingI'mMIssing('for %i in (*jp2) do gdal_translate -of GTiff -co TFW=YES %i %~ni_tiff.tif') continue with python processing