আর-তে জিওরিফারেন্সযুক্ত তথ্য সহ একটি রাস্টার তৈরি করুন

Question

আমার পূর্ববর্তী প্রশ্নের উল্লেখ করে , আমি কীভাবে "জিওকোড" / আর - ব্যবহার করে একটি উত্পন্ন রাস্টার চিত্রকে একটি বিন্যাসে রূপান্তর করব যা মানচিত্রে একটি স্তর হিসাবে প্রদর্শিত হতে পারে? সংক্ষেপে, আমি বাণিজ্যিক স্প্যাটিয়ালকি ( স্ক্রিনশট ) সফ্টওয়্যারটির ফলাফলগুলি পুনরায় তৈরি করার লক্ষ্য নিয়ে যাচ্ছি ।

আমি মানচিত্রের জন্য ডেটা হোস্ট করার জন্য একটি জিও সার্ভার উদাহরণ ব্যবহার করছি, তবে আমি কীভাবে প্রয়োজনীয় প্রযোজনাটি অন্তর্ভুক্ত করতে পারি এবং তথ্যটিকে চিত্রের সাথে সমন্বয় করতে পারি তাই এটির সত্যিকারের ভৌগলিক অবস্থানের সাথে সঠিকভাবে সারিবদ্ধ হবে। স্থানিক ডেটা ব্যবহার করে ঘনত্বের হিটম্যাপ তৈরি করা হয়েছিল।

Comments

যেহেতু এটি কেবলমাত্র একটি প্রোগ্রামিং প্রশ্ন এটি আসলে কোনও বিষয় নয়। আপনি কি স্প্যাটস্যাট প্যাকেজ সম্পর্কিত তথ্য এবং শর্ট কোর্স পিডিএফটি দেখেছেন যে এটি সেখানে আচ্ছাদিত রয়েছে? আমি নিশ্চিত যে এটি এর আগে এসেছিল এবং একটি দ্রুত গুগল অনুসন্ধান এই প্রাসঙ্গিক আর-সহায়তা উত্তর নিয়ে আসে । আরও কিছু অনুসন্ধানের পাশাপাশি আরও সরাসরি সমাধান উপস্থিত থাকলে আমি অবাক হব না।

— অ্যান্ডি ডব্লিউ

আপনি উত্পন্ন রাস্টার চিত্রটি কীভাবে তৈরি করেছেন? এটি কি অন্য কোনও অ্যাপ্লিকেশন থেকে এসেছে, বা আর-এর মধ্যে তৈরি হয়েছে?

— সিম্বামাঙ্গু

এটি চিত্রের () ফাংশনটি ব্যবহার করে আর এর মধ্যে উত্পন্ন হয়।

— মিমো

এক্স এর অসম অনুভূমিক এবং উল্লম্ব রেজোলিউশন রয়েছে। এই জাতীয় ডেটা আর্ক-এস্কি ফর্ম্যাটে সংরক্ষণ করা যাবে না আমার এক্সটির 4684, 1950 মাত্রা রয়েছে, এটি রেজোলিউশনের বিষয়ে অভিযোগ করে। আমি যদিও এটি তৈরি করেছিলাম। বাগি

Answer 1

রাস্টার প্যাকেজ আপনাকে স্বেচ্ছাচারী রাস্টারগুলি তৈরি করতে দেয় (আকার, অভিক্ষেপ) বা ম্যাট্রিকের মতো বিদ্যমান বস্তুগুলি ব্যবহার করতে পারে; তারপরে আপনি অনুমান এবং এক্সটেন্টস নির্ধারণ করতে পারেন। আপনার image()ফাংশনটি ম্যাট্রিক্স যুক্তি নেয় যা সরাসরি ব্যবহার করা যেতে পারে:

library(raster)    
## Create a matrix with random data & use image()
xy <- matrix(rnorm(400),20,20)
image(xy)

# Turn the matrix into a raster
rast <- raster(xy)
# Give it lat/lon coords for 36-37°E, 3-2°S
extent(rast) <- c(36,37,-3,-2)
# ... and assign a projection
projection(rast) <- CRS("+proj=longlat +datum=WGS84")
plot(rast)

তারপরে আপনি যে writeRaster()কোনও সংখ্যক ফর্ম্যাট ব্যবহার করতে পারেন

writeRaster(rast, "~/myraster.asc", format = "ascii")

Answer 2

বিভিন্ন ধরণের বিভিন্ন উপায় রয়েছে তবে আপনি উইন্ডোজে ওপেন-সোর্স সরঞ্জামগুলি ব্যবহার করতে চান এমন ছুরিকাঘাতে নেওয়া 3 টি বিকল্প প্রস্তাব করতে পারে:

1. জিডিএল (ওপেন-সোর্স জিও-রাস্টার ম্যানিপুলেশন লাইব্রেরি) এ gdal_translate সরঞ্জাম ব্যবহার করুন । আপনার সম্ভবত -gcp pixel line easting northing elevationপ্যারামিটারের দরকার আছে যেখানে পিক্সেল এবং লাইনটি রাস্টার ( http://www.gdal.org/gdal_datamodel.html ) এবং এক্সপ্লোর পরিচালনা / ভৌগলিক স্থানাঙ্কের নূরিংয়ের উপর x / y । আপনার সাথে একটি প্রক্ষেপণও নির্দিষ্ট করতে হবে-a_srs srs_def
2. জিডিএল- এ আর-বাইন্ডিং ব্যবহার করুন এবং একই কাজ করুন
3. একটি সহজ "জিইউআই" উপায় হতে পারে এই টিউটোরিয়ালটি অনুসরণ করা এবং জিডিএল সহ কিউজিআইএস ব্যবহার করা। এর অসুবিধা হ'ল যদি আপনি ইতিমধ্যে না করে থাকেন তবে আপনাকে QGIS এবং GDAL উভয়ই সেট আপ করতে হবে।

আপনি যে কোনও রুট বেছে নেওয়ার আগে আপনার প্রয়োজনীয় প্রতিটি বিট সফটওয়্যারটির নির্ভরতাগুলি আমি মনোযোগ সহকারে লক্ষ্য করব।