Artur Bercik | 24 Oct 08:26 2014
Picon

(no subject)

Dear Proj list members:

Being a very new in this field, could someone answer me about the following problem:

I have a pair of Geodetic Latitude and Geodetic Longitude (50.89524, 117.69043).
How can I convert it into Latitude and Longitude with WGS84 Ellepsoid?

Thanks in the advance.

Artur
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Charles Clercq | 23 Oct 10:13 2014
Picon

Convertion Lat/Long WGS84 -> UTM-WGS84 problem?

Dear all,

I have written a simple basic code inspired from the doc to convert data in lat/long WGS84 in meters (UTM/WGS84)
My problem is that when I tested it with lat: 51.0242 / long: 3.96338 and obtained  x: 7627600.282316 / y: 742130.588781 which I compared with online converter which gave me : 5652958 / 567564

I would like to understand the error. I that because proj can't convert from WGS84 lat/long to UTM WGS84 ?

I tested with command line:
echo 51.0242 3.96338 | cs2cs.exe +proj=latlong +ellps=WGS84 +to +proj=utm +ellps=WGS84 +datum=WGS84 +units=m
and obtained the exact same result.

Any advice?

Thanks

Clercq

The code
***********************************************************************************************************
#include <iostream>

#include <proj_api.h>

int _tmain(int argc, _TCHAR* argv[])
{
projPJ pj_merc, pj_latlong;
    double x, y;

    if (!(pj_merc = pj_init_plus("+proj=utm +ellps=WGS84 +datum=WGS84 +units=m")) )
        exit(1);
    if (!(pj_latlong = pj_init_plus("+proj=latlong +ellps=WGS84")) )
        exit(1);
std::cout << "Lat: ";
double latitude;
std::cin >> latitude;
std::cout << std::endl;
std::cout << "Long: ";
double longitude;
std::cin >> longitude;

    x = latitude*DEG_TO_RAD;
y = longitude*DEG_TO_RAD;

    int p = pj_transform(pj_latlong, pj_merc, 1, 1, &x, &y, NULL );

//std::cout << "The coordinate in meters are: " << x << "/" << y << std::endl;
printf("Lat/Long: %f/%f - in meters ->: %f/%f\n", latitude, longitude, x, y);

return 0;
}
***********************************************************************************************************
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Mikael Rittri | 22 Oct 11:34 2014

Map projection used for CADRG polar areas?

Hello,

I am reading the specification for CADRG,

http://earth-info.nga.mil/publications/specs/printed/89038/89038_CADRG.pdf .

But I am confused about the polar areas.

 

Section 3.5.5 says: “CADRG shall use a polar stereographic projection“.

 

And Section 3.5.6 says “Polar frames are defined in polar stereographic

projection and non-polar frames are in the ARC projection.

 

But Section 3.6.2 says: “Pixel coordinates in the polar zones are

proportional to rectangular coordinates of the Azimuthal

Equidistant projection, polar aspect, spherical form (as

defined in Map Projections–A Working Manual, page 191).

 

And section 30.1 says: “Polar zone equations are based on

the Azimuthal Equidistant projection, polar aspect,

spherical form.

 

Since Polar Stereographic is not the same thing as Azimuthal Equidistant

in polar aspect, I assume that one of the two alternatives is just a typo

(unless I am missing something). But which one is the typo of  

“Stereographic” and “Azimuthal Equidistant”? Does anyone know?

 

Best regards,

 

Mikael Rittri

Carmenta

Sweden

http://www.carmenta.com

 

 

_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Olivier PENA | 21 Oct 11:34 2014
Picon

pj_set_searchpath support for wince

Hi,

Small patchs are attached in order to have minimal support for opening files under wince using pj_set_searchpath.

I hope I’m in the right place to send this patch.



>>>>-------------------------------------------------------------Ce message a été traité contre les virus par quatre outils différents (Kaspersky, McAfee, Symantec et ThreatSeeker).
This message has been scanned for viruses (by Kaspersky, McAfee, Symantec and ThreatSeeker).
>>>>-------------------------------------------------------------

Attachment (pj_open_lib.c.patch): application/octet-stream, 1098 bytes
Attachment (projects.h.patch): application/octet-stream, 692 bytes
_______________________________________________
Proj mailing list
Proj <at> lists.maptools.org
http://lists.maptools.org/mailman/listinfo/proj
Ed Weber | 29 Sep 18:43 2014
Picon

spamming the wiki

Hi Proj4 list,
There is a problem with the Proj4 wiki. A user named "domtheos" has been
editing pages and adding completely unrelated links (maybe the
account has been hacked?).  I changed all pages back to there correct form
a few days ago but domtheos has done it again. I looked around the Proj4
web site but I don't see a place to report problems.  Could someone advise
about how to report the problem and flag the account (or let the
user know there is a problem)?

Thanks for your help,
Ed Weber

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

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

Gmane