[gdal-dev] GDAL Utilities: Keeping output format identical to input format
Zoltan Szecsei
zoltans at geograph.co.za
Mon Dec 15 03:27:34 PST 2014
On 2014/12/15 12:53, Even Rouault wrote:
> See http://www.gdal.org/gdal_edit.html
Hi Even,
I had a look at gdal_edit but could not get it to work on the server as
it is Ubuntu 12.04, so first I need to work out how to fix the gdal
installation on there.
*zls at gs0:/mnt/geo_lvm0/gisdata/raster/orthos$* gdalinfo --version
GDAL 1.9.1, released 2012/05/15
*zls at gs0:/mnt/geo_lvm0/gisdata/raster/orthos$* gdal_edit.py --version
Traceback (most recent call last):
File "/usr/bin/gdal_edit.py", line 33, in <module>
from osgeo import gdal
File "/usr/lib/python2.6/dist-packages/osgeo/__init__.py", line
21, in <module>
_gdal = swig_import_helper()
File "/usr/lib/python2.6/dist-packages/osgeo/__init__.py", line
17, in swig_import_helper
_mod = imp.load_module('_gdal', fp, pathname, description)
ImportError: /usr/local/lib/libgdal.so.1: undefined symbol: GEOSSimplify
*zls at gs0:/mnt/geo_lvm0/gisdata/raster/orthos$*
and:
*root at gs0:/usr/share/proj#*apt-get update | grep -i gis
*Ign* http://ppa.launchpad.net/ubuntugis/ubuntugis-unstable/ubuntu/
lucid/main Translation-en_ZA
*root at gs0:/usr/share/proj#*uname -a
Linux gs0 2.6.32-41-generic #94-Ubuntu SMP Fri Jul 6 18:00:34 UTC
2012 x86_64 GNU/Linux
*root at gs0:/usr/share/proj#*
To get back to this thread topic, then "Bad" TIFF is:
zls at gs0:/mnt/geo_lvm0/gisdata/raster/orthos$ gdalinfo
ngi_20141205/3121DA/3121DA_09_2011_407_RGB_RECT.TIF
Driver: GTiff/GeoTIFF
Files: ngi_20141205/3121DA/3121DA_09_2011_407_RGB_RECT.TIF
ngi_20141205/3121DA/3121DA_09_2011_407_RGB_RECT.TFW
Size is 10881, 11781
Coordinate System is:
LOCAL_CS["IMAGINE GeoTIFF Support ERDAS IMAGINE 2013 13.0.0.92
Geocoding information not available Projection Name = Unknown Units
= meters GeoTIFF Units = meters",
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (61356.750000000000000,-3491979.750000000000000)
Pixel Size = (0.500000000000000,-0.500000000000000)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
TIFFTAG_XRESOLUTION=1
TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
COMPRESSION=YCbCr JPEG
INTERLEAVE=PIXEL
SOURCE_COLOR_SPACE=YCbCr
Corner Coordinates:
Upper Left ( 61356.750,-3491979.750)
Lower Left ( 61356.750,-3497870.250)
Upper Right ( 66797.250,-3491979.750)
Lower Right ( 66797.250,-3497870.250)
Center ( 64077.000,-3494925.000)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
Overviews: 5440x5890, 2720x2945, 1360x1472, 680x736, 340x368,
170x184, 85x92, 42x46, 21x23
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
Overviews: 5440x5890, 2720x2945, 1360x1472, 680x736, 340x368,
170x184, 85x92, 42x46, 21x23
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue
Overviews: 5440x5890, 2720x2945, 1360x1472, 680x736, 340x368,
170x184, 85x92, 42x46, 21x23
zls at gs0:/mnt/geo_lvm0/gisdata/raster/orthos$
and the good version (tile next to this tile) is:
zls at gs0:/mnt/geo_lvm0/gisdata/raster/orthos$ gdalinfo
ngi_20141205/3121DA/3121DA_10_2011_407_RGB_RECT.TIF
Driver: GTiff/GeoTIFF
Files: ngi_20141205/3121DA/3121DA_10_2011_407_RGB_RECT.TIF
ngi_20141205/3121DA/3121DA_10_2011_407_RGB_RECT.TFW
Size is 10881, 11781
Coordinate System is:
PROJCS["WG21",
GEOGCS["GCS_Hartebeesthoek_1994",
DATUM["Hartebeesthoek94",
SPHEROID["WGS_1984",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6148"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",21],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (66102.750000000000000,-3492009.750000000000000)
Pixel Size = (0.500000000000000,-0.500000000000000)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
TIFFTAG_XRESOLUTION=1
TIFFTAG_YRESOLUTION=1
Image Structure Metadata:
COMPRESSION=YCbCr JPEG
INTERLEAVE=PIXEL
SOURCE_COLOR_SPACE=YCbCr
Corner Coordinates:
Upper Left ( 66102.750,-3492009.750) ( 21d41'46.15"E, 31d32'54.95"S)
Lower Left ( 66102.750,-3497900.250) ( 21d41'47.58"E, 31d36' 6.19"S)
Upper Right ( 71543.250,-3492009.750) ( 21d45'12.41"E, 31d32'53.78"S)
Lower Right ( 71543.250,-3497900.250) ( 21d45'13.94"E, 31d36' 5.02"S)
Center ( 68823.000,-3494955.000) ( 21d43'30.02"E, 31d34'30.00"S)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
Overviews: 5440x5890, 2720x2945, 1360x1472, 680x736, 340x368,
170x184, 85x92, 42x46, 21x23
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
Overviews: 5440x5890, 2720x2945, 1360x1472, 680x736, 340x368,
170x184, 85x92, 42x46, 21x23
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue
Overviews: 5440x5890, 2720x2945, 1360x1472, 680x736, 340x368,
170x184, 85x92, 42x46, 21x23
zls at gs0:/mnt/geo_lvm0/gisdata/raster/orthos$
I am using the following (instead of gdal_edit) but was wondering if my
band pixels were being altered.
gdal_translate -a_srs +proj="tmerc +lat_0=0 +lon_0=19 +k=1 +x_0=0
+y_0=0 +axis=wsu +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m
+no_defs" -co "COMPRESS=JPEG" -co "PHOTOMETRIC=YCBCR"
3121DA_09_2011_407_RGB_RECT.tif 3121da09.tif
Then because the gdal_translate strips out the overviews:
gdaladdo --config COMPRESS_OVERVIEW JPEG --config
PHOTOMETRIC_OVERVIEW YCBCR --config INTERLEAVE_OVERVIEW PIXEL
3121da09.tif 2 4 8 16 32 64 128 256 512
So, I was looking at a 1-step way of doing the above, ideally without
having to manually peek into an adjoining sheet to get the correct proj.
*This is what I am aiming at:*
I have about 40000 of these images, and I have written a BASH with awk
script to go through all the images and pull out certain info to build
a CSV file so that I can create attributed polygons to have info about
the underlying raster image. When this awk script fails to find embedded
proj info, it would be nice to make it peek into the adjoining ortho,
extract the relevant info, and edit it into the 'bad' ortho.
Thanks again for your help.
Kind regards,
Zoltan
--
===========================================
Zoltan Szecsei PrGISc [PGP0031]
Geograph (Pty) Ltd.
GIS and Photogrammetric Services
P.O. Box 7, Muizenberg 7950, South Africa.
Mobile: +27-83-6004028
Fax: +27-86-6115323 www.geograph.co.za
===========================================
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20141215/43773468/attachment.html>
More information about the gdal-dev
mailing list