Wiki » History » Version 237
« Previous -
Version 237/317
(diff) -
Next » -
Current version
Giulio Di Anastasio, 03/05/2021 10:46
Wiki¶
CSR Geomatics Team¶
Collaborations¶
Data_analysis¶
CSR server¶
Our workflow¶
Surveys¶
Most field surveys are executed by our team of surveyors, using advanced DGPS equipment.
Other surveyors might also be contracted.
CAD¶
The survey data are imported to a CAD software (Aurocad/Civil 3D).
Editing shapefiles in Autocad¶
FEATURES CREATION¶
1. Assign a CRS to the drawing (TM-AUSPOS) (MAPCSLIBRARY command)
2. Create features in CAD (Points, lines, polygons)
3. Export shapefile (a) from CAD (Output > DWG to SDF) (Convert to LL84 – 3D)
FEATURES IMPORT INTO DB¶
4. Create zip file of the shapefile
5. Upload into the GISAF Shapefiles Basket
6. Import the shapefile into DB
7. Delete the shapefile from Local Machine
FEATURES EDITING¶
8. Open the table in QGis
9. Save as a shapefile (b) in TM AUSPOS CRS
10. In CAD, open a new drawing and assign ASUPOS CRS
11. Import the shapefile (b) (MapImport) with all Object Data
12. Edit features
13. Export shapefile (a) from CAD (Output > DWG to SDF) with ONLY the id selected (Data Tab > Select Attributes > Object Data > Filename > id) (Convert to LL84 – 3D)
FEATURES IMPORT INTO DB¶
14. Create zip file of the shapefile
15. Upload into the GISAF Shapefiles Basket
16. Import the shapefile into DB
17. Delete the shapefile from Local Machine
QGis¶
Conventions¶
Shapefiles¶
We work with QGis , a widely used open source, free software for working on maps and geographical data.
Shapefile is a standard file format for mapping, that Gisaf can import. QGis can open and save shapefiles.
We have defined some simple norms about these shapefiles for integration with Gisaf.
Coordinate system¶
We use CRS SRID 32644.
Column (attributes)¶
- All objects in a shapefile (layer) have a unique identifier named "id": numerical value.
Field (attribute) names¶
- All fields are lower case (
UpperCase, use:lowercase
)
- They do not contain spaces, but underscores (
not this, use:but_that
)
- The field names cannot have more than 8 characters (
this_is_too_long, use:this_short
)
Attribute types¶
- The dates are written in ISO format:
YYYY-MM-DD
(eg.1968-02-25
)
Fields to remove¶
- Eventually, remove the fields containing the coordinates: northing, easting, elevation, latitude, longitude, area, length, etc (these might be present when the data is imported from speadsheet)
Foreign keys¶
We often deal with labels or categories of objects in a layer.
A common use case to explain: Alice creates a layer of stones, and wants to tag each stone with a rating: these are picked from a list of choices, like: Beautiful, Interesting, Pretty, etc.
For these kind of attribute: define a column like rating_id
(something ending with _id
) as a numerical value, and assign values 1, 2, 3, etc. The text is defined in another table (typically a CSV file), that looks like:
id | name |
---|---|
1 | Beautiful |
2 | Interesting |
3 | Pretty |
Code¶
We have defined a standard set of codes, that defines the type of data. They can be found here (TODO: add link).
Add a column code_name
, matching with one the code, eg. V25
for TDEF.
Surveyor¶
We keep a record of the people who realized the surveys (the surveyors).
The shapefiles must contain an attribute srvyr_id
, which refers to this table (TODO: add link).
Accuracy¶
We keep a record of the accuracy of the surveys.
The shapefiles must contain an attribute accur_id
, which refers to this table (TODO: add link).
Date of survey¶
As nothing is absolutely permanent, it's also important to keep track of the date of the surveys: the shapefiles must contain an attribute date
.
Working with Gisaf¶
Survey data¶
Raw survey data are contained in CSV files, typically downloaded from surveying instruments.
See more information of the process for the survey data (including a flow diagram)
Survey data¶
Workflow summary¶
Import basket¶
Gisaf provides an "import basket" dedicated for raw survey data, which is generated by land survey equipment (Leica's Total Station and RTK).
These are CSV files, like:
100081,370633.969,1327742.157,51.187,,, 100083,370628.876,1327702.913,51.565,T52,, 100082,370628.729,1327720.019,51.261,T52,, 100081,370633.969,1327742.154,51.179,,, 100083,370628.876,1327702.913,51.565,T52,, 20800,370633.969,1327742.154,51.180,,, 20801,370618.795,1327713.172,52.817,E30,, 20802,370623.674,1327711.436,51.283,B35,, 20803,370619.314,1327713.407,51.383,B35,,
Each category (5th column) must be defined in the Category table (see Categories).
Organization of the raw survey data basket¶
The basket should be organized in a directory structure:
- Project name (these can be themselves put in a hierarchy of (sub)directories)
- Surveyor's organization
- Equipment (eg. TS, RTK)
- Survey files (eg. Our_project-Some_comment-2018-02-23.txt
)
Format of the survey file names¶
Our_project-Some_comment-2018-02-23.txt
The date of the survey follows the ISO date standard: YYYY-MM-DD
.
Import to the database¶
When importing raw survey data files to the database, Gisaf does 2 steps as described below.
It's worth noting that, in this process, no reprojection is done.
Feed the raw_survey table¶
Each point of the imported raw survey data file is inserted to the raw_survey table:
- Creation of a Point geometry: the raw_survey table has a geometry column for a single point (
geom
) with x,y and z coordinates - Save the
id
of the original point to theorig_id
column - A unique
id
is computed from the following fields:id
,project
,equipment
,date
- The project is saved in the
project_id
column - The surveyor identification in
srvyr_id
- The date of survey is saved in the
date
column - The accuracy is tagged in the
accur_id
, according to a mapping defined in theaccuracyequimentsurveyormapping
table, which depends on the surveyor and equipment - The category of the point
Feed the RAW_V_*
tables¶
From the raw_survey
table, each point is then copied to its respective RAW_V_
table, with basically the same information.
These tables (which should be created manually or with the admin notebook called create_tables
, as of today's writing), only contain points.
The project is saved along: see below.
Import the points¶
For categories that define points (opposite to lines and polygons, which require line work carried over in CAD or in a GIS software, see Line work), the points can be imported automatically to their final destination: the V_*
tables.
Note: in this process, the geometries are reprojected.
Auto import of the points¶
The points found in the RAW_V_*
tables can be imported automatically, project per project, the project page of the admin interface.
Import of the line work (lines and polygons)¶
See Line work - to be written¶
The shapefiles generated manually (line work) should be put in the project's basket, and imported from it.
Categories¶
The categories define the types of the geographical features and they are mapped according to ISO standard layer naming conventions: see https://www.nationalcadstandard.org/ncs5/pdfs/ncs5_clg_lnf.pdf
Gisaf uses:
- a table
category
where the layers are defined - a table per category
Fields for the categories¶
to be written - TODO
Creation of the RAW_* tables¶
This step must be done manually (as of today's date of writing).
QGis: work on shapefiles¶
Go to shapefiles - to be written¶
Backup and restoration of the database¶
Database¶
Troubleshooting¶
Layers missing in the map's tree¶
Gisaf relies on counting features through Postgres statistics collector subsystem.
In case the server is restarted dirty (eg. without clean shutdown), then the count of the tables might be wrong or just 0, leaving the layers apparently empty and thus not even appearing.
The fix is as easy as:
sudo -u postgres psql avgis -c VACUUM
Installation¶
This documentation assumes that the Postgis package has been installed (see CSR_server).
Configure the server¶
Allow connections from other hosts in the local network¶
Set the server to listen to addresses, set listen_addresses to *
in /etc/postgresql/9.6/main/postgresql.conf
.
Allow the connections, add in /etc/postgresql/9.6/main/pg_hba.conf
:
host all all 192.168.0.0/24 md5
Creation of the database¶
As postgres
user:
createdb -E utf8 -T template0 avgis
Backups¶
Primary¶
The database is backed up every day at midnight. The dump file is located in /var/backups/postgres/
.
Secondary¶
There are other backups (daily, weekly, monthly) thanks to Debian package autopostgresqlbackup
), located (default) in /var/lib/autopostgresqlbackup
.
Tertiary (dom0)¶
The whole virtual machine is backed up by BackupNinja on the "dom0" controller, using:
- rdiff backups every day
- tar files on Saturdays.
See /etc/backups.d
on the dom0 (192.168.0.12).
Remote¶
TODO: remote backup to be written
Restoration¶
If the VM is not shutdown properly, there's a chance that the database is corrupt, and needs to be restored from one of the backups.
After the restoration, restart gisaf:
systemctl restart uwsgi.service
From primary backup¶
Note: the roles aren't restored with this method.
With user postgres
:
# Optionally, rename the corrupt database (selecting a name for a database like "avgis_c2")... psql -c "ALTER DATABASE avgis RENAME TO avgis_c2;" # ... or drop the existing database psql -c "drop database avgis;" # Create a new database: createdb -E utf8 -T template0 avgis # Restore the database pg_restore -d avgis /var/backups/postgres/avgis.pg_dump
From secondary backup¶
autopostgresqlbackup
backs up the roles in postgres_globals
.
zcat /var/lib/autopostgresqlbackup/daily/postgres_globals/postgres_globals_2018-10-24_06h25m.Wednesday.sql.gz | psql zcat /var/lib/autopostgresqlbackup/daily/avgis/avgis_2018-10-24_06h25m.Wednesday.sql.gz | psql
Gear¶
Survey equipment¶
See survey equipment - to be written¶
Weather station¶
Ambient Weather weather station¶
We have purchased a WS2902A weather station (https://www.ambientweather.com/amws2902.html).
Firmware version: 4.0.2.
Manual¶
The operating manual of the weather station can be found at https://p10.secure.hostingprod.com/@site.ambientweatherstore.com/ssl/Manuals/WS-2902C.pdf
Connection¶
Wifi¶
Set up the wifi of the console using the "Ambient Tool" phone application. IP address given by DHCP on the router is: 192.168.1.101
Local communication¶
Fail so far: the only exposed port is TCP/45000. Telnet doesn't show any activity. Nothing found on Internet on this protocol.
One interesting project may come, hijacking the connection to cloud services: https://www.wxforum.net/index.php?topic=35033.0
Cloud connection¶
We'll create an account on AmbientWeather.net (and eventually on WUnderground.net and/or weathercloud.net), and:
- have the console upload data to there
- gisaf to retrieve our WS data from there
Notes on Console¶
The daily rainfall data displayed in the console resets at 00.30 every night
Plan for future¶
Beside living well,
Plan¶
Some interesting projects that might be integrated:
- https://github.com/Oslandia/albion : Build 3D geological model from wells information
Other¶
GDAL (OGR) tools¶
GDAL is a translator library for raster and vector geospatial data formats. It is used by many software (including QGIS and many other open source ones, including Gisaf). Some command line utilities are supplied, like:
ogr2ogr
can easily convert one data format to anotherogrinfo
displays information about files.
Using Windows¶
On a computer with Windows and GQIS installed:
1. Open a command line console (eg. <Windows Key> to display the Start menu, then just type cmd
and <Enter>)
2. In the console window, type (adjust with the QGIS version and location, this seems to be the standard one):
"c:\Program Files\QGis 3.10\OSGeo4W.bat"
3. GDAL utilities can be used: ogr2ogr
, etc.
Example: convert Geopackage to Shapefiles¶
Output the content of the geopackage 9wdoogfr_2019-11-13_12_26_07.gpkg
to the folder shapefiles
:
ogr2ogr -progress -f "ESRI Shapefile" shapefiles 9wdoogfr_2019-11-13_12_26_07.gpkg
To output the content of the geopackage 9wdoogfr_2019-11-13_12_26_07.gpkg
to the root
folder:
ogr2ogr -progress -f "ESRI Shapefile" c:\shapefiles 9wdoogfr_2019-11-13_12_26_07.gpkg
With reprojection¶
Same as above, reprojecting to UTM44N:
ogr2ogr -progress -f "ESRI Shapefile" -t_srs EPSG:32644 c:\shapefiles 9wdoogfr_2019-11-13_12_26_07.gpkg
Links and references¶
Links¶
Water management¶
Modflow¶
The reference software for underground water modelling and simulation.
In conjunction with flopy (https://water.usgs.gov/ogw/flopy/) and Jupyter (https://jupyter.org/), it provides a relatively easy to use interface.
Freewat¶
This project is partly based on modflow, and integrates with QGis.
QGIS¶
- Tools for Geology
Construction of geological cross sections in QGIS - http://www.geokincern.com/?p=1452
Autocad¶
- Overview of Converting Geospatial Data to Drawing Objects:
http://docs.autodesk.com/CIV3D/2013/ENU/index.html?url=filesMAPC3D/GUID-C38FD485-3CC2-4B52-8264-0D8C0F45422B.htm,topicNumber=MAPC3Dd30e41809
- CAD-DB connection:
https://knowledge.autodesk.com/support/autocad-civil-3d/learn-explore/caas/video/youtube/watch-v-AQoB--nyUJA.html
Orfeo¶
Remote sensing
Old docs¶
Shapefiles - to be written
Data (measurements auxiliary tables)¶
Besides the importation of shapefiles, Gisaf can import non-geophical information: auxiliary data (typically categories like the list of locations names, well types, etc), and temporal informations (well levels, etc).
Command line¶
The import_to_db.py
script imports files, fetched from a set of URLs (typically, in the Redmine Files section of this project), formats and pre-process, and imports to the database.
import_to_db.py
is a support tool, that is planned to be integrated with the web interface.
Import all with:
phil@phil-mbp:~/BlueLight/gisaf_src/gisaf$ python import_to_db.py
The script currently accepts an argument for filtering the URLs to import.
Pavneet's docs (imported from gisaf's wiki)¶
Basic Rules of Map making¶
Regardless of the cartographic style or content, most maps have the following common elements.
TITLE¶
The title should be in a large font, easily identifiable as the title of the map and should include descriptive text as to the location and purpose of the map. If the map is thematic, the theme should be included in the title. For example: Corn Production in Washington, 1990. The title is usually the largest font size of all lettering on the layout, however, it should not dominate the map graphic itself. The title may or may not be in a box and does not need to be at the top of the page (though it often is). For published materials (e.g., books or articles) the title may be included in a figure caption instead.
SCALE INDICATOR¶
The scale of the map is typically indicated by a graphic bar scale, a representative fraction or a verbal scale. The reader must be able to determine the relationship between a unit of measure on the map and a unit of measure in the real world.
ORIENTATION¶
A map should indicate which way is north (and/or south, east and west). Commonly this is done by a north arrow or compass rose. Orientation may also be shown by graticule or grid marks (e.g. lines of latitude and longitude). By convention north is towards the top of the page (thus some maps do not have north arrows), but the orientation must still be given for a 'proper' map. North does not have to be at the top of the page and a north arrow is essential in maps where it is not.
BORDER¶
A border identifies exactly where the mapped area stops. The border is often the thickest line on the map and should be close to the edges of the mapped area. The distance between the map and the border should be the same on all sides (balanced).
There can also be a border around the entire map layout (enclosing and grouping the title, legend, text boxes, etc.).
Both of these borders are sometimes referred to as a 'neatline.' In addition, there is sometimes a thin additional line just outside of a border (accentuating it and ideally making it more visually appealing) that may also be referred to as a neatline.
LEGEND¶
A legend defines the symbols or colors (including shades of gray and patterns) used on the map. Maps do not need legends if the symbology is so common or simple as to be easily understood by the reader. However, it must be clear what each marker or line type, weight and pattern represents. The legend does not need to be labeled "Legend." The more complicated the symbology on a map the more important the legend becomes.
MAP CREDITS¶
- SOURCE OF DATA (especially on thematic maps)
- NAME of the cartographer
- DATE of the map creation/publication
- DATE of the map data
- PROJECTION of the map (especially small-scale maps)
LOCATOR MAP (INSET)¶
A locator map is needed if the area of the map is not easily recognizable or is of large scale. For example, if you map Whatcom County, there should be an inset map of Washington, showing the location of Whatcom County. Inset DETAIL map(s) may also be used to show an area of the map in greater detail (larger scale).
EFFECTIVE GRAPHICAL DESIGN¶
The layout design is as important as effective sentence structure is to written text. Layout design refers to the planning and decision making processes involved in the visual display of the spatial data. You can achieve balance by rearranging the map elements (north arrow, legend, scale, title, etc.) and changing size of the text, border. etc. The map and map elements should be:
- Neatly drawn
- Appropriately and consistently generalized
- Symmetrically balanced (avoid crowding or large blank areas)
- Without unnecessary clutter (keep it simple, be wary of 'artistic' details)
VISUAL HIERARCHY¶
A hierarchy of symbology should be used for the lettering, line weights and shading. More important features are typically larger and/or darker, less important/background information should be smaller and/or lighter. At the same time, do not "over weight" or "under weight" features.
PURPOSE¶
All maps have a purpose which should influence every element of the map and the map layout. A cartographer should be able to clearly articulate the purpose of their map and should keep the audience (who the map is going to be used by) and the client (who the maps is being produced for) in mind.
NOTE: Any, or all, of the above 'rules' can be (and frequently have been) violated at the discretion of the cartographer IF doing so produces a better map (better serving its purpose and audience).¶
In general, with cartography, less is more (avoid excessive clutter).
GIS - Survey_Database¶
DB- Database
C3D- Civil 3D
How to create Survey database in Civil 3D¶
1. Setting up the Working Folder¶
A working folder needs to be created where the survey databases gets stored. In the Toolbox > Right click on Survey databases > Set the working folder - Save the folder in the desired location. We are setting the Working Folder for the Civil 3D databases in C:
2. Creating a database¶
Right click on Survey databases > New local Survey database - Enter the name> click ok. A survey DB is created (It is in bold which signifies it is the current DB under use). Multiple DB's can be created. A survey DB can be opened for editing by right click > open for edit. It can be closed by right click > close for edits. A survey DB can also be opened as read-only by same procedures.
3. Setting up the database¶
Survey Database(abc) > Right click + Survey Database Settings > Specify co-ordinate Zone > Distance- Metre > Temperature - Celsius > Pressure - Millibars > Distance type - Horizontal > Vertical type - Vertical Distance.
Then the next step is to create the Networks under which the point data gets uploaded. For example in our case we have - TS and RTK
Components of Survey Database¶
Survey database in civil 3D has the following components when expanded.- Import Events
- Survey Queries
- Networks
- Network Groups
- Figures
- Figure Groups
- Survey Points
- Survey point Groups
Survey- Field to Finish¶
Steps from field work¶
- ?? To be added by Raj and Ram
- Copy the .txt files from the equipment- Controller/Total station, to system using a pendrive (Storing the data- "D:"> "AVSM"> "Water Projects"> Respective project folder)> "TS"/"RTK")
- Cleaning up the .txt files/removing errors (If any)
Feeding the survey data into Civil 3D- by Surveyor¶
- After cleaning up of files, the files are stored in a desired location, to be used as import events into the Civil 3D file
- Open up project in Civil 3D
- Importing of Events (explain the step)
Processing- by Surveyor¶
Initial line work is generated from survey points in the drawing. Points for reconciliation (Changing the point codes) are reported, if any. It is important that there is a fair knowledge of Autodesk Civil 3D for this work.
Post-Processing¶
Final linework is generated on top of the processed linework. At this stage, a drawing is cleaned and prepared for sharing. problems like - overlaps in line-work, proper assigning of layers etc. Points for reconciliation are reported, if any. Post-Processing is crucial to generate correct and standardised survey drawings (in .dwg format) keeping in mind that those drawings will be used to generate shape files (.shp format). Therefore, it becomes important to follow a certain workflow based on compatibility of .dwg elements and .shp elements - for example for creating a shapefile of "Polygon" type, the elements in .dwg should be all "Polylines" of "closed" nature. All features should be in zero " 0 " elevation and the shapefile generated should be 3d type. This has been elaborated here < insert > At this stage, a good knowledge of Autodesk Civil 3D and GIS (to a certain extent) is a must.
Creating shapefiles in AutoCAD¶
?? add by pavneet - to be written
Sharing on WebGIS - "GISAF- gis.auroville.org.in"¶
After creating the shp files from AutoCAD, the shp files are ready to be shared on WebGIS platform. Following are the steps to do so.
- Upload and import Raw survey points to GISAF - Before uploading and importing the shapefiles, it is necessary to upload and import the corresponding Raw survey points into GISAF.
- Auto-import of point shapefiles - After the import of point files (RTK and TS), the next step is to Auto-import the points which are point type shp files. Go to Admin> Others> Projects> Select the corresponding project > With selected > Auto- import to GIS Database ( GISAF)
- Upload shp to basket - Zip together all the types of files obtained after saving a shapefile on the system (shp, prj, dbf etc). Follow the same naming standards when creating this. Upload it to the basket in GISAF- Admin> Basket> Shapefiles> click on relevant project> click on relevant sub-project (if any)> upload file> select the zipped file from system> select the category from directory >save.
- Import to GISAF (Only with Authorization) - Click on the import arrow.
Adding new codes¶
ADD !! - to be written
Reconciliation of points¶
The points are noted with their codes from and codes to reconcile into. ADD!! - to be written
Survey Data Post-Processing¶
Softwares used: Civil 3D, QGIS
AutoCAD Civil 3D - C3D,
C3D is being used for post processing of survey data. C3D offers a BIM solution for Land Surveyors. Basic knowledge about C3D can be obtained by following tutorials online though a pre-acquired knowledge of Autodesk AutoCAD is a plus point and sufficient enough to start working with C3D.
QGIS is an Open source GIS.
Getting Started:
Standards and Workflow for C3D
The layers (Nomenclature of layers is coming from standards- U.S. National CAD Standard Version-AIA) are pre-assigned a geometry type (Point, Line and Polygon) in their layers description (Use layer manager to check). This is done keeping in mind inter-operability with GIS (In shapefiles format), geometry type Shape files are of three kinds- points, lines and polygons.
Geometry types-
Based on the description of Layers, if the layer Geometry is specified as 'Point', the raw survey data can be auto imported to webGIS - GISAF, no post processing is needed for 'Points'. For lines and polygons, We use Polylines (for 2D lines/Polygons) and 3Dpolylines (for 3D lines/polygons). In case of curved lines in 2D, Polylines command (PLINE/PL) is enough but its not possible to create curved lines in 3D using the 3Dpolylines, in such cases feature lines come into the picture. (Important: 2D polylines, splines, ellipses and circles are not exportable to shapefile format, they have to be always converted into Polylines).
*We are at the moment generating shapefiles (lines and polygons) in 2D (z=0) but keeping the format of files as 3D.
Using Feature lines for creating Curved geometries in 3D-
There are a series of steps involved when creating curved lines using feature lines. You can follow two methods depending on the situation :-
Method 1(For curved objects in 3D using elevation of points):- Create the object using 3DPOLY. Now use the create feature lines from objects and select the objects to be converted to feature lines. You will not notice any difference in the geometry on the screen as such but in properties you will see the object type description of selected object as feature lines.
Method 2(For curved objects in 3D by draping the lines/curves on a surface):- Create the object using PLINE. You will notice all the lines are straight. Now use the create feature lines from objects and select the objects to be converted to feature lines. You will not notice any difference in the geometry on the screen as such but in properties you will see the object type description of selected object as feature lines.
Now modify tab is used to convert the straight lines to curved lines. There are many ways to do it and different methods can be adopted in different situations. The easiest and the most commonly used method is by using the smooth command in the modify tab. After using the smooth command, a curve would be visible. This is a curve in 3D. To view it, object viewer option can be used by selecting the object and right clicking.
More complex methods have to be studied and each case has to be taken into consideration separatly when using - FITCURVEFEATURE In some cases to obtain a desired curve(Meaning more fragments) more PI's (Point of insertion) need to be added into the feature line. PI are the points where the feature lines gets fragmented when exploded.This step often becomes difficult to manage due to complexities involved in mathematical functions behind creation of feature lines.But it should be explored by all means.
Because the feature lines are not exportable to shapefiles, the feature line has to converted back into a 3D polyline. This can be done simply by explolding the object. Once exploded into 3D polylines, you will see that the curve gets fragmented into smaller segments of straight lines. You can repeat the process of creating the 3Dpolyline object into feature line and then smoothening it till you get a desired geometry resembling the curve.
Situations of sharing of points with different layers and proceeding with linework in such circumstances-
There is always a best way to take the survey points taken in the field by the surveyor. For instance, in case of a road and a curb adjacent to it, the surveyor takes the points only once. It is only during the post processing that the lines are generated (You can create the lines in field while taking the survey points as well, but that is not the most feasible method majorly due to time constraints) and in this case, two line on the same position shall be generated, one under layer of roads and one under layer of curbs.
The line work becomes complex in case of generating 3d lines and 2d lines on the same place. This can arise due to having some layers as 2D (for example building outlines) and some as 3D (paved surfaces- <taking into consideration the future use of the layers. In this case, for the purpose of water management, it becomes essential to have survey of such surfaces as 3D. Survey of building outlines doesn't need to be 3D because the elevation points are not recorded while surveying them due to feasibility).
So, in case of generating 3D linework from survey points with no elevation, Feature lines and SURFACES come into picture. The surfaces are generated by an interpolation method in C3D by giving a set of Point Groups. We have selected the triangulation method as the appropriate method due to availability of a dense set of points.
While creating the feature lines from objects (as explained above, check the box "assign elevation" and a dialog box appears to select the surface you want to select). Surfaces will have to be generated prior to using this. (add wiki for generating surface). Surface should be generated including the points of the layer in the point group of it.
Purpose of Survey
- The survey is conducted for generating a Base Map which is representative of Topography, Infrastructure, Drainage, Buildings etc, with a focus on Water Management in Auroville. Furthermore this survey can be used as base map to build upon more detailed surveys for the various purposes such as Town planning, Land Surveys, Transportation Planning etc.
Use of Layers
- V-BLDG-SHED (Shed) layer is used when the structure is not entirely enclosed, if At least one side is open (No walls/partition etc). One thumb rule is to answer the question if the building can be locked safely or not. If not, then it is a shed.
- V-BLDG-HUT (Hut) is used when the building (For human activity/living space) is single storey with roof made up of perishable material such as keet, straw etc.
- Use of V-BLDG-ROOF (Roof) is for buildings with permanent roofs extending out about more than 1 meter from the outline at ground level. It can also be used for complex roof structures for example Matrimandir Petals.
- Use of V-BLDG-OTLN (Building outline) is for depicting the outer line of walls of a structure on ground level. This is used when the roof of the building is more or less of the same profile.
- V-BLDG-RTWL (Retaining Wall) is for walls with the function of retaining either earth/water and the side of the wall should be visible. In case of very thin walls and no side visible, it comes as V-BLDG-RWLL (Retaining wall line). This is depicted as a single line in drawing.
- Use of V-ROAD-CYCP is for designated cycle path as a single line.
- Use of V-ROAD-FPAT is for pedestrian paths as a single line. It is used when there is a path of very small width. If the need is to represent a path/road which is unpaved in nature, V-ROAD-UPVD is used, which is polygon type. To represent a path/road of paved nature V-ROAD-PAVD is used, which is also a polygon type.
- V-WMNG-PIPP for pipe points, V-WMNG-PIPL- for pipe lines.
- V-BLDG-BMRK is for unidentified bench marks and V-BLDG-MHOL is for unidentified manholes.
- V-BLDG-PLTF is for impervious surfaces. for example sand and gravel will not come under this layer. cemented paths etc can come.
- V-WMNG-POOL is for water bodies made for leisure activities- like swimming pool etc.
- V-WMNG-SUMP is collectively for artificial water retention structures.
- V-WATR-POND is for top of natural water retention.
- V-WATR-PONB is for bottom of natural water retention.
- V-WMNG-DRBA is for bottom of artificial drains.
- V-WMNG-DRTA is for top of artificial drains.
- V-WATR-DRBN is for bottom of Natural drains.
- V-WATR-DRTN is for top of Natural drains.
- V-WMNG-DRNC is for drain covers.
- V-WATR-DRAS is for indicating slope of drainage, these are arrows.
Wells Documentation¶
Documentation of the wells in Auroville started afresh in September 2017 by Bala working with the Water Group. He used a mobile GPS to record co-ordinates and took pictures along with other data related to a well. An effort was made to reconcile data by identification of the wells with existing data from Auroville Water Harvest which ceased to exist around 2007/2009. In some cases the codes on the well on location (sometimes there is a code on the pump, on the casing or on a nearby wall) helped in reconciliation but mostly by spatial mapping using QGIS. Some wells could not be reconciled due to absence of any nearby well in old data. However, we are publishing all the wells we have documented so far.
Terms of reference
Well status
- In use Well is equipped with functioning a pump
- Not in use Well is not equipped with a functioning pump. (Special case- well is equipped with a non functioning pump)
- Closed Well is closed/ sealed/ abandoned, (Historical reference)
Note: Wells with pump under repair (temporary measure) at the time of survey are treated as In use.
Civil 3D useful commands¶
Making Feature line to polyline-¶
Select the feature line> go to elevation editor (under feature line tab> edit elevations)> select all the points in table, give elevation(this will give same elevation to all the points thereby making it possible to retain the curves in the polyline) > explode the feature line (use X enter)
Converting Circle to polyline-¶
BR > enter > break the circle at two points to obtain a part of circle.
PEDIT> enter> select the part of leftover circle> J > enter> Close> enter.
Converting 2DPOLYLINE to POLYLINE¶
Explode the 2D polyline, and use the PEDIT command to convert the segments to polylines. Then join the polylines.
Viewing only the used layer in ACAD-¶
Set the value of SHOWLAYERUSAGE from 0 to 1
Convert 3d polyline to polyline-¶
COVERT3DPOLYS (Change the elevation of the polyline using properties manager to 0 in our case)
Generating contours from a Surface in Civil 3D¶
- In prospector, Create a point group of points needed to generate the contours from for example- TOPO elevation points or BLDG floor levels. If the point group is already present, then proceed to next step.
- Create a surface in the prospector> when creating the new surface, give the name and styling in the dialog box- styling used by us is "Triangles and surface- 0.1 and 0.5"> Assign the point groups in the surface as created in previous step.
- Go to the surface created in prospector and right click > Edit surface style> Turn on the Major and Minor contour.
- Turn on the layers for Major and Minor contours in Layers Manager (LA > enter)
Shortcuts:¶
• Properties manager - ctrl + 1
• To view 3D- Select and right click> object Viewer.
• To copy the Line - Copy and select the line enter
• To convert spline to polyline >Splinedit
• To export shapefile - Mapexport > follow the process.
• To assign coordinate system to the drawing - MAPCSLIBRARY
• To check coordinate system of the drawing- > TOOLSPACE> Settings> Right click on the drawing name > Edit settings
Online references for Civil 3D¶
Description Key Sets
https://www.youtube.com/watch?v=Von8oCwYcTk
Autodesk civil 3D Geotechnical module- for borehole data
https://www.youtube.com/watch?v=0Be9kShBou0
Civil 3D for Surveyors
https://www.youtube.com/playlist?list=PL1EzH8XlwSxuxCMdeLvzqfgsAlmeeHMv2
Civil 3D: Survey - Survey Database
https://www.youtube.com/watch?v=VPWdAfYJt5Y
Connections in QGIS- Using browser panel and Add postGIS¶
Working using QGIS as interface using PostGIS connections.
- Server:
gisdb.csr.av
- Database:
avgis
- Leave
Service
empty
Adding tables (With geometries) in form of shape files from database using Browser panel-¶
Go to QGIS> view> Panels> Browser panel> PostGIS> Expand the connection> give credentials> add the desired file by double clicking on it.
Adding tables (For non-geometry type) using PostGIS connections-¶
Add PostGIS Layers > Give credentials > select "Also list tables with no geometry" > expand public > click on the desired table> add.
Joining tables¶
Right click/double click on file> go to Join > perform the desired joins - add/subtract the joins.
Reconcilation of Raw survey data using pgAdmin¶
For changing the layer codes, to be done using Pgadmin.
Enter Schema> AVSM RAW Survey > Tables > select the table> Right click - View edit data > All rows. Apply filter to the original id and note the Database id's to be reconciled.¶
GISAf Admin> Other > Reconciliation. Create > Add the database (point) id and give the new target corresponding to the database id to be changed. (Use the Layer name in target not Raw layer name).¶
Other > Project > select the project > with selected> reconcile RAW survey points.¶
Importing point data (TS and RTK) to GISAF¶
- Gisaf Admin > Basket > Survey data > "Project" > "surveyor" > TS/ RTK > upload > import.
- Auto import of Raw points data (changes from Raw to Shapefiles, the point files)
In GISAF Admin > Other> Project > select the project for which you want to import the data > with selected > Auto import to GIS...
Editing Z value of features in Shapefiles in QGIS¶
using vertex editor tool - https://www.youtube.com/watch?v=8V8i1AtcA74&t=256s¶
Miscellaneous- Civil 3D¶
Autodesk civil 3D Geotechnical module- for borehole data :- To analyse borehole data, To make profiles and calculate volumes
https://www.youtube.com/watch?v=0Be9kShBou0
Reference videos
Exploring:
https://www.youtube.com/watch?v=gr-ISPzLcU0
Modelling:
https://www.youtube.com/watch?v=Pz0-HOoiBrs
solids
https://www.youtube.com/watch?v=mQ0Yeh6tZA8
https://www.youtube.com/watch?v=4eMsUiYBhuE
CIVIL 3D Survey
Getting started- Always open a new drawing with a template.
Description key sets- (till 18:30) https://www.youtube.com/watch?v=mmwkkRyBkS0
As when the points come into the drawing, they are going to be filtered with description key sets and the description key sets will assign properties to the points.
Tool space > settings > points > description key sets
“Description key sets name”> edit key- to view the points in a list:
Automatic linework (18:30- till end)
Survey > Linework code sets- for automatic linework
Using Master view
Uses-
copy styles from one drawing to other
copy drawing data from one drawing to other (data referencing)
Civil 3d surface model- from points- break lines
https://www.youtube.com/watch?v=wC40rdmDEMo
moving for example- point groups under point groups-
Civil 3D Planning and Analysis
Workspace: Planning and Analysis
1. Working with Object data (GIS Attributes)
Map Setup> Define object data> New table> Define new object data table start defining fields
Documentation- Rain Gauge¶
Manual Rain Gauge¶
Why it is important to comply to standards?
1. Consistency for comparing rainfall data in different places within Auroville since rainfall varies in different parts of Auroville.
2. Making it possible to use for any kind of scientific analysis by contribution to the primary data in sustainable water management for Auroville.
Proposed standards and ethics:
1. The time of taking the reading - 08:30 am.
The rain is recorded over a period of 24 hours- 8:30 of previous day till 8:30 of the present day, and the date is put as the present day. For example, if the rain is measured at 8:30 am on 10 dec 2017, then it is logged in as on 10 dec 2017.
2. Manual rain gauge typically used- green cylindrical and a collection jar.
3. Measuring jar: 10 mm corresponding jar
Area of the rim of manual rain gauge= 200 square cm
Diameter of rim= approximately 16 cm
NOTE: - If another kind of rain gauge is used, it should be used with the corresponding measuring jar as type of jar and calibration on jar depends on the area of rim. Any cases of different rain gauge other than above should be reported before use for recording purposes.
4. The minimum record-able unit is 0.2 mm (Least measure on the 10 mm cylinder). T is marked for below 1 mm.
5. The reading should be preferably submitted in an excel sheet (the data sheet will be provided). The frequency of sending data can be daily or weekly.
6. The rain gauge should be kept in a safe and open to sky area clear from tree cover or objects in surrounding so that direct rain falls into it and not from trees/ objects etc. Preferably roof top or a higher place on ground.
7. A check on the rain gauge should be kept when there is long period of gap in rainy days to make sure it was emptied out and there is no blockage.
8. In case of doubt in reading or incorrect measure, "incorrect measure" should be marked
Ethics
1. Readings should not be missed. Time to time check on rain gauge is required during long periods of no rain (especially when the rain season is approaching).
This is important because firstly, if it rained in the night and it wasn’t recorded, reading for the day is lost; Secondly, if the jar was not emptied, it gives incorrect subsequent reading.
2. If the person is not going to be available temporarily where the rain gauge is kept, the responsibility should be passed on to someone they can rely on after explaining the standards.
Important-
Data won’t be published in case of non-compliance to standards. This is to maintain the sanctity of proper scientific data collection and to keep it reliable for sharing by publishing. Discussions regarding standards are welcomed.
CSR Geomatics Team is placed on first floor in CSR, Auroshilpam.
AV rain data publication: (http://gis.auroville.org.in/measures/raingauge_av).
Adding Rain gauge to our web portal (GISAF), following information is needed:
1. GPS co-ordinates of location of placing the rain gauge (can be obtained using mobile phones easily)
2. Name, place of residence & contact number- Mobile & Landline
Automatic Rain Gauge¶
There is a possibility of publishing rain data coming from Automatic rain gauges like (id 15, auro orchard) and Weather stations.
For Automatic rain gauges, we can upload the files coming from the rain gauge directly into GISAF.
Thank you for contribution towards a sustainable management of Water in Auroville through Data collection.
Documentation- Wells Monitoring (Manual and Automatic / Piezometer by Bala)¶
Manual- Using tape with a sensor¶
Timings: The person assigned the job of monitoring collects the readings in three slots.
1. Between 6 am to 7 am till about 9 am to 10 am.
2. Between 11 am and 1 pm
3. Between 2 pm and 5 pm
Things required:
• Notebook
• Pen
• Measuring tape
• Vehicle for movement
Ethics-
• The monitoring should happen in coordination with the community members/ care taker (assigned by the person in charge from within the community). There should be a clear communication from monitor side as to which days the monitoring of a well has to happen and at around what time so that it can be made sure that the pump is not turned on before monitoring. In case a pump was turned on in a well, the monitor should have a gap of about 5 hours on the same day before going again for monitoring.
• The monitor should have contact numbers of the person in charge/ care taker for any communication.
• If the monitoring is stopped for any reason at any point, the monitor should communicate the same to the person concerned.
• The monitor shall take responsibility to inform any kind of changes in a well in terms of its functioning etc.
• The monitor should report to the CSR Geomatics Team who has responsibility to publish data.
• The monitor published the data to the website and works with the geomatics team.
Definitions:
A master file is maintained with the records of the wells from the field. It is updated when a new well is located. The following set of information are filled out in the Masterfile.
Well type
1. Open
2. Borewell
3. Dug cum borewell
Well status
1. In use
2. Not in use
3. Closed
Pump Status
1. Functioning
2. Not functioning
Pump Automation
1. Manual
2. Automatic
- Non-accessibility factors*
1. Dogs
2. Heavy slab
3. Narrow casing
4. enclosed/locked
5. permission
Data Matching Accuracy (reconnecting with harvest wells data)
1. High
2. Medium
3. Low
Well coordinates
For a new well, coordinates of the well position are taken on a mobile GPS.
Person In charge
For communication purpose, contact the person in charge as recorded in directory.
Automatic- Piezometer (by bala, to be edited)¶
The calibrations and setting are already done by Azha
The piezometer is taken to the site of unused well. First the depth and Water level is checked manually using water meter
Then the sensor of the piezometer is inserted into the well 1m above the bottom of the well, so the sensor does not get affected from mud or water inside
The sensor sends the data to an electronic board which translates the signal into the proper output, and then it sends the signal to the transmitting device above, on the ground surface. The transmitter then sends it to the receiver at CSR, the signal is then routed to Talam office through Interneet, enters the software and the reading is processed.
The piezometer sometimes does not have proper signal and so it would not be able to send it
Documentation- Flow meter by Bala¶
Flow meter – for checking the flow of water in pipe
CSR bought ultrasonic flow meter from Chennai and the company (company name? )trained Bala and Vijai(CSR) on how to set it up.
we tested in the west water system pump in csr. then we test in many more place .
ami, aurodam and buddha garden borewell pumps . we got request from the water service
Cross check they flow meters .
so before that we want to know how our meter works.
so we tested in our tank . the pump pumping from our sump tank so the flow will be stranded. when you pump in the bore well the flow goes up and down.
we tested two times in 15 minutes.and one time 30 minutes,so we got variation between this three.
we find 1,5 % error but the flow meret they said 1 % error only.
we cross checked the flow meter of the water service we find some error in they meter also.
we tried they bore well flowmeter and they sump tank flow meter also.and we find some error also in they meters.
then we find some error in (AVWS) meter also.
It is used to check the flow of water in pipe
tools used
First the outer diameter of the pipe has to be entered in the device. It is measured using Vernier caliper
Then the thickness of the pipe is set which is also measured using vernier
Also the device asks for the material of the pipe. If the pipe material is know it can be set and if it not know then there is an option which is other (mostly PVC, HDPE, and iron)
After entering these details the device gives the spacing for the sensors. There are 2 sensors up and down which has to be fixed accordingly
When the sensors are fixed the motor is turned on and the sensors send reading to the display device
This is noted down once every minute and taken for 15 to 30 minutes. Then the average of this is determined. This is done because of the variation in the flow. This gives the flow rate in 1 hour.
This was done in different places to check to flow rate
Documentation- DST- Vegetation Indexing¶
Steps for Dzetsaka Classification tool for Vegetation indexing in QGIS¶
1. Install the plugin Dzetsaka classfication tool.
2. Open the Raster from the Survey.
3. Create a polygon shapefile for index sampling. Mark polygons and give the ID's (1,2,3 for Tree, grass, bare land etc) Cover the variations in samples as much as possible.
More the samples, better the indexing.
4. Apply Dzetsaka Classification tool, Select the base raster and the sample- index polygon shapefile created in step 3.
5. The result is a Raster with DN numbers specified in the Shapefile in step 3.
6. Apply the Sieve raster command (Raster> Analysis> Sieve)- Try different threshold numbers and view the results till the noise is removed from the Raster.
7. Polygonise the Raster to Vector (From processing)
8. Run the v.generalise tool on the shapefile. This tool removes the pixelated boundaries of the polygons in the Vector.
Documentation- DST- Interpolation (Processing toolbox)¶
The following 4 tools have been mostly explored and the results were compared. The ones used for quick analysis are 1. Cubic Spline and 2. V.surf.spline . The rest of the tools are for further exploration and used depending on the need of the project. In some tools, the elevation values of points should be stored in the attribute table (Using field calculator and giving command - Z($Geometry) )
- Interpolate (Cubic spline) - SAGA
- V.surf.bspline - GRASS. Parameters to set - cell size. Set this parameter above 0.001 ( 0.00001, 0.000001 etc) and check the results.
- V.surf.rst - GRASS
- Krigging - SAGA
Documentation- DST- Survey- Office workflow - to be written
From CAD to GIS by Giulio¶
FEATURES CREATION
1. Assign a CRS to the drawing (TM-AUSPOS) (MAPCSLIBRARY command)
2. Create features in CAD (Points, lines, polygons)
3. Export shapefile (a) from CAD (Output > DWG to SDF) (Convert to LL84 – 3D)
FEATURES IMPORT INTO DB FIRST TIME
4. Create zip file of the shapefile
5. Upload into the GISAF Shapefiles Basket
6. Import the shapefile into DB
7. Save the shapefile on Local Machine
FEATURES IMPORT INTO DB EVERYTIME
8. Combine the new features to corresponding last shape files (Insert the process here).
9. Follow step 4-8 again
FEATURES EDITING IN QGIS
10. Open the table in QGis
11. Save as a shapefile (b) in TM AUSPOS CRS
12. In CAD, open a new drawing and assign AUSPOS CRS
13. Import the shapefile (b) (MapImport) with Object Data (Data tab > Create Object Data > OK), tick “Import polygons as closed polylines”, then press OK
14. Edit features
15. Change workspace into “Planning and analysis”
16. Export shapefile (a) from CAD (Output > DWG to SDF) with ONLY the id selected (Data Tab > Select Attributes > Object Data > Filename > id) (Convert to LL84 – 3D)
FEATURES IMPORT INTO DB
17. Create zip file of the shapefile
18. Upload into the GISAF Shapefiles Basket
19. Import the shapefile into DB
20. Delete the shapefile from Local Machine
QGIS- Miscellaneous¶
QGIS Introduction:
https://www.birdseyeviewgis.com/blog/2018/2/22/my-favorite-features-of-qgis-30to-date
3D visualization of raster DEM- https://www.youtube.com/watch?v=2KrCsbP0kUs
Spatial Query is selection of features that satisfies a certain condition which relates to other features in a space.
Using plugin- Spatial query
http://www.geodose.com/2018/03/spatial-query-in-qgis-3.html
Labelling:
Labelling with more than one field names and in different lines
'\n' |
Hierarchy of extensions
File levels and their uses.
https://nathanw.net/2014/03/22/all-the-q-files/
The Project file (.qgs)
It contains: Layer source pointer + Style information + Composers + a whole heap of other stuff
The Layer Definition file (.qlr)
It contains: Layer source pointer + Style information
The QML file (.qml)
It contains: Style information
Giulio's documentation¶
Documentation - Reconciliation of points using Gisaf¶
Reconciliation of points is a procedure used when a point is stored in a wrong table, because its category was either wrongly recorded in the field by the surveyor, or it has ben reviewed later by the surveyor or the data validator and found to be wrong.
Definition: Raw points are all points coming from the field survey. Raw points can be points referring to a Point feature (e.g. trees, or elevation points, or floor level), or points measured in the field to draw a line (e.g. the vertices of a fence) or a polygon (e.g. the corners of a building outline).¶
In the overall workflow, reconciliation affects only raw points. It takes place after the field textfile is uploaded into the basket and its points imported (raw points stored in the raw survey tables). Here you can visualize the Workflow diagram: Survey_data.
If a raw point refers to a point feature, reconciling it means moving it to another category/table meant for point features, not for lines/polygons. So, only raw points of point features can be reconciled. Raw points pertaining to line features and/or polygon features cannot be reconciled, so these raw points will remain in their original wrong table.
PLEASE NOTE:
The attributes of a line/polygon, which are derived from their defining raw points, will not be modified by any reconciliation, because only raw points of point features can be reconciled.
How to perform Reconciliation¶
To perform reconciliation of points: Login to Gisaf -> click on the G icon on the upper-left corner of the website page -> Manage -> Reconciliation by orig.ID
On the right end side of the screen, clicking on the field "Destination" a list of all categories will appear: these categories refer not to the Raw survey tables, but to the V_ tables of the database (points, lines, polygons).
Under it, in the field "Original ID", the original point number of the point to be reconciled is to be entered.
Clicking on the "Search points" button, the result shows the database unique id of the point, its survey category, its survey date, its geometry type (point, line, polygon), and the Project the point belongs to. In case of multiple points with the same original id (in case of different Projects, the field number of points might be not unique if the numbering of points in the field has restarted) all points having that original id are displayed: thanks to their date or Project or type, it is easy to identify the correct point to be reconciled.
Once the point to be reconciled is identified, clicking on the button "Reconcile" will run the reconciliation, and a message will appear stating that it has been done successfully. An error message can appear if a reconciliation of a raw point of a line/polygon feature has been attempted: this type of points cannot in fact be reconciled.
Once a raw point has been reconciled, it cannot be reconciled again. In case a raw point has been wrongly reconciled, it cannot be reconciled again through the above procedure, but it has to be reconciled manually through QGis or pgadmin software.
Documentation - Status and Status Changes¶
Status have been created to keep track og changes in surveyed features.
It is an additional (though provided for in AIA standards), single digit value, at the end of the Gisaf Category/Cad layer name.
Each Status need to have a corresponding CAD layer/Gisaf Category, with a short code associatedto it, so that field entries can be done easily.
Status have been defined as follows:
N - New Work
E - Existing o remain
D - Existing to demolish, Demolished or Changed
F- Future work, Proposed feature
T - Temporary work
M - Item to be moved
X - Not in contract
By default, Status is defined as E (Existing). Status can anyhow be changed later using gAdmin, or QGis (through the PostGis connection).
It needs to be done manually, one feature (point, line, polygon) at a time.
In the future Status changes might be incorporated in the Admin panel.
Documentation - Tags retained after re-import of same geometry¶
Documentation - Tags retained after re-import of same geometry
The linework for infrastructure survey carried out by Eric Chacra in May 2020 was originally imported with a problem of ambiguity in the "Accuracy" table and in the "Accuracy" table.
The result was that lines did not inherit the attributes survey date, accuracy, equipment, surveyor.
Nevertheless lines were displayed on the Gisaf map, without these attributes, and tags were given to some of these lines.
The values for the two tables ("Accuracy", "Accuracy") have been corrected, ambiguity resolved.
The lines in the layers V-ELEC-UGND------E and V-COMM-CABL------E have been then reimported, the attributes have been properly assigned, and the tags have been retained.
4 August 2020
Access to data¶
Connection to server directly from CSR¶
To connect to the server directly without going through Aurinoco server, the correct url is
http://gis.csr.av
Connection to Gisaf via QGis through WFS / OGC API¶
This works only on QGis from version 3.14.15 onward
In the browser, click on WFS/OGC API, then right-click to create a new connection
Give a name (e.g. OGC API Qgis Gisaf)
Give the url https://gis.auroville.org.in/ogcapi
Under the WFS Options box, on Version dropdown, the default option "Maximum" works just fine
Click on OK
The list of layers will appear in the Browser under WFS/OGC API.
How to create a new projection in QGis¶
To create a new projection in QGis, go to menu "Settings", and click on "Custom Projections".
A pop-up window appears with a list of all projections defined in QGis projects used by the user so far.
Click on the green "+" sign on the right top part of the window to create a new projection.
In the "Name" box, type "TM CSRAUSPOS SF1" (which means TM = Transverse Mercator projection; CSRAUSPOS = theparameters for this projection are derived from the processing of DGPS raw data by AUSPOS - Online GPS Processing Service - https://www.ga.gov.au/scientific-topics/positioning-navigation/geodesy/auspos; SF1 = Scale Factor is 1).
In the "Format" dropdown list, select "Proj String (legacy - Not Recommended)"
In the "Parameters" box, paste the following "+proj=tmerc +lat_0=12.01605433+lon_0=79.80998934 +k=1 +x_0=370455.630 +y_0=1328608.994 +ellps=WGS84+towgs84=0,0,0,0,0,0,0 +units=m +no_defs".
Finally, click on OK.
In a more explicit way, the parameters mean the following:
Map Projection: TransverseMercator (TM)
False Easting: 370455.6300
False Northing: 1328608.9940
Latitude of Origin: 12°00'57.79560" (DMS) 12.01605433 (DD)
Central Meridian: 79°48'35.96164" (DMS) 79.80998934 (DD)
Scale Factor: 1.00000000
Zone Width: 6.0°
Elimination of Duplicate points – General criteria¶
It might happen that the same physical feature (e.g. a tree, or a pole) is surveyed more than once: this can happen because there are many physical features in an area, and the survey needs more than one station. So, for example a tree is surveyed from a station, and gets a serial number on that date. When the station is then changed, it might happen that the same tree is resurveyed: another serial number is given, and possibly a different date, if the survey from the second station happened on a different day.
It is clear that the same tree is then represented with two different points, which means that two different trees exist: but only one tree really exist in the physical reality.
It is clear that one of the two points is redundant and needs to be removed. If this is noted by the surveyor directly in the field, then the issue is solved by the surveyor himself during processing time.
If instead, due to various reasons, it was not noted by the surveyor in the field, it will need to be cleaned after the processing, possibly by post-processing staff.
How to identify duplicate points?
The following criteria can be used:
1. The distance between the two points is less than 30 cm (trees are surveyed if their trunk diameter is at least about 20 cm, so in 30 cm cannot exist two of them)
2. The orig_id (serial number) of the points are not in series
3. The survey date is not the same
4. In case of trees, the species of trees is the same
5. 5. In case of trees, the tree type is not TDEF (because TDEF are mapped irrespective of their diameter, so they can actually have a small trunk, and two of them might exist in 30 cm), not OT (many TDEF species are surveyed as OT if not otherwise indicated by a botanist)
6. The context needs to be evaluated: if one tree is deleted in an area where many trees exist in a limited space, then loosing one in the map is not a big error. If instead one tree is deleted where there are very few trees, then it might be a big loss.
Linework for the Survey Area¶
1. Creation of Initial Linework in QGIS using Survey points import - (Ram, System 4)¶
Initial Linework in QGIS is started by surveyor with the knowledge from the Field. For this step, points are simply imported into the QGIS from the field text file (.csv or .txt). CRS needs to be TM-AUSPOS. The box of “First record has field names” shall not be ticked. In Point Coordinates, select the correct field for x, for y and for z (usually “field_2” for x, “field_3” for y and “field_4” for z). Points can be styled using the “Categorized” style in “Symbology”, using “Field_5” as value, or using a Rule-based symbology using the category (field 5) as filter.
Linework is created by connecting points having same description and belonging to the same physical feature. All line and polygon features are created as lines.
The Initial Linework for the Survey Area is also stored temporarily in
D: > AVSM > Zone-Survey number (eg RZ-01) > Survey Area (eg J) > Temporary WD
Note: The line shapefiles / Geopackages shall be in CRS: TM AUSPOS¶
2. Creation of final working drawing Shapefiles / Geopackages - (Selvarani, System 1)¶
Final working drawing Shapefiles / Geopackages are created from the Initial Linework of Survey Area.
As the Surveyor draws all features as lines (both for lines and polygons features), the following actions shall be done:
1. If features are lines:
• Export the shapefile / geopackage into the final working drawing folder (Final WD), in separate folders according to its type (e.g. BLDG, FENC, ROAD, etc).
The CRS for the export shall be EPSG:4326 - WGS 84¶
2. If features are polygons:
• Lines shall be converted into polygons:
to do it, first click on the layer to be converted to make it active (e.g. WD-CZ-01-F-LL84_V-BLDG-MHOL------E), then go to “Vector” Menu, click on Geometry Tools, click on Line to Polygons:
The new window for “Lines to Polygons” conversion will appear:
• Always cross check the input layer, to make sure that the input layer is the active one
• Save the output in a temporary layer
• The temporary layer will be listed in the list of layers, it shall be exported to the saving location as D: > Survey > Zone-Survey Number > Final WD > Survey Area SHP (eg . D: > Survey > GB-01 > Final WD > A-Shp)
The CRS for the export shall be EPSG:4326 - WGS 84¶
Once all the shapefiles / geopackages are exported in Final WD, for each of the newly exported layers the Topology Checker Tool shall be used.
Linework for the whole Survey Zone¶
1. Merging Shapefiles / Geopackages - (Selvarani, System 1)¶
A copy of the Zone Master shapefiles / geopackages are taken from System 4 and stored in Temp Folder on Desktop in System 1.
Master shapefiles / geopackages are merged with the Survey Area shapefiles / geopackages:
• To do it, go to “Vector” Menu, click on Geoprocessing Tools, then click on Union:
The new window for “Union” will appear:
• To make sure that the right geometry is generated by this process (“line” type, not “Multiline”, and similarly “Polygon” type, not “Multipolygon), we need to always keep the Master shapefile (e.g. Final-CZ-01-2021-02-05-LL84_V-BLDG-MHOL------E) as Input layer, and the Survey Area shapefile as Overlay Layer (e.g. WD-CZ-01-F-LL84_V-BLDG-MHOL------E).
• (The output can be saved to a file, as the CRS should already be EPSG4326 – WGS84.)
2. Storing Shapefiles / Geopackages - (Selvarani, System 1)¶
Save the merged shape file in the correct location in Final folder as D: > Survey > Zone-Survey Number > Final (eg . D: > Survey > GB-01 > Final)
Date in the name of Final Shapefile / Geopackage needs to be updated.
Once the merging operation is completed, the copy of Master shapefile / geopackage is deleted from the Temp folder.
3. Topology check of merged shapefiles¶
The topology checker is applied again on the merged shapefiles / geopackages.
The “id_field” shall be removed from the attribute table.
4. Archive and replace the Master Shapefiles / Geopackages (Ram, System 4)¶
Archive the previous master shapefiles / geopackages on system 4, and copy the new merged shapefiles / geopackages in its place.
*Then delete the Merged Shapefile / Geopackage folder from System 1.
*
5. Note about Shapefiles and Geopackages¶
All the above works are usually done using shapefile format, in QGIS latest version (3.16.3).
The Geopackage export is done in QGis versions older than 3.12 (e.g. 3.4, 3.6, 3.8, 3.10) so that the lines are not saved as “Multilines” but as “Lines”and polygons are not saved as “Multipolygons” but as “Polygons”. This is very important to be noted, as Gisaf database does not accept the Multipolygon and Multiline geometry types.
A different way to create "Polygons" is to use the command Vector -> Geometry Tools -> Multipart to Single Parts and apply it to the layer: from "Multipolygon" it will become "Polygon" (check in layer Properties).
As on 13 March 2021, Gisaf can accept Multipolygon layers, because the command "Multipart to Single parts" has been integrated into the importing command (see Redmine ticket #11691)
Creating 3D Shapefile/Geopackage¶
While creating a Shapefile/Geopackage, the File name/Database-table name, Geometry type and CRS have to be entered.
In order to create a 3D Shapefile/Geopackage, the additional dimensions "Z(+M values)"/ "Include Z dimension" has to be ticked: this way the 3D Shapefile/Geopackage is accepted by Gisaf without errors, otherwise the Shapefile/Geopackage can't be imported in Gisaf because the Z dimension is missing (the geometries in the database are all 3D).
Exporting from QGis (shapefiles and/or geopackages) to CAD dxf format¶
An algorithm has been created by Selvarani, to see the whole process click here:
https://redmine.auroville.org.in/projects/gis/wiki/Automatic_export_from_QGis_to_dxf