This tutorial covers the complete procedure to create a land cover change raster from a comparison of generated vegetation index (NDVI) rasters by the use of Python and the Numpy and Rasterio libraries. Results of the NDVI for given years and NDVI change are plotted on Jupyter Lab as color grid and contour grid.
Read MoreGeopandas is an amazing library for spatial analysis since it combines the spatial tools from Shapely and Fiona with the versatility of Pandas Dataframes. Rasterio is one of the most complete libraries to deal with raster files in Python.
We are aware that most geoscientists, water resources specialists and related professionals work on Windows, therefore we are always in the search of new ways to get Python working with all its geospatial capabilities in every computer.
Based on the new Anaconda releases now we can install Geopandas and Rasterio in a easy way.
Read MoreWe are allways on the research of new visualization tools for groundwater models and that would imply the review of different apis, codes and backends. This time we have done an applied example for the interactive generation of cross section from a Modflow 6 model with a voronoi mesh on a Jupyter notebook. The cross section shows the model grid, head distribution as a color grid and contours and its generated from the line delineation over the 2d view of the model grid. The code can be easily adapted to any Modflow model by editing the Python and Flopy code.
Read MoreWe have developed a version of our Angascancha basin case study that runs on MODFLOW 6 with discretized by vertices (DISV). The model implements a Voronoi mesh that is generated from the basin boundary and river network and the refinement levels are defined by a minimum and maximum cell size together with determined refinement steps. The tutorial is complete in all the steps involving the model discretization, construction, simulation and visualization.
Read MoreIf you have a set of drone imagery and you want to know the location and direction of the camera, this tutorial might be interesting for you. We have done an applied example that retrieves the geospatial metadata from the drone camera for a group of images and makes a geostapatial repretacion on a map with the image name available as popup. The tutorial is done under a Jupyter notebook with Python and Folium.
Read MoreAnalysis of groundwater chemistry is a difficult task for the limited set of monitoring samples, the limited samples and the limited analized components. In order to assess the actual extension of a contamination plume or the efficiency of remediation techniques we need new and innovative methods to plot and analyze water chemistry data with open source data. We have done an applied case of interactive VOCs and PFAs representation on a Jupyter notebook with Python, Folium and Ipywidgets. The dataset has more than 3300 samples of 127 points over a period of 30 years and corresponds to a contaminated site of a former airfield.
Read MoreOpen lithology datasets are scarce and most times they come in certain exchange formats that can´t be easily coupled with other softwares. Python has the tools to read and extract data from those files and provide them in more friendly formats as csv or xlsx. With the use of spatial libraries as Fiona or 3D visualization libraries as Pyvistas we can go one step further and process our data as shapefiles or vtks.
Read MoreGempy is an open-source library for modeling geology written in Python. The library is capable of creating complex 3D geological models including structures, fault networks, and unconformities and it can be coupled with uncertainty analysis.
Due to the required packages and specific configuration the installation of Gempy is a challenge on Windows, we have developed an alternative installation process by the use of a Debian kernel under a Windows Subsystem for Linux. This tutorial shows the complete procedure to install Gempy from the Debian kernel setup to the review of the Gempy in Python.
Read MoreStorage, management and analysis of geospatial vector data as an ESRI shapefile is a common procedure of GIS and related professionals. The generation of these spatial files can be done not only on a desktop software but also by Python commands. We have created an applied example that shows the procedure in Python to create point, line, and polygon shapefiles from a csv file by the use of the Fiona library.
Read MorePython is a great tool for spatial analysis and geomachine learning, however sometimes the Windows operating system presents some difficulties to install and run the bunch of Python libraries such as Gdal, Fiona, Geopandas, Rasterio. We are aware that most geoscientists, water resources specialists and related professionals work on Windows, therefore we are always in the search of new ways to get Python working with all its geospatial capabilities in every computer. We have created a tutorial that shows the installation process of the Python geospatial libraries in Windows by the use of a Conda environment; the process is simple on its steps, however the sequence and factors related to the package compatibility are important on the installation.
Read MoreWe did a simple procedure in QGIS but unknown to us that extracts the centroids of lines and polygons with Field Calculator commands inside the attribute table. The procedure is straightforward and does not involve the use of any intermediate layer.
Read MoreOver the last months we have researched the different tools in Python and QGIS available to recognize crops and vegetation as geospatial vector files. We have used a variety of techniques that range from machine learning algorithms with Scikit Learn and Scikit Image to just more innovative band combinations and reclassifications in QGIS. This article shows the summary of the tutorials produced so far that we are sure will be very helpful for GIS professionals and geoscientists.
Read MoreSpatial analysis and machine learning sometimes require massive coding in order to achieve decent results such as identifying plants from a drone orthophoto. We wanted to create a simple workflow for beginner and intermediate Python users to work with these libraries without much pain or frustration. This tutorial has the complete procedure to use a Python class that recognizes plants from an orthophoto based on sample points and creates intermediate plots and identifies plants as point shapefiles.
Read MoreWhile doing research on machine learning algorithms for drone orthophotos we found that crops can be delineated with just standard QGIS with excellent performance. Based on the addition of the blue and red band divided by the green band we can have a new vegetation index where the most healthy crops have low index values (0-1.5) and the dry crops / barren soil have high values (more than 1.9). This tutorial shows the complete procedure on QGIS to perform the delineation from bean plants on a drone orthophoto with a resolution of 5cm.
Read MorePython is a programming language capable of performing calculations for hydrological studies and water resources evaluations. We have done a tutorial for the volume-elevation curve determination of the lake Patillas in Puerto Rico with Python and numerical / spatial libraries as Numpy and Rasterio. Finally, results were compared to the volume-elevation curve form a USGS survey. The procedure was done this time for a lake, but can be easily applied to any reservoir or water body when the bottom elevation is available as a raster file.
Read MorePython is a useful tool for data analysis but also for data representation and as a graphic tool. This is an applied tutorial for the representation of a wind rose with Python from wind speed and direction stored on an Excel spreadsheet. The tutorial explores the options of the library to represent windroses as bars, boxes, polygons or contours.
Read MoreTutorial for the calculation of the solubility and thermodynamic stability of gypsum and anhydrite. The example simulates the dissolution of two minerals in a beaker at equilibrium, and the beaker is heated step-wise from 25ºC to 75ºC. Concentrations and saturation indices for the initial solution and the batch reaction are shown as Pandas dataframes and plotted as bar diagrams on a Jupyter Notebook. Finally a plot of the saturation index for anhydrite and gypsum with temperature is generated from an iteration over the batch reaction steps.
Read MoreThere are three essential things for the representation of geological data on QGIS: the spatial geological information, an appropriate symbology and the software knowledge. Once these three things are available, the potential of QGIS to represent geological maps is unlimited.
We have done a tutorial for the representation of synclines, anticlines, overturned synclines and overturned anticlines on the regional scale with QGIS.
Read MoreThe speciation modeling allows to calculate the distribution of aqueous species in a solution. Phreeqc is capable to simulate this speciation calculation and we are going to demonstrate this capability on a study case of aqueous species in seawater.
We have done a tutorial for the speciation modeling of seawater with Phreeqc that runs under Python in a Jupyter Lab enviroment. The code can run the Phreeqc executable, define the databases and stablish the output files. Results from simulation are available as Pandas dataframes and plots are made for the main components and the distribution of saturation indices.
Read MoreIn our perspective, 3D visualization of geospatial data has been a long desired feature that has been covered in some features from SAGA GIS or in some plugins from QGIS. This time we developed a Python script that converts point / line / polygon ESRI shapefiles (or any vector file) to unstructured grid Vtk format type (Vtu) by the use of the Python libraries Geopandas and Pyvista. The tutorial has files, scripts, and videos that show the whole procedure with some remarks on the software and spatial files and a discussion about the nature of the spatial files that presents some challenges in the data conversion.
Read MoreHatarilabs presents its educational program designed for mastering Python in real professional and academic environments. The program has an extensive practical work that goes from the basic concepts of Python, Numpy and Pandas to specific applications in water resources and geosciences coupled with geospatial analysis and machine learning.
We have designed a complete program on MODFLOW with Model Muse that ranges from the essential topics in flow and transport modeling to applied cases of regional modeling, infrastructure projects, and seawater intrusion. MODFLOW with Model Muse is a recent open-source software that provides various tools for simulating the groundwater flow regime in different environments.
This program is designed to provide systematical learning from the most basic and general topics to intermediate and advanced topics, allowing the student to understand the various steps in developing hydrological models using the SWAT+ interface for QGIS. This program has a series of theoretical and practical concepts that allow the student to process spatial data, manage meteorological databases, build hydrological models, develop sensitivity analyses, and perform manual and automatic calibrations using versatile tools such as SWAT+ Toolbox, SQLite, and Python.
ONLINE COURSES:
Phreeqc is a software developed by the USGS written in C++ capable of modeling a variety of aqueous geochemical processes such as mixing of waters, modeling equilibrium between solid and aqueous phase, modeling impact of temperature, calculation of element concentration among others.
Hatarilabs has developed an applied and unique course of Phreeqc integrated with Python undes a class (object type) that can run, parse and provide Phreeqc output as Pandas dataframes and integrate it to plots and further analysis with other Python libraries on a Jupyter notebook.
HEC-RAS is a software designed to perform hydraulic 1D and 2D calculations on river channels either natural or man constructed. It allows steady flow, unsteady flow, sediment transport and water quality modelling. It is developed by the US army corps of engineers, US Department of defence and it is openly available software. The HEC-RAS steady flow is based on one dimensional energy equation, while for the unsteady flow, HEC-RAS uses an implicit, finite difference method to solve the fully dynamic Saint-Venant equations.
The objective of HEC-RAS is to solve river dynamics and river engineering problems. Applications are delimitation of flood areas, embankment and protection work design, river restoration, emergency plans for dam break, optimization of hydraulic work and risk assessment and management.
This course develops the main functions and applications of the latest version of the MODFLOW 6 groundwater modeling code through the ModelMuse interface, both developed by the United States Geological Survey (USGS). This version includes innovative tools for the construction and simulation of hydrogeological models, mainly highlighting the incorporation of the discretization option for discretized by vertices grids.
The development of machine learning and geospatial libraries in Python as Scikit Image, Rasterio and Fiona gave us a new range of tools to analyze vegetation from traditional sources as satellite imagery and new sources as drone orthophotos. The appliance of Python programming for crop identification depends on many factors such as the image resolution, the vegetation stage and even the algorithm setup.
We have developed a course crop identification / delineation with geospatial and machine learning tools of Python, the course is aimed to GIS or related professionals that have basic Python knowledge. The course covers the introductory concepts of the geospatial libraries and a series of machine learning applications for crop indentification in olive trees, palms, agave and corn fields.
Note:
This course is for professionals with previous knowledge of Python. Basic concepts of Python programming will not be covered in the lessons and course material. In case you don’t have any knowledge of Python, we recommend you to take this online course.
The Hydrologic Modeling System (HEC-HMS) is a free program developed by the US Corps of Engineers that implements a series of hydrological methods to represent different physical processes of the water cycle.
HEC-HMS also includes procedures necessary for continuous simulation including evapotranspiration, infiltration, hydrological routing and soil moisture accounting. The tools and options of HEC HMS make it a very versatile and powerful software for the hydrological simulation of different scenarios such as extreme events in arid regions, or water balances in wet climates.
Hatarilabs offers the online course of Hydrological Modeling with HEC-HMS for students and professionals of any country. The course is focused on theory and practical to know the different HEC-HMS tools and its application to the hydrological modeling.
Hatarilabs offers the online course of Hydrological Modeling with QSWAT for students and professionals of any country. The course is focused on the acquisition and preprocessing of data for SWAT+ using QGIS 3.4 for develop hydrological watershed modeling.
This course covers the whole procedure of developing regional groundwater models through the ModelMuse interface, including the model construction, boundary conditions, parameter assignment, model simulation, results analysis and visualization, and model calibration. Flopy is also introduced in the course to import, create, refine and simulate the model as in ModelMuse. Two basins with different aquifer setting and boundary conditions are used as case studies to demonstrate the use of these tools so that the participants learn how to build these models and analyze the results for future decision-making.
Evaluation of long term impacts of the mining activity is a key topic on the public concern and regulatory and sustainability policies. Due to the scale of the mining projects and the different related activities, impacts could be made to several parts of the physical environment. Groundwater modeling is a useful tool for the simulation of mine closure impacts to the groundwater flow system and groundwater quality. Hatarilabs has developed an applied course about the simulation, calibration, representation and impact analysis of a waste dump closure on the environment and water bodies with a long term scope (closure and post closure). The course will be given on Model Muse for flow and contaminant transport modeling and scripts in Flopy for model calibration, water balance analysis and evaluation of remediation techniques performance.
This course is focused on how to use the different packages related to numerical tools, data analysis, representation of graphs, and others. Many of these tools are applicable to the field of Hydrology in the preparation of data for hydrological models, statistical analysis of hydrological parameters and model result representation.
Graphical user interfaces (GUIs) are commonly used to construct and post-process numerical groundwater flow. Flopy is the Python library that builds and executes MODFLOW models; this library has tools and options for the complete geo-referencing of a MODFLOW model.
Python is a simple and powerful programming language; its simplicity is remarkable compared to other programming languages and its power is based on the number of tools available for different areas of study.
Use of Python with the available scientific packages and FloPy, facilitates data exploration, alternative model evaluations and model analyses that can be difficult to perform with GUIs. One advantage of Flopy and Python is that there are many packages available to facilitate the model development process, including packages for plotting, array manipulation, optimization, and data analysis.
This course allows acquired knowledge of the open source programming language to construct model input files, run the model, and read and plot simulation results.
LATEST POSTS
With the upcoming release of mf6Voronoi v0.0.38, optimizing mesh generation workflows remains a top priority. Constructing complex unstructured Voronoi grids for MODFLOW 6 can become computationally intensive as spatial resolution increases. To address this challenge, parallel execution capabilities were introduced into mf6Voronoi leveraging dask-geopandas.
This tutorial explores a powerful interactive GeoPandas tool that facilitates the visual analysis and representation of water chemistry data. Through this Python workflow, you can display sampling points over satellite imagery and various basemaps, map concentrations using automatic color scales, and adjust marker sizes based on analyzed parameters. This direct integration allows you to customize the variables displayed on hover, streamlining data evaluation for baseline studies or water pollution diagnostics.
Interpolating large volumes of spatial data through the QGIS graphical interface is often a headache: the process takes too long or, even worse, the program freezes while consuming all system resources. In this tutorial, we will show you how to interpolate nearly 2 million points in less than two minutes using the OSGeo4W Shell and GDAL commands with the IDW method. By creating a virtual layer file (.vrt), you will learn a highly efficient, reliable, and lightweight workflow that enables you to work with millions of data points without being constrained by the limits of the GUI.
This tutorial demonstrates how to accurately and efficiently extract the water surface elevation at a specific point from a HEC-RAS 2D unsteady flow model using Python in Jupyter Lab. Since water levels vary across the mesh and native software tools can limit automation, the author teaches how to open HDF5 output files (.hdf), load a point shapefile using GeoPandas, and apply a geospatial index K-d tree to locate the nearest computational cell. Finally, the time series data is processed to generate a water level hydrograph, which is essential for hydrogeological model calibration.
estimations of crop water requirements. This applied case offers a deep dive into modeling crop evapotranspiration using the standardized FAO-56 Dual Crop Coefficient methodology, implemented programmatically through Python.
Using a real-world case study of olive cultivation (Olea europaea) in the hyper-arid coastal desert of Bella Unión (Arequipa, Peru), we will demonstrate how to build an end-to-end hydrological data pipeline. Participants will learn how to integrate spatial data, process satellite-derived and weather station climate parameters, define soil-water-crop constraints, and simulate complex deficit drip-irrigation schedules.
SEARCH
 |
 |
 |
 |
 |
 |
FloPy is the Python library that builds and executes MODFLOW models; this library has been enhanced to provide full support of MODFLOW 6 with most of its recent development is related to functionality for MODFLOW 6, tools to use vector and raster spatial data and common plotting and export functionality.