moving raster maps across UTM zone boundaries
Bruce K. Wylie (303) 490 8336
wylie at gpsrv1.gpsr.colostate.edu
Thu Apr 15 15:13:24 EDT 1993
Look at the SCS contrib commands m.setproj and v.proj to move a
vector file from one projection (MapsetX with UTM zone x) to
another projections (MapsetY with UTM zone z). Note that MapsetX
and MapsetZ must lie under the same Location.
P.S. Note: I have not done this but this is how I would approach
> Hello fellow GRASS users.
> I've been fortunate in the past to have data sets well within UTM zones.
> Unfortunately, that happy period is at an end. I am dealing with county soils
> and land use data for Oklahoma which is split between zones 14 and 15. What I
> want to do is bring the raster map for a county in zone 15, and patch it with
> a raster map in zone 14. I've tried this procedure:
> 1) used r.poly to convert the zone 15 map to a vector map.
> 2) used v.out.ascii to convert the grass vector map to ascii.
> 3) used m.u2ll to convert the ascii vector UTM map to lat.-lon. coordinates
> 4) used m.ll2u to convert lat-lon. ascii vector map to zone 14 utm coordinates
> 5) used v.in.ascii to convert new zone14 ascii vector map to grass vector map.
> Now, I'm stuck. The manual hints that I should apply v.digit to re-assign
> attributes to the vector map before using v.to.rast. BUT, I'm dealing with
> _thousands_ of cells. Is there a better way to do this whole process? Or
> does anyone have any suggestions? I suspect this question may be posted
> frequently, but I only recently began monitoring this discussion. Thanks in
> advance for any and all replies.
> Fred L. Ogden
> Iowa Institute of Hydraulic Research
> The University of Iowa
> Iowa City, IA 52242
> flogden at hml.uiowa.edu
More information about the grass-user