[GRASSLIST:9072] Create new location with shellscript

GRASS GRASS user
1

dear list,

i want to create a new temporary location (and mapset) with a shellscript in batch mode.

i know you can do it with r.in.gdal if you have already one existing location and a raster file.
but i have no location and only points in ascii format and the projection can change.

lg
michaela

--
Institut fuer Geographie und Regionalforschung
Universitaet Wien
Kartografie und Geoinformation

Departement of Geography and Regional Research
University of Vienna
Cartography and GIS

Universitaetstr. 7, A-1010 Wien, AUSTRIA

Tel: (+43 1) 4277 48646
Fax: (+43 1) 4277 48649
E-mail: kinb@atlas.gis.univie.ac.at

WWW: http://www.gis.univie.ac.at/karto

2

the only thing, you have to do, is to create directory within grass
database, with one subdirectory -- PERMANENT with all necessary files
(DEFAULT_WIND, PROJ_INFO, ...) files in it

jachym

P.S. Ich werde am Wochenende im Wien sein -- vielleicht könnte ich dir
dort helfen

On Tue, Nov 15, 2005 at 09:46:38AM +0100, Michaela Kinberger wrote:

dear list,

i want to create a new temporary location (and mapset) with a
shellscript in batch mode.

i know you can do it with r.in.gdal if you have already one existing
location and a raster file.
but i have no location and only points in ascii format and the
projection can change.

lg
michaela

--
Institut fuer Geographie und Regionalforschung
Universitaet Wien
Kartografie und Geoinformation

Departement of Geography and Regional Research
University of Vienna
Cartography and GIS

Universitaetstr. 7, A-1010 Wien, AUSTRIA

Tel: (+43 1) 4277 48646
Fax: (+43 1) 4277 48649
E-mail: kinb@atlas.gis.univie.ac.at

WWW: http://www.gis.univie.ac.at/karto

--
Jachym Cepicky
e-mail: jachym.cepicky@centrum.cz
URL: http://les-ejk.cz
GPG: http://les-ejk.cz/gnupg_public_key/
-----------------------------------------
OFFICE:
Department of Geoinformation Technologies
LDF MZLU v Brnì
Zemìdìlská 3
613 00 Brno
e-mail: xcepicky@node.mendelu.cz
URL: http://mapserver.mendelu.cz
Tel.: +420 545 134 514

3

is there a documentation about the parameters necessary for the fiffrent projections known in GRASS.

in /usr/local/grass6.0.1/etc i found the file projections but in PROJ_INFO the is an number for the projection.

lg
michaela

Jachym Cepicky wrote:

the only thing, you have to do, is to create directory within grass
database, with one subdirectory -- PERMANENT with all necessary files
(DEFAULT_WIND, PROJ_INFO, ...) files in it

jachym

P.S. Ich werde am Wochenende im Wien sein -- vielleicht könnte ich dir
dort helfen

On Tue, Nov 15, 2005 at 09:46:38AM +0100, Michaela Kinberger wrote:

dear list,

i want to create a new temporary location (and mapset) with a shellscript in batch mode.

i know you can do it with r.in.gdal if you have already one existing location and a raster file.
but i have no location and only points in ascii format and the projection can change.

lg
michaela

--
Institut fuer Geographie und Regionalforschung
Universitaet Wien
Kartografie und Geoinformation

Departement of Geography and Regional Research
University of Vienna
Cartography and GIS

Universitaetstr. 7, A-1010 Wien, AUSTRIA

Tel: (+43 1) 4277 48646
Fax: (+43 1) 4277 48649
E-mail: kinb@atlas.gis.univie.ac.at

WWW: http://www.gis.univie.ac.at/karto

--
Institut fuer Geographie und Regionalforschung
Universitaet Wien
Kartografie und Geoinformation

Departement of Geography and Regional Research
University of Vienna
Cartography and GIS

Universitaetstr. 7, A-1010 Wien, AUSTRIA

Tel: (+43 1) 4277 48646
Fax: (+43 1) 4277 48649
E-mail: kinb@atlas.gis.univie.ac.at

WWW: http://www.gis.univie.ac.at/karto

4

is there a documentation about the parameters necessary for the
diffrent projections known in GRASS.

Generally it follows PROJ4 parameters:
  http://www.remotesensing.org/proj/
  http://www.remotesensing.org/proj/gen_parms.html

in /usr/local/grass6.0.1/etc i found the file projections but in
PROJ_INFO the is an number for the projection.

That is class of projection, the number is not of much practical
interest to the user, but here they are anyway:

include/gis.h

#define PROJECTION_XY 0
#define PROJECTION_UTM 1
#define PROJECTION_SP 2
#define PROJECTION_LL 3
#define PROJECTION_OTHER 99

Usually you just need to run g.setproj and type "list" if you don't
know what to choose. see also g.proj.

Hamish

5

Hello all

On Thu, 17 Nov 2005, Hamish wrote:

is there a documentation about the parameters necessary for the
diffrent projections known in GRASS.

Generally it follows PROJ4 parameters:
http://www.remotesensing.org/proj/
http://www.remotesensing.org/proj/gen_parms.html

Something Frank said on the PROJ list recently kind of hinted that there might be some new programmatic way of finding out about the parameters needed for different projections:
http://lists.maptools.org/pipermail/proj/2005-November/001876.html
But I haven't had time to investigate it.

Usually you just need to run g.setproj and type "list" if you don't
know what to choose. see also g.proj.

Yes g.setproj is interactive only so can't be used from a shell script. g.proj has become quite versatile in creating new locations from various sources of projection information; have you looked through the examples in the g.proj man page?

Paul