ich versuche bei meinem Raster Layer das Gebiet zu selektieren mit einem bestimmten Wert. Der raster Layer hat nur 4 Werte und ich möchte nur den 4. ansprechen um einen neuen Layer daraus zu erstellen.
das ist mein Rasterlayer
Code: Alles auswählen
str(geology)
Formal class 'RasterLayer' [package "raster"] with 12 slots
..@ file :Formal class '.RasterFile' [package "raster"] with 13 slots
.. .. ..@ name : chr "/Users/......./geologie.asc"
.. .. ..@ datanotation: chr "FLT4S"
.. .. ..@ byteorder : chr "little"
.. .. ..@ nodatavalue : num -9999
.. .. ..@ NAchanged : logi FALSE
.. .. ..@ nbands : int 1
.. .. ..@ bandorder : chr "BIL"
.. .. ..@ offset : int 6
.. .. ..@ toptobottom : logi TRUE
.. .. ..@ blockrows : int 0
.. .. ..@ blockcols : int 0
.. .. ..@ driver : chr "ascii"
.. .. ..@ open : logi FALSE
..@ data :Formal class '.SingleLayerData' [package "raster"] with 13 slots
.. .. ..@ values : logi(0)
.. .. ..@ offset : num 0
.. .. ..@ gain : num 1
.. .. ..@ inmemory : logi FALSE
.. .. ..@ fromdisk : logi TRUE
.. .. ..@ isfactor : logi FALSE
.. .. ..@ attributes: list()
.. .. ..@ haveminmax: logi FALSE
.. .. ..@ min : num Inf
.. .. ..@ max : num -Inf
.. .. ..@ band : int 1
.. .. ..@ unit : chr ""
.. .. ..@ names : chr ""
..@ legend :Formal class '.RasterLegend' [package "raster"] with 5 slots
.. .. ..@ type : chr(0)
.. .. ..@ values : logi(0)
.. .. ..@ color : logi(0)
.. .. ..@ names : logi(0)
.. .. ..@ colortable: logi(0)
..@ title : chr(0)
..@ extent :Formal class 'Extent' [package "raster"] with 4 slots
.. .. ..@ xmin: num 45806
.. .. ..@ xmax: num 71906
.. .. ..@ ymin: num 85296
.. .. ..@ ymax: num 101996
..@ rotated : logi FALSE
..@ rotation:Formal class '.Rotation' [package "raster"] with 2 slots
.. .. ..@ geotrans: num(0)
.. .. ..@ transfun:function ()
..@ ncols : int 522
..@ nrows : int 334
..@ crs :Formal class 'CRS' [package "sp"] with 1 slot
.. .. ..@ projargs: chr NA
..@ history : list()
..@ z : list()
Code: Alles auswählen
geology_MA<- geology[geology == 4]
Code: Alles auswählen
geology_MA<- geology[attributes==4]
Code: Alles auswählen
plot(geology_MA)
is.factor(geology)
x<-as.factor(geology)
factorValues(x , 4)
zzzz<<- geology[level=4]
plot(zzzz)
Hier ist der gesamte Layer und ich möchte wie gesagt nur den Teil haben, der den Wert 4 hat.
Ich danke im voraus für eine Antwort