Re: [gdal-dev] gdal.AutoCreateWarpedVRT returning different results than osr.CoordinateTransformation

2012-05-06 Thread Etienne Tourigny
Billy, You might consider creating a warped image with a larger pixel size (a multiple of your real target size, say 3000m) to calculate the correct extents and then throwing it away. Etienne On Sat, May 5, 2012 at 6:22 PM, Billy Newman newman...@gmail.com wrote: Still very confused. I have

Re: [gdal-dev] gdal.AutoCreateWarpedVRT returning different results than osr.CoordinateTransformation

2012-05-05 Thread Billy Newman
Frank, Thanks so mush for the quick response. One more quick question. I am grabbing the metadata from my image files to store in the database as an index of where my imagery lives. I server out my images in EPSG:4326 so that is how I would like to store the metadata (bounding box) of the

Re: [gdal-dev] gdal.AutoCreateWarpedVRT returning different results than osr.CoordinateTransformation

2012-05-05 Thread Frank Warmerdam
On 12-05-05 01:19 PM, Billy Newman wrote: Can you think of a way that I can get the same metadata as gdalwarp would give for a warped image without actually warping the image? Billy, Just open the .vrt and read the corners (that is compute them from the size and geotransform). There will not

Re: [gdal-dev] gdal.AutoCreateWarpedVRT returning different results than osr.CoordinateTransformation

2012-05-05 Thread Billy Newman
Still very confused. I have a tiff file. I want to get the metadata for that tiff but in EPSG:4326, without actually warping the image. I can use CoordinateTransformation to do this, but again it gives me different results than gdalwarp. I will use gdalwarp later (realtime) to serve out an

Re: [gdal-dev] gdal.AutoCreateWarpedVRT returning different results than osr.CoordinateTransformation

2012-05-05 Thread Billy Newman
Ok I think I get exactly what you are saying. Dataset dataset = // image file dataset in wrong projection Dataset epsg4326Dataset = gdal.AutoCreateWarpedVrt(); double[] geoTransform = epsg4326Dataset.GetGeoTransform(); // grab out the corner point in epsg:4326 So doing that won't actually

[gdal-dev] gdal.AutoCreateWarpedVRT returning different results than osr.CoordinateTransformation

2012-05-04 Thread Billy Newman
I ran a quick test and confirmed that gdal.AutoCreateWarpedVRT is returning different results than osr.CoordinateTransformation. I.E. I opened a file using: Dataset dataset = gdal.Open(/images/image.tif); I grabbed the corner point and used osr.CoordinateTransformation.TransformPoint(), to

Re: [gdal-dev] gdal.AutoCreateWarpedVRT returning different results than osr.CoordinateTransformation

2012-05-04 Thread Frank Warmerdam
On 12-05-04 03:39 PM, Billy Newman wrote: I ran a quick test and confirmed that gdal.AutoCreateWarpedVRT is returning different results than osr.CoordinateTransformation. I.E. I opened a file using: Dataset dataset = gdal.Open(/images/image.tif); I grabbed the corner point and used