আমি একটি অনিয়মিত ব্যবধানে পয়েন্ট ডেটাবেস থেকে রাস্টার ইমেজ তৈরি করার চেষ্টা করছিলাম। ডেটা দেখে মনে হচ্ছে-
> head(s100_ras)
x y z
1 267573.9 2633781 213.29545
2 262224.4 2633781 69.78261
3 263742.7 2633781 51.21951
4 259328.4 2633781 301.98413
5 264109.8 2633781 141.72414
6 255094.8 2633781 88.90244
আমি জাজের মধ্যে এই 'জেড' মানগুলি তৈরি করতে চাই যা আমি তৈরি করেছি
# Create a fine mesh grid
my_mesh=expand.grid(seq(min(s100_ras$Y),max(s100_ras$Y),l=100),
seq(min(s100_ras$X),max(s100_ras$X),l=100))
আমি আরও চাই যে জেড-মানগুলি সেই জাতীয় জাল পয়েন্টগুলির জন্য 'এনএ' হিসাবে অর্পণ করা হোক যা ডেটা পয়েন্টের বাইরে থাকে। জাল উপরের পয়েন্টগুলি এর মতো দেখাচ্ছে: https://drive.google.com/file/d/0B6GUNg-8d30vYzlwTkhvaHBFTnc/edit?usp=sarsing যখন আমি প্লট করি
plot(my_mesh)
points(s100_ras$Y, s100_ras$X, pch="*", col='blue')
সমস্যাটি হ'ল আমি নিশ্চিত যে এটি কীভাবে তৈরি করবেন, নিম্নলিখিত পদক্ষেপগুলি কার্যকর হয় না কারণ আমার জাল গ্রিড এবং ডেটা পয়েন্ট একই স্কেলের নয় !!
library(rgdal)
library(raster)
xyz<-cbind(my_mesh, s100_ras)
r <- rasterFromXYZ(xyz)
image(r)
আমি যদি ডেটা পয়েন্টগুলি (কোনও জাল ছাড়াই) ব্যবহার করে কোনও রাস্টার তৈরি করার চেষ্টা করি, তবে আমার ডেটা অনিয়মিতভাবে ব্যবধানযুক্ত হওয়ায় আর একটি ত্রুটি ছুঁড়েছে!
library(sp)
s100_ras <- data.frame(expand.grid(x = s100_ras$Y, y = s100_ras$X),
z = as.vector(s100_ras$mean))
coordinates(s100_ras) <- ~x+y
proj4string(s100_ras) <- CRS("+proj=utm +zone=46 +datum=WGS84")
gridded(s100_ras) = TRUE
suggested tolerance minimum: 0.916421
Error in points2grid(points, tolerance, round) :
dimension 1 : coordinate intervals are not constant
তদুপরি, আমি 'রাস্টার প্যাকেজ' এর 'রাস্টারাইজ' ফাংশন (অনিয়মিত গ্রিডগুলির জন্য) দিয়ে খেলার চেষ্টা করছিলাম, তবে এটির সাথে কোনও উপায় পেলাম না :( আমি কীভাবে ফাঁকে ফাঁকে ফাঁকে ফাঁকে নিয়মিত গ্রিড বানাতে পারি তা জানার জন্য, তবে প্রয়োজনে মৌলিকত্বের, আমি ইন্টারপোলেশনটি এড়াতে চাই id আইডিডাব্লু বা ক্রিগিং পদ্ধতি ছাড়াই অনিয়মিতভাবে ব্যবধানযুক্ত ডেটা পয়েন্টগুলির একটি রাস্টার তৈরি করা সম্ভব?
SpatialPixelsDataFrame
প্রস্তাবিত tolerance
যুক্তি (আপনার ক্ষেত্রে 0.916421) দিয়ে ব্যবহার করা।