[gdal-dev] transforming epsg:900913 to Calif. Teale Albers y-shift problem

[Jennifer Strahan]

Hello All,

I'm using OpenLayers and Featureserver to draw shapes on top of a Google 
base map and save into my PostGIS database.  The geometry column is 
using ESPG: 900913.  Now, I'd like to export out a shapefile for use in 
ArcGIS, using the California Teale Albers projection. I tried the 
following and am ending up with an ~18,000 meter Y shift and a little 
bit of an X shift.
$ret = shell_exec ("ogr2ogr -f 'ESRI Shapefile' 
/var/tmp/shapefile/parkinfo_edits.shp 'PG:dbname=parkinfo_edit 
host=localhost user=pgsql password=*****' -s_srs EPSG:900913 -sql 
'Select * from holdingsfeb08_new' -t_srs '+proj=aea +lat_1=34.00 
+lat_2=40.50 +lat_0=0.00 +lon_0=-120.00 +x_0=0.000 +y_0=-4000000.000 
+ellps=GRS80 +units=m +datum=NAD83'");

I went to spatialreference.org and have been trying various 
projections.  In each case I'm getting the same results.
Do you have any suggestions?

Thanks,
Jennifer