Frank Warmerdam | 14 Sep 01:33 2014
Picon

PROJ 4.9.0RC1

Folks,

I have discovered that PROJ.4 was left with a hanging beta
again some time ago without my following up on the release. 
Grr, I'm lame!  Anyways, Howard Butler and I have made a 
pass through open bugs and fixed some.  Howard has also
introduced CMake support for PROJ.4 (though the automake
way is also still supported). 

Source is at:


I would appreciate feedback before finalizing the release in the coming week.

Best regards,
Frank


--
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam <at> pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Software Developer
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Frank Warmerdam | 13 Sep 19:40 2014
Picon

EPSG 8.5 Upgrade

Folks,

I have upgraded PROJ.4 in trunk to EPSG 8.5.  We are hoping to prepare a new release candidate today at the FOSS4G 2014 code sprint.

Best regards,  
--
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up   | Frank Warmerdam, warmerdam <at> pobox.com
light and sound - activate the windows | http://pobox.com/~warmerdam
and watch the world go round - Rush    | Geospatial Software Developer
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Yutaka KURASHINA | 11 Sep 08:01 2014
Picon

Problem about egm96 geoid hight calculation

Dear all,

I use proj.4(4.8.0) and egm96 geoid data below.

   http://download.osgeo.org/proj/vdatum/egm96_15/egm96_15.gtx
     Last modified 15-Feb-2014 08:15

I'm trying to calculate egm96 geoid height with cs2cs.
But I could not get the height at boundary of the geoid data.

At first I ran a command below.

$ cs2cs +proj=latlong +geoidgrids=./egm96_15.gtx +to +proj=latlong <<EOF
heredoc> 0 0 0
heredoc> EOF
0dE     0dN 17.162

This works well. But, the next command didn't.

$ cs2cs +proj=latlong +geoidgrids=./egm96_15.gtx +to +proj=latlong <<EOF
heredoc>     0  90 0
heredoc>     0 -90 0
heredoc>   180   0 0
heredoc>  -180   0 0
heredoc> EOF
Rel. 4.8.0, 6 March 2012
<cs2cs>: while processing file: <stdin>, line 1
pj_transform(): point not within available datum shift grids
*       * 0.000
Rel. 4.8.0, 6 March 2012
<cs2cs>: while processing file: <stdin>, line 2
pj_transform(): point not within available datum shift grids
*       * 0.000
Rel. 4.8.0, 6 March 2012
<cs2cs>: while processing file: <stdin>, line 3
pj_transform(): point not within available datum shift grids
*       * 0.000
Rel. 4.8.0, 6 March 2012
<cs2cs>: while processing file: <stdin>, line 4
pj_transform(): point not within available datum shift grids
*       * 0.000

Is there anything wrong in the commands?
I would like to ask for advice.

Thank you in advance.

Yutaka Kurashina

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

Roger Bivand | 30 Aug 14:04 2014
Picon

NAD_1983_To_WGS_1984_5

Hi,

I'm aware of Frank's reply:

http://lists.maptools.org/pipermail/proj/2008-September/003833.html

but would value clarification of why this happens (further details on:

https://stat.ethz.ch/pipermail/r-sig-geo/2014-August/021611.html

from WGS84 geographical coordinates to NAD27 via WGS84_5).

(code attached to get around line breaks)

The first test using cs2cs is wrong because it doesn't use +towgs84= to 
use *_5:

(the test is against benchmark coordinates)

The works in the first sub-step (geographical to WGS84_5 State Plane):

So locs.step1 are OK, but the next step gets the same errors as ignoring 
*_5

Turning on -v shows:

cs2cs -v -f "%16.9f" +proj=lcc +lat_1=38.43333333333333 
+lat_2=37.06666666666667 +lat_0=36.5 +lon_0=-120.5 +x_0=2000000.0 
+y_0=500000.0 +ellps=GRS80 +datum=NAD83 +units=us-ft +no_defs 
+towgs84=-0.991,1.9072,0.5129,0.025789908,0.0096501,0.0116599,0.0 +to 
+init=epsg:26743 < locs.step1

