[GRASSLIST:3062] (no subject)

Tue Mar 30 13:41:53 EST 2004

>From liboe at linux  Sun Jun 11 19:48:24 2000
Return-Path: <owner-GRASSLIST at baylor.edu>
Received: from [129.62.4.29] (HELO ccis01.baylor.edu)
  by frontend01.cg.ifxnetworks.com (CommuniGate Pro SMTP 4.1.8)
  with ESMTP id 31892749 for gfmb at tutopia.com; Mon, 08 Mar 2004 06:07:46 -0500
Received: from localhost (ccis01.baylor.edu [129.62.4.29])
	by ccis01.baylor.edu (8.9.3p2/8.9.1) with SMTP id EAA1404298;
	Mon, 8 Mar 2004 04:58:00 -0600 (CST)
Received: from waylon.baylor.edu (waylon.baylor.edu [129.62.3.228])
	by ccis01.baylor.edu (8.9.3p2/8.9.1) with ESMTP id EAA1368474
	for <GRASSLIST at listproc.baylor.edu>; Mon, 8 Mar 2004 04:57:04 -0600 (CST)
Received: from willie.baylor.edu ([129.62.8.12]) by waylon.baylor.edu with Microsoft SMTPSVC(6.0.3790.0);
	 Mon, 8 Mar 2004 04:57:04 -0600
Received: from camelot.itc.it ([195.223.171.5])
 by willie.baylor.edu (SAVSMTP 3.1.1.32) with SMTP id M2004030804570104960
 for <grasslist at baylor.edu>; Mon, 08 Mar 2004 04:57:03 -0600
Received: from ntvir.itc.it (ntvir [10.0.20.5])
	by camelot.itc.it (8.12.8/8.12.8) with SMTP id i28Av0Pg008664;
	Mon, 8 Mar 2004 11:57:00 +0100 (MET)
Received: from orchestra.itc.it ([10.0.10.11])
 by ntvir.itc.it (SAVSMTP 3.1.2.35) with SMTP id M2004030811570026415
 ; Mon, 08 Mar 2004 11:57:00 +0100
Received: from thuille.itc.it. (thuille [10.40.0.110])
	by orchestra.itc.it (8.11.6/8.11.6) with ESMTP id i28AuvW21341;
	Mon, 8 Mar 2004 11:56:57 +0100
Received: from thuille.itc.it. (localhost.localdomain [127.0.0.1])
	by thuille.itc.it. (8.12.8/8.12.8) with ESMTP id i28AulIs022261;
	Mon, 8 Mar 2004 11:56:57 +0100
Received: (from neteler at localhost)
	by thuille.itc.it. (8.12.8/8.12.8/Submit) id i28AuV45022259;
	Mon, 8 Mar 2004 11:56:31 +0100
Date: Mon, 8 Mar 2004 11:56:31 +0100
From: Markus Neteler <neteler at itc.it>
To: Paul Kelly <paul-grass at stjohnspoint.co.uk>
Cc: Morten Hulden <morten at untamo.net>,
 werchowyna at pf.pl,
 grasslist at baylor.edu
Subject: [GRASSLIST:2871] Re: Roussilhe Oblique Stereographic
Message-ID: <20040308105631.GG20834 at thuille.itc.it>
Mail-Followup-To: Paul Kelly <paul-grass at stjohnspoint.co.uk>,
	Morten Hulden <morten at untamo.net>, werchowyna at pf.pl,
	grasslist at baylor.edu
References: <404B3FE8.E7EC46A2 at stjohnspoint.co.uk>
Mime-Version: 1.0
Content-Type: text/plain;
  charset=us-ascii
Content-Disposition: inline
In-Reply-To: <404B3FE8.E7EC46A2 at stjohnspoint.co.uk>
User-Agent: Mutt/1.4.1i
X-OriginalArrivalTime: 08 Mar 2004 10:57:04.0658 (UTC) FILETIME=[178CA720:01C404FC]
Sender: owner-GRASSLIST at baylor.edu
X-Listprocessor-Version: 8.2.10/991025/16:55 -- ListProc(tm) by CREN

On Sun, Mar 07, 2004 at 03:29:44PM +0000, Paul Kelly wrote:
> Hello
> 
> Morten Hulden wrote:
> > 
> > On Sat, 6 Mar 2004 werchowyna at pf.pl wrote:
> > 
> > > I need to make a mapset in a projection called Roussilhe Projection (alias
> > > Quasi-Stereographic) on the Krassovsky ellipsoid, modified Pulkovo 1942
> > > datum different than those included in libproj but the datum os not an issue
> > > here, I solved it myself
> > >
> > > the complete system (ellipsoid, datum and projection) is an old, the most
> > > popular and widely used system in Poland (called "System 1965" ("Uklad 1965"
> > > in Polish))
> > >
> > > the issue is the projection only - there is no projection called like it in
> > > Grass
> > >
> > > notice: the projection name is my literal translation from Polish; I'm not
> > > sure what is the proper name in English though I have met some sound similar
> > >
> > > any ideas how to apply this projection in Grass?
> > 
> > I have:
> > 
> > PROJ_INFO
> > 
> > name: Transverse Mercator
> > datum: S-42
> > towgs84: 33.4,-146.6,-76.3,-0.359,-0.053,0.844,-0.84
> > proj: tmerc
> 
> tmerc is only used in zone V in Poland; for the other 4 zones it should
> be
> proj: stere
> Plenty of information at
> http://crs.ifag.de/country_select.php?country=PL
> or ask on the PROJ.4 mailing list http://remotesensing.org/proj/
> 
> If you have GRASS 5.7 it is much easier to set up projections: try 
> g.proj -c proj4='+init=epsg:2171' location=poland1
> g.proj -c proj4='+init=epsg:2172' location=poland2
> 
> etc. to create new locations with these projections. Should work if you
> have PROJ.4 and GDAL installed properly.
> 
> Paul

Based in this nice feature (thanks, Paul!) I have made a new
script:

Script to create a new LOCATION from EPSG code
Usage:
   make_location_epsg_g57.sh epsg newlocation_name [GISDBASE]

       epsg: EPSG code number of projection (see /usr/local/share/proj/epsg)
       newlocation_name: new location to be created
       GISDBASE: full path to GRASS database directory (optional)
                 e.g. /ssi0/ssi/neteler/grassdata

You have to run it from *outside* GRASS.

Example:

  make_location_epsg_g57.sh  2171 poland1

does the same outside GRASS 5.7 as above 'g.proj' command line
from inside a GRASS location.

Fetch it here:

 http://mpa.itc.it/grasstutor/scripts/make_location_epsg_g57.sh

Cheers

 Markus

-- 
Markus Neteler     <neteler itc it>       http://mpa.itc.it
ITC-irst -  Centro per la Ricerca Scientifica e Tecnologica
MPBA - Predictive Models for Biol. & Environ. Data Analysis
Via Sommarive, 18        -       38050 Povo (Trento), Italy