[gdal-dev] Re: Reprojecting w/ gdalwarp .. cant figure it out

[fork]

Frank Warmerdam <warmerdam <at> pobox.com> writes:

> Try:
> gdalwarp -s_srs utm10n.prj -t_srs ESRI::spwasth.prj \
>    -co "TFW=YES" ortho_imagery/ortho_1-1_1n_s_wa033_2009_1.tif foo.tif
> 
> The ESRI:: prefix before the filename tells GDAL/OGR to convert from
> ESRI format into "normal" WKT.

So, as I said in my other email, this seems to work nicely, and lines up with a
shapefile of interest when I use ArcMap.

However, I am trying to get the raster image to line up with a PostGIS geometry.
 I would think that the bounding box on the geometry and the tfw would be really
close, but I get this from PostGIS:


tiger2010=# select id, getbbox(the_geom) from annexpoly;
      id      |                      getbbox
--------------+----------------------------------------------------
 my_id | BOX(1206737.375 866782.5625,1234465.5 882399.875)

While the new foo.tfw is:

10.7547912261
0.0000000000
0.0000000000
-10.7547912261
358134.0684926199
272894.5606149968

Just for interest, here is my geometry_columns of interest:

tiger2010=# select f_table_name, f_geometry_column, srid from geometry_columns;
   f_table_name   | f_geometry_column |  srid
------------------+-------------------+---------
 annexpoly        | the_geom          | 9102349

Any ideas ....?? It is very possible (likely) I made a stupid mistake, so I TIA
for patience in working out projection subtleties with me....