আমি {এসপি} প্যাকেজটিতে 'রাস্টারাইজ' ফাংশনটি ব্যবহার করে একটি বিশ্ব বাথাইমেট্রি রাস্টার স্তর থেকে আমার আগ্রহের ক্ষেত্রের স্নানের মানচিত্রগুলি বের করার চেষ্টা করছি।
* সম্পাদনা: আমি 'এক্সট্রাক্ট' ফাংশনটি পেয়েছি যা মনে হচ্ছে যা আমি খুঁজছি তার চেয়ে বেশি।
আমি এ পর্যন্ত এটিই করেছি:
> class(subarea0) #This is my area of interest (Eastern Canadian Arctic Sea)
[1] "SpatialPolygons"
attr(,"package")
[1] "sp"
> extent(subarea0)
class : Extent
xmin : -82.21997
xmax : -57.21667
ymin : 60.2
ymax : 78.16666
library(marelac)
data("Bathymetry")#World bathymetric data in library (marelac)
names(Bathymetry);class(Bathymetry);str(Bathymetry)
[1] "x" "y" "z"
[1] "list"
List of 3
$ x: num [1:359] -180 -179 -178 -177 -176 ...
$ y: num [1:180] -89.5 -88.5 -87.5 -86.5 -85.5 ...
$ z: num [1:359, 1:180] 2853 2873 2873 2873 2873 ...
raster_bath<-raster(Bathymetry)#Transformed into a raster layer
extent(raster_bath) <- extent(subarea0)#Transform the extend of my raster to the extend of my SpatialPolygons
>ras_sub0<-rasterize(subarea0,raster_bath)#rasterize my SpatialPolygons (*Edits: not the function that I need here, but I am still interested to learn what results mean)
Found 2 region(s) and 10 polygon(s)
> plot(ras_sub0)
> plot(subarea0, add=TRUE)
> ras_sub0
class : RasterLayer
dimensions : 180, 359, 64620 (nrow, ncol, ncell)
resolution : 0.06964709, 0.0998148 (x, y)
extent : -82.21997, -57.21667, 60.2, 78.16666 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +ellps=WGS84
values : in memory
min value : 1
max value : 2
layer name : layer
আমি ফলাফল বুঝতে পারি না। আমি কেন আমার প্রতিটি বহুভুজের জন্য 2 টি রঙ পাচ্ছি? তাঁরা কি বোঝাতে চাইছেন?
আমি কীভাবে অবশেষে স্নানের গভীরতা কনট্যুর পেতে পারি? এটি কি আমার রেজোলিউশন বা মাত্রা পরিবর্তন করার সাথে কিছু করার?
* সম্পাদনা: ঠিক আছে, আমি এখন নিম্নলিখিতগুলি সম্পাদন করেছি:
v <- extract(raster_bath, subarea0)#Extract data from my Raster Layer for the locations of my SpatialPolygons
ভি একটি তালিকা এবং আমি কীভাবে / কী আকারে আমার স্থানীয় বহুভুজের সাথে এই তথ্যটি ফেরত দিতে চাই তা নিশ্চিত নই ...
ধন্যবাদ!