Re: [JPP-Devel] [Fwd: Re: Projection with Geotools2 CRS]

2007-09-17 Thread Michaël Michaud 

>it is used to transform from one ellipsoid (e.g. WGS84) to another 
>ellipsoid (e.g. Clarke 1866). (actually on the geotools page they speak 
>of spheroid instead ellipsoid, where the latter is a specialisation)
>
>and - i am not really sure, but there is second method that can be used 
>for such transformation. But anyway.. as far as i know the use of 
>Bursa-Wolf is quite common.
>
>if you stay with the ellipsoid you should not need to define params (or 
>maybe they are zero?)
>  
>
Bursa-Wolf transform is used to change from one datum to another.
It is a 3-dimensional affine transformation between geocentric 
coordinates and it is not direcly related to ellipsoid, except that 
often, different datum use different ellipsoids, and you have also to 
transform from/to geocentric to/from geographic coordinates using 
appropriate ellipsoids.
There are other methods like coordinate frame rotation which is about 
the same but with inverse sign convention (which is quite vicious), or 
even a direct transformation to go from datum1/ellipsoid1/geographic 
coordinate to datum2/ellipsoid2/geographic coordinate, but I think this 
one accept only translation between datums, no rotation...

Michael

>stefan
>
>  
>
>>That is one of our most common questions:
>>- http://docs.codehaus.org/display/GEOTDOC/11+Referencing+FAQ
>>
>>My impression is that this is scary math I never quite understood. The 
>>javadocs describe it all detail (and have links to papers etc..).
>>
>>On a pragmatic side: Bursa Wolf Parameters defines a transform; and you 
>>can include the definition of a transform with your crs. At the very 
>>least it is nice to have a transform defined to take you to 
>>DefaultGeographic.WSG84.
>>
>>Jody
>>
>>-
>>This SF.net email is sponsored by: Microsoft
>>Defy all challenges. Microsoft(R) Visual Studio 2005.
>>http://clk.atdmt.com/MRT/go/vse012070mrt/direct/01/
>>___
>>Jump-pilot-devel mailing list
>>Jump-pilot-devel@lists.sourceforge.net
>>https://lists.sourceforge.net/lists/listinfo/jump-pilot-devel
>>
>>
>>
>>
>
>-
>This SF.net email is sponsored by: Microsoft
>Defy all challenges. Microsoft(R) Visual Studio 2005.
>http://clk.atdmt.com/MRT/go/vse012070mrt/direct/01/
>___
>Jump-pilot-devel mailing list
>Jump-pilot-devel@lists.sourceforge.net
>https://lists.sourceforge.net/lists/listinfo/jump-pilot-devel
>
>
>  
>


-
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2005.
http://clk.atdmt.com/MRT/go/vse012070mrt/direct/01/
___
Jump-pilot-devel mailing list
Jump-pilot-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/jump-pilot-devel

Re: [JPP-Devel] [Fwd: Re: Projection with Geotools2 CRS]

2007-09-17 Thread Stefan Steiniger 
>> just one question .. what is this bursa wolf parameter option?

it is used to transform from one ellipsoid (e.g. WGS84) to another 
ellipsoid (e.g. Clarke 1866). (actually on the geotools page they speak 
of spheroid instead ellipsoid, where the latter is a specialisation)

and - i am not really sure, but there is second method that can be used 
for such transformation. But anyway.. as far as i know the use of 
Bursa-Wolf is quite common.

if you stay with the ellipsoid you should not need to define params (or 
maybe they are zero?)

stefan

