Hi there,

Yes a nube with very little experience so I am seeking help.  Here's some background info:

I have a python script 'gdalcopyproj.py' (I did not write it), and it works great.  Basically it takes a tiff geo-referenced image file and transfers that information into ANOTHER tiff image that does not have geo-referenced data.

The script is written to apply this one file at a time "gdalcopyproj.py source_file dest_file".  One uses it inside the command prompt of a program called OSGeo4W to implement the script.  For example, it looks like so in the command window.

“C:\OSGeo4W>python gdalcopyproj.py SOCAL_CZ62.tif SOCAL_CZ62_LM.tif”

My goal is to process a whole folder of files at once.  I do have a batch file that has been created for another user.  Appears to be somewhat custom for his needs. Here is that batch file:

ECHO OFF
SETLOCAL EnableDelayedExpansion
for /f "delims=" %%i in ('2^>nul dir/a-d/b .tif^| findstr/i "^......\.tif$"') do (
for /f "delims=" %%j in ('2^>nul dir/a-d/b "%%~ni_*.tif"') do (
gdalcopyproj.py "%%i" "%%j"
)
)

What I do not understand is how to use/setup this batch file.  It appears to be seeking a certain type of file which has six characters and converts similarly named files which have custom filenames after the “_”.  Is this correct?  Can I customized that to fit my needs?

Second, where does one put this batch file.  I have tried in the same folder as the files. I have tried in the directory of OSGeo4W but no luck.  I have tried running the batch directly or from inside a command window.  Still no luck.  Either I get an error message or simply nothing happens.  This where I am stuck.

Can someone explain how to use this batch file somewhat step-by-step?

Clutch

P.S. – not sure if you need the python script to help explain but here she is:


try:
    from osgeo import gdal
except ImportError:
    import gdal
import sys
if len(sys.argv) < 3:
    print("USAGE: gdalcopyproj.py source_file dest_file")
    sys.exit(1)
input = sys.argv[1]
dataset = gdal.Open( input )
if dataset is None:
    print('Unable to open', input, 'for reading')
    sys.exit(1)
projection   = dataset.GetProjection()
geotransform = dataset.GetGeoTransform()
if projection is None and geotransform is None:
    print('No projection or geotransform found on file' + input)
    sys.exit(1)
output = sys.argv[2]
dataset2 = gdal.Open( output, gdal.GA_Update )
if dataset2 is None:
    print('Unable to open', output, 'for writing')
    sys.exit(1)
if geotransform is not None and geotransform != (0,1,0,0,0,1):
    dataset2.SetGeoTransform( geotransform )
if projection is not None and projection != '':
    dataset2.SetProjection( projection )
gcp_count = dataset.GetGCPs()
if gcp_count != 0:
    dataset2.SetGCPs( gcp_count, dataset.GetGCPProjection() )
dataset = None
dataset2 = None


P.S.S -  Here is even another batch file for the same purpose if this helps but still no luck... sorry coding is not my best suit.  Appears this one would take any file name up to the point of a "_"?  am I understanding this correct?

echo off
for /f "delims=" %%a in ('dir /b *.tif ^| findstr /v /c:"_.tif"') do gdalcopyproj.py "%%~a" "%%~na_.tif"Not sure what your first batch file is doing, but try modifying your second one as such:
Code: [Select]echo off
for /f "delims=" %%a in ('dir /a-d /b *.tif ^| find /v "_." ') do python gdalcopyproj.py "%%~a" "%%~na_.tif"

the segment 'for /f "delims=" %%A in () do' tells the program to loop through the list in the parenthesis.  In this case that is 'dir /a-d /b *.tif ^| find /v "_."' which is a list of everything in the current directory, excluding files whose names contain the string "_." and folders.  For each of these, it will run the command 'python gdalcopyproj.py "%%~a" "%%~na_.tif"' where %%~a is replaced with each thing in the list, and %%~na is replaced with only that file's name.

If you run into the problem where nothing shows up, try removing the "echo off" at the beginning and post the screen.  (you can do this by right clicking and using the 'mark' function)Hey Lemonilla,

Thank you so much for taking a stab at this.  Tried a few times without success so let me give you a play-by-play:

Created a folder on my S: drive called Geotest.  In there I placed two files: ALBQ_T20.tif (CONTAINS geo info), and a second file ALBQ_T20_SH.tiff (which contains no geo info).  I placed your bat file there and ran it.  I got a quick "blip" of the command window but then nothing.  Same thing happens if I remove the ECHO off text.  No difference.

2nd test.  I moved the tiff files and the bat file into the C:\OSGeo4W directory.  Same thing happens there as well.  Are my filenames different than what you expected them to be?


Adding a 'pause' command at the end of a batch file is a good way to debug it.  If you wouldn't mind, remove "echo off" and put the pause at the end and paste the error message it gives.  This can also be accomplished by running the batch file from a cmd.exe window (command prompt window).

When I run the program (with a slight edit) it looks like it should run.  Try running the below outside of the batch file and see if there is a syntax error in calling the python script.
Code: [Select]python gdalcopyproj.py "ALBQ_T20.tif" "ALBQ_T20_.tif"


Code: (my code) [Select]echo off
for /f "delims=" %%a in ('dir /a-d /b *.tif ^| find /v "_." ') do (
REM by changing ".tif" to "%%~xa" we are maintaining the exact extension of the previous file.
REM e.g. if we call it on a .tiff file instead of .tif, it will maintain the double f at the end
REM the 'echo' is there for debugging purposes only, it shows me how the python script will be called, but bypasses the actual execution
echo python gdalcopyproj.py "%%~a" "%%~na_%%~xa"
)
pause
Code: (output) [Select]python gdalcopyproj.py "ALBQ_T20.tif" "ALBQ_T20_.tif"
python gdalcopyproj.py "ALBQ_T20_SH.tiff" "ALBQ_T20_SH_.tiff"
Press any key to continue . . .



EDIT: You might try removing the quotes around the file names, that might be messing with your python script.  Just know then that your file names cannot have spaces or they will be interpreted as new arguments.Oh look!  Somebody else with exactly the same question!  What are the odds??

http://www.dostips.com/forum/viewtopic.php?f=3&t=6409 Quote from: foxidrive on April 18, 2015, 02:15:27 AM

Oh look!  Somebody else with exactly the same question!  What are the odds??

http://www.dostips.com/forum/viewtopic.php?f=3&t=6409

Hey, I am back.  Got tied up with a few other pressing issues.  Yes, I did post this issue at a few different forums.  Looking for options.

• to sys.argv[#-1]