Wie visualisieren Sie eine Karte aus einer netcdf-Datei?

Habe ich eine netcdf-Datei möchte ich nur die Visualisierung der bodentiefe Karte

   [1] "file C:\\Users\\SoilDepth-gswp.nc has 3 dimensions:"
     [1] "x   Size: 360"
     [1] "y   Size: 150"
     [1] "land   Size: 15238"
     [1] "------------------------"
     [1] "file C:\\SoilDepth-gswp.nc has 3 variables:"
     [1] "float nav_lon[x,y]  Longname:Longitude Missval:1e+30"
     [1] "float nav_lat[x,y]  Longname:Latitude Missval:1e+30"
     [1] "float SoilDepth[land]  Longname:Soil depth Missval:1.00000002004088e+20"

Es scheint, dass ich haben, um die Verbindung der breiten mit Längen-sowie der land Punkte um eine Karte von der bodentiefe.Ich bin wirklich verwirrt .Kann mir jemand helfen mit dieser Art von Daten.

Die Größe des Rasters ist (360 * 150 = 54e3), während die Größe Ihrer land variable 15238, die nicht ein Vielfaches von Ihr gridsize. Haben Sie eine Erklärung dafür?
Wo hast du diese Daten? Oder hast du ihn selbst erstellen? Der Unterschied kann verursacht werden durch die Anwesenheit von NA-Werte, ich.e es gibt keinen Boden, die Tiefe in den Ozean.
land ist nicht eine variable, sondern eine andere dimension, so dass es nicht zu der Größe x * y
du hast Recht. Ich bin mir nicht ganz sicher, was in die Datei, aber ich denke, das ist eine Datei in einem bestimmten format, die von einigen Forscher, und, dass Sie müssen einige sein/Ihr code in der Lage sein zu interpretieren, den Inhalt der Datei korrekt. Ich würde versuchen, Kontakt mit dem Ersteller der Datei.
Die Abstimmung zu schließen, als auch lokalisiert, dies ist ein problem, das sehr spezifisch für diese bestimmte Datei.

InformationsquelleAutor sacvf | 2012-09-05

10
```
library(ncdf)
# I'm assuming this is the netcdf file you are working with:
download.file("http://dods.ipsl.jussieu.fr/gswp/Fixed/SoilDepth.nc", destfile="SoilDepth.nc")
soil <- open.ncdf("SoilDepth.nc")
#The way to extract a variable is as following:
soil$var[[3]] -> var3 # here, as shown in your question, SoilDepth is the 3rd variable
get.var.ncdf(soil, var3) -> SoilDepth

dim(SoilDepth)
[1] 15238
```
Wie gesagt wurde in der Zusammenfassung für Ihre netcdf-Datei, die variable SoilDepth hängt dimension land nur und nicht auf x und y also ich bin mir nicht sicher, wo bleiben Sie, wenn es um das Plotten dieses dataset.

Bearbeiten

