Hi,

As suggested, I am transforming my data from Sapper Hill 1943 /
UTM_Zone_20S to WGS84  /UTM_Zone_20S to run variograms of my plots that
have the x axis in meters.  I get the following error message and I need
some hint of what is going on. The SpatialPointsDataFrame has 6816
points. Thanks.

> tuc1pro3.50.utm <- spTransform(tuc1pro3.50frm, CRS("+proj=utm +zone=20
+south +ellps=WGS84 +datum=WGS84 +units=m"))

non finite transformation detected:
        x         y
-65.29769 -26.39063       Inf       Inf
Error in spTransform(xSP, CRSobj) : failure in points 1
In addition: Warning message:
In spTransform(xSP, CRSobj) : 6816 projected point(s) not finite



2011/11/16 Ben Madin <lists@remoteinformation.com.au>

> You will need to transform (re-project) your data into a suitable
> coordinate system - use the sptransform() function in the rgdal package, or
> look at p 64 in the ASDAR book.
>
> depending on where you data is and how large an area it covers will
> influence your choice of coordinate reference system, but you might get
> away with a UTM grid 19 North for Puerto Rico - have a look here :
>
> http://spatialreference.org/ref/epsg/32619/
>
> cheers
>
> Ben
>
>
> On 16/11/2011, at 11:42 PM, carolina monmany wrote:
>
> > Sorry. My problem is that I want the units in the x axis to be in meters
> > but they are in degrees. How can I make R read the projected coordinate
> > system, which is in meters, and not the geographic one?
> >
> > Thanks!
> >
> > 2011/11/16 Paul Hiemstra <paul.hiemstra@knmi.nl>
> >
> >> Hi,
> >>
> >> This is mail is a bit on the short side :), so please try and be more
> >> complete next time. I think what you want to know is what the units are
> >> in a semi-variogram. For the y-axis, they are the unit of the variable
> >> of interest squared. So if the variable is m/s, than the semi-variogram
> >> has units (m/s)^2. For the x-axis they are most probably the distance
> >> unit in which the coordinate system of your data set is, e.g. kilometer.
> >>
> >> hope this helps,
> >> Paul
> >>
> >> On 11/16/2011 03:08 PM, carolina monmany wrote:
> >>
> >> --
> >> Paul Hiemstra, Ph.D.
> >> Global Climate Division
> >> Royal Netherlands Meteorological Institute (KNMI)
> >> Wilhelminalaan 10 | 3732 GK | De Bilt | Kamer B 3.39
> >> P.O. Box 201 | 3730 AE | De Bilt
> >> tel: +31 30 2206 494
> >>
> >> http://intamap.geo.uu.nl/~paul
> >> http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770
> >>
> >>
> >
> >
> > --
> > ---------------------------
> > A. CAROLINA MONMANY
> > Universidad de Puerto Rico
> > Departamento de Biologia - CN 235
> > POBOX 23360
> > San Juan, Puerto Rico 00931-3360
> > Tel: +1 787 764 0000 x2847
> > Fax: +1 787 764 2610
> >
> >       [[alternative HTML version deleted]]
> >
> > _______________________________________________
> > R-sig-Geo mailing list
> > R-sig-Geo@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>


-- 
---------------------------
A. CAROLINA MONMANY
Universidad de Puerto Rico
Departamento de Biologia - CN 235
POBOX 23360
San Juan, Puerto Rico 00931-3360
Tel: +1 787 764 0000 x2847
Fax: +1 787 764 2610

	[[alternative HTML version deleted]]

