[GRASS-user] partial vector copy to smaller territory

Jón Eiríksson jeir at hi.is
Wed Oct 31 04:35:48 PDT 2018

The following partial script may be helpful (having the source map in shape format):

assume folder in home directory containing source map - in this case contours.shp
v_w,v_s,v_e,v_n are string variables with destination region
the commands create a subset defined by the destination region: contours.shp
note that the identical naming of input and output is perhaps confusing but harmless
as long as the home directory does not contain files named contours.*
importing required modules

import os
import sys
import subprocess
from array import array
import string
import math
from grass.script import core as grass
from grass.script import db as grass


cmd = 'ogr2ogr -clipsrc' + ' ' + v_w + ' ' + v_s + ' ' + v_e + ' ' + v_n
cmd = cmd + ' ' + '~/contours.shp' + ' ~/folder_name/contours.shp'
subprocess.call(cmd, shell = True)
grass.run_command('v.in.ogr', input='contours.shp', layer='contours', output='contours')


On 31 Oct 2018, at 10:12, Frank David wrote:


What is the best and fast way to make a copy of a part of vector map (national) to get a smaller territory ? (partial copy).
Thank you for your advice !

grass-user mailing list
grass-user at lists.osgeo.org
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20181031/6d9a2fa1/attachment.html>

More information about the grass-user mailing list