<div dir="ltr"><br><br>On Sat, Mar 16, 2019 at 5:20 PM Ken Mankoff <<a href="mailto:mankoff@gmail.com">mankoff@gmail.com</a>> wrote:<br>><br>> Hi Moritz [and Markus]<br>><br>> On 2019-03-16 at 08:24 -0700, Moritz Lennert <<a href="mailto:mlennert@club.worldonline.be">mlennert@club.worldonline.be</a>> wrote...<br>> > Try running v.out.ogr three times, changing layer each time and using<br>> > the -u flag for updating the existing gpkg file:<br>> ><br>> > v.out.ogr -c input=pal output_layer=pal1 output=pal.gpkg layer=1<br>> > v.out.ogr -c input=pal output_layer=pal2 output=pal.gpkg layer=2 -u --o<br>> > v.out.ogr -c input=pal output_layer=pal3 output=pal.gpkg layer=3 -u --o<br>><br>> That seems to work. When I then read in the GeoPackage and diff, things look better.<br>><br>> v.in.ogr input=pal.gpkg output=pal2 --o<br>> diff <(v.db.connect -g pal) <(v.db.connect -g pal2)<br>> # no major diff, but table name is not preserved.<br>><br><div>> But the vector data is also included 3x now.</div><div><br></div><div>Why do you think the vector data are included 3x in the output? Did you check with ogrinfo on the different layers in the output GPKG?<br></div><div><br></div><div>> This may be due to v.patch? I'm adding Markus because he suggested the method to patch different layers. After I patch the three layers:</div>><br>> rm -fR ./nc_spm_08_grass7/gpkg<br>> grass ./nc_spm_08_grass7/PERMANENT<br>> g.mapset -c gpgk<br>> g.copy vector=poi_names_wake@PERMANENT,points<br>> g.copy vector=lakes@PERMANENT,areas<br>> g.copy vector=railroads@PERMANENT,lines<br>> v.category in=areas out=areas_2 op=chlayer layer=1,2 --o<br>> v.category in=lines out=lines_3 op=chlayer layer=1,3 --o<br>> v.patch in=points,areas_2,lines_3 output=pal --o<br>><br>> I would expect that <a href="http://v.info">v.info</a> reports different information if I use the "layer" option. But all three layers report IDENTICAL info. This diff command is silent:<br>><br>> diff3 -A <(<a href="http://v.info">v.info</a> pal layer=1) <(<a href="http://v.info">v.info</a> pal layer=2) <(<a href="http://v.info">v.info</a> pal layer=3)<br><div><br></div><div>The layer option of <a href="http://v.info">v.info</a> has only effect if used together with the -c flag (list columns for given layer). By default <a href="http://v.info">v.info</a> reports info on the whole vector. You might want to check instead with v.category option=report to see how many features are assigned to each layer. <br></div><div><br></div><div>></div>> Also, the final export commands suggested by Moritz:<br>><br>> v.out.ogr -c input=pal output_layer=tbl_pts output=pal.gpkg layer=1 --o<br>> v.out.ogr -c input=pal output_layer=tbl_areas output=pal.gpkg layer=2 -u --o<br>> v.out.ogr -c input=pal output_layer=tbl_lines output=pal.gpkg layer=3 -u --o<br>><br><div>> seem to export everything 3x. This is based on 1) the message "Exporting 15279 areas (may take some time)..." repeated each time 2) the file size increasing by 12 MB each time, and 3) reading in, displaying, and clicking brings up everything 3x in the "Query Results" window.</div><div><br></div><div>I suggest to check first with v.category input=pal option=report for the assignment of features to layers, then with ogrinfo on the output pal.gpkg.</div><div><br></div><div>Markus M<br></div>><br>> Should I just keep my vectors and databases each as a separate vector and not use layers?<br>><br>>   -k.<br>><br>><br></div>