আপনি এটি অনলাইন এপিআইয়ের পরিবর্তে স্থানীয়ভাবে করতে পারতেন। আর এর সাথে একটি সমাধান: মার্কিন যুক্তরাষ্ট্রের জন্য আদমশুমারির সীমানা তথ্য টাইগার এর মাধ্যমে আদমশুমারিতে পাওয়া যায় ; আপনি যদি মার্কিন যুক্তরাষ্ট্র সম্পর্কে চিন্তাভাবনা করেন তবে আপনি মার্কিন যুক্তরাষ্ট্রের রাজ্যগুলি ডাউনলোড করতে পারেন এবং ট্র্যাক্টগুলি পেতে কোনও ফাংশন ব্যবহার করতে পারেন। উদাহরণস্বরূপ, ক্যালিফোর্নিয়া এবং একটি এলোমেলো (বা র্যান্ডম-ইশ) পয়েন্টটি ব্যবহার করে:
library(maps)
library(maptools)
tractLookup <- function(x, y, state) {
pt <- SpatialPoints(data.frame(x = x, y = y))
overlay.pt <- overlay(pt, state) # what index number does pt fall inside?
return(census$TRACT[overlay.pt]) # give the Tract number from the census layer
}
চেষ্টা কর:
california <- readShapePoly("~/Downloads/US_2000_Census_Tracts/tr06_d00_shp/tr06_d00.shp")
tractLookup(-123.123, 40.789, california)
0004 দেয়, যা সঠিক।
# Look at the map
plot(census)
map('state', c('California'), lwd = 2, col = 'green', add = F) # optional
points(-123.123, 40.789, col = 'red', lwd = 2)
এটি কাজ করে, তবে আমার ম্যাকে 5 সেকেন্ড সময় লাগে, সম্ভবত আপনি যদি এটি অনেক কিছু করতে যাচ্ছেন তবে খুব বেশি হবে; আমার সন্দেহ হয় যে কেউ পোস্টজিআইএস সমাধানের সাথে শীঘ্রই ওজন করবে যা মিলিয়নগুণ দ্রুত হবে ...