Re: [GRASS-user] problem creating secant Lambert Conformal Conic projection

[Patton, Eric (NRCan/RNCan)]

Markus,

Responding to a months-old post (sorry!)

Upgrading to proj 4.9.2 and gdal 2.1.2 seemed to fix the problem.

Both secant parallels in my Lambert Conformal projection are now being reported:

> g.proj -we
PROJCS["Lambert_Conformal_Conic",
GEOGCS["GCS_WGS_1984",
DATUM["D_WGS_1984",
SPHEROID["WGS_1984",6378137,298.257223563]],
PRIMEM["Greenwich",0],
UNIT["Degree",0.017453292519943295]],
PROJECTION["Lambert_Conformal_Conic"],
PARAMETER["standard_parallel_1",65],
PARAMETER["standard_parallel_2",75],
PARAMETER["latitude_of_origin",65],
PARAMETER["central_meridian",-60],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["Meter",1]]

Sorry for the noise!

--
Eric 


From: neteler.os...@gmail.com <neteler.os...@gmail.com> on behalf of Markus 
Neteler <nete...@osgeo.org>
Sent: Tuesday, May 17, 2016 5:33 PM
To: Patton, Eric (NRCan/RNCan)
Cc: grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] problem creating secant Lambert Conformal Conic 
projection

On Tue, May 17, 2016 at 7:28 PM, Patton, Eric (NRCan/RNCan)
<eric.pat...@canada.ca> wrote:
> Hi Markus,
>
> Sorry for the delay in responding. I am using Grass 7.0.5svn, proj 4.8.0,
> gdal 1.10.1.
>
> I really don’t understand what could be the cause of this. I tried a fresh
> checkout and install, and encountered the same behavior. I'm using the proj
> and gdal packages provided by my package manager, so there's no custom
> installations installed in parallel.

proj-4.8.0 is from roughly 13-Mar-2012

I use proj 4.9.1, 04 March 2015

Perhaps some fixes in between?

Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] problem creating secant Lambert Conformal Conic projection

[Markus Neteler]

On Tue, May 17, 2016 at 7:28 PM, Patton, Eric (NRCan/RNCan)
 wrote:
> Hi Markus,
>
> Sorry for the delay in responding. I am using Grass 7.0.5svn, proj 4.8.0,
> gdal 1.10.1.
>
> I really don’t understand what could be the cause of this. I tried a fresh
> checkout and install, and encountered the same behavior. I'm using the proj
> and gdal packages provided by my package manager, so there's no custom
> installations installed in parallel.

proj-4.8.0 is from roughly 13-Mar-2012

I use proj 4.9.1, 04 March 2015

Perhaps some fixes in between?

Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] problem creating secant Lambert Conformal Conic projection

[Patton, Eric (NRCan/RNCan)]

??Hi Markus,

Sorry for the delay in responding. I am using Grass 7.0.5svn, proj 4.8.0, gdal 
1.10.1.

I really don't understand what could be the cause of this. I tried a fresh 
checkout and install, and encountered the same behavior. I'm using the proj and 
gdal packages provided by my package manager, so there's no custom 
installations installed in parallel.

~ Eric.


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] problem creating secant Lambert Conformal Conic projection

[Moritz Lennert]

On 14/05/16 16:14, Markus Neteler wrote:

Hi Eric,

On Fri, May 13, 2016 at 6:16 PM, Eric Patton  wrote:

Hi guys,

I am trying to make a new Grass location with the secant Lambert Conformal
Conic projection.


To be sure, that's this one, right?
http://www.remotesensing.org/geotiff/proj_list/lambert_conic_conformal_2sp.html


I have noticed that if I am selecting the projection parameters through the
Grass GUI dialog, and lat_0 = lat_1, the value of lat_2 doesn't seem to be
included in the projection definition (see attached screenshot
lat_2_is_gone.png).


Strange. I have replicated it here and it looks fine:

GRASS 7.0.3 (lcc):~ > g.proj -w
PROJCS["Lambert Conformal Conic",
 GEOGCS["wgs84",
 DATUM["WGS_1984",
 SPHEROID["WGS_1984",6378137,298.257223563],
 TOWGS84[0.000,0.000,0.000]],
 PRIMEM["Greenwich",0],
 UNIT["degree",0.0174532925199433]],
 PROJECTION["Lambert_Conformal_Conic_2SP"],
 PARAMETER["standard_parallel_1",65],
 PARAMETER["standard_parallel_2",75],
 PARAMETER["latitude_of_origin",60],
 PARAMETER["central_meridian",60],
 PARAMETER["false_easting",0],
 PARAMETER["false_northing",0],
 UNIT["Meter",1]]




In your case let_0 != lat_1, so this corresponds to the case where it 
works also for Eric. It's the case when lat_0 = lat_1 that doesn't work 
for him.


Testing here with grass64, grass70, grass71 it works, though (see 
attached screenshot). No idea what is going wrong for you, there, Eric.


Moritz
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] problem creating secant Lambert Conformal Conic projection

[Markus Neteler]

Hi Eric,

On Fri, May 13, 2016 at 6:16 PM, Eric Patton  wrote:
> Hi guys,
>
> I am trying to make a new Grass location with the secant Lambert Conformal
> Conic projection.

To be sure, that's this one, right?
http://www.remotesensing.org/geotiff/proj_list/lambert_conic_conformal_2sp.html

> I have noticed that if I am selecting the projection parameters through the
> Grass GUI dialog, and lat_0 = lat_1, the value of lat_2 doesn't seem to be
> included in the projection definition (see attached screenshot
> lat_2_is_gone.png).

Strange. I have replicated it here and it looks fine:

GRASS 7.0.3 (lcc):~ > g.proj -w
PROJCS["Lambert Conformal Conic",
GEOGCS["wgs84",
DATUM["WGS_1984",
SPHEROID["WGS_1984",6378137,298.257223563],
TOWGS84[0.000,0.000,0.000]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",65],
PARAMETER["standard_parallel_2",75],
PARAMETER["latitude_of_origin",60],
PARAMETER["central_meridian",60],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["Meter",1]]


> If however, lat_0 != lat_1, lat_2 *does* get included in the projection
> parameters (see attached lat_2_exists.png).
>
> Is this an error? Or is it incorrect usage on my part to have lat_0 and
> lat_1 equal?

Which GRASS GIS version do you use?

Markus
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user