def __init__(self,iface):
QtGui.QDialog.__init__(self)
# Set up the user interface from QTDesigner.
self.ui = Ui_spatial_autocorrelation()
self.ui.setupUi(self)
self.iface = iface
self.dir = os.path.realpath(os.path.curdir)
self.layers = []
for i in range(self.iface.mapCanvas().layerCount()): #this for loop adds current layers; layercount: current map layer
layer = self.iface.mapCanvas().layer(i) #to dropdown menu #get the currently active layer
self.layers += [layer] #can not use i+=1, due to i is a string, not a list. In order to do calculations, set a empty list first.
#set a filter to select vector layers into combobox. It is necessary to set a filter to exclude the raster layers because whole pysal can not include any raster, and pysal can handle data by columns only in .dbf files.
if layer.type() ==layer.VectorLayer:
self.ui.activecombobox.addItem(layer.name()) #set dynamic labels in the combobox
elif layer.type == layer.RasterLayer:
pass
else:
pass
class spatial_autocorrelationDialog(QtGui.QDialog):
def __init__(self,iface):
QtGui.QDialog.__init__(self)
# Set up the user interface from QTDesigner.
self.ui = Ui_spatial_autocorrelation()
self.ui.setupUi(self)
self.iface = iface
self.dir = os.path.realpath(os.path.curdir)
self.layers = []
for i in range(self.iface.mapCanvas().layerCount()): #this for loop adds current layers; layercount: current map layer
layer = self.iface.mapCanvas().layer(i) #to dropdown menu #get the currently active layer
self.layers += [layer] #can not use i+=1, due to i is a string, not a list. In order to do calculations, set a empty list first.
#set a filter to select vector layers into combobox. It is necessary to set a filter to exclude the raster layers because whole pysal can not include any raster, and pysal can handle data by columns only in .dbf files.
if layer.type() ==layer.VectorLayer:
self.ui.activecombobox.addItem(layer.name()) #set dynamic labels in the combobox
elif layer.type == layer.RasterLayer:
pass
else:
pass
#for revised records: if not right operation, raise: QMessageBox.information(None,"Raster Scale","Layer is ok")
@pyqtSignature('') #prevents actions being handled twice
def on_inputshpbutton_clicked(self):
myFile = QFileDialog.getOpenFileName (self, "Select a shapefile","","*.shp")
self.ui.inputshpline.setText(myFile)
#create a new combobox(1)use pysal to open file(2) read in
#for saved shapefile to show columns
openfile=str(self.ui.inputshpline.text())
f=pysal.open(openfile)
opendbf=openfile[:-3] + "dbf" #open the same file only with dbf
f_dbf = pysal.open(opendbf)
self.fileheader=f_dbf.header #find columns, already in a list
for i in self.fileheader: #i is in a string
self.ui.selectcombobox.addItem(i)
@pyqtSignature('') #prevents actions being handled twice
def on_activecombobox_currentIndexChanged(self): #when selecting any shapefile from active layers, but do not know how to get the file path from active layers?
#for active layer to show columns
self.ui.activecombobox.currentIndexChanged(int)
QMessageBox.information(self,"Vector file","Layer is ok")
#self.ui.activecombobox.text()
openfile=str(self.ui.activecombobox.getPath().Text())
f=pysal.open(openfile)
opendbf=openfile[:-3] + "dbf" #open the same file only with dbf
f_dbf = pysal.open(opendbf)
self.fileheader=f_dbf.header #find columns, already in a list
for i in self.fileheader: #i is in a string
self.ui.selectcombobox.addItem(i)
@pyqtSignature('') #prevents actions being handled twice
def on_Inputweightsbutton_clicked(self):
myFile2 = QFileDialog.getOpenFileName (self, "Select a weights file","","*.gal;;*.gwt;;*.mat")
self.ui.Inputweightsline.setText(myFile2)
@pyqtSignature('') #prevents actions being handled twice
def on_Inputweightscreate_clicked(self):
dlg = WeightsDialog(self.iface)
dlg.show()
results = dlg.exec_()
if self.ui.Inputweightsline != None: #how to call variable from other files?
#self.ci=WeightsDialog(WeightsDialog.accept)
#myfile4=WeightsDialog.accept(savefile)
#self.ui.Inputweightsline.setText(myfile4)
pass
else:
pass
@pyqtSignature('') #prevents actions being handled twice
def on_outputbutton_clicked(self):
dlg = QFileDialog()
myFile3 = dlg.getSaveFileName(self, "Save Output", "", "comma_separatedfile(*.csv)")
myFile3 += dlg.selectedNameFilter()[0] #?
self.ui.outputline.setText(myFile3[0:-1])
###############################################################################################
#### ####
#### ####
#### This is the method we need to implement. When they click OK this method runs ####
#### ####
#### ####
###############################################################################################
def accept(self):
if self.ui.savedshpradio.isChecked(): #when selecting saved shp
openfile=str(self.ui.inputshpline.text()) #make a string of saved file
savefile = str(self.ui.outputline.text()) #this will be a string like "c:\output.(.csv)"
weightsfile=str(self.ui.Inputweightsline.text())
if self.ui.MoranIcheck.checkState(): #run moran's I value
np.random.seed(10) #? Is this step necessary?
f=pysal.open(openfile) #read a shp file, not need to read
w=pysal.open(weightsfile).read() #read a weights file
opendbf=openfile[:-3] + "dbf" #open the same file only with dbf
f_dbf = pysal.open(opendbf) #read the dbf attribute file
fileheader=f_dbf.header
#select a column and let it function
columnindex=self.ui.selectcombobox.currentText() #when select a column
y=np.array(f_dbf.by_col[columnindex]) #change into array, by_col function is only for dbf file
mi=pysal.Moran(y,w) #value of Moran's I
MI=mi.I #list
savestring=str(MI)
#savestring=','.join(str(n) for n in MI) #change from list to string for saving
output=pysal.open(savefile, 'w')
output.write(savestring)
output.close()
#should rewrite the following part, because they can not do multi-selection?
if self.ui.normalradioButton.isChecked(): #under the assumption of normal distribution
if self.ui.expectedcheckbox.checkState():
NE=mi.EI
savestring1=columnindex+'\n'+'Moran\'s I'+','+savestring+'\n'+'\n'+'Normality Assumption'+'\n'+'Expected Value'+','+str(NE)
output=pysal.open(savefile, 'w')
output.write(savestring1)
output.close()
#show in the screen?
#if self.ui.showcheck.checkState():
# dlg = WeightsDialog(self.iface)
# dlg.show()
# results = dlg.exec_()
#else:
# pass
else:
pass
if self.ui.variancecheckbox.checkState():
NV=mi.VI_norm
savestring2=columnindex+'\n'+'Moran\'s I'+','+savestring+'\n'+'\n'+'Normality Assumption'+'\n'+'Variance'+','+str(NV)
output=pysal.open(savefile, 'w')
output.write(savestring2)
output.close()
else:
pass
if self.ui.standardcheckbox.checkState():
NS=mi.seI_norm
savestring3=columnindex+'\n'+'Moran\'s I'+','+savestring+'\n'+'\n'+'Normality Assumption'+'\n'+'Standard Deviation'+','+str(NS)
output=pysal.open(savefile, 'w')
output.write(savestring3)
output.close()
else:
pass
if self.ui.Zcheckbox.checkState():
Nz=mi.z_norm
savestring4=columnindex+'\n'+'Moran\'s I'+','+savestring+'\n'+'\n'+'Normality Assumption'+'\n'+'z-value'+','+str(Nz)
output=pysal.open(savefile, 'w')
output.write(savestring4)
output.close()
else:
pass
if self.ui.Pcheckbox.checkState():
Np=mi.p_norm
savestring5=columnindex+'\n'+'Moran\'s I'+','+savestring+'\n'+'\n'+'Normality Assumption'+'\n'+'p-value'+','+str(Np)
output=pysal.open(savefile, 'w')
output.write(savestring5)
output.close()
else:
pass
elif self.ui.randomradiobutton.isChecked(): #under the assumption of random distribution
if self.ui.variancecheckbox.checkState():
RV=mi.VI_rand
savestring6=columnindex+'\n'+'Moran\'s I'+','+savestring+'\n'+'\n'+'Randomization Assumption'+'\n'+'Variance'+','+str(RV)
output=pysal.open(savefile, 'w')
output.write(savestring6)
output.close()
else:
pass
if self.ui.standardcheckbox.checkState():
RS=mi.seI_rand
savestring7=columnindex+'\n'+'Moran\'s I'+','+savestring+'\n'+'\n'+'Randomization Assumption'+'\n'+'Standard Deviation'+','+str(RS)
output=pysal.open(savefile, 'w')
output.write(savestring7)
output.close()
else:
pass
if self.ui.Zcheckbox.checkState():
Rz=mi.z_rand
savestring8=columnindex+'\n'+'Moran\'s I'+','+savestring+'\n'+'\n'+'Randomization Assumption'+'\n'+'z-value'+','+str(Rz)
output=pysal.open(savefile, 'w')
output.write(savestring8)
output.close()
else:
pass
if self.ui.Pcheckbox.checkState():
Rp=mi.p_rand
savestring9=columnindex+'\n'+'Moran\'s I'+','+savestring+'\n'+'\n'+'Randomization Assumption'+'\n'+'p-value'+','+str(Rp)
output=pysal.open(savefile, 'w')
output.write(savestring9)
output.close()
else:
pass
else:
pass
else:
return #without checking Moran's I, close directly
#differences between pass and return: "pass" just skip the code, while "return" will terminate the program
elif self.ui.activecombobox.isChecked(): #when selecting active shp and then import pysal
layer = self.layers[self.ui.activecombobox.currentIndex()] #select a shp layer
savefile = str(self.ui.outputline.text())
weightsfile=str(self.ui.Inputweightsline.text())
pass
if self.ui.MoranIcheck.checkState():
np.random.seed(10)
#f=pysal.open() #calculate Moran's I and other value, but do not know how to get the file path from active layers?
else:
return
#how to show all results in the screen by creating a new dialogue?
self.close() #close the dialog window
