On Sat, May 7, 2011 at 11:41 AM, Andrea Aime
<[email protected]>wrote:
> Hi,
> nowadays when generating .prj files for shapefile and world files
> GeoTools produces something like
> the following:
>
> PROJCS["WGS 84 / UTM zone 35S",
> GEOGCS["WGS 84",
> DATUM["World Geodetic System 1984",
> SPHEROID["WGS 84", 6378137.0, 298.257223563,
> AUTHORITY["EPSG","7030"]],
> AUTHORITY["EPSG","6326"]],
> PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
> UNIT["degree", 0.017453292519943295],
> AXIS["Geodetic latitude", NORTH],
> AXIS["Geodetic longitude", EAST],
> AUTHORITY["EPSG","4326"]],
> PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]],
> PARAMETER["central_meridian", 27.0],
> PARAMETER["latitude_of_origin", 0.0],
> PARAMETER["scale_factor", 0.9996],
> PARAMETER["false_easting", 500000.0],
> PARAMETER["false_northing", 10000000.0],
> UNIT["m", 1.0],
> AXIS["Easting", EAST],
> AXIS["Northing", NORTH],
> AUTHORITY["EPSG","32735"]]
>
> This is fine and standard, but unfortunately ESRI ArcMap does not
> understand it.
> The ESRI equivalent is:
>
>
> PROJCS["WGS_1984_UTM_Zone_35S",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],
>
> PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["Central_Meridian",27],
>
> PARAMETER["Latitude_Of_Origin",0],PARAMETER["Scale_Factor",0.9996],PARAMETER["False_Easting",500000],PARAMETER["False_Northing",10000000],
> UNIT["Meter",1]
>
> Notable differences are the way the datum, unit of measure, projection
> and projection parameters are spelled out,
> plus the fact that it's all in one line.
>
> Now, I've just committed a few changes on trunk that make the following
> code:
>
> Formattable f = (Formattable) CRS.decode("EPSG:32735", true);
> System.out.println(f.toWKT(Citations.ESRI, 0));
>
> generate this instead:
>
> PROJCS["WGS 84 / UTM zone 35S", GEOGCS["WGS 84", DATUM["D_WGS_1984",
> SPHEROID["D_WGS_1984", 6378137.0, 298.257223563,
> AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]],
> PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["Degree",
> 0.017453292519943295], AXIS["Geodetic longitude", EAST],
> AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4326"]],
> PROJECTION["Transverse_Mercator", AUTHORITY["EPSG","9807"]],
> PARAMETER["Longitude_Of_Center", 27.0],
> PARAMETER["Latitude_Of_Origin", 0.0], PARAMETER["Scale_Factor",
> 0.9996], PARAMETER["false_easting", 500000.0],
> PARAMETER["False_Northing", 10000000.0], UNIT["Meter", 1.0],
> AXIS["Easting", EAST], AXIS["Northing", NORTH],
> AUTHORITY["EPSG","32735"]]
>
> While not 1-1 equal, the datum, projection, projection parameters and
> units are spelled out the same
> way. I don't have an ArcMap to test that, but my guess would be that
> it should read it.
>
> Ah, what was the code above? First I casted the CRS into a Formattable
> object, that gives
> the developer more control on how things are converted in WKT, and
> then asked to generate
> the WKT using the ESRI aliases and 0 indentation when generating the output
> (the default is CITATIONS.ESPG, 2 instead).
>
> Now, I was wondering, shall we make this wkt format the default for
> the shapefile
> and image world stores?
>
Might be a little nicer to do a regular deprecation cycle like we normally
do for api changes since you never really know who is relying on the current
structure.
2c.
>
> Shapefile wise it would also be possible to add a store parameter or a
> store property
> (I would go for the latter, since prj file creation is a programmatic
> business), whilst
> for WorldImageWriter it might be either a creation Hint or a write
> time GeneralParameterValue.
>
> Opinions, experiences? Anyone with an ArcMap that can try this out? :-)
>
> Cheers
> Andrea
>
> --
> -------------------------------------------------------
> Ing. Andrea Aime
> GeoSolutions S.A.S.
> Tech lead
>
> Via Poggio alle Viti 1187
> 55054 Massarosa (LU)
> Italy
>
> phone: +39 0584 962313
> fax: +39 0584 962313
>
> http://www.geo-solutions.it
> http://geo-solutions.blogspot.com/
> http://www.youtube.com/user/GeoSolutionsIT
> http://www.linkedin.com/in/andreaaime
> http://twitter.com/geowolf
>
> -------------------------------------------------------
>
>
> ------------------------------------------------------------------------------
> WhatsUp Gold - Download Free Network Management Software
> The most intuitive, comprehensive, and cost-effective network
> management toolset available today. Delivers lowest initial
> acquisition cost and overall TCO of any competing solution.
> http://p.sf.net/sfu/whatsupgold-sd
> _______________________________________________
> Geotools-devel mailing list
> [email protected]
> https://lists.sourceforge.net/lists/listinfo/geotools-devel
>
--
Justin Deoliveira
OpenGeo - http://opengeo.org
Enterprise support for open source geospatial.
------------------------------------------------------------------------------
WhatsUp Gold - Download Free Network Management Software
The most intuitive, comprehensive, and cost-effective network
management toolset available today. Delivers lowest initial
acquisition cost and overall TCO of any competing solution.
http://p.sf.net/sfu/whatsupgold-sd
_______________________________________________
Geotools-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/geotools-devel
