Paul Chavent | 23 Mar 23:14 2015
Picon

How to offset height above ellipsoid.

Hi.

I try to convert lat/long coordinates to eqc. The origin of the eqc projection plane is at (43.315959
1.402896 320.000). My eqc specification take into account the altitude of the projection plane, so the
plane should be tangeant to a sphere with a radius of 6368437.553703493.

The specification of the translation is then :

+proj=longlat +datum=WGS84 +no_defs +to +proj=eqc +a=6368437.553703493 +lat_ts=43.315959
+lat_0=43.315959 +lon_0=1.402896 +no_defs

At this point, if i give an altitude in the LLA input, it isn't translated since no datum have been specified
for my destination projection.

For datum conversion, +towgs84 and +nadgrids seems to be the preferred choice for horizontal datum
shifting, although +geoidgrid seems the preferred way for vertical datum shift.

In order to add a constant offset to the LLA altitude, i could provide a "gtx" file with only one value for the
whole lat/long space.

I wonder if it is a good practice, or if there is other way of doing it ?

Wouldn't be possible to specify the offset on the description string for such a simple case (that also would
be easier in the case of the usage of the c API with frequent altitude offset changes) ?

The expected result is

cs2cs -rs \
     +proj=longlat \
     +datum=WGS84 \
(Continue reading)

Richard Greenwood | 18 Mar 02:59 2015
Picon

why doesn't the epsg file reference the hpgn grids?

The epsg file included with proj has definitions for hpgn coordinate systems that do not use the grid shift files, for example:
# NAD83(HARN) / Wyoming West (ftUS)
<3758> +proj=tmerc +lat_0=40.5 +lon_0=-110.0833333333333 +k=0.9999375 +x_0=800000 +y_0=100000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs  <>

So there is no datum transformation applied. But proj has the appropriate harn grid shift files here:

And changing the above definition to include the grid shift files like below seems more appropriate:
# NAD83(HARN) / Wyoming West (ftUS)
<3758> +proj=tmerc +lat_0=40.5 +lon_0=-110.0833333333333 +k=0.9999375 +x_0=800000 +y_0=100000 +ellps=GRS80 +nadgrids=wyhpgn.gsb +units=us-ft +no_defs  <>

Does the proj epsg file get created from data at www.epsg.org? And is that data lacking the reference to the grid shift files, or is it "lost in translation"? It seems like there is the wonderful resource in the harn grids that isn't getting used, but maybe I'm overlooking something.

Rich

--
Richard W. Greenwood, PLS
www.greenwoodmap.com
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Howard Butler | 13 Mar 19:27 2015
Picon

Migrate to github?

All,

In my role as self-appointed maintainer, I became quickly aware that my brain had been poisoned by github's
ticket and development workflow. I found myself rather hamstrung without it in the context of doing the
proj 4.9.1 release.

I would like to migrate the proj Trac and Subversion instances, including all tickets, revision control,
branches, and tags, to a github repository at http://github.com/OSGeo/proj and decommission them once
there's sign-off that the migration looks good.

Unless I hear significant pushback *against* doing this, consider this message your notice that the
effort is ongoing, and I will notify once a test migration is complete for people to look around. At that
time I will then motion to switch over once folks are satisfied.

Note that this same migration was done in the past by Thomas Bonfort for MapServer, with a much larger Trac
instance (which itself was preserved from a Bugzilla instance). It's annoying to have software survive
multiple generations of bug tracking and revision control software, but I guess it's a good thing that the
software is useful enough to warrant the migration in the first place.

Howard

_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj

Adrien Barbaste | 12 Mar 14:29 2015
Picon

How to create a new EPSG for a local system

Hello,

 

I’m doing an internship for an archaeology school in Greece and I have a problem in creating of a new EPSG for local coordinates.

 

I use Proj4js to work on a WebGIS and I try to determine a new EPSG for the layers that are in local coordinates.

I've all the layers in a local system built by the archeologists and some of this points in the greek coordinates system too (EGSA87).

 

First, I can calculate a 2D transformation between the local coordinates and the greek coordinates. So I obtained the 5 parameters of transformation : the 3 translations (Tx Ty and Tz), the scale factor (K) and the angle of rotation (Rz).

 

