[GRASS-user] Python Script not working now

[Phillip Allen]

Hi all,

I wrote a python script ~6 years ago to create sample catchment polygons
for stream sediment samples and all worked well. I created 10,000s of
catchements successfully!

Today I reconfigured the script for some new data and I am now getting this
error:

(Mon Aug  7 12:57:09 2023)

D:\DEM_12_5\mcc_area\basin_code04_mcc.py

Traceback (most recent call last):
  File "D:\DEM_12_5\mcc_area\basin_code04_mcc.py", line 43,
in <module>
    gsetup.init(gisbase, gisdb, location, mapset)
  File "C:\OSGeo4W\apps\grass\grass83\etc\python\grass\scrip
t\setup.py", line 329, in init
    raise ValueError(
ValueError: Mapset D:\DEM_12_5\grassdata\dem12_5_a is not
valid: <D:/DEM_12_5/grassdata> is not a valid GRASS Location
because PERMANENT Mapset is missing
(Mon Aug  7 12:57:10 2023) Command ended with non-zero return code 1 (0
sec)

I know the mapset PERMANENT exists and is valid since I just opened it in
GRASS.

Why is this python code below not running properly now?
Has something in GRASS changed 'recently' ?

regards,

Phil
Consulting Geochemist
(someday I really need to become a better python coder!)

#####
#####  basin python code
#! C:\OSGeo4W64\bin\python.exe

# IMPORT CSV & LOOP THROUGH COORDINATES FOR r.water.outlet
import os
import sys
import psycopg2


#set up GRASS environment variables
sys.path.append(os.path.join(os.environ['GISBASE'], 'etc', 'python'))
import grass.script as g
import grass.script.setup as gsetup
gisbase = os.environ['GISBASE']

gisdb = 'D:/DEM_12_5/grassdata'
location = 'dem12_5_a'
mapset = 'PERMANENT'
rast_draindir = 'mcc_d1_draindir150'


#     first join the new snap table's cat column to the original attribute
table
# v.out.ascii input=th_ss_aem_snap_r5 at PERMANENT type=point
output=D:/DEM_12_5/mcc_area/mcc_ss_aga_pt_temp.csv columns=sample_number
format=point separator=comma precision=20
txtfile = 'D:/DEM_12_5/mcc_area/mcc_ss_aga_pt_temp.csv'
txtfileWKT = 'D:/DEM_12_5/mcc_area/mcc_ss_WKT.txt'
pgsrid = '32618'

pg_user = 'postgres'
pg_user_pwd = 'notmyrealpassword'
pg_schema = 'basin'
pg_tbl_name = 'mcc_samp_basin'
pg_host = 'localhost'
pg_dbname = 'ss_basin_test'

txt_pg_tbl_create = 'CREATE TABLE ' + pg_schema + "." + pg_tbl_name + ' (x
double precision, y double precision, samp_id character varying(254), geom
geometry(POLYGON, ' + pgsrid + ') );'
txt_pg_db =  "host='" + pg_host + "' dbname='" + pg_dbname + "' user='" +
pg_user + "' password='" + pg_user_pwd + "'"

txt_gpg_db = "PG:dbname=" + pg_dbname + " host='" + pg_host + "' user='" +
pg_user + "' password='" + pg_user_pwd + "'"


#print "Connecting to database\n ->%s" % (txt_pg_db)
#print "Grass PG Connection String\n ->%s" % (txt_gpg_db)

gsetup.init(gisbase, gisdb, location, mapset)
g.run_command('g.region', flags='a', rast = rast_draindir)


# Create PostGIS Polygon Table
conn = psycopg2.connect(txt_pg_db)
cur = conn.cursor()
cur.execute(txt_pg_tbl_create)
conn.commit()

#csv file reading and importation
import csv
rows = list(open(txtfile))
totalrows = len(rows) - 1

#loop through the csv(coordinates) file in r.water.outlet module
f = open(txtfile, 'r')
element = list(csv.reader(f))
i = 0
j = 0
while True:
if i <= totalrows:
g.run_command('r.water.outlet', overwrite = True , input = rast_draindir,
output = 'b' , coordinates = element[i][j] + ',' + element[i][j + 1] )
g.run_command('r.null' , map='b', setnull = 0 )
g.run_command('r.to.vect' , overwrite = True , input='b', output = 'b_rc',
type = 'area' )
g.run_command('v.out.ascii', overwrite = True, input='b_rc',
type='boundary', output=txtfileWKT, format='wkt')
g.run_command('v.out.postgis', overwrite = True, input='b_rc',
type='boundary', output=txt_gpg_db , output_layer = pg_schema + ".b_rc" )
sql_txt = "INSERT INTO " + pg_schema + "." + pg_tbl_name + " SELECT " +
element[i][j] + " AS x, " + element[i][j + 1] + " AS y, '" + element[i][j +
3] + "' AS samp_id, b_rc.geom AS geom FROM " + pg_schema + ".b_rc;"
cur.execute(sql_txt)
conn.commit()
i = i + 1
else:
break

cur.close()
conn.close()
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20230807/c87efcf6/attachment.htm>