Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Difficulties running the vic_auto script #15

Open
WalterVIC opened this issue Jun 17, 2022 · 1 comment
Open

Difficulties running the vic_auto script #15

WalterVIC opened this issue Jun 17, 2022 · 1 comment

Comments

@WalterVIC
Copy link

Hello there, thanks for the script, it might help me a lot in preprocessing the VIC-data :) However, when running the vic_auto script, I get the following warning in the soil part and therefore a bad soil parameter file:

WARNING 000916: The input feature class does not appear to contain projected data.

Any suggestions on how to fix this? The input features seem to me to have a projection, all GCS WGS1984. The temp_soil csv contains all the data correctly but when it writes it to the soil_parameter_new file, it goes wrong.

Next to that I have a small second issue: in the vegetation parameter file preparation, I get the following error:

ValueError: labels ['OID'] not contained in axis in line 187

I am sure that the OID column is present, I can see that in the temp_lulc csv file. When replacing the OID label with the column index in the script, it does remove the column, but then the remainder of the script does not function anymore because the index values change. Thank already for helping me out!

@LSRathore
Copy link
Owner

Your first error seems to be because of the projection of the input files. Make sure that none of your input files have a projection (or they should have a geographic coordinate system). If the temp_soil.csv file is generated without any error then the final file should be generated too. If it does not work, you can modify the temp_soil.csv file following lines 162-171. It is just basic pandas and python.

Regarding the second issue, can you check the file named temp_lulc.csv in your workspace? Check if the column OID exists there. If not, then there must be the ID column with a different name. Replace that column name with the OID in the code (line 185 temp_l=temp_l.drop(["OID"],axis=1))

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants