[R-sig-Geo] Define projection using Raster package

[Robert Hijmans]

Tim,

"proj" is not a valid argument in raster(x, ... ) when x (the first
argument) is of class character (i.e. interpreted as a filename); but
I can add that argument. For now, what should work is:

r<-raster(paste(asc.in.dir, "\\", asc.in.files[i], sep=""))
projection(r) <- "+proj=tmerc +lat_0=0 +lon_0=0 +k=9996 +x_0=500000
+y_0=10000000"

Robert

On Fri, May 8, 2009 at 10:26 AM, Tim Sippel <tsippel at gmail.com> wrote:
> I need to take an ESRI grid (.asc) which is in geographic coordinates
> (lat/lon) and project it to UTM coordinates.  The following command is
> giving me an error saying I'm giving it an undefined argument.  I can read
> the raster without the projs arguement, but I'm stuck when adding that on.
>
>> r<-raster(paste(asc.in.dir, "\\", asc.in.files[i], sep=""),
> projs="+proj=tmerc +lat_0=0 +lon_0=0 +k=9996 +x_0=500000 +y_0=10000000")
> Error in .local(x, ...) :
>  unused argument(s) (projs = "+proj=tmerc +lat_0=0 +lon_0=0 +k=9996
> +x_0=500000 +y_0=10000000")
>
> Thanks,
>
> Tim
>
>        [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>