PROSAGA码农传奇-时空索引-数据帧到R中的ArcGIS ASCII

数据帧到R中的ArcGIS ASCII

作者: 筱梨
发布时间: 2024-08-16 09:12:20 (3天前)
转自：

<div class =“post-text”itemprop =“text”>
  <P>
    检查包裹
     <code>
      raster
    </code>
    。
  </p>
  <P>
    首先需要将数据帧转换为栅格文件（重新排列列以匹配xyz）：
  </p>
   <pre>
    <code>
      r <- rasterFromXYZ(predHB.Lv1Hard[,c(2,1,3)])

</code>
  </pre>
  <P>
    然后使用：
  </p>
   <pre>
    <code>
      writeRaster(r,'filename.asc', format='ascii')

</code>
  </pre>
  <P>
    如果这些点不在常规网格中，您可以尝试该功能
     <code>
      rasterize
    </code>
    。我可以用该文件的头部做很多事情，但一个例子可能是：
  </p>
   <pre>
    <code>
      # Convert to SpatialPoinDataFrane
coordinates(predHB.Lv1Hard) <- ~x+y

# Calculate the extent of the observations
ext <- extent(predHB.Lv1Hard)

# Get the distance between points
table(diff(predHB.Lv1Hard@coords[,1]))
#0 0.00010000000000332 
#3                   2
# With the limited info I have I could assume that the resolution is 0.00010000000000332 degrees.
# If you know the original resolution, much better.

# Create the template
r <- raster(ext=ext,res=0.00010000000000332)
# Final raster
to_ascii <- rasterize(predHB.Lv1Hard, r, 'predLv1Hard')

# Write raster to ascii
writeRaster(to_ascii,'filename.asc', format='ascii')

</code>
  </pre>
</DIV>