আমার কাছে বিভিন্ন সিআরএসে বেশিরভাগ শেফফিল রয়েছে (বেশিরভাগ ডাব্লুজিএস ৮৮ ল্যাট / দীর্ঘ) যা আমি একটি সাধারণ প্রজেকশনে রূপান্তর করতে চাই (সম্ভবত অ্যালবার্স ইক্যুয়াল এরিয়া কনিক, তবে আমার সমস্যাটি আরও ভাল হয়ে গেলে আমি অন্য প্রশ্নে বাছাইয়ের জন্য সাহায্য চাইতে পারি) -defined)।
আমি আরে স্থানিক পরিসংখ্যান স্টাফ করতে কয়েক মাস ব্যয় করেছি, তবে এটি 5 বছর আগে। আমার জীবনের জন্য, আমি মনে করতে পারি না কীভাবে কোনও sp
বস্তুর (যেমন SpatialPolygonsDataFrame
) একটি প্রজেকশন থেকে অন্য প্রজেক্টে রূপান্তর করা যায়।
উদাহরণ কোড:
P4S.latlon <- CRS("+proj=longlat +datum=WGS84")
hrr.shp <- readShapePoly("HRR_Bdry"), verbose=TRUE, proj4string=P4S.latlon)
# Shapefile available at
# http://www.dartmouthatlas.org/downloads/geography/hrr_bdry.zip
# but you must rename all the filenames to have the same
# capitalization for it to work in R
এখন আমার কাছে SpatialPolygonsDataFrame
উপযুক্ত প্রজেকশন তথ্য রয়েছে তবে আমি এটি পছন্দসই অভিক্ষেত্রে রূপান্তর করতে চাই। আমি মনে করি এটির জন্য কিছুটা অনিচ্ছাকৃত-নামক কার্য রয়েছে, তবে এটি কী তা মনে করতে পারছি না।
মনে রাখবেন যে আমি কেবল সিআরএস পরিবর্তন করতে চাই না তবে সমন্বয়গুলি মেলে ("পুনরায় প্রজেক্ট", "রূপান্তর", ইত্যাদি) পরিবর্তন করতে চাই।
সম্পাদন করা
একে / এইচআই বাদ দিচ্ছেন যা এই শেফফিলের জন্য বিরক্তিকরভাবে মেক্সিকোয় স্থাপন করা হয়েছে:
library(taRifx.geo)
hrr.shp <-
subset(hrr.shp, !(grepl( "AK-" , hrr.shp@data$HRRCITY ) |
grepl( "HI-" , hrr.shp@data$HRRCITY )) )
proj4string(hrr.shp) <- P4S.latlon