Balazs, If you want to programatically perform the actions of gdalwarp, it's best to go through the code of grapwarp.cpp http://trac.osgeo.org/gdal/browser/trunk/gdal/apps/gdalwarp.cpp
There are some more test scripts you can refer, if you want to do this through command line. http://trac.osgeo.org/gdal/browser/trunk/autotest/warp/warp.py http://trac.osgeo.org/gdal/browser/trunk/autotest/utilities/test_gdalwarp.py On Thu, Aug 4, 2011 at 5:26 PM, Balazs Szabó <[email protected]> wrote: > Dear list! > > How can I use this two bash command with gdal-python bindings? > > gdal_translate -of GTiff -a_srs EPSG:xxxx -gcp 775 560 20.58329 48.12650 > -gcp 9518 539 20.8329 48.12650 -gcp 9525 7136 20.83329 48.00150 -gcp 758 > 7134 20.58329 48.00150 o.tif s.tif > gdalwarp -tps s.tif t.tif > > [I can done this task by os.popen('gdalwarp -tps s.tif t.tif') ] > I found some code here: > http://svn.osgeo.org/gdal/trunk/autotest/gcore/transformer.py > But after this line ( tr = gdal.Transformer( ds, None, [ 'METHOD=GCP_TPS' ] > )) I dont know how can I warping my hole raster. > > Here is my python code: > > import osgeo.gdal > import osgeo.osr > import os > > osgeo.gdal.AllRegister() > driver=osgeo.gdal.GetDriverByName('GTiff') > > src='path/to/orig/o.tif' > dst='path/to/t.tif' > gcps=give_me_gcps() # I got a list like [gcp1,gcp2...] > > t_srs=osgeo.osr.SpatialReference() > t_srs.ImportFromEPSG(xxxx) > > s_img=osgeo.gdal.Open(src) > s_img_cols=s_img.RasterXSize > s_img_rows=s_img.RasterYSize > > > t_img=driver.Create(dst,s_img_cols,s_img_rows,3,osgeo.gdal.GDT_Byte) # > Create a new image > t_img.SetGCPs(gcps,t_srs.ExportToWkt()) > > for band_num in range(s_img.RasterCount): > # copying the original tif bands > band=s_img.GetRasterBand(band_num+1) > inData=band.ReadAsArray(0,0,cols,rows) > outBand=t_img.GetRasterBand(band_num+1) > outBand.WriteArray(inData,0,0) > outBand.FlushCache() > > > # ???? > > s_img=None > t_img=None > > > Thanks in advance > Best regards, > > Balazs > > _______________________________________________ > gdal-dev mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev > -- Best regards, Chaitanya kumar CH. +91-9494447584 17.2416N 80.1426E
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