> That is one of our most common questions:
> - http://docs.codehaus.org/display/GEOTDOC/11+Referencing+FAQ
> 
> My impression is that this is scary math I never quite understood. The 
> javadocs describe it all detail (and have links to papers etc..).
> 
> On a pragmatic side: Bursa Wolf Parameters defines a transform; and you 
> can include the definition of a transform with your crs. At the very 
> least it is nice to have a transform defined to take you to 
> DefaultGeographic.WSG84.
> 
> Jody
> 
> -
> This SF.net email is sponsored by: Microsoft
> Defy all challenges. Microsoft(R) Visual Studio 2005.
> http://clk.atdmt.com/MRT/go/vse012070mrt/direct/01/
> ___
> Jump-pilot-devel mailing list
> Jump-pilot-devel@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/jump-pilot-devel
> 
> 

-
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2005.
http://clk.atdmt.com/MRT/go/vse012070mrt/direct/01/
___
Jump-pilot-devel mailing list
Jump-pilot-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/jump-pilot-devel

Re: [JPP-Devel] [Fwd: Re: Projection with Geotools2 CRS]

2007-09-17 Thread Jody Garnett 
Edgar Soldin wrote:
> Not really ;) ...
>
> after posting the snippet mentioned, i also had a look in your docs 
> and found the utility class, even posted about it in the list (see 
> attachment).
> once again .. respect to you geotools folks, you really made a nice 
> feature even more handy by creating the utility class.
>
> just one question .. what is this bursa wolf parameter option?
That is one of our most common questions:
- http://docs.codehaus.org/display/GEOTDOC/11+Referencing+FAQ

My impression is that this is scary math I never quite understood. The 
javadocs describe it all detail (and have links to papers etc..).

On a pragmatic side: Bursa Wolf Parameters defines a transform; and you 
can include the definition of a transform with your crs. At the very 
least it is nice to have a transform defined to take you to 
DefaultGeographic.WSG84.

Jody

-
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2005.
http://clk.atdmt.com/MRT/go/vse012070mrt/direct/01/
___
Jump-pilot-devel mailing list
Jump-pilot-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/jump-pilot-devel

Re: [JPP-Devel] [Fwd: Re: Projection with Geotools2 CRS]

2007-09-17 Thread Edgar Soldin 

Not really ;) ...

after posting the snippet mentioned, i also had a look in your docs and 
found the utility class, even posted about it in the list (see attachment).
once again .. respect to you geotools folks, you really made a nice 
feature even more handy by creating the utility class.


just one question .. what is this bursa wolf parameter option?

Also I stumbled across a problem converting DHDN / Gauss Kruger 3 into 
WGS84 with the GT2.3 EPSG. It was impossible due to Bursa Wolf (?*) 
parameters missing. Setting a GT2.3 parameter named 
LENIENT_DATUM_SHIFT 
 

enabled the conversion but gave inexact results. Maybe the dialogue 
should enable the user to set options like that to modify the way 
GT2.3 works?
Finally I compared the wkt strings from GT2.3 and my original  
cts.0.1a  and found some missing parameters  in the DHDN definition. 
Using the old WKT string it worked perfectly and gave nearly identical 
results. 


cheers too .. ede
--

Stefan Steiniger wrote:

 Original-Nachricht 
Betreff: Re: [jump-devel] Projection with Geotools2 CRS
Datum: Sat, 15 Sep 2007 18:22:55 -0700
Von: Jody Garnett <[EMAIL PROTECTED]>
Antwort an: JUMP Developers <[EMAIL PROTECTED]>
An: JUMP Developers <[EMAIL PROTECTED]>
Referenzen: <[EMAIL PROTECTED]>

Hi Edgar:

We have a different idea of simplicity ;-) I tend to use the utility
class to get the job done.

CoordinateReferenceSystem src = CRS.decode("EPSG:31467");
CoordinateReferenceSystem trg  = CRS.decode("EPSG:4326");
MathTransform transform = CRS.findMathTransform( src, trg );

Geometry targetGeometry = JTS.transform( sourceGeometry, transform);

This example and others are in the user guide:
- http://docs.codehaus.org/display/GEOTDOC/01+CRS+Helper+Class

Cheers,
Jody
  

Hello All,

