<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
</head>
<body bgcolor="#FFFFFF" text="#000000">
Hi,<br>
I've been struggling to get some 3200 orthos mosaicked and tiled.<br>
They are very old and have really bad projection definitions, pixel
sizes and metadata funnies.<br>
I'm down to the last 261, and have hit a brick wall. <br>
Given that gdalwarp works on this image when I leave out s_srs and
t_srs (see below), I'm not sure if this is a bug related hiccough or
a proj issue.<br>
<br>
I have compiled the latest "GDAL 2.1.1, released 2016/07/07" myself,
on Ubuntu 16.04.1 Desktop: Linux gl05 4.4.0-36-generic #55-Ubuntu
SMP Thu Aug 11 18:01:55 UTC 2016 x86_64 x86_64 x86_64 GNU/Linux<br>
<br>
<br>
<br>
I get the following error from gdalwarp:<br>
<br>
<blockquote>ERROR 6: Unable to load PROJ.4 library (libproj.so),
creation of<br>
OGRCoordinateTransformation failed.<br>
</blockquote>
I have tried various different gdalwarp parameters - even leaving
gdalwarp to sense in input projection and me only specifying the
output projection.<br>
If I leave out both s_srs and t_srs then gdalwarp does work but
obliviously does not reproject to Lo19. <br>
I have even tried feeding the output of this redundant gdalwarp into
a gdalwarp with t_srs, but no success.<br>
<br>
Specifying only t_srs (as above) on many other orthos on Lo17 &
Lo21 does work so it is likely that the way this projection is
embedded in this ortho is the problem.<br>
I have tried stripping it out with gdal_edit.py but haven't gotten
that right yet.<br>
<br>
Anyone with pointers or ideas on this?<br>
<br>
My last effort before email this list is below, as is the gdalinfo
output.<br>
<br>
Thanks in advance,<br>
Zoltan<br>
<br>
<tt>gdalwarp -s_srs '+proj=tmerc +lat_0=0 +lon_0=21 +k=1 +x_0=0
+y_0=0 +datum=Hartebeesthoek94 +units=m +no_defs' \</tt><tt><br>
</tt><tt> -t_srs '+proj=tmerc +lat_0=0 +lon_0=19 +k=1
+x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs' \</tt><tt><br>
</tt><tt> -r cubicspline -multi -overwrite -of GTiff -co
TFW=YES -co COMPRESS=JPEG -co JPEG_QUALITY=75 -tr 0.75 -0.75 \</tt><tt><br>
</tt><tt> /mnt/sdb1/used_images/dwa_2000/3220/3220CC01.tif \</tt><tt><br>
</tt><tt>
/mnt/sda1/gisdata/raster/used_2016/dwa_2000/3220/3220CC01.tif</tt><tt><br>
</tt><tt>#</tt><tt><br>
</tt><tt>#</tt><tt><br>
</tt><tt>gdalinfo /mnt/sdb1/used_images/dwa_2000/3220/3220CC01.tif</tt><tt><br>
</tt><tt>#</tt><tt><br>
</tt><tt>Driver: GTiff/GeoTIFF</tt><tt><br>
</tt><tt>Files: /mnt/sdb1/used_images/dwa_2000/3220/3220CC01.tif</tt><tt><br>
</tt><tt>Size is 6729, 7720</tt><tt><br>
</tt><tt>Coordinate System is:</tt><tt><br>
</tt><tt>LOCAL_CS[" Geocoding information not available Projection
Name = WGS 1984 Lo21 Units = meters GeoTIFF Units = meters",</tt><tt><br>
</tt><tt> UNIT["metre",1,</tt><tt><br>
</tt><tt> AUTHORITY["EPSG","9001"]]]</tt><tt><br>
</tt><tt>Origin = (-93872.189308887900552,-3625321.216973003931344)</tt><tt><br>
</tt><tt>Pixel Size = (0.750000000000000,-0.750000000000000)</tt><tt><br>
</tt><tt>Metadata:</tt><tt><br>
</tt><tt> AREA_OR_POINT=Area</tt><tt><br>
</tt><tt> TIFFTAG_DOCUMENTNAME=hetotiff</tt><tt><br>
</tt><tt> TIFFTAG_IMAGEDESCRIPTION=ER Mapper 6.4 hetotiff</tt><tt><br>
</tt><tt> TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)</tt><tt><br>
</tt><tt> TIFFTAG_SOFTWARE=IMAGINE TIFF Support</tt><tt><br>
</tt><tt>Copyright 1991 - 1999 by ERDAS, Inc. All Rights Reserved</tt><tt><br>
</tt><tt>@(#)$RCSfile: etif.c $ $Revision: 1.9.1.3 $ $Date:
2002/07/29 15:39:06EDT $</tt><tt><br>
</tt><tt> TIFFTAG_XRESOLUTION=1</tt><tt><br>
</tt><tt> TIFFTAG_YRESOLUTION=1</tt><tt><br>
</tt><tt>Image Structure Metadata:</tt><tt><br>
</tt><tt> INTERLEAVE=PIXEL</tt><tt><br>
</tt><tt>Corner Coordinates:</tt><tt><br>
</tt><tt>Upper Left ( -93872.189,-3625321.217) </tt><tt><br>
</tt><tt>Lower Left ( -93872.189,-3631111.217) </tt><tt><br>
</tt><tt>Upper Right ( -88825.439,-3625321.217) </tt><tt><br>
</tt><tt>Lower Right ( -88825.439,-3631111.217) </tt><tt><br>
</tt><tt>Center ( -91348.814,-3628216.217) </tt><tt><br>
</tt><tt>Band 1 Block=6729x1 Type=Byte, ColorInterp=Red</tt><tt><br>
</tt><tt>Band 2 Block=6729x1 Type=Byte, ColorInterp=Green</tt><tt><br>
</tt><tt>Band 3 Block=6729x1 Type=Byte, ColorInterp=Blue</tt>
</body>
</html>