On Mon, Apr 6, 2015 at 1:59 PM, Peter Borissow <peter.boris...@yahoo.com>
wrote:
> Hello Andrea-
> I have attached an updated copy of the CylindricalEqualArea class
> along with a test script. The CylindricalEqualArea class supports 2 "alias"
> projections:
>
> - World Behrmann (ESRI:54017)
> - Lambert Cylindrical Equal Area
>
>
> You will need to add 3 lines to the
> org.geotools.referencing.operation.MathTransformProvider for it to work:
>
>
> org.geotools.referencing.operation.projection.CylindricalEqualArea$Provider
>
> org.geotools.referencing.operation.projection.CylindricalEqualArea$BehrmannProvider
>
> org.geotools.referencing.operation.projection.CylindricalEqualArea$LambertCylindricalEqualAreaProvider
>
> The CylindricalEqualArea class is supposed to support both spheriods and
> ellipsoids. I have verified that the ellipsoidal calculations are correct
> using the Behrmann projection (ESRI:54017) which uses the WGS84 ellipsoid.
>
Odd, I've applied the projection and tests you provided (on GeoTools
master, we don't accept new features on 10.x), and the Berhman case is
failing:
TransformException at line 33:
Expected GeneralDirectPosition[964862.8, 1269436.74]
but got GeneralDirectPosition[964862.802508965, 1269436.7435378402]
If I reduce the precision for X to 0.01 then I have the following:
TransformException at line 45:
Expected GeneralDirectPosition[1.254321643E7, 6351420.0] but got
GeneralDirectPosition[1.2543216432616547E7, 6351419.997326233]
because a few lines above you setup a 0.001 tolerance for Y... oh well,
I've reduced it all to 0.01 and
this part works.
Then it gets to the autalic sphere, and does not even start working:
OperationNotFoundException at line 82: Bursa wolf parameters required.
I've added TOWGS84[0,0,0,0,0,0,0], to the autalic sphere datum (not sure
if correct, but it gets things going) and
finally see:
TransformException at line 85:
Expected GeneralDirectPosition[963010.77, 1277505.4] but got
GeneralDirectPosition[963010.7746492706, 1269160.2484886926]
Indeed quite off. This is happening because cs2cs is not applying any
transformation to lon and lat,
it is just taking them as valid values on the autalic sphere without
compensating for the difference
in flattening:
cs2cs -v +datum=WGS84 +proj=lonlat +to +proj=lonlat +a=6371228 +b=6371228
+units=m +no_defs# ---- From Coordinate System ----
#Lat/long (Geodetic)
#
# +datum=WGS84 +proj=lonlat +ellps=WGS84 +towgs84=0,0,0
# ---- To Coordinate System ----
#Lat/long (Geodetic)
#
# +proj=lonlat +a=6371228 +b=6371228 +units=m +no_defs
10 10
10dE 10dN 0.000
(the lat should have been something like 9.93 instead)
And I'm not sure why it's doing that... theoretically it should compensate
for the different size
of the elliposoid/spheroid, the behavior we're seeing should be something
we get only
with +@nadgrids=null instead:
http://proj.maptools.org/faq.html
I don't have an answer here, so I've asked on the proj.4 mailing list
Cheers
Andrea
--
==
GeoServer Professional Services from the experts! Visit
http://goo.gl/NWWaa2 for more information.
==
Ing. Andrea Aime
@geowolf
Technical Lead
GeoSolutions S.A.S.
Via Poggio alle Viti 1187
55054 Massarosa (LU)
Italy
phone: +39 0584 962313
fax: +39 0584 1660272
mob: +39 339 8844549
http://www.geo-solutions.it
http://twitter.com/geosolutions_it
*AVVERTENZE AI SENSI DEL D.Lgs. 196/2003*
Le informazioni contenute in questo messaggio di posta elettronica e/o
nel/i file/s allegato/i sono da considerarsi strettamente riservate. Il
loro utilizzo è consentito esclusivamente al destinatario del messaggio,
per le finalità indicate nel messaggio stesso. Qualora riceviate questo
messaggio senza esserne il destinatario, Vi preghiamo cortesemente di
darcene notizia via e-mail e di procedere alla distruzione del messaggio
stesso, cancellandolo dal Vostro sistema. Conservare il messaggio stesso,
divulgarlo anche in parte, distribuirlo ad altri soggetti, copiarlo, od
utilizzarlo per finalità diverse, costituisce comportamento contrario ai
principi dettati dal D.Lgs. 196/2003.
The information in this message and/or attachments, is intended solely for
the attention and use of the named addressee(s) and may be confidential or
proprietary in nature or covered by the provisions of privacy act
(Legislative Decree June, 30 2003, no.196 - Italy's New Data Protection
Code).Any use not in accord with its purpose, any disclosure, reproduction,
copying, distribution, or either dissemination, either whole or partial, is
strictly forbidden except previous formal approval of the named
addressee(s). If you are not the intended recipient, please contact
immediately the sender by telephone, fax or e-mail and delete the
information in this message that has been received in error. The sender
does not give any warranty or accept liability as the content, accuracy or
completeness of sent messages and accepts no responsibility for changes
made after they were sent or for other risks which arise as a result of
e-mail transmission, viruses, etc.
-------------------------------------------------------
------------------------------------------------------------------------------
BPM Camp - Free Virtual Workshop May 6th at 10am PDT/1PM EDT
Develop your own process in accordance with the BPMN 2 standard
Learn Process modeling best practices with Bonita BPM through live exercises
http://www.bonitasoft.com/be-part-of-it/events/bpm-camp-virtual- event?utm_
source=Sourceforge_BPM_Camp_5_6_15&utm_medium=email&utm_campaign=VA_SF
_______________________________________________
GeoTools-Devel mailing list
GeoTools-Devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geotools-devel