Es stellt sich heraus, es ist ein Schlüssel, der links x, y und land:
```
download.file("http://dods.ipsl.jussieu.fr/gswp/Fixed/landmask_gswp.nc", destfile="landmask.nc")
landmask <- open.ncdf("landmask.nc")
landmask$var[[3]] -> varland
get.var.ncdf(landmask, varland) -> land
sum(land==1)
[1] 15238
```
So, in Reihenfolge zu zeichnen:
```
# The values where stored in an expected order, hence the transpose
land = t(land)
land[land==1] <- SoilDepth
land[land==0] <- NA
land = t(land)
image(land)
```
- Nur zum telefonieren Handlung auf das Ergebnis des get.var.ncdf(nc, "SoilDepth") auch nicht die Ausbeute eines schönen timeseries. Seltsame Datei...
- Es ist wahrscheinlich eine implizite Beziehung zwischen den Dimensionen land und Abmessungen x und y aber ich kann nicht scheinen, um es zu finden...
- Ich habe auch eine harte Zeit der Suche eine Beziehung. Die OP muss, bitten Sie den Ersteller der Daten...
- Ok, vielleicht ist diese Datei hier: dods.ipsl.jussieu.fr/gswp/Fixed/landmask_gswp.nc ist der Schlüssel.
- +1 (was ich bereits Tat), die für das Graben, um die Lösung zu erhalten.
- Die Daten, die gelesen wurde fälschlicherweise die Reihenfolge, in SoilDepth falsch war. Erste Umsetzung, als die Zuordnung SoilDepth und der Umsetzung wieder das problem gelöst.
- Danke, hatte nicht gesehen, dass. Ich bearbeitet habe, das Grundstück entsprechend.
- Wenn Sie haben ein neues problem, poste bitte eine neue Frage.
- die OP geschrieben, dass das Ergebnis falsch war, wie ein Kommentar auf die Antwort von AllenR.
- Millionen Dank für Euch alle
- Ich fügte hinzu, diese codes zu extrahieren Sie die Datei und schreiben in eine binäre Datei mit"""""" :schreiben =file(paste("C:\\Users\\aalyaari\\Desktop\\landmask_gswp.bin",sep=""),"wb")""""""" writeBin(als.Doppelzimmer(land),zu.schreiben,size=4)""""" veranschaulichen, wurde die Karte umgedreht!!! – jeder Vorschlag, warum. wie die Vergangenheit einer Karte zusammen mit der Frage hier?Ich konnte nicht
- Wenn die Eingabe einer Frage oder einer Antwort ist ein kleines icon (für eine polaroid oder so) müssen Sie klicken, um ein Bild (Sie benötigen mindestens 10 Rufpunkte zu dürfen, was der Fall ist). In Anbetracht der neuen Ausgabe, die Sie gerade haben, scheint zu sein, sehr Verschieden von Ihrer ersten Ausgabe an, sollten Sie veröffentlichen es als neue Frage statt der Bearbeitung deiner ursprünglichen Frage.
InformationsquelleAutor plannapus
11

Ich lieber mit den ggplot2 - Paket für die Visualisierung. Mit der ausgezeichneten Lösung von @plannapus:
```
require(reshape)
require(ggplot2); theme_set(theme_bw())
land_df = melt(land)
ggplot(aes(x = X1, y = X2, fill = value), data = land_df) + 
  geom_raster() + coord_equal() + 
  scale_fill_continuous(na.value = "transparent")
```
Wenn Sie möchten, ändern Sie den Titel einer Achse, tun nicht ändern Sie den Namen der Variablen in aes. Diese Werte beziehen sich auf die Spalten der Daten, und ändern Sie führt zu der Fehlermeldung, die Sie bekommen, gibt es keine Achse mit Namen X im land_df. Wenn Sie möchten, ändern Sie den Namen auf der Achse platziert:
```
ggplot(aes(x = X1, y = X2, fill = value), data = land_df) + 
  geom_raster() + coord_equal() + 
  scale_fill_continuous(na.value = "transparent") + 
  scale_x_continuous("X")
```
- Lesen ?scale_fill_continuous oder had.co.nz/ggplot2/scale_continuous.html
- Wenn Sie sich die Mühe, Sie sollte in der Lage sein, um Sie zu entfernen, indem Sie tun: +theme(axis.title=element_blank())
InformationsquelleAutor Paul Hiemstra
5

Tun, die Sie visualisieren möchten, um es mit R ?

Wenn es nicht ein problem zu visualisieren, die mit einer anderen software, die Sie verwenden können, ncBrowse, verfügbar hier, oder Palette, eine komplexere, zur Verfügung gestellt von der NASA, die können Sie donwload hier.

Wenn Sie wollen, um die Arbeit mit R, die Sie verwenden können ncdf Paket. Sie werden in der Lage, extrahieren Sie Ihre Daten Dank der get.var.ncdfFunktion. Sie können zeichnen Sie es Dank der sp Paket und spplot - Funktion, oder verwenden Sie die rglPaket (oder sonst scatterplot).
- Ich habe einige code, um dieses problem zu beheben, auf die Antwort von @plannapus
InformationsquelleAutor Pop
4

Für die schnelle Suche in Dateien, die Sie verwenden können, ncview. Die Karten sind nicht besonders schön, aber sehr funktionell für das erhalten einer Idee von, was eine bestimmte Datei aussieht. Auch dies funktioniert problemlos über remote-Server.

Siehe hier: http://meteora.ucsd.edu/~pierce/ncview_home_page.html

InformationsquelleAutor tbc

Schreibe einen Kommentar

Du musst angemeldet sein, um einen Kommentar abzugeben.