Then, I don’t know how to create a new EPSG, because I tried with a transverse mercator projection as the greek projection, but I think that it’s impossible to insert a rotation in a proj4js (only with "+axis=", but I cannot choose the angle)

So, I tried to use an oblique mercator projection with "+alpha=" (the rotation angle) and "+gamma=" (the meridian convergence)

 

The results aren't good. I've an error of less than 3cm on an archaeological site when the rotation is very low, but when I try this method on another site where the rotation angle is about 33 degrees, I have an error that can reach 15-20 centimeters.

 

I’d like to know how I can build a good EPSG for the local coordinates system ?

 

Thank you for your help

 

_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Roger Bivand | 11 Mar 14:58 2015
Picon

SO names in 4.9.1 and current trunk

Is anyone else seeing a change in the naming of shared objects?

In PROJ trunk of January 21, 2015, and in the abandoned 4.9.0RC1 I see:

$ ls -l lib*
-rw-rw-r--. 1 rsb rsb 3193442 Jan 21 09:54 libproj.a
lrwxrwxrwx. 1 rsb rsb      13 Jan 21 09:54 libproj.la -> ../libproj.la
-rw-rw-r--. 1 rsb rsb     953 Jan 21 09:54 libproj.lai
lrwxrwxrwx. 1 rsb rsb      16 Jan 21 09:54 libproj.so -> libproj.so.0.8.0
lrwxrwxrwx. 1 rsb rsb      16 Jan 21 09:54 libproj.so.0 ->
   libproj.so.0.8.0
-rwxrwxr-x. 1 rsb rsb 1584463 Jan 21 09:54 libproj.so.0.8.0

in current trunk:

$ ls -l lib*
-rw-rw-r--. 1 rsb rsb 3205082 Mar 11 14:49 libproj.a
lrwxrwxrwx. 1 rsb rsb      13 Mar 11 14:49 libproj.la -> ../libproj.la
-rw-rw-r--. 1 rsb rsb     955 Mar 11 14:49 libproj.lai
lrwxrwxrwx. 1 rsb rsb      16 Mar 11 14:49 libproj.so -> libproj.so.9.0.0
lrwxrwxrwx. 1 rsb rsb      16 Mar 11 14:49 libproj.so.9 ->
  libproj.so.9.0.0
-rwxrwxr-x. 1 rsb rsb 1593661 Mar 11 14:49 libproj.so.9.0.0

and in the 4.9.1 tarball:

$ ls -l lib*
-rw-rw-r--. 1 rsb rsb 3206058 Mar 11 12:36 libproj.a
lrwxrwxrwx. 1 rsb rsb      13 Mar 11 12:36 libproj.la -> ../libproj.la
-rw-rw-r--. 1 rsb rsb     955 Mar 11 12:36 libproj.lai
lrwxrwxrwx. 1 rsb rsb      16 Mar 11 12:36 libproj.so -> libproj.so.9.0.0
lrwxrwxrwx. 1 rsb rsb      16 Mar 11 12:36 libproj.so.9 ->
  libproj.so.9.0.0
-rwxrwxr-x. 1 rsb rsb 1593661 Mar 11 12:36 libproj.so.9.0.0

The consequence is that downstream software linking to libproj.so.0 gets 
lost. Should libproj.so.9.0.0 be libproj.so.0.9.0? (All cases Fedora 21, 
compiled locally using autotools, not Cmake)

Roger

--

-- 
Roger Bivand
Department of Economics, Norwegian School of Economics,
Helleveien 30, N-5045 Bergen, Norway.
voice: +47 55 95 93 55; fax +47 55 95 91 00
e-mail: Roger.Bivand <at> nhh.no

_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj

Howard Butler | 10 Mar 20:44 2015
Picon

proj 4.9.1 released

All,

proj 4.9.1 is released as approved by the MetaCRS list. It cumulates a number of bug fixes and small updates,
and it went through five release candidates before being declared final. 

http://download.osgeo.org/proj/proj-4.9.1.tar.gz

Thanks,

Howard

