Roger,
projInfo("datum") works fine and outputs 10 datum,
ED50 and ED87 not listed.
Anyway, looking for an operational solution, is this correct for Norway?
CRS("+proj=utm +zone=35 +ellps=intl +units=m +no_defs
+towgs84=-82.981,-99.719,-110.709,-0.5076,0.1503,0.3898")
(I'm assuming that the Kol
On Mon, 11 Oct 2010, Agustin Lobo wrote:
Casper,
Thanks, this is very useful, although does not solve my problem.
The list of datums known to proj4 is given by projInfo("datum") - untried,
please try yourself. If they are not given, it is because they are
ambiguous, so proj4 sensibly errs o
I'm really not an expert, but I know that in GRASS exist four different sets of +towgs84 parameters
for the potsdam datum (due to regional differences that seem to go back to inaccuracies in the
original land survey by Karl Friedrich Gauß in Germany). If one of these is specified, it is not
nece
Thanks,
I see the problem: have 11 datums and ED50 or ED87 are not
among them. Is there any way to add datums? or is this something that
comes with proj4 and it is just as it is?
Agus
2010/10/11 Tom Gottfried :
> Hi Agustin,
>
> you can use `cs2cs -ld' directly in your shell to list all datums
Hi Agustin,
you can use `cs2cs -ld' directly in your shell to list all datums known to
proj.4
Tom
--
Technische Universität München
Department für Pflanzenwissenschaften
Lehrstuhl für Grünlandlehre
Alte Akademie 12
85350 Freising / Germany
Phone: ++49 (0)8161 715324
Fax: ++49 (0)8161 713243
e
Casper,
Thanks, this is very useful, although does not solve my problem.
What I would need is to define:
proj4string(lakes) <- CRS("+proj=utm +zone=35 +datum=ED87
+towgs84=-82.981,-99.719,-110.709,-0.5076,0.1503,0.3898")
according to
http://www.asprs.org/resources/grids/10-99-norway.pdf
My ques
HI Agus,
see
A=make_EPSG()
A
in rgdal
Caspar
On Sun, Oct 10, 2010 at 10:36 PM, Agustin Lobo wrote:
> HI!
>
> While
>> CRS("+proj=utm +zone=32 +datum=WGS84")
> CRS arguments:
> +proj=utm +zone=32 +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
>
> is correct, the same with ED50 is not:
>
>> CRS("
HI!
While
> CRS("+proj=utm +zone=32 +datum=WGS84")
CRS arguments:
+proj=utm +zone=32 +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
is correct, the same with ED50 is not:
> CRS("+proj=utm +zone=32 +datum=ED50")
Error in CRS("+proj=utm +zone=32 +datum=ED50") :
unknown elliptical parameter name
Is t
