2012-09-26 20 views
5

मेरे पास एक डेटा सेट है जिसमें एक .csv फ़ाइल में 3 कॉलम होते हैं। पहले 2 कॉलम नक्शा समन्वय होते हैं और तीसरा संबंधित नक्शा समन्वय में बोरहेल में पाए जाने वाले जस्ता का प्रतिशत होता है। मैं दूरी के साथ जेएन एकाग्रता परिवर्तन दिखाने के लिए एक समोच्च मानचित्र बनाना चाहता हूं। कोड के सभी उदाहरण मैं मैट्रिक्स के रूप में उपयोग डेटा ढूंढने में सक्षम हूं, जबकि मेरा सूचियों में है। मैंने इसे साकार करने के कई अलग-अलग तरीकों की कोशिश की है जो मैंने नीचे दिखाया है, अधिकांश तरीकों से मुझे "ऑब्जेक्ट एक्स नहीं मिला" की रेखाओं के साथ त्रुटि संदेश मिलते हैं, जो मुझे लगता है कि मेरे डेटा के लेआउट के साथ करना है। क्या किसी को भी यह करना आता है? मैंने नीचे मेरा एक समान डेटा सेट जोड़ा है। किसी भी सहायता के लिए अग्रिम धन्यवाद। होलीसूचियों में डेटा का उपयोग कर एक भरे समोच्च साजिश का निर्माण

डेटा सेट:

Statsrep <- structure(list(X = c(156000L, 156010L, 156421L, 156450L, 156500L, 156700L, 158420L, 158646L, 158970L, 159050L, 159050L, 159130L, 159155L), Y = c(143630, 143980, 147260, 145000, 146000, 142800, 146700, 145207, 147170, 145200, 144800, 147815, 145890), Zn = c(2, 8, 4, 0, 3, 0, 2, 7, 12, 0, 4, 19, 0)), .Names = c("X", "Y", "Zn"), row.names = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L), class = "data.frame") 
Statsrep 

कोड:

library(ggplot2) 
Grade <- read.csv(file="filename.csv", header=TRUE, sep=",") 
ggplot(Grade, aes(x$x="X", y$y="Y", z$z="Zn")) + 
stat_contour() 

library(lattice) 
Grade <- read.csv(file="filename.csv", header=TRUE, sep=",") 
levelplot(Grade ~x*y, data = Zn, 
xlab = "Eastings", ylab = "Northings", 
col.regions = terrain.colours) 

Grade <- read.csv(file="filename.csv", header=TRUE, sep=",") 
x$x <- X 
y$y <- Y 
z$z <- Zn 
filled.contour(x$x, y$y, z$z, color = terrain.colours, 
xlab = "Eastings", ylab = "Northings"), 
plot.axes = {axis(1, seq(156000, 165000, by=1000)); axis(2, seq(142000, 150000, by=1000))}, 
key.title = title(main="Zn content\n(percent)"), 
key.axes= axis(4, seq(0, 20, by = 2))) 

उत्तर

8

ggplot2 के साथ कार्य करना है, तो आप का उपयोग कर सेट अपने उदाहरण डेटा के साथ एक समोच्च साजिश बना सकते हैं:

ggplot(Statsrep, aes(x=X, y=Y, z=Zn)) + 
    geom_density2d() 

को

दें

enter image description here

आपको ggplot2 कोड के साथ कुछ समस्याएं थीं। विशेष रूप से, जहां आपने सौंदर्यशास्त्र सेट किया था:

aes(x$x="X", y$y="Y", z$z="Zn") 
+0

आह इसलिए मैं गलत कार्य का उपयोग कर रहा था। यह काम करता है लेकिन रंगों के साथ लाइनों के बीच भरना संभव है? –

+0

आज़माएं: 'g + stat_density2d (aes (fill = .. घनत्व ..), समोच्च = गलत, geom =" टाइल ")' – csgillespie

+0

धन्यवाद! जब दोनों संयुक्त होते हैं जो निश्चित रूप से साजिश का प्रकार है जिसे मैं ढूंढ रहा हूं। हालांकि, घनत्व की कुंजी जस्ता का मूल्य नहीं देती है, यह ई -08 –

 संबंधित मुद्दे

  • कोई संबंधित समस्या नहीं^_^