ল্যাটের / লম্বা থেকে ইউটিএম-তে আর-তে পুনরায় প্রসেসকোটিং করা হচ্ছে?


13

বাফারকে কার্যকর করতে আমাকে এটি একটি ইউটিএম রূপান্তর করতে হবে।

wets<-readOGR(dsn=".",layer="shapefile")
r.raster <- raster()
extent(r.raster) <- extent(wets)
res(r.raster) <- 100 

wets.r <- rasterize(wet,r.raster)
plot(wets.r)
wetsbuf<-buffer(wets.r,width=500)

কোডের শেষ লাইনটি বাফার তৈরির সময় এটি এই সতর্কতা দেয়:

Warning message:  
In couldBeLonLat(x) :
  raster has a longitude/latitude CRS, but coordinates do not match that

এখানে তথ্য

  summary(wets.r)
          layer
 Min.        1
 1st Qu.     1
 Median      2
 3rd Qu.     9
 Max.       11
 NA's    52629

summary(wets)

  Object of class SpatialPolygonsDataFrame
Coordinates:
      min       max
 x  683705  714088.8
 y 4326266 4343768.0
 Is projected: TRUE 
 proj4string :
 [+proj=tmerc +lat_0=0 +lon_0=24 +k=0.9996 +x_0=500000 +y_0=0 +datum=GGRS87
 +units=m +no_defs +ellps=GRS80 +towgs84=-199.87,74.79,246.62]
 Data attributes:
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
    0.0     2.5     5.0     5.0     7.5    10.0 






 wets.r

class       : RasterLayer 
dimensions  : 175, 304, 53200  (nrow, ncol, ncell)
resolution  : 100, 100  (x, y)
extent      : 683705, 714105, 4326268, 4343768  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
data source : in memory
names       : layer 
values      : 1, 11  (min, max)
attributes  :
   ID FID
 from:  1   0
 to  : 11  10

বাফারটি সম্ভব হওয়ার জন্য আমাকে প্রক্ষেপণটি পরিবর্তন করতে হবে।


আপনি কি কোনও অনুমানিত সমন্বয় ব্যবস্থা বা কোনও ভৌগলিক সিএসে ডেটা করছেন?
— হারুন

প্রাথমিক ভেক্টর ডেটা আমার মনে হয় এমন একটি প্রজেক্টড কোঅর্ডিনেট সিস্টেমে।
— gsa

প্রস্তাবিত (যেমন ইউটিএম) বা ভৌগলিক (ল্যাট / লম্বা)?
— হারুন

আমি কীভাবে এটি চেক করতে জানি না আমি মনে করি ইউটিএম
— বুটি সম্পর্কে

স্থানাঙ্কগুলি কী এবং আপনি কোন অঞ্চলে (রাজ্য, শহর)?
— এড.হঙ্ক

উত্তর:


16

এইভাবে আপনি রাস্টার প্যাকেজটি ব্যবহার করে আর-তে কোনও রস্টারকে পুনঃপ্রজনিত করতে পারেন । এই উদাহরণস্বরূপ, ইনপুট জিওটিফটি একটি NAD83 ভৌগলিক স্থানাঙ্ক সিস্টেমে ছিল এবং আমি একটি এনএডি 83 ইউটিএম 15 প্রজেক্টড কোঅর্ডিনেট সিস্টেমটিতে পুনরায় প্রজেক্ট করি। আরজিডিএল দ্বারা ব্যবহৃত প্রজ 4 ফর্ম্যাট প্রজেকশনগুলির জন্য একটি ভাল রেফারেন্স স্পটিয়ালারিফারেন্স . org এ পাওয়া যাবে ।

library(raster)

# Create RasterLayer object
r <- raster('C:/temp/binary_nad83.tif')

# Define the Proj.4 spatial reference 
# http://spatialreference.org/ref/epsg/26915/proj4/
sr <- "+proj=utm +zone=15 +ellps=GRS80 +datum=NAD83 +units=m +no_defs" 

# Project Raster
projected_raster <- projectRaster(r, crs = sr)

# Write the RasterLayer to disk (See datatype documentation for other formats)
writeRaster(projected_raster, filename="C:/temp/binary_utm15.tif", datatype='INT1U', overwrite=TRUE)

উত্তরের জন্য ধন্যবাদ, যদিও এটি ফেরত দেয়: প্রজেক্টে ত্রুটি (প্রজেক্ট থেকে): এই রূপান্তরটি করতে পারে না এছাড়াও: সতর্কতা বার্তা: আরজিডাল :: কাঁচা ট্রান্সফর্ম (প্রজফ্রোম, প্রজেক্ট, নরো (এক্সআই), এক্সআই [, 1], xy [,: 218 প্রস্তাবিত পয়েন্ট (গুলি) সীমাবদ্ধ নয়
— gsa

আপনি কি একটি একক ব্যান্ড রাস্টার বা একটি মাল্টিব্যান্ড ব্যবহার করছেন? এই উদাহরণটি একটি একক ব্যান্ড রাস্টারের জন্য।
— হারুন

এটি কোনও ভেক্টর (জলাভূমি) স্যাটেলাইট চিত্র নয়, তবে আমি একক ব্যান্ডের সাথে যাব a
— gsa

আপনার রাস্টার / রাগাল প্যাকেজগুলি আপডেট করার চেষ্টা করুন: r-sig-geo.2731867.n2.nabble.com/…
— হারুন

@ জিএসএ, যদি এটি কাজ করে তবে উত্তরের উত্তরটি গ্রহণ করা এবং গ্রহণ করা ভাল ধারণা! অন্যথায়, আপনার মূল প্রশ্নটি সম্পাদনা করুন এবং স্পষ্ট করুন।
— সিম্বামাঙ্গু
আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.