> o 4.9.0RC2 release was abandoned because it was not promoted in a 
>   timely fashion. Subsequent maintenance of tickets has continued, 
>   and a new 4.9.1 release was issued in its place.
> 
> o Implement inverse solution for Winkel Tripel from Drazan Tutic #250
> 
> o More CMake configuration tweaks. The CMake configuration is probably 
>   not at feature parity with the autotools builds at this point but it 
>   is converging #256
> 
> o Tweak initialization ordering around setlocal which may have caused
>   issues #237
> 
> o Support out-of-tree autoconf builds more completely #247
> 
> o Fix NaN handling by geod_inverse and geod_polygon_addedge #251 & #253
> 
> o Update config.sub and config.guess #257
> 
> o Adapt Charles Karney's CMake patches for smoother build #258
> 
> o Define default PROJ_LIB location for CMake compilation #261
> 
> o Fix Windows compilation on PJ_aitoff.c
> 
> o Align CMake SOVERSION with autotools #263
> 
> o Regenerate nad/epsg with GDAL r28536 to avoid precision loss in TOWGS84
>   parameters, e.g. on Amersfoort / RD EPSG:4289 (#260)
> 
> o Add CMake project-config.cmake scripts (#264 from Charles Karney)
> 
> o Dial back test sensitivity #255
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj

Roger Bivand | 9 Mar 10:06 2015
Picon

RC5 - SOVERSION

Hi,

Should I be seeing on Unix (Fedora 21) using autotools:

$ ls -l /usr/local/lib/*proj*
-rw-r--r--. 1 root root 3206058 Mar  8 20:00 /usr/local/lib/libproj.a
-rwxr-xr-x. 1 root root     955 Mar  8 20:00 /usr/local/lib/libproj.la
lrwxrwxrwx. 1 root root      16 Mar  8 20:00 /usr/local/lib/libproj.so ->
libproj.so.9.0.0
lrwxrwxrwx. 1 root root      16 Mar  8 20:00 /usr/local/lib/libproj.so.9 ->
libproj.so.9.0.0
-rwxr-xr-x. 1 root root 1593661 Mar  8 20:00 /usr/local/lib/libproj.so.9.0.0

while with the 4.9.0 RC I see:

-rw-r--r--. 1 root root 3193442 Jan 21 15:49 /usr/local/lib/libproj.a
-rwxr-xr-x. 1 root root     953 Jan 21 15:49 /usr/local/lib/libproj.la
lrwxrwxrwx. 1 root root      16 Jan 21 15:49 /usr/local/lib/libproj.so ->
libproj.so.0.8.0
lrwxrwxrwx. 1 root root      16 Jan 21 15:49 /usr/local/lib/libproj.so.0 ->
libproj.so.0.8.0
-rwxr-xr-x. 1 root root 1584463 Jan 21 15:49 /usr/local/lib/libproj.so.0.8.0

that is, was the resolution to #263 (r2629) the right one? The consequence
is, of course, that the original libproj.so.0 is now libproj.so.9, and at
least GDAL and Spatialite complain (2 of a sample of 2). This means that all
downstream use of libproj.so.0 has to be rebuilt, if I understand correctly.

I hope I'm just confused, sorry for only reporting now, I've been travelling.

Roger

_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj

Howard Butler | 4 Mar 18:36 2015
Picon

proj 4.9.1RC5 released

All,

I have issued a RC5 which moved the problematic test into its own, manual, test script, and Charles has made
some CMake tweaks. Unless significant issues are identified, I will motion this for release by the
MetaCRS committee on Friday.

http://download.osgeo.org/proj/proj-4.9.1RC5.tar.gz

Howard

_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj

Luís de Sousa | 26 Feb 10:53 2015
Picon

How to use proj from a Ruby application

Dear all,

I am trying to use Proj from a Ruby application. There are various
Ruby libraries out there [0, 1] typically called "proj4rb"; none seems
to have had active development for multiple years and so far I was not
able to install any of them [2].

What other alternatives are there to use proj with Ruby?

Thank you,

Luís

[0] https://github.com/cfis/proj4rb

[1] https://github.com/Caged/proj4rb

[2] http://gis.stackexchange.com/questions/136550/how-to-install-proj4rb-on-ubuntu-14-04
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Bernhard Jenny | 24 Feb 05:03 2015
Picon

Patterson and Compact Miller projections

All,

Attached is code for the Patterson and the Compact Miller projections for addition to PROJ.4.

The Patterson projection is a cylindrical described here: http://shadedrelief.com/patterson/
Also see Ticket #252: http://trac.osgeo.org/proj/ticket/252 

The Compact Miller is another cylindrical described in this paper: http://cartography.oregonstate.edu/pdf/2015_Jenny_etal_ACompromiseAspect-adaptiveCylindricalProjectionForWorldMaps.pdf

Below are a few test points for both projections.

Thank you,
Bernie Jenny

Assistant Professor
Cartography and Geovisualization
Oregon State University
jennyb <at> geo.oregonstate.edu
http://cartography.oregonstate.edu/

Patterson projection (lon, lat, 0, X, Y):

0 0.0 0 0.0 0.0
0 22.5 0 0.0 2551415.729669344
0 45.0 0 0.0 5366413.421153781
0 67.5 0 0.0 8729502.054111844
0 90.0 0 0.0 1.1409566822831295E7
45 0.0 0 5003778.588046594 0.0
45 22.5 0 5003778.588046594 2551415.729669344
45 45.0 0 5003778.588046594 5366413.421153781
45 67.5 0 5003778.588046594 8729502.054111844
45 90.0 0 5003778.588046594 1.1409566822831295E7
90 0.0 0 1.0007557176093187E7 0.0
90 22.5 0 1.0007557176093187E7 2551415.729669344
90 45.0 0 1.0007557176093187E7 5366413.421153781
90 67.5 0 1.0007557176093187E7 8729502.054111844
90 90.0 0 1.0007557176093187E7 1.1409566822831295E7
135 0.0 0 1.501133576413978E7 0.0
135 22.5 0 1.501133576413978E7 2551415.729669344
135 45.0 0 1.501133576413978E7 5366413.421153781
135 67.5 0 1.501133576413978E7 8729502.054111844
135 90.0 0 1.501133576413978E7 1.1409566822831295E7
180 0.0 0 2.0015114352186374E7 0.0
180 22.5 0 2.0015114352186374E7 2551415.729669344
180 45.0 0 2.0015114352186374E7 5366413.421153781
180 67.5 0 2.0015114352186374E7 8729502.054111844
180 90.0 0 2.0015114352186374E7 1.1409566822831295E7

Compact Miller (lon, lat, 0, X, Y):

0 0.0 0 0.0 0.0
0 22.5 0 0.0 2537439.6610749415
0 45.0 0 0.0 5391682.432264133
0 67.5 0 0.0 8661480.510260897
0 90.0 0 0.0 1.2009484264916677E7
45 0.0 0 5003778.588046594 0.0
45 22.5 0 5003778.588046594 2537439.6610749415
45 45.0 0 5003778.588046594 5391682.432264133
45 67.5 0 5003778.588046594 8661480.510260897
45 90.0 0 5003778.588046594 1.2009484264916677E7
90 0.0 0 1.0007557176093187E7 0.0
90 22.5 0 1.0007557176093187E7 2537439.6610749415
90 45.0 0 1.0007557176093187E7 5391682.432264133
90 67.5 0 1.0007557176093187E7 8661480.510260897
90 90.0 0 1.0007557176093187E7 1.2009484264916677E7
135 0.0 0 1.501133576413978E7 0.0
135 22.5 0 1.501133576413978E7 2537439.6610749415
135 45.0 0 1.501133576413978E7 5391682.432264133
135 67.5 0 1.501133576413978E7 8661480.510260897
135 90.0 0 1.501133576413978E7 1.2009484264916677E7
180 0.0 0 2.0015114352186374E7 0.0
180 22.5 0 2.0015114352186374E7 2537439.6610749415
180 45.0 0 2.0015114352186374E7 5391682.432264133
180 67.5 0 2.0015114352186374E7 8661480.510260897
180 90.0 0 2.0015114352186374E7 1.2009484264916677E7

Attachment (PJ_comill.c): application/octet-stream, 1345 bytes
Attachment (PJ_patterson.c): application/octet-stream, 1443 bytes

_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
ed | 23 Feb 18:51 2015
Picon

Proj 4.9.1RC3 Released

Hi Howard and all,
I have a small suggested addition to the 4.9.1 Release Notes. It would be nice to add a bullet:
o Added the CalCOFI pseudo-projection, #135I'm very happy that this will be re leased soon. Thanks for your help,EdEdward D. Weber, Ph.D. Research Fisheries Biologist NOAA Southwest Fisheries Science Center 8601 La Jolla Shores Drive La Jolla, CA 92037, U.S.A. (858) 546-5676 ed.weber <at> noaa.gov
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj

Gmane