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
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

Gmane