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


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.