Hi Frank:
You mentioned that I have not provided a datum definition for the source
coordinate system, only an ellipsoid definition. Could you tell me how
I can provide the datum definition to get the same result from PROJ4.5?
The command I used is:
ogr2ogr -t_srs "+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96
+x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs " -s_srs
"+proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=33 +lat_2=45 +lat_0=40
+lon_0=-97" wrf12km_nc_g152.shp wrf12km_nc.shp
Thank you very much,
Limei
Limei,
Generally speaking, when going between "real" datums and these spherical
based projections, the right thing to do is to apply no shift which is
why
the default behavior was changed in PROJ 4.6.
But to get consistent overlay, you really need to know quite a bit about
your data.
Best regards,
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
