On 19/11/2010 18:04, Frank Warmerdam wrote:
geographika wrote:
Hi,
I have upgraded from GDAL 1.6 (32bit Windows) to 1.7 (64 bit Windows)
and the following command no longer works:
C:\mapserver\bin\gdal\apps\gdalwarp C:\Data\Rasters\MiscSuit.tif
C:\RasterClips\mytest.tif -cutline c:\RasterClips\hello.json -te
118008.672141 177232.164284 138695.761666 206164.398565 -dstnodata -9999
I get the following message in the command prompt:
Creating output file that is 414P x 579L.
Processing input file C:\Data\Rasters\MiscSuit.tif.
for band 1, destination nodata value has been clamped to 0, the
original value being out of range.
ERROR 1: Failed to parse CUTLINE geometry wkt.
The CUTLINE is valid GeoJSON (but would clearly be invalid WKT) and
produces the correct results in 1.6.
The ogr formats lists GeoJSON as read/write.
I am using 64-bit builds of GDAL taken from http://vbkto.dyndns.org/sdk/
I have also tried using the development version of GDAL 1.8dev but
get the same message. If I do not use a CUTLINE the gdalwarp
completes successfully but the new image does not contain any data
from my original image (all cells are NoData).
Seth,
Can you make the .json file available? Does it work with OGR? Try
ogrinfo on it.
Hmm, tracking through the code the LoadCutline() function in gdalwarp.cpp
will read from the OGR datasource and convert the geometry to WKT
which is
attached to the cutline property of the warp. Later GDALWarpOperation::
Initialize() turns that into a polygon object and that is where the
message
is coming from. So I am *suspecting* improper WKT of some type is
getting
produced somehow. We will really need to see the WKT that is causing the
error to know more.
If you are comfortable rebuilding things then try changing:
const char *pszCutlineWKT =
CSLFetchNameValue( psOptions->papszWarpOptions, "CUTLINE" );
if( pszCutlineWKT )
{
in the file gdal/alg/gdalwarpoperation.cpp to something like:
const char *pszCutlineWKT =
CSLFetchNameValue( psOptions->papszWarpOptions, "CUTLINE" );
if( pszCutlineWKT )
{
printf( "WKT = %s\n", pszCutlineWKT );
Best regards,
Thanks very much for your reply Frank.
I have uploaded the GeoJSON to http://geographika.co.uk/downloads/test.json
Running ogrinfo C:\test.json on the file produces the following results:
ERROR 4: GeoJSON Driver doesn't support update.
Had to open data source read-only.
INFO: Open of `C:\RasterClips\hello.json'
using driver `GeoJSON' successful.
1: OGRGeoJSON (Polygon)
Is it possible to use the ogr2ogr to convert GeoJSON to WKT to check if
this is successful? WKT is not in the list of formats at
http://www.gdal.org/ogr/ogr_formats.html
Compiling a 64bit version of GDAL is on my list of things to learn, so
I will try your second suggestion then.
The gdalwarp with the same GeoJSON works in 1.6 but it does throw up the
warning that:
"the source raster dataset and the input vector layer do not have the
same SRS. Results will be probably incorrect."
Regards,
Seth
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
