Skip to content

Commit fca156e

Browse files
new version
Add labels to elements in results file, add relative distance between points
1 parent 15e7ca4 commit fca156e

1 file changed

Lines changed: 32 additions & 0 deletions

File tree

Lines changed: 32 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,32 @@
1+
library(rstudioapi)
2+
library(svDialogs)
3+
library(stringr)
4+
library(stringi)
5+
setwd(selectDirectory(caption="Select 0000 QNT folder"))
6+
NumberofSamples<-as.numeric(dlgInput("Enter number of samples")$res)
7+
NumberofElements<-as.numeric(dlgInput("Enter number of Elements")$res)
8+
StoreKraws<-data.frame(matrix(NA,NumberofSamples,NumberofElements))
9+
StorePositions<-data.frame(matrix(NA,NumberofSamples,4))
10+
for(i in 1:NumberofSamples){
11+
QuantResult<-read.table(paste(getwd(),"/Pos_",stri_pad(i,4,pad="0"),"/1.wt",sep=""),header=FALSE,sep="",fill=TRUE)
12+
FolderName<-QuantResult[which(QuantResult=='Project',arr.ind=T)[1,1],which(QuantResult=='Project',arr.ind=T)[1,2]+2]
13+
elementlabels<-"El No."
14+
for(ii in 1:NumberofElements){
15+
El1<-which(QuantResult[,1]==ii)
16+
elementlabels<-c(elementlabels,(QuantResult[El1,2]))
17+
rowdecimals<-str_extract(QuantResult[El1,],".*[0-9]\\...[0-9]")
18+
rowdecimals<-rowdecimals[!is.na(rowdecimals)]
19+
ExtractedKraw<-rowdecimals[length(rowdecimals)]
20+
StoreKraws[i,ii]<-ExtractedKraw
21+
}
22+
StorePositions[i,1]<-QuantResult[which(QuantResult=='X=',arr.ind=T)[1,1],which(QuantResult=='X=',arr.ind=T)[1,2]+1]
23+
StorePositions[i,2]<-QuantResult[which(QuantResult=='Y=',arr.ind=T)[1,1],which(QuantResult=='Y=',arr.ind=T)[1,2]+1]
24+
StorePositions[i,3]<-QuantResult[which(QuantResult=='Z=',arr.ind=T)[1,1],which(QuantResult=='Z=',arr.ind=T)[1,2]+1]
25+
xmove<-as.numeric(StorePositions[1,1])-as.numeric(StorePositions[i,1])
26+
ymove<-as.numeric(StorePositions[1,2])-as.numeric(StorePositions[i,2])
27+
StorePositions[i,4]<-sqrt((xmove^2)+(ymove^2))
28+
}
29+
colnames(StoreKraws)<-elementlabels[2:length(elementlabels)]
30+
colnames(StorePositions)<-c("X","Y","Z","Rel Dist")
31+
#write.csv(StoreKraws,file=file.path(getwd(),paste(FolderName,"RStudio_Extract_Kraws_Points.csv",sep="_")))
32+
write.csv(cbind(StoreKraws,StorePositions),file=file.path(getwd(),paste(FolderName,"RStudio_Extract_Kraws_Points_RelDist.csv",sep="_")))

0 commit comments

Comments
 (0)