<div dir="ltr"><br><br>On Wed, Nov 7, 2018 at 11:18 AM Frank David <<a href="mailto:frank.david@geophom.fr">frank.david@geophom.fr</a>> wrote:<br>><br>> Le 06/11/2018 à 21:36, Markus Metz a écrit :<br>><br>><br>><br>> On Tue, Nov 6, 2018 at 4:03 PM Frank David <<a href="mailto:frank.david@geophom.fr">frank.david@geophom.fr</a>> wrote:<br>> ><br>> ><br>> > Le 05/11/2018 à 14:23, Frank David a écrit :<br>> > ><br>> > > Le 02/11/2018 à 22:27, Markus Metz a écrit :<br>> > >><br>> > >><br>> > >> On Fri, Nov 2, 2018 at 6:32 PM Frank David <<a href="mailto:frank.david@geophom.fr">frank.david@geophom.fr</a><br>> > >> <mailto:<a href="mailto:frank.david@geophom.fr">frank.david@geophom.fr</a>>> wrote:<br>> > >> ><br>> > >> > Le 02/11/2018 à 14:43, Markus Metz a écrit :<br>> > >> ><br>> > >> ><br>> > >> ><br>> > >> > On Fri, Nov 2, 2018 at 11:58 AM Frank David <<a href="mailto:frank.david@geophom.fr">frank.david@geophom.fr</a><br>> > >> <mailto:<a href="mailto:frank.david@geophom.fr">frank.david@geophom.fr</a>>> wrote:<br>> > >> > ><br>> > >> > > Le 31/10/2018 à 20:57, Markus Metz a écrit :<br>> > >> > ><br>> > >> > > Hi Frank,<br>> > >> > ><br>> > >> > > On Wed, Oct 31, 2018 at 1:39 PM Frank David<br>> > >> <<a href="mailto:frank.david@geophom.fr">frank.david@geophom.fr</a> <mailto:<a href="mailto:frank.david@geophom.fr">frank.david@geophom.fr</a>>> wrote:<br>> > >> > > ><br>> > >> > > > Hello Vero,<br>> > >> > > ><br>> > >> > > > Thanks for your reply. It took almost one hour to extract my<br>> > >> area (I guess because of topology). I do the same with Qgis from a<br>> > >> shape in less than 10 minutes (but with several operations)... it's<br>> > >> may be a bad idea to have a so heavy vector map in Grass...<br>> > >> > ><br>> > >> > > this kind of operation could be much faster in latest GRASS 7.4<br>> > >> than in any previous versions. Which GRASS version are you using?<br>> > >> > ><br>> > >> > > Hello Markus,<br>> > >> > ><br>> > >> > > I'm using 7.4.1 on debian9<br>> > >> > ><br>> > >> > > But my grass db folder  is on a small server (Synology NAS), so<br>> > >> may be the network makes slow down the operation ?<br>> > >> ><br>> > >> > This is possible, but it should not slow it down that much: 10 min<br>> > >> vs 1 hour<br>> > >> ><br>> > >> > Can you provide a link to the vector from which you want to extract<br>> > >> a small area, and the command you used to extract a small area? Maybe<br>> > >> there is scope for optimization.<br>> > >> ><br>> > >> > Markus,<br>> > >> ><br>> > >> > The data is CorineLandCover for france. You can download 400mo here<br>> > >> > <a href="https://we.tl/t-tvldlmkVw3">https://we.tl/t-tvldlmkVw3</a><br>> > >> ><br>> > >> > First I have imported the shape in PERMANENT<br>> > >> ><br>> > >> > and I do : v.extract --verbose input=CLC12_FR_RGF@PERMANENT<br>> > >> where=CODE_12 IN (111,112,131,311,312,313,324) output=clc12<br>> > >> ><br>> > >> > Hope it helps you...<br>> > >><br>> > >> import of the shapefile with v.in.ogr takes here 5 minutes<br>> > >> v.extract takes here 27 seconds instead of 1 hour<br>> > >><br>> > >> "here" means a laptop with Intel Core i7 4th generation and standard<br>> > >> SATA SSD, i.e. no fancy new high-speed system at all.<br>> > >><br>> > >> I am afraid this does not help you. What could take long is writing<br>> > >> attributes, depending on the database connection. I am using the<br>> > >> default sqlite database.<br>> > >><br>> > >> You could use top to check if v.extract is running at 100%.<br>> > ><br>> > > Hello Markus,<br>> > ><br>> > > Thank you for your help.<br>> > ><br>> > > I'm trying to do again the operation on local machine. I want fisrt,<br>> > > create a vector with only selected data (on full France area). But<br>> > > I've met a problem with v.extract.<br>> > ><br>> > > v.extract --overwrite --verbose input=CLC12_FR_RGF@PERMANENT type=area<br>> > > where=CODE_12 IN (111,112,121,311,312,313,324)<br>> > > output=CLC12_FR_RGF_BOISBATI@PERMANENT<br>> > > Loading categories from table <CLC12_FR_RGF>...<br>> > > DBMI-DBF erreur de pilote :<br>> > > SQL parser error (syntax error, unexpected NAME processing<br>> > > 'IN') in statement:<br>> > > SELECT cat FROM CLC12_FR_RGF WHERE CODE_12 IN<br>> > > (111,112,121,311,312,313,324)<br>> > > Unable to open cursor.<br>> > > DBMI-DBF erreur de pilote :<br>> > > SQL parser error (syntax error, unexpected NAME processing<br>> > > 'IN') in statement:<br>> > > SELECT cat FROM CLC12_FR_RGF WHERE CODE_12 IN<br>> > > (111,112,121,311,312,313,324)<br>> > > Unable to open cursor.<br>> > > ERROR: Unable select records from table <CLC12_FR_RGF><br>> > ><br>> > > It looks that grass does not recognized "IN" sql condition. If I force<br>> > > layer=-1, it runs bur the output vector is the same as input one.<br>> > ><br>> > > Am I wrong somewhere ? I use Grass 7.4.2 on debian9.<br>> > ><br>> > > Thank you for your help.<br>> > ><br>> > > Frank<br>> > ><br>> > Hi,<br>> ><br>> > I've repeated the operation v.in.ogr and v.extract with success on<br>> > different grass mapset, but I get always the same error (on v.extract :<br>> > see above) only with one of my local grass mapset (PERMANENT). I do not<br>> > understand what happens...<br>><br>> In the current mapset, the default database connection must not be DBF. You can reset the default database connection to sqlite with<br>><br>> db.connect -d<br>> After that verify with<br>> db.connect -p<br>> that sqlite is indeed the default database connection.<br>><br>> Now you can repeat with v.in.ogr + v.extract<br>><br>> Markus M<br>><br>> I'm afraid it's not the problem...<br>><br>> db.connect -d                                                                  <br>> Pilote par défaut / base de données définie à :<br>> pilote : sqlite<br>> base : /media/hdd1/grassdata/France_L93/PERMANENT/sqlite/sqlite.db<br>> (Wed Nov  7 10:46:01 2018) La commande s'est terminée (0 sec)                  <br>> (Wed Nov  7 10:46:29 2018)                                                     <br>> db.connect -p                                                                  <br>> driver: sqlite<br>> database: /media/hdd1/grassdata/France_L93/PERMANENT/sqlite/sqlite.db<br>> schema:<br>> group:<br>><br>> and when execute v.extract I get the same error...<br><div><br></div><div>Have you re-imported the vector?</div><div>What is the output of v.db.connect map=CLC12_FR_RGF -p</div><div>The "IN" SQL syntax works only if the vector CLC12_FR_RGF is not using dbf. You must make sure with "v.db.connect map=CLC12_FR_RGF -p" that CLC12_FR_RGF is really using sqlite as database driver.</div><div><br></div><div>HTH,</div><div><br></div><div>Markus M</div><div><br></div><div>></div>> I have tested again the same operation in other mapset than PERMANENT, and it works. I don't understand... if any idea ?...<br>><br>> Cheers,<br>> Frank<br>><br>><br>><br></div>