[Gdal-dev] Image Reprojection

[Matt Lynch]

I ran gdalinfo on out.tif and saw that I was missing gdal.dll.  I
grabbed that off the gdall5a release.  Now I get the following - again
it must be that I need to feed proj4 the specifics of a 10TM projection.

C:\gdal>gdalwarp -t_srs "+proj=latlong +datum=WGS84" jeffpan2000_10m.tif
out.tif

ERROR 1: Can't find requested entry point: pj_strerrno

ERROR 1: Failed to import coordinate system `'.
:0...10...20...30...40...50...60...70...80...90...100 - done.

When I run gdalinfo on out.tif now I get ...

C:\gdal>gdalinfo out.tif
Driver: GTiff/GeoTIFF
Size is 2000, 2000
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.2572235630016,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin = (644995.948217,5518006.066294)
Pixel Size = (10.000000,-10.000000)
Corner Coordinates:
Upper Left  (  644995.948, 5518006.066) ERROR 1: Can't find requested
entry poin
t: pj_strerrno

(644995d56'4294967349.58"E,5518006d 3'21474836538.66"N)
Lower Left  (  644995.948, 5498006.066)
(644995d56'4294967349.58"E,5498006d 3'21
474836538.66"N)
Upper Right (  664995.948, 5518006.066)
(664995d56'4294967349.58"E,5518006d 3'21
474836538.66"N)
Lower Right (  664995.948, 5498006.066)
(664995d56'4294967349.58"E,5498006d 3'21
474836538.66"N)
Center      (  654995.948, 5508006.066)
(654995d56'4294967349.58"E,5508006d 3'21
474836538.66"N)
Band 1 Block=2000x4 Type=Byte, ColorInterp=Gray


Thanks,

Matt

-----Original Message-----
From: gdal-dev-admin at remotesensing.org
[mailto:gdal-dev-admin at remotesensing.org] On Behalf Of Matt Lynch
Sent: Thursday, April 03, 2003 11:24 AM
To: gdal-dev at remotesensing.org
Subject: RE: [Gdal-dev] Image Reprojection


Hi Frank,

>I believe this is a commandline parsing issue.  The '+proj=latlong
+datum=WGS84' would be parsed as a single argument by a Unix shell, but
the DOS shell doesn't work the same.  

That makes sense.  I added a forth option and used double quotes (I can
recall using them in the past within dos to bypass the whitespace).

Now I get
------------------------------------------------------------------------
--------
C:\gdal>gdalwarp -t_srs "+proj=latlong +datum=WGS84" jeffpan2000_10m.tif
out.tif

ERROR 1: Failed to import coordinate system `'.
Creating output file is that 2000P x 2000L.
ERROR 1: Failed to import coordinate system `'.
:0...10...20...30...40...50...60...70...80...90...100 - done.
------------------------------------------------------------------------
--------

The output file is created, but it doesn't look like the conversion
happens.  The line 'FAILED to import coordinate system' is probably the
source of my trouble now (image is in 10TM).

How can I specify a input coordinate system like 10TM

Thanks,

Matt

---
Outgoing mail is certified Virus Free.
Checked by AVG anti-virus system (http://www.grisoft.com).
Version: 6.0.467 / Virus Database: 266 - Release Date: 4/1/2003