# ---- From Coordinate System ----
#Lambert Conformal Conic
#	Conic, Sph&Ell
#	lat_1= and lat_2= or lat_0
# +proj=lcc +lat_1=38.43333333333333 +lat_2=37.06666666666667 +lat_0=36.5
# +lon_0=-120.5 +x_0=2000000.0 +y_0=500000.0 +ellps=GRS80 +datum=NAD83
# +units=us-ft +no_defs
# +towgs84=-0.991,1.9072,0.5129,0.025789908,0.0096501,0.0116599,0.0
#--- following specified but NOT used
# +ellps=GRS80 +towgs84=0,0,0
# ---- To Coordinate System ----
#Lambert Conformal Conic
#	Conic, Sph&Ell
#	lat_1= and lat_2= or lat_0
# +init=epsg:26743 +proj=lcc +lat_1=38.43333333333333
# +lat_2=37.06666666666667 +lat_0=36.5 +lon_0=-120.5 
+x_0=609601.2192024384
# +y_0=0 +datum=NAD27 +units=us-ft +no_defs +ellps=clrk66
# +nadgrids= <at> conus, <at> alaska, <at> ntv2_0.gsb, <at> ntv1_can.dat

Puzzled,

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
cat > locs.xy << EOF
-121.524764291826 37.6600366036405
-121.523480804667 37.6543604613483
EOF

cs2cs -f "%16.9f" +proj=longlat +datum=WGS84 +to +init=epsg:26743 locs.xy | awk '{if (NR == 1)
 printf("%16.9f %16.9f\n", 1703671.30566227-$1, 424014.398045834-$2)
 else
 printf("%16.9f %16.9f\n", 1704020.20113366-$1, 421943.708664294-$2)
}'

cs2cs -f "%16.9f" +proj=longlat +datum=WGS84 +to +proj=lcc +lat_1=38.43333333333333
+lat_2=37.06666666666667 +lat_0=36.5 +lon_0=-120.5 +x_0=2000000.0 +y_0=500000.0 +ellps=GRS80
+datum=NAD83 +units=us-ft +no_defs
+towgs84=-0.991,1.9072,0.5129,0.025789908,0.0096501,0.0116599,0.0 locs.xy > locs.step1
awk '{
if (NR == 1)
 printf("%16.9f %16.9f\n", 6265039.1378244-$1, 2064418.92932968-$2)
 else
 printf("%16.9f %16.9f\n", 6265388.04257557-$1, 2062348.22239488-$2)
}' locs.step1

cs2cs -f "%16.9f" +proj=lcc +lat_1=38.43333333333333 +lat_2=37.06666666666667 +lat_0=36.5
+lon_0=-120.5 +x_0=2000000.0 +y_0=500000.0 +ellps=GRS80 +datum=NAD83 +units=us-ft +no_defs
+towgs84=-0.991,1.9072,0.5129,0.025789908,0.0096501,0.0116599,0.0 +to +init=epsg:26743 <
locs.step1 | awk '{
if (NR == 1)
 printf("%16.9f %16.9f\n", 1703671.30566227-$1, 424014.398045834-$2)
 else
 printf("%16.9f %16.9f\n", 1704020.20113366-$1, 421943.708664294-$2)
}'

_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Paul Boettinger | 30 Aug 10:26 2014
Picon
Picon

Datum shift not working

Hello,
i want to do a transformation with the EPSG Code 3440
In my epsg file 3440 references the following parameters:
+proj=utm +zone=40 +ellps=clrk80
+towgs84=-180.624,-225.516,173.919,-0.81,-1.898,8.336,16.7101 +units=m +no_defs

Executing proj 4.8.0_0 from the MacPorts Repository on OSX 10.9.4 makes no differences using the
transformation with or without the datum shift.
Any help welcome,

Paul

proj +proj=utm +zone=40 +ellps=clrk80
+towgs84=-180.624,-225.516,173.919,-0.81,-1.898,8.336,16.7101 +units=m +no_defs -f "%.6f"
59.28632619 24.27334106
732098.167370	2686171.824604

proj +proj=utm +zone=40 +ellps=clrk80 +units=m +no_defs -f "%.6f"
59.28632619 24.27334106
732098.167370	2686171.824604
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj

Cleo Drakos | 23 Aug 08:38 2014
Picon

Convert the array into required coordinates

Dear guys,

