R / में UTM को lat / lon से रेखापुंज को निरस्त कर रहा है?

बफर को क्रियाशील बनाने के लिए मुझे इसे UTM में बदलना होगा।

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

अनुमानित (जैसे UTM) या भौगोलिक (lat / lon)?

— हारून

मुझे नहीं पता कि यह कैसे जांचना है, मुझे लगता है कि UTM bu t i'm कुछ निश्चित नहीं है

— gsa

निर्देशांक और आप किस क्षेत्र (राज्य, शहर) में हैं?

— ed.hank

यह है कि आप R का उपयोग कैसे कर सकते हैं एक रेखापुंज में रेखापुंज पैकेज का उपयोग करके । इस उदाहरण में, इनपुट जियोटीफ़ एक एनएडी83 भौगोलिक समन्वय प्रणाली में था और मैं एक एनएडी 83 यूटीएम 15 अनुमानित समन्वय प्रणाली को फिर से प्रस्तुत करता हूं। Proj4 प्रारूप अनुमानों के लिए एक अच्छा संदर्भ, जिसका उपयोग RGDAL द्वारा किया जाता है, spatialreference.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)

— हारून
स्रोत

इस उत्तर के लिए धन्यवाद, हालांकि यह इसे लौटाता है: ProjectExtent में त्रुटि (से, प्रोजो): यह परिवर्तन नहीं कर सकता इसके अलावा: चेतावनी संदेश: rgdal में :: rawTransform (projfrom, projto, nrow (xy), xy [, 1], 1] xy [:: 218 अनुमानित बिंदु (s) नहीं परिमित

— gsa

क्या आप सिंगल बैंड रास्टर या मल्टीबैंड का उपयोग कर रहे हैं? यह उदाहरण एकल बैंड रेखापुंज के लिए है।

— हारून

यह एक सदिश (वेटलैंड्स) की एक रेखापुंज प्रक्रिया से आता है न कि एक उपग्रह चित्र के साथ ताकि मैं एकल बैंड के साथ जाऊं।

— जीएसए

: अपने रेखापुंज / rgal संकुल को अपडेट करके देखें r-sig-geo.2731867.n2.nabble.com/...

— हारून

@ जीएसए, अगर यह काम करता है तो उत्तर देना और जवाब स्वीकार करना एक अच्छा विचार है! अन्यथा, अपने मूल प्रश्न को संपादित करें और स्पष्ट करें।

— Simbamangu