[Gdal-dev] gdalwarp does not use my +towgs84 parameter
Uwe Raudszus
URaudszus at STADTDO.DE
Tue Oct 2 05:20:18 EDT 2007
Hallo List,
I'm trying to reproject a tif image using gdalwarp. Specially,
from GAUSS-KRUEGER Zone 2 to UTM Zone 32 (epsg:25832). For more precision i
try to use local towgs84
parameter. But it seems so, that gdalwarp does not use my +towgs84
parameter.
Is there any way to include these parameters in my gdalwarp code?
Uwe
C:\work\rasterfactory\work>gdalwarp -s_srs "+proj=tmerc +lat_0=0 +lon_0=6
+k=1.0
00000 +x_0=2500000 +y_0=0 +ellps=bessel +datum=potsdam
+towgs84=600.75,20,413 +
units=m" -t_srs EPSG:25832 -multi -dstnodata 255 -of Gtiff -co
"INTERLEAVE=PIXEL
" -co COMPRESS=PACKBITS -co TFW=true
c:\work\rasterfactory\output\8\8_31466_wgs.
tif c:\work\rasterfactory\output\8\8_25832.tif
Output dataset c:\work\rasterfactory\output\8\8_25832.tif exists,
but some commandline options were provided indicating a new dataset
should be created. Please delete existing dataset and run again.
C:\work\rasterfactory\work>gdalwarp -s_srs "+proj=tmerc +lat_0=0 +lon_0=6
+k=1.0
00000 +x_0=2500000 +y_0=0 +ellps=bessel +datum=potsdam
+towgs84=600.75,20,413 +
units=m" -t_srs EPSG:25832 -multi -dstnodata 255 -of Gtiff -co
"INTERLEAVE=PIXEL
" -co COMPRESS=PACKBITS -co TFW=true
c:\work\rasterfactory\output\8\8_31466_wgs.
tif c:\work\rasterfactory\output\8\8_25832.tif
Copying color table from c:\work\rasterfactory\output\8\8_31466_wgs.tif to
new f
ile.
Creating output file that is 8861P x 9557L.
Processing input file c:\work\rasterfactory\output\8\8_31466_wgs.tif.
:0...10...20...30...40...50...60...70...80...90...100 - done.
C:\work\rasterfactory\work>listgeo
c:\work\rasterfactory\output\8\8_25832.tif
TIFFReadDirectory: Warning, Unknown field with tag 42113 (0xa481)
encountered.
Geotiff_Information:
Version: 1
Key_Revision: 1.0
Tagged_Information:
ModelTiepointTag (2,3):
0 0 0
387002.011 5710514.84 0
ModelPixelScaleTag (1,3):
0.0564260761 0.0564260761 0
End_Of_Tags.
Keyed_Information:
GTModelTypeGeoKey (Short,1): ModelTypeProjected
GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
GTCitationGeoKey (Ascii,22): "ETRS89 / UTM zone 32N"
GeogCitationGeoKey (Ascii,7): "ETRS89"
GeogAngularUnitsGeoKey (Short,1): Angular_Degree
ProjectedCSTypeGeoKey (Short,1): Unknown-25832
ProjLinearUnitsGeoKey (Short,1): Linear_Meter
End_Of_Keys.
End_Of_Geotiff.
PCS = 25832 (ETRS89 / UTM zone 32N)
Projection = 16032 (UTM zone 32N)
Projection Method: CT_TransverseMercator
ProjNatOriginLatGeoKey: 0.000000 ( 0d 0' 0.00"N)
ProjNatOriginLongGeoKey: 9.000000 ( 9d 0' 0.00"E)
ProjScaleAtNatOriginGeoKey: 0.999600
ProjFalseEastingGeoKey: 500000.000000 m
ProjFalseNorthingGeoKey: 0.000000 m
GCS: 4258/ETRS89
Datum: 6258/European Terrestrial Reference System 1989
Ellipsoid: 7019/GRS 1980 (6378137.00,6356752.31)
Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)
Projection Linear Units: 9001/metre (1.000000m)
Corner Coordinates:
Upper Left ( 387002.011, 5710514.843)
Lower Left ( 387002.011, 5709975.579)
Upper Right ( 387502.002, 5710514.843)
Lower Right ( 387502.002, 5709975.579)
Center ( 387252.007, 5710245.211)
C:\work\rasterfactory\work>gdalwarp -s_srs "+proj=tmerc +lat_0=0 +lon_0=6
+k=1.0
00000 +x_0=2500000 +y_0=0 +ellps=bessel +datum=potsdam +units=m" -t_srs
EPSG:25
832 -multi -dstnodata 255 -of Gtiff -co "INTERLEAVE=PIXEL" -co
COMPRESS=PACKBITS
-co TFW=true c:\work\rasterfactory\output\8\8_31466_wgs.tif
c:\work\rasterfacto
ry\output\8\8_25832_2.tif
Copying color table from c:\work\rasterfactory\output\8\8_31466_wgs.tif to
new f
ile.
Creating output file that is 8861P x 9557L.
Processing input file c:\work\rasterfactory\output\8\8_31466_wgs.tif.
:0...10...20...30...40...50...60...70...80...90...100 - done.
C:\work\rasterfactory\work>listgeo
c:\work\rasterfactory\output\8\8_25832_2.tif
TIFFReadDirectory: Warning, Unknown field with tag 42113 (0xa481)
encountered.
Geotiff_Information:
Version: 1
Key_Revision: 1.0
Tagged_Information:
ModelTiepointTag (2,3):
0 0 0
387002.011 5710514.84 0
ModelPixelScaleTag (1,3):
0.0564260761 0.0564260761 0
End_Of_Tags.
Keyed_Information:
GTModelTypeGeoKey (Short,1): ModelTypeProjected
GTRasterTypeGeoKey (Short,1): RasterPixelIsArea
GTCitationGeoKey (Ascii,22): "ETRS89 / UTM zone 32N"
GeogCitationGeoKey (Ascii,7): "ETRS89"
GeogAngularUnitsGeoKey (Short,1): Angular_Degree
ProjectedCSTypeGeoKey (Short,1): Unknown-25832
ProjLinearUnitsGeoKey (Short,1): Linear_Meter
End_Of_Keys.
End_Of_Geotiff.
PCS = 25832 (ETRS89 / UTM zone 32N)
Projection = 16032 (UTM zone 32N)
Projection Method: CT_TransverseMercator
ProjNatOriginLatGeoKey: 0.000000 ( 0d 0' 0.00"N)
ProjNatOriginLongGeoKey: 9.000000 ( 9d 0' 0.00"E)
ProjScaleAtNatOriginGeoKey: 0.999600
ProjFalseEastingGeoKey: 500000.000000 m
ProjFalseNorthingGeoKey: 0.000000 m
GCS: 4258/ETRS89
Datum: 6258/European Terrestrial Reference System 1989
Ellipsoid: 7019/GRS 1980 (6378137.00,6356752.31)
Prime Meridian: 8901/Greenwich (0.000000/ 0d 0' 0.00"E)
Projection Linear Units: 9001/metre (1.000000m)
Corner Coordinates:
Upper Left ( 387002.011, 5710514.843)
Lower Left ( 387002.011, 5709975.579)
Upper Right ( 387502.002, 5710514.843)
Lower Right ( 387502.002, 5709975.579)
Center ( 387252.007, 5710245.211)
--
View this message in context: http://www.nabble.com/gdalwarp-does-not-use-my-%2Btowgs84-parameter-tf4553338.html#a12993966
Sent from the GDAL - Dev mailing list archive at Nabble.com.
More information about the Gdal-dev
mailing list