Re: [gdal-dev] gdalwrap and crossing 180W/180E

2010-10-04 Thread Manuel Massing
Hi, To follow myself up, I just had a look at the documentation (i know, next time I'll look before posting :-)) - the bounding window is sampled in destination coordinates, so my remark viz choosing SAMPLE_STEPS as the divisor of x-resolution is not correct - but increasing the number of SAMPLE_

Re: [gdal-dev] gdalwrap and crossing 180W/180E

2010-10-04 Thread Manuel Massing
Hi Ghislain, you could try increasing the "SAMPLE_STEPS" warp option, e.g. gdalwarp -wo "SAMPLE_STEPS=210" which should improve the localization of the source bounding window. Choosing SAMPLE_STEPS as a divisor of the source datasets x-resolution should help avoid this discretization a

Re: [gdal-dev] 4 point, projective warping

2010-08-02 Thread Manuel Massing
Hi, > I'd be very interested in a native GDAL implementation of perspective > warping. I've done something similar to the work Manuel has done on > Qgis, having written an in-house application that uses GSL/CBLAS. So I'd > be happy to help where I can, but as Manuel says, writing robust SVD > rout

Re: [gdal-dev] GDALSuggestedWarpOutput returns negative infinity in padfGeoTransformOut

2010-08-02 Thread Manuel Massing
Hi Lance, I suspect the most probable cause for the transform failure is a non- invertible geotransform matrix in CSTrans (i.e., due to incorrectly set scale values). What are m_srcGT / m_dstGT set to (the 6 actual values, not the code fragment used to set them)? m_srcGT should be set to the sou

Re: [gdal-dev] 4 point, projective warping

2010-07-30 Thread Manuel Massing
Hi Jeffrey, > I've thought about implementing a projective transform in the qgis > georeferencer, so if you're interested I'll let you know when I get around > to it. If this interest you, I have just commited an initial implementation of projective warping (see https://trac.osgeo.org/qgis/change

Re: [gdal-dev] 4 point, projective warping

2010-07-28 Thread Manuel Massing
Hi Jeffrey, > Hi, all - I need 4-point (GCP) projective warping ( > http://docs.bentley.com/en/I-RASB/irasbhelp171.html) for my very basic > rubbersheeting app at http://cartagen.org/maps (getting full-res TMS output > working). These are for aerial photographs, not scanned images. to implement t

Re: [gdal-dev] Re: Add GCP's to a file

2010-02-19 Thread Manuel Massing
Hi, > Gdal_translation has on option > [-gcp pixel line easting northing [elevation]]* this should also work using the VRT output format, which has the advantage that the raster data doesn't need to be written to an intermediate file (VRT "wraps" the original data). E.g.: gdal_translate [-gcp