Kirk McKelvey wrote:
When I try to warp the entire world (WGS84) into a UTM zone using
gdalwarp, I get a black image.
This is not surprising given the valid domain for such a transformation.
The first problem is that the bounding box gets reprojected into outer
space (again, b/c it is outside the valid domain). But once I “found”
the data, I was able to get something quasi-sensible out of GDAL (image
attached). South America shows twice, which, although pretty cool, is
another reflection of the domain problem.
So my question is: How can I tell programmatically whether (1) any given
warp operation makes sense and (2) what its valid input domain is?
Kirk,
I do not have a good answer for this. I have not been able to do it
in a general way.
I was surprised to find that the success vector of
GDALReprojectionTransform() does not return failures for out-of-domain
inputs (at least it did not in this case). Is it intended to?
This depends in part on what PROJ.4 does. Traditionally the transverse
mercator project did not consider any values out of range and so all
transformation succeeded even if the results were crazy. I recently
modified the PROJ.4 tm projection to consider anything 90 degrees off
the central meridian as a failure which should get propagated up through
the transformer as FALSE in the success vector.
Some projections with distinct horizons do already fail for over the horizon
transformations.
Best regards,
--
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up | Frank Warmerdam, warmer...@pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush | Geospatial Programmer for Rent
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev
