[Qgis-user] Problem Running v.sample in QGIS 3

Terry Morse Tue, 17 Sep 2019 02:33:14 -0700

I am trying to run the GRASS tool v.sample in 64 bit QGIS 3.8.2 and 3.4.11 with GRASS 7.6.1, and 3.2.3 with GRASS 7.4.1 on a 64-bit Windows 10 laptop, following directions in Learn QGIS, Fourth Edition: Your step-by-step guide to the fundamental of QGIS 3.4, by Andrew Cutts and Anita Glaser (Packt Publishing, 2018).

The data layers are from the QGIS Alaska sample data set: airports.shp and landcover.img, both EPSG 2964.. I've attached a screen capture of the dialog box showing the main settings. All advanced settings are the defaults.

When I run v.sample, I get the following errors:

1) WARNING: Concurrent mapset locking is not supported on Windows. [I've read on the internet that this error isn't significant]

2) WARNING: Unable to initialise PROJ with the following parameter list: +proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 +x_0=0 +y_0=0 +no_defs +a=6378206.4 +rf=294.9786982 +nadgrids=C:\PROGRA~1\QGIS3~1.4\share\proj\@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat
3) WARNING: Unable parse GRASS PROJ_INFO file
4) ERROR: Unable to create OGR spatial reference

5) The following layers were not correctly generated.<ul><li>C:/Users/Terry/Documents/My Books/Learn QGIS 4th ed/Data/samples_2.shp</li></ul>You can check the 'Log Messages Panel' in QGIS main window to find more information about the execution of the algorithm.

Perplexingly, the algorithm ran successfully back in January, but not this week. I ran it from the processing history from January to ensure the settings were all the same. The only difference I can think of between then and now is that I was running QGIS under 64-bit Windows 7 at the time. I installed Windows 10 a couple weeks ago. I'm not sure which version of QGIS 3 I had installed at the time, which is why I tried running it with 3.2, 3.4, and 3.6 this week. It didn't work in any of them.

Can anyone help me figure this out? I've copied the full log file below.

Thank you,

Terry Morse

Full log file:

QGIS version: 3.8.2-Zanzibar

QGIS code revision: 4470baa1a3

Qt version: 5.11.2

GDAL version: 2.4.1

GEOS version: 3.7.2-CAPI-1.11.0 b55d2125

PROJ version: Rel. 5.2.0, September 15th, 2018

Processing algorithm…

Algorithm 'v.sample' starting…

Input parameters:

{ 'GRASS_MIN_AREA_PARAMETER' : 0.0001, 'GRASS_OUTPUT_TYPE_PARAMETER' : 0, 'GRASS_REGION_CELLSIZE_PARAMETER' : 0, 'GRASS_REGION_PARAMETER' : None, 'GRASS_SNAP_TOLERANCE_PARAMETER' : -1, 'GRASS_VECTOR_DSCO' : '', 'GRASS_VECTOR_EXPORT_NOCAT' : False, 'GRASS_VECTOR_LCO' : '', 'column' : 'ID', 'input' : 'C:/Users/Terry/Documents/GIS/QGIS Sample Data/qgis_sample_data/shapefiles/airports.shp', 'method' : 0, 'output' : 'C:/Users/Terry/Documents/My Books/Learn QGIS 4th ed/Data/sampled_points.shp', 'raster' : 'C:/Users/Terry/Documents/GIS/QGIS Sample Data/qgis_sample_data/raster/landcover.img', 'zscale' : 1 }

g.proj -c proj4="+proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 +x_0=0 +y_0=0 +datum=NAD27 +units=us-ft +no_defs"

v.in.ogr min_area=0.0001 snap=-1.0 input="C:\Users\Terry\Documents\GIS\QGIS Sample Data\qgis_sample_data\shapefiles\airports.shp" output="vector_5d809faeb5d0d2" --overwrite -o

r.external input="C:\Users\Terry\Documents\GIS\QGIS Sample Data\qgis_sample_data\raster\landcover.img" band=1 output="rast_5d809faebafab3" --overwrite -o

g.region n=7809680.0 s=1367760.0 e=4897040.0 w=-7117600.0 res=3280.0

v.sample input=vector_5d809faeb5d0d2 column=ID raster=rast_5d809faebafab3 zscale=1 method="nearest" output=output3138c3949b014923adb9b8d44c37e6ba --overwrite

v.out.ogr type="point" input="output3138c3949b014923adb9b8d44c37e6ba" output="C:\Users\Terry\Documents\My Books\Learn QGIS 4th ed\Data\sampled_points.shp" format="ESRI_Shapefile" --overwrite

Starting GRASS GIS...

WARNING: Concurrent mapset locking is not supported on Windows

Executing <C:\Users\Terry\AppData\Local\Temp\processing_2b250d50fb6d452a811e6d4d9b4d009b\grassdata\grass_batch_job.cmd> ...

C:\Users\Terry\Documents>chcp 1252 1>NUL

C:\Users\Terry\Documents>g.proj -c proj4="+proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 +x_0=0 +y_0=0 +datum=NAD27 +units=us-ft +no_defs"

Default region was updated to the new projection, but if you have multiple mapsets `g.region -d` should be run in each to update the region from the default

Projection information updated

C:\Users\Terry\Documents>v.in.ogr min_area=0.0001 snap=-1.0 input="C:\Users\Terry\Documents\GIS\QGIS Sample Data\qgis_sample_data\shapefiles\airports.shp" output="vector_5d809faeb5d0d2" --overwrite -o

Over-riding projection check

Check if OGR layer <airports> contains polygons...

0..2..5..7..9..11..13..15..17..19..21..23..25..27..30..32..34..36..38..40..42..44..46..48..50..52..55..57..59..61..63..65..67..69..71..73..75..77..80..82..84..86..88..90..92..94..96..98..100

Creating attribute table for layer <airports>...

Importing 76 features (OGR layer <airports>)...

0..2..5..7..9..11..13..15..17..19..21..23..25..27..30..32..34..36..38..40..42..44..46..48..50..52..55..57..59..61..63..65..67..69..71..73..75..77..80..82..84..86..88..90..92..94..96..98..100

-----------------------------------------------------

Building topology for vector map <vector_5d809faeb5d0d2@PERMANENT>...

Registering primitives...

C:\Users\Terry\Documents>r.external input="C:\Users\Terry\Documents\GIS\QGIS Sample Data\qgis_sample_data\raster\landcover.img" band=1 output="rast_5d809faebafab3" --overwrite -o

Over-riding projection check

Reading band 1 of 1...

Link to raster map <rast_5d809faebafab3> created.

C:\Users\Terry\Documents>g.region n=7809680.0 s=1367760.0 e=4897040.0 w=-7117600.0 res=3280.0

C:\Users\Terry\Documents>v.sample input=vector_5d809faeb5d0d2 column=ID raster=rast_5d809faebafab3 zscale=1 method="nearest" output=output3138c3949b014923adb9b8d44c37e6ba --overwrite

Reading points...

2..5..9..13..17..21..25..28..31..34..38..42..46..50..53..56..59..63..67..71..75..78..81..84..88..92..96..100

Building topology for vector map <output3138c3949b014923adb9b8d44c37e6ba@PERMANENT>...

Registering primitives...

C:\Users\Terry\Documents>v.out.ogr type="point" input="output3138c3949b014923adb9b8d44c37e6ba" output="C:\Users\Terry\Documents\My Books\Learn QGIS 4th ed\Data\sampled_points.shp" format="ESRI_Shapefile" --overwrite

WARNING: Unable to initialise PROJ with the following parameter list: +proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 +x_0=0 +y_0=0 +no_defs +a=6378206.4 +rf=294.9786982 +nadgrids=C:\PROGRA~1\QGIS3~1.8\share\proj\@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat

WARNING: Unable parse GRASS PROJ_INFO file

ERROR: Unable to create OGR spatial reference

C:\Users\Terry\Documents>exit

Execution of <C:\Users\Terry\AppData\Local\Temp\processing_2b250d50fb6d452a811e6d4d9b4d009b\grassdata\grass_batch_job.cmd> finished.

Cleaning up default sqlite database ...

Cleaning up temporary files...

Press any key to continue . . .

Execution completed in 17.65 seconds

Results:

{'output': <QgsProcessingOutputLayerDefinition {'sink':C:/Users/Terry/Documents/My Books/Learn QGIS 4th ed/Data/sampled_points.shp, 'createOptions': {'fileEncoding': 'System'}}>}

Loading resulting layers

The following layers were not correctly generated.<ul><li>C:/Users/Terry/Documents/My Books/Learn QGIS 4th ed/Data/sampled_points.shp</li></ul>You can check the 'Log Messages Panel' in QGIS main window to find more information about the execution of the algorithm.

_______________________________________________
Qgis-user mailing list
Qgis-user@lists.osgeo.org
List info: https://lists.osgeo.org/mailman/listinfo/qgis-user
Unsubscribe: https://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Problem Running v.sample in QGIS 3

Reply via email to