ভৌগলিক সমন্বয় ব্যবস্থা আর তে রূপান্তর করা

Question

আমার ভৌগলিক স্থানাঙ্ক ব্যবস্থায় পয়েন্ট রয়েছে এবং আমি সেগুলি সুইস গ্রিডে (CH1903 +) রূপান্তর করতে চেয়েছিলাম।

নমুনা তথ্য:

 id       lon      lat
  2 7.173500 45.86880
  3 7.172540 45.86887
  4 7.171636 45.86924
  5 7.180180 45.87158
  6 7.178070 45.87014
  7 7.177229 45.86923  
  8 7.175240 45.86808
  9 7.181409 45.87177
  10 7.179299 45.87020

সম্মানিত ফলাফল:

id       E              N      
2     2579408.2431  1079721.1499
3     2579333.7158  1079729.1852
4     2579263.6502  1079770.1125
5     2579928.0358  1080028.4605
6     2579763.6471  1079868.9218
7     2579698.0756  1079767.9762
8     2579543.1019  1079640.6512
9     2580023.6226  1080049.2672
10    2579859.1889  1079875.2740 

Comments

অ্যারন, আমি একই লোক আমি সেখানে সঠিক উত্তর পাই না। তুমি আমাকে সাহায্য করবে? আমি সত্যিই অবাক হয়েছি আপনি কতটা বাছাই করছেন।

— টপডম্বিলি

@ টপ এটি চটজলদি নয়, এটি স্ট্যাক এক্সচেঞ্জ নীতি। ক্রস পোস্টিং সমস্ত ধরণের অসঙ্গতি এবং সমস্যা তৈরি করে। (আপনিও লক্ষ্য করেছেন যে ভুল ফোরামে পোস্ট করা দরকারী উত্তরগুলির চেয়ে কিছুটা কম পেতে পারে)) আপনার পোস্ট করা এসও সংস্করণ মুছুন।

— whuber

@ টপডম্বিলি, আমি কেবল উল্লেখ করতে চেয়েছিলাম যে হোবারের উত্তর অনুসারে ইনপুট মানগুলি ডাব্লুজিএস ৮৪-তে রয়েছে এবং ডেটাম ট্রান্সফর্মেশন প্লাস সিএইচ 903 + / এলভি 95 গ্রিডে প্রজেকশন সহ চলছে।

— মকনেডে

@ মেকনেডি সেই পর্যবেক্ষণের জন্য আপনাকে ধন্যবাদ। আমি এটি ব্যাখ্যা না করেই ছিলাম যে আমি আসল (ল্যাট, লম্বা) স্থানাঙ্কগুলি ডাব্লুজিএস ৮৪ (যে অনুমানটি কোডের একটি মন্তব্যে সমাহিত করা হয়েছে) ধরে নিয়েছি । যদি তা না হয় তবে proj4string(d)সেই অনুযায়ী মান পরিবর্তন করুন । আমার দৃষ্টি আকর্ষণটি প্রাথমিকভাবে ভ্রান্ত পূর্বদিকে এবং নোথিংয়ের পরামিতিগুলির প্রতি আকৃষ্ট হয়েছিল x0এবং y0কারণ ওয়েবে কিছু জনপ্রিয় উল্লেখ (যেমন কোডের প্রথম মন্তব্যে) তাদের সর্বাধিক উল্লেখযোগ্য অঙ্কগুলি ফেলেছে, যার ফলে কয়েক হাজার কিলোমিটার দ্বারা সমস্ত পয়েন্টগুলি স্থানচ্যুত করে :-)।

— হোয়াইট

@ শুভ, আউট রে: তথ্যসূত্র! আমি WGS 84 এ সেট করা ইনপুট সম্পর্কে আপনার মন্তব্য দেখেছি, তবে টপডম্বিলি এটি দেখেছে তা নিশ্চিত করতে চেয়েছিলাম।

— মেকনেডি

Answer 1

