-
Notifications
You must be signed in to change notification settings - Fork 0
/
transectsurveydialog.py
149 lines (120 loc) · 7.24 KB
/
transectsurveydialog.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
# -*- coding: utf-8 -*-
import os
from PyQt5 import uic
from PyQt5 import QtWidgets
from PyQt5.QtCore import *
from PyQt5.QtWidgets import QDialog
from PyQt5.QtWidgets import QFileDialog
from qgis.core import *
from qgis.gui import *
from .surveyutils import *
from .transectsample import *
from .surveyproperties import *
from .xlsxwriter import *
FORM_CLASS = uic.loadUiType(os.path.join( os.path.dirname(__file__), 'transectsurveydialogbase.ui'))[0]
class TransectSurveyDialog( QtWidgets.QDialog, FORM_CLASS ):
def __init__(self, iface, parent=None):
super(TransectSurveyDialog, self).__init__(parent)
self.setupUi(self)
self.iface = iface
fillLayerComboBox( self.mSurveyAreaLayerComboBox, QgsWkbTypes.PolygonGeometry, True )
fillLayerComboBox( self.mStrataLayerComboBox, QgsWkbTypes.PolygonGeometry, False )
fillLayerComboBox( self.mSurveyBaselineLayerComboBox, QgsWkbTypes.LineGeometry, False )
self.strataLayerComboBoxChanged()
self.mCreateSampleButton.clicked.connect( self.createSample )
def strataLayerComboBoxChanged( self ):
comboIndex = self.mStrataLayerComboBox.currentIndex()
if comboIndex is None:
return
layerId = self.mStrataLayerComboBox.itemData( comboIndex )
layer = QgsProject.instance().mapLayer( layerId )
if layer is None:
return
fillAttributeComboBox( self.mMinimumDistanceAttributeComboBox, layer )
fillAttributeComboBox( self.mNSamplePointsComboBox, layer )
fillAttributeComboBox( self.mStrataIdAttributeComboBox, layer )
def createSample( self ):
#survey properties
surveyProps = SurveyProperties( self )
if surveyProps.exec_() == QDialog.Rejected:
return
#save dir
fileDialog = QFileDialog( self, QCoreApplication.translate( 'SurveyDesignDialog', 'Select output directory for result files' ) )
fileDialog.setFileMode( QFileDialog.Directory )
fileDialog.setOption( QFileDialog.ShowDirsOnly )
if fileDialog.exec_() != QDialog.Accepted:
return
saveDir = fileDialog.selectedFiles()[0]
outputPointShape = saveDir + '/transect_points.shp'
outputLineShape = saveDir + '/transect_lines.shp'
usedBaselineShape = saveDir + '/used_baselines.shp'
#strata map layer
strataMapLayer = self.stratumLayer()
if strataMapLayer is None:
return
strataMinDistance = self.mMinimumDistanceAttributeComboBox.currentText()
strataNSamplePoints = self.mNSamplePointsComboBox.currentText()
strataId = self.mStrataIdAttributeComboBox.currentText()
#baseline map layer
comboIndex = self.mSurveyBaselineLayerComboBox.currentIndex()
baselineLayerId = self.mSurveyBaselineLayerComboBox.itemData( comboIndex )
baselineMapLayer = QgsProject.instance().mapLayer( baselineLayerId )
#assume everything is in stratum units
minDistanceUnits = TransectDistanceUnits.StrataUnits
transectSample = TransectSample( strataMapLayer, strataId , strataMinDistance, strataNSamplePoints, minDistanceUnits, baselineMapLayer, True,
'', outputPointShape, outputLineShape, usedBaselineShape, self.mMinimumTransectLengthSpinBox.value(), -1.0, -1.0 )
transectSample.createSample()
#open line shape and create a dict station_id / bearing
#first find out 'station_co' attribute name
transectLineLayer = QgsVectorLayer( outputLineShape, 'transectLineLayer', 'ogr' )
fields = transectLineLayer.fields()
stationCodeFieldName = "station_co"
for field in fields:
if field.name().startswith( "station_co" ):
stationCodeFieldName = field.name()
break
bearingDict = {}
transectLineIt = transectLineLayer.getFeatures()
for transectLineFeature in transectLineIt:
station_co = transectLineFeature.attribute( stationCodeFieldName )
bearing = transectLineFeature.attribute( 'bearing' )
bearingDict[station_co] = bearing
#write bearing into output point shape
transectPointLayer = QgsVectorLayer( outputPointShape, 'transectPointLayer', 'ogr' )
transectPointLayer.startEditing()
transectPointLayer.addAttribute( QgsField( 'bearing', QVariant.Double, "Double" ) )
bearingIndex = transectPointLayer.fields().indexFromName( 'bearing' )
transectPointIt = transectPointLayer.getFeatures()
for transectPointFeature in transectPointIt:
station_co = transectPointFeature.attribute( stationCodeFieldName )
transectPointLayer.changeAttributeValue( transectPointFeature.id(), bearingIndex, bearingDict[ station_co] )
transectPointLayer.commitChanges()
#write gpx file
gpxFileInfo = QFileInfo( outputPointShape )
gpxFileName = gpxFileInfo.path() + '/' + gpxFileInfo.baseName() + '.gpx'
writePointShapeAsGPX( outputPointShape, stationCodeFieldName, 'bearing', gpxFileName )
transectLayer = QgsVectorLayer( outputLineShape, "transect", "ogr" )
#write XLSX output
workbook = Workbook( saveDir + '/survey.xlsx')
writeSurveyXLSX( workbook, surveyProps.survey(), surveyProps.projectCode(), surveyProps.date_s() , surveyProps.date_f(), surveyProps.contactName(), surveyProps.areas(), surveyProps.mainspp(), surveyProps.comments() )
writeStationXLSX( workbook, transectLayer, "stratum_id", "station_id", surveyProps.survey() )
writeStratumXLSX( workbook, self.stratumLayer(), self.mStrataIdAttributeComboBox.currentText(), surveyProps.survey() )
writeStratumBoundaryXLSX( workbook, self.stratumLayer(), self.mStrataIdAttributeComboBox.currentText(), surveyProps.survey() )
writeCatchXLSX( workbook )
writeLengthXLSX( workbook )
workbook.close()
#write csv files
#Survey.csv
writeSurveyCSV( fileDialog.selectedFiles()[0], surveyProps.survey(), surveyProps.projectCode(), surveyProps.date_s() , surveyProps.date_f(), surveyProps.contactName(), surveyProps.areas(), surveyProps.mainspp(), surveyProps.comments() )
writeStationCSV( saveDir, transectLayer, "stratum_id", "station_id", surveyProps.survey() )
writeStratumCSV( fileDialog.selectedFiles()[0], self.stratumLayer(), self.mStrataIdAttributeComboBox.currentText(), surveyProps.survey() )
writeStratumBoundaryCSV( fileDialog.selectedFiles()[0], self.stratumLayer(), self.mStrataIdAttributeComboBox.currentText(), surveyProps.survey() )
writeCatchCSV( saveDir )
writeLengthCSV( saveDir )
self.iface.addVectorLayer( outputLineShape, 'transects', 'ogr' )
self.iface.addVectorLayer( outputPointShape, 'transect_stations', 'ogr' )
def stratumLayer(self):
comboIndex = self.mStrataLayerComboBox.currentIndex()
strataLayerId = self.mStrataLayerComboBox.itemData( comboIndex )
strataMapLayer = QgsProject.instance().mapLayer( strataLayerId )
return strataMapLayer