I have 2d numpy array of 480 rows and 1440 columns as named by 'data' below:

The first element belongs to (49.875S,179.875W), the second element belongs to (49.625S,179.625W), and the last element belongs to (49.875N,179.875E). import os, glob, gdal, numpy as np fname = '3B42RT.2014010606.7.bin' with open(fname, 'rb') as fi: fi.seek(2880,0) data = np.fromfile(fi,dtype=np.uint16,count=480*1440) data = data.byteswap() data = data.reshape(1440,480)

How can I convert this numpy array so that its first element belongs to (49.875N,179.625W), i.e., upper left latitude and longitude respectively; and the last element belong to (49.625S,179.875E), i.e., lower right latitute and longitude respectively.

I tried to rotate it, but I do not think it is correct.

data = np.rot90(data,1)

Have some of you experienced with this type of problem? The binary file I am using is here:ftp://trmmopen.gsfc.nasa.gov/pub/merged/3B42RT/3B42RT.2014010606.7.bin.gz

cleo
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Aleksander Yanovskiy | 21 Aug 09:35 2014
Picon

direct geodesic problem

Greetings!

Maybe somebody can clarify how to solve the following direct geodesic 
problem: find the end point of a geodesic given its starting point and 
initial azimuth and the angle between the tangent planes at the 
beginning and the end points (or, equivalently, the angel between the 
normals at the points).
It seems the most easy way is to find the corresponding arc length on 
the auxiliary sphere, but I haven't found anywhere the explicit formula 
for it as the function of the point coordinates, the azimuth and the 
angle between the tangent planes at the beginning and the end points. 
And I'm not sure that such formula would be valid for the the beginning 
and the end points laying in different hemispheres.

And why is the geodesic problem in such formulation not being used in 
practice ?

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

Chris HUDDART | 18 Aug 19:04 2014

Does anyone know of any Ebola mapping activities?

Dear All,

 

Does anyone know of any West Africa Ebola mapping activities?

 

Many thanks & regards,

 

 

 

Chris.

The information contained in this electronic message and any attachments is intended for specific individuals or entities, and may be confidential, proprietary or privileged. If you are not the intended recipient, please notify the sender immediately, delete this message and do not disclose, distribute or copy it to any third party or otherwise use this message. The content of this message does not necessarily reflect the official position of the World Food Programme. Electronic messages are not secure or error free and may contain viruses or may be delayed, and the sender is not liable for any of these occurrences. The sender reserves the right to monitor, record and retain electronic messages.
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Cleo Drakos | 18 Aug 16:13 2014
Picon

Can I convert the ellipsoids using Proj?

Dear guys,

I have the data (latitudes, longitudes, and elevations) given in ICESat/GLAS ellipsoid. This equatorial radius and  reciprocal flattening (1/f) on this ellipsoid are 6378136.30 m and is 298.257 respectively.

Now, I want to convert my data into WGS84 ellepsoid with equatorial radius and reciprocal flattening (1/f) of 6378137.000000 and 298.25722356 respectively.

Can I convert the ellipsoids using Proj?

I hope some of you can help me.

cleo
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Roger Bivand | 18 Aug 11:32 2014
Picon

#229 - trunk init still broken

Please could someone patch pj_init.c with the bugfix for the missing 
found_def after line 263? The fix has been known for months, but no action 
has been taken.

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

Charles Karney | 10 Aug 23:49 2014

rhumb lines + arbitrary precision arithmetic

I've released version 1.37 of GeographicLib.  See

   http://geographiclib.sf.net/html/changes.html

A couple of items might be of interest to the readers of this list.

(1) Classes and a utility for computing rhumb lines.  For details, see

   http://geographiclib.sf.net/html/rhumb.html

(2) Support for high precision arithmetic.  A compilation option is
provided to allow linking against Boost (for quad precision) or MPFR
(for arbitrary precision).  Thanks to C++'s operator overloading this
entailed minimal changes to the code itself.  For details, see

   http://geographiclib.sf.net/html/highprec.html

--

-- 
Charles Karney <charles.karney <at> sri.com>
SRI International, Princeton, NJ 08543-5300

Tel: +1 609 734 2312
Fax: +1 609 734 2662
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj


Gmane