আরজিডাল প্যাকেজটি ব্যবহার করুন । কোন সিআরএস ব্যবহার করা উচিত তা নিয়ে একটি সমস্যা রয়েছে; আরজিডাল ইপিএসজি কোডটি স্বীকৃতি দেয় না। আপনাকে এখানে স্পষ্টভাবে প্যারামিটার সরবরাহ করতে হবে provide (আপাতদৃষ্টিতে এগুলি প্রায় অনুমান, তবে সেগুলি বেশ ভাল বলে মনে হচ্ছে They এগুলি মনে হয় মানগুলি প্রায় 0.1 মিটারের মধ্যে রয়েছে))

# References:
# http://lists.maptools.org/pipermail/proj/2001-September/000248.html (has typos)
# http://www.remotesensing.org/geotiff/proj_list/swiss_oblique_cylindrical.html
#
# Input coordinates.
#
x <- c(7.173500, 7.172540, 7.171636, 7.180180, 7.178070, 7.177229, 7.175240, 7.181409, 7.179299)
y <- c(45.86880, 45.86887, 45.86924, 45.87158, 45.87014, 45.86923, 45.86808, 45.87177, 45.87020)
#
# Define the coordinate systems.
#
library(rgdal)
d <- data.frame(lon=x, lat=y)
coordinates(d) <- c("lon", "lat")
proj4string(d) <- CRS("+init=epsg:4326") # WGS 84
CRS.new <- CRS("+proj=somerc +lat_0=46.9524056 +lon_0=7.43958333 +ellps=bessel +x_0=2600000 +y_0=1200000 +towgs84=674.374,15.056,405.346 +units=m +k_0=1 +no_defs")
# (@mdsumner points out that
#    CRS.new <- CRS("+init=epsg:2056")
# will work, and indeed it does. See http://spatialreference.org/ref/epsg/2056/proj4/.)
d.ch1903 <- spTransform(d, CRS.new)
#
# Plot the results.
#
par(mfrow=c(1,3))
plot.default(x,y, main="Raw data", cex.axis=.95)
plot(d, axes=TRUE, main="Original lat-lon", cex.axis=.95)
plot(d.ch1903, axes=TRUE, main="Projected", cex.axis=.95)
unclass(d.ch1903)

        lon        lat  

[1,] 2579408,24 1079721,15
[2] 2579333,69 1079729,18
[3] 2579263,65 1079770,55
[4] 2579928,04 1080028,46
[5] 2579763,65 1079868,92
[6] 2579698,00 1079767,97
[7] 2579543,10 1079640,65
[8] 2580023,55 1080049,26
[9 ,] 2579859.11 1079875.27

Comments

আমি জানতে চেয়েছিলাম যে দশমিক মান না থাকাকালীন রূপান্তর নির্ভুলতার ফলাফল কম হতে পারে যখন শ্রদ্ধেয় ফলাফলগুলিতে উপলব্ধ স্থানাঙ্কগুলি নিকটতম 10 ডিগ্রিতে থাকে, আমার অর্থ 2 দশমিক দশমিক।

— টপডম্বিলি

আমি আপনার মন্তব্য বুঝতে পারি না। আপনি কি জিজ্ঞাসা করছেন যে যখন মূল (ল্যাট, লম্বা) মানগুলি সীমিত নির্ভুলতার সাথে নির্দিষ্ট করা হয় তখন অনুমানিত স্থানাঙ্কগুলি কতটা নির্ভুল হয়? যদি তা হয় তবে আপনি এই উত্তরটি দরকারী বলে মনে করতে পারেন ।

— whuber

rgdal EPSG: 2056,

— এফডাব্লুআইডাব্লু

@ এমডি আপনাকে ধন্যবাদ! আমি একটি রেফারেন্স পেয়েছি যা জানিয়েছে এটি ইপিএসজি: 9814, কিন্তু আরজিডাল এটি স্বীকৃতি দেয়নি।

— whuber