# work on a copy of the original geodetic points map g.copy vector=geodetic_pts,mygeodetic_pts # data preparation: convert z-values from string to double format v.db.addcolumn map=mygeodetic_pts columns="Z_VALUE_D double precision" v.db.update map=mygeodetic_pts column=Z_VALUE_D qcolumn=Z_VALUE # verify: should show identical z-values v.db.select map=mygeodetic_pts columns=cat,Z_VALUE,Z_VALUE_D # convert 2D vector point map to 3D based on attribute v.to.3d input=mygeodetic_pts output=mygeodetic_pts_3d column=Z_VALUE_D # create TIN v.delaunay input=mygeodetic_pts_3d output=mygeodetic_pts_3d_delaunay # rasterize TIN to 500m resolution raster map g.region vector=mygeodetic_pts_3d_delaunay res=500 -p v.tin.to.rast input=mygeodetic_pts_3d_delaunay output=mygeodetic_pts_3d_delaunay r.colors mygeodetic_pts_3d_delaunay color=srtm_plus
Available at: v.tin.to.rast source code (history)
Latest change: Sunday May 01 15:34:14 2022 in commit: ba4fb54ae0451966bc8767ffe8dee4dd049334ba
Main index | Vector index | Topics index | Keywords index | Graphical index | Full index
© 2003-2024 GRASS Development Team, GRASS GIS 8.4.1dev Reference Manual