Hi geotoolers,
I wanted to layout an issue I encountered and am not sure about.
I have a shp file created with QGIS in 4326, which results into having a
prj like:
GEOGCS["GCS_WGS_1984", DATUM["D_WGS_1984", SPHEROID["WGS_1984", 6378137.0,
298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree",
0.017453292519943295], AXIS["Longitude", EAST], AXIS["Latitude", NORTH]]
If I use the crs object from the featurestore's schema (when reading the
shp) to create a geotiff, I results into having a crs (from gdalinfo):
GEOGCRS["Greenwich", DATUM["unnamed", ELLIPSOID["unnamed",6378137,
298.257223563, LENGTHUNIT["metre",1, ID["EPSG",9001]]]], PRIMEM["Greenwich",
0, ANGLEUNIT["degree",0.0174532925199433, ID["EPSG",9122]]], CS[ellipsoidal,
2], AXIS["latitude",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433,
ID["EPSG",9122]]], AXIS["longitude",east, ORDER[2], ANGLEUNIT["degree",
0.0174532925199433, ID["EPSG",9122]]]]
Funny enough, even if the parameters seem to be right, QGIS then is not
able to understand the crs.
I dived into the Geotiff writer and found (gosh I can't find where right
now, sorry) that the test done to get the epsg code is:
*final* Set< ? *extends* Identifier> identifiers = crs
.getIdentifiers();
*final* Iterator< ? *extends* Identifier> it = identifiers
.iterator();
String code = "";
*while* (it.hasNext()) {
*final* Identifier identifier = it.next();
*final* Citation cite = identifier.getAuthority();
*if* (Citations.*identifierMatches*(cite, "EPSG")) {
code = identifier.getCode();
*break*;
}
}
which breaks in case there is no epsg definition.
I am currentl double-checking this with:
Integer epsg = CRS.*lookupEpsgCode*(crs, *false*);
which gives the right epsg.
So all in all things do work, but for some cases it might not?
I probably forgot what my original question was at this point, so I will
just go with: did anyone experience something like that?
Thanks for any input,
Cheers,
Andrea
_______________________________________________
GeoTools-GT2-Users mailing list
GeoTools-GT2-Users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