because of all the comments related to projection (interfaces) , here 
is an example of the simplicity in using GT2 CRS for 
projection/transformation


SNIPPET-->
   CoordinateOperationFactory coFactory = 
FactoryFinder.getCoordinateOperationFactory(null);

   MathTransform transform;
   src = 
FactoryFinder.getCRSFactory(null).createFromWKT("PROJCS[\"DHDN / 
Gauss-Kruger Zone 3\", GEOGCS[\"DHDN\", DATUM[\"Deutsches 
Hauptdreiecksnetz\", SPHEROID[\"Bessel 1841\", 6377397.155, 
299.1528128, AUTHORITY[\"EPSG\",\"7004\"]], TOWGS84[598.1, 73.7, 
418.2, 0.2018, 0.04495, -2.45499974, 
6.7], AUTHORITY[\"EPSG\",\"6314\"]], PRIMEM[\"Greenwich\", 0.0, 
AUTHORITY[\"EPSG\",\"8901\"]], UNIT[\"degree of 
angle\",0.017453292519943295], AXIS[\"Geodetic longitude\",EAST], 
AXIS[\"Geodetic latitude\",NORTH], AUTHORITY[\"EPSG\",\"4314\"]], 
PROJECTION[\"Transverse_Mercator\"], PARAMETER[\"semi_major\", 
6377397.155], PARAMETER[\"semi_minor\", 6356078.962818189], 
PARAMETER[\"central_meridian\", 8.991], 
PARAMETER[\"latitude_of_origin\", 0.0], PARAMETER[\"scale_factor\", 
1.0], PARAMETER[\"false_easting\", 350.0], 
PARAMETER[\"false_northing\", 0.0], UNIT[\"metre\",1.0], 
AXIS[\"Easting\",EAST], AXIS[\"Northing\",NORTH], 
AUTHORITY[\"EPSG\",\"31467\"]]");
//src = FactoryFinder.getCRSAuthorityFactory("EPSG", 
null).createCoordinateReferenceSystem("EPSG:31467");
   trg = FactoryFinder.getCRSAuthorityFactory("EPSG", 
null).createCoordinateReferenceSystem("EPSG:4326");

   CoordinateOperation co = coFactory.createOperation(src, trg);
 transform = co.getMathTransform();
<--SNIPPET

the snippet is only rough cut out of my tests. It will not compile in 
this state, but still it shows how easy it is with GT2 to get a 
transformation (from WKT or EPSG ID, as you choose) for projecting 
coordinates. The EPSG entries of the factory are of course 
request-/listable by name(description).


just trying to lobby for the GT2 library, as I get the impression the 
list is a bit afraid of it ;).


regards ede




___
jump-devel mailing list
[EMAIL PROTECTED]
http://lists.refractions.net/mailman/listinfo/jump-devel



-
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2005.
http://clk.atdmt.com/MRT/go/vse012070mrt/direct/01/
___
Jump-pilot-devel mailing list
Jump-pilot-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/jump-pilot-devel
  



--
public class WhoDidIt{ // A comment. I love comments 
 private static Person sender;


 public static void main (String[] foo){

 sender = new Person();
 sender.setName(new String[]{"Edgar", "Soldin"});

 Address address = new Address();
 address.setStreet("Stadtweg 119");
 address.setZip(39116);
 address.setCity("Magdeburg");
 address.setCountry("Germany");

 sender.setAddress(address);

 sender.setMobilePhone(" +49(0)171-2782880 ");
 sender.setWebSiteUrl(" http://www.soldin

[JPP-Devel] [Fwd: Re: Projection with Geotools2 CRS]

2007-09-16 Thread Stefan Steiniger 
 Original-Nachricht 
Betreff: Re: [jump-devel] Projection with Geotools2 CRS
Datum: Sat, 15 Sep 2007 18:22:55 -0700
Von: Jody Garnett <[EMAIL PROTECTED]>
Antwort an: JUMP Developers <[EMAIL PROTECTED]>
An: JUMP Developers <[EMAIL PROTECTED]>
Referenzen: <[EMAIL PROTECTED]>

Hi Edgar:

We have a different idea of simplicity ;-) I tend to use the utility
class to get the job done.

CoordinateReferenceSystem src = CRS.decode("EPSG:31467");
CoordinateReferenceSystem trg  = CRS.decode("EPSG:4326");
MathTransform transform = CRS.findMathTransform( src, trg );

Geometry targetGeometry = JTS.transform( sourceGeometry, transform);

This example and others are in the user guide:
- http://docs.codehaus.org/display/GEOTDOC/01+CRS+Helper+Class

Cheers,
Jody
> Hello All,
>
> because of all the comments related to projection (interfaces) , here 
> is an example of the simplicity in using GT2 CRS for 
> projection/transformation
>
> SNIPPET-->
>CoordinateOperationFactory coFactory = 
> FactoryFinder.getCoordinateOperationFactory(null);
>MathTransform transform;
>src = 
> FactoryFinder.getCRSFactory(null).createFromWKT("PROJCS[\"DHDN / 
> Gauss-Kruger Zone 3\", GEOGCS[\"DHDN\", DATUM[\"Deutsches 
> Hauptdreiecksnetz\", SPHEROID[\"Bessel 1841\", 6377397.155, 
> 299.1528128, AUTHORITY[\"EPSG\",\"7004\"]], TOWGS84[598.1, 73.7, 
> 418.2, 0.2018, 0.04495, -2.45499974, 
> 6.7], AUTHORITY[\"EPSG\",\"6314\"]], PRIMEM[\"Greenwich\", 0.0, 
> AUTHORITY[\"EPSG\",\"8901\"]], UNIT[\"degree of 
> angle\",0.017453292519943295], AXIS[\"Geodetic longitude\",EAST], 
> AXIS[\"Geodetic latitude\",NORTH], AUTHORITY[\"EPSG\",\"4314\"]], 
> PROJECTION[\"Transverse_Mercator\"], PARAMETER[\"semi_major\", 
> 6377397.155], PARAMETER[\"semi_minor\", 6356078.962818189], 
> PARAMETER[\"central_meridian\", 8.991], 
> PARAMETER[\"latitude_of_origin\", 0.0], PARAMETER[\"scale_factor\", 
> 1.0], PARAMETER[\"false_easting\", 350.0], 
> PARAMETER[\"false_northing\", 0.0], UNIT[\"metre\",1.0], 
> AXIS[\"Easting\",EAST], AXIS[\"Northing\",NORTH], 
> AUTHORITY[\"EPSG\",\"31467\"]]");
> //src = FactoryFinder.getCRSAuthorityFactory("EPSG", 
> null).createCoordinateReferenceSystem("EPSG:31467");
>trg = FactoryFinder.getCRSAuthorityFactory("EPSG", 
> null).createCoordinateReferenceSystem("EPSG:4326");
>CoordinateOperation co = coFactory.createOperation(src, trg);
>  transform = co.getMathTransform();
> <--SNIPPET
>
> the snippet is only rough cut out of my tests. It will not compile in 
> this state, but still it shows how easy it is with GT2 to get a 
> transformation (from WKT or EPSG ID, as you choose) for projecting 
> coordinates. The EPSG entries of the factory are of course 
> request-/listable by name(description).
>
> just trying to lobby for the GT2 library, as I get the impression the 
> list is a bit afraid of it ;).
>
> regards ede
>

___
jump-devel mailing list
[EMAIL PROTECTED]
http://lists.refractions.net/mailman/listinfo/jump-devel



-
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2005.
http://clk.atdmt.com/MRT/go/vse012070mrt/direct/01/
___
Jump-pilot-devel mailing list
Jump-pilot-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/jump-pilot-devel