>> Are you sure this transform is in the EPSG tables?
Probably not. The grid was calculated specifically for the state of
Schleswig Holstein in Germany. In my head it is like this: I register a
user-defined transformation in geotools of 31467 to 4647. This is always
done taking into account the grid I provide (the specific gsb file). Am I
wrong?
Benjamin
Von: "Ian Turton" <ijtur...@gmail.com>
An: "Benjamin Sokolowski" <benjamin.sokolow...@data-experts.de>
Kopie: "geotools-users" <geotools-gt2-users@lists.sourceforge.net>
Datum: 08.04.2020 14:33
Betreff: Re: Re: Re: Re: [Geotools-gt2-users] Coordinate
Transformation
On further investigation at epgs.io I can get a gsb transform for
epsg:31467->epsg:4326 (WGS84) using BETA2007.gsb but not the epsg:4258
(ETRS) which is the base for EPSG:4647 (ETRS89 / UTM zone 32N (zE-N)).
I can't find any reference to the need for a grid file for this transform
so that may be why it isn't picked up by the code? Are you sure this
transform is in the EPSG tables?
Ian
On Wed, 8 Apr 2020 at 05:55, Benjamin Sokolowski <
benjamin.sokolow...@data-experts.de> wrote:
With and without gsb file the same output:
CONCAT_MT[INVERSE_MT[PARAM_MT["Transverse_Mercator",
PARAMETER["semi_major", 6377397.155],
PARAMETER["semi_minor", 6356078.962818189],
PARAMETER["central_meridian", 9.0],
PARAMETER["latitude_of_origin", 0.0],
PARAMETER["scale_factor", 1.0],
PARAMETER["false_easting", 3500000.0],
PARAMETER["false_northing", 0.0]]],
PARAM_MT["Ellipsoid_To_Geocentric",
PARAMETER["dim", 2],
PARAMETER["semi_major", 6377397.155],
PARAMETER["semi_minor", 6356078.962818189]],
PARAM_MT["Affine",
PARAMETER["num_row", 4],
PARAMETER["num_col", 4],
PARAMETER["elt_0_0", 1.00000255],
PARAMETER["elt_0_1", 0.00001356027323924231],
PARAMETER["elt_0_2", 0.00000027634450290912],
PARAMETER["elt_0_3", 612.4],
PARAMETER["elt_1_0", -0.00001356027323924231],
PARAMETER["elt_1_1", 1.00000255],
PARAMETER["elt_1_2", 0.00000026180005538759],
PARAMETER["elt_1_3", 77.0],
PARAMETER["elt_2_0", -0.00000027634450290912],
PARAMETER["elt_2_1", -0.00000026180005538759],
PARAMETER["elt_2_2", 1.00000255],
PARAMETER["elt_2_3", 440.2]],
PARAM_MT["Geocentric_To_Ellipsoid",
PARAMETER["dim", 2],
PARAMETER["semi_major", 6378137.0],
PARAMETER["semi_minor", 6356752.314140356]],
PARAM_MT["Transverse_Mercator",
PARAMETER["semi_major", 6378137.0],
PARAMETER["semi_minor", 6356752.314140356],
PARAMETER["central_meridian", 9.0],
PARAMETER["latitude_of_origin", 0.0],
PARAMETER["scale_factor", 0.9996],
PARAMETER["false_easting", 32500000.0],
PARAMETER["false_northing", 0.0]]]
which i think, is the default transformation of geotools. Here's the gsb
file
Benjamin
Von: "Ian Turton" <ijtur...@gmail.com>
An: "Benjamin Sokolowski" <benjamin.sokolow...@data-experts.de>
Kopie: "geotools-users" <geotools-gt2-users@lists.sourceforge.net>
Datum: 07.04.2020 15:53
Betreff: Re: Re: Re: [Geotools-gt2-users] Coordinate Transformation
Can you try this program with and without the GSB file
import org.geotools.referencing.CRS;
import org.opengis.referencing.FactoryException;
import org.opengis.referencing.NoSuchAuthorityCodeException;
import org.opengis.referencing.crs.CoordinateReferenceSystem;
import org.opengis.referencing.operation.MathTransform;
public class gridShiftTest {
public static void main(String[] args) throws
NoSuchAuthorityCodeException, FactoryException {
CoordinateReferenceSystem sourceCRS = CRS.decode( "EPSG:31467" );
CoordinateReferenceSystem targetCRS = CRS.decode( "EPSG:4647" );
MathTransform transform = CRS.findMathTransform(sourceCRS, targetCRS);
System.out.println(transform);
}
}
Or mail me the gsb file if you can and I'll try - I don't seem to have a
grid file for this transformation on my system because I've mainly used UK
and USA ones in the past.
Ian
On Tue, 7 Apr 2020 at 13:21, Benjamin Sokolowski <
benjamin.sokolow...@data-experts.de> wrote:
Hi Ian,
even with the .gsb in the new folder with your given name, no changes at
all. Same wrong results. Do you have another idea what i'm missing? The
"epsg_operations.properties" is not mandatory right? I don't have that
right now.
Benjamin
Von: "Ian Turton" <ijtur...@gmail.com>
An: "Benjamin Sokolowski" <benjamin.sokolow...@data-experts.de>
Kopie: "geotools-users" <geotools-gt2-users@lists.sourceforge.net>
Datum: 07.04.2020 11:33
Betreff: Re: Re: [Geotools-gt2-users] Coordinate Transformation
I've never tried adding a gsb file directly to the jar file, I usually
create the folder
src/main/resources/org/geotools/referencing/factory/gridshift in my
project and store the file in there.
Ian
On Tue, 7 Apr 2020 at 10:14, Benjamin Sokolowski <
benjamin.sokolow...@data-experts.de> wrote:
Hi Ian,
thanks for the reply. Found the discussion at stackexchange already, but
it didn't help me. Maybe i didn't understand it correctly.
>> where exactly did you put it?
Do you mean in which jar-File? gt-referencing-16.0.jar
I only have geotools in the build path, not the sources.
Benjamin
Von: "Ian Turton" <ijtur...@gmail.com>
An: "Benjamin Sokolowski" <benjamin.sokolow...@data-experts.de>
Kopie: "geotools-users" <geotools-gt2-users@lists.sourceforge.net>
Datum: 07.04.2020 11:06
Betreff: Re: [Geotools-gt2-users] Coordinate Transformation
Have a look at
https://gis.stackexchange.com/questions/313512/specifying-epsg-transformation-method-in-geotools/313523#313523
for an example.
When you say you added the .gsb to
'org.geotools.referencing.factory.gridshift' where exactly did you put it?
Ian
On Tue, 7 Apr 2020 at 09:33, Benjamin Sokolowski <
benjamin.sokolow...@data-experts.de> wrote:
Hello,
first time posting. Hope you guys can help.
I want to transform coordinates (point geometry) from Gauß Krueger Zone 3
to ETRS:89 UTM32 using geotools and a grid file, which i already have. How
exactly do i do this?
I already got the codelines:
CoordinateReferenceSystem sourceCRS = CRS.decode( "EPSG:31467" );
CoordinateReferenceSystem targetCRS = CRS.decode( "EPSG:4647" );
MathTransform transformation = CRS.findMathTransform( sourceCRS,
targetCRS );
Geometry result = JTS.transform( geometry, transformation );
and put the .gsb to 'org.geotools.referencing.factory.gridshift', but that
doesn't seem to work. Geotools is still using the default Transformation,
not the more precise with the grid file.
Best regards
Benjamin Sokolowski_______________________________________________
GeoTools-GT2-Users mailing list
GeoTools-GT2-Users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
--
Ian Turton
--
Ian Turton
--
Ian Turton
--
Ian Turton
_______________________________________________
GeoTools-GT2-Users mailing list
GeoTools-GT2-Users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geotools-gt2-users
