<div dir="ltr"><br><br><div class="gmail_quote"><div dir="ltr">On Wed, 2 Dec 2015 at 06:34 Jonathan Greenberg <<a href="mailto:jgrn@illinois.edu">jgrn@illinois.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">GDAL Developers:<div><br></div><div>How do I go about requesting a new driver?  R's "raster" format is a pretty straightforward flat binary/header format (very similar to ENVI's).  The specifications are laid out in:</div><div><br></div><div><a href="https://cran.r-project.org/web/packages/raster/vignettes/rasterfile.pdf" target="_blank">https://cran.r-project.org/web/packages/raster/vignettes/rasterfile.pdf</a><br></div><div><br></div></div></blockquote><div><br></div><div><br></div><div><br></div><div><div>Hi Jonathan, FWIW, you can already do this via VRT. Can I suggest you flesh out examples using the hdr() pathway? </div><div><br></div><div><div>library(raster)</div><div>## example basic raster from matrix</div><div>r <- raster(volcano)</div><div>## two layers (bands), out to native grd/gri format</div><div>b <- writeRaster(brick(r, r/2), "b.grd")</div><div><br></div><div>## round trip to VRT on the .gri binary</div><div>## and read in with rgdal</div><div>v <- hdr(b, "VRT")</div><div>if (v) {</div><div><span class="Apple-tab-span" style="white-space:pre">       </span>library(rgdal)</div><div><span class="Apple-tab-span" style="white-space:pre">       </span>summary(readGDAL("b.vrt"))</div><div>}</div></div><div><div>b.vrt has GDAL driver VRT </div><div>and has 87 rows and 61 columns</div><div>Object of class SpatialGridDataFrame</div><div>Coordinates:</div><div>      min       max</div><div>x 0.0e+00 0.9999998</div><div>y 2.5e-07 1.0000000</div><div>Is projected: NA </div><div>proj4string : [NA]</div><div>Grid attributes:</div></div><div>...</div><div><br></div><div>I know this only complicates the issue in some ways - since it won't work for a lot of existing grd/gri files out there, but also does provide a workaround as well as a neat way to build test cases. </div><div><br></div><div>Cheers, Mike. </div></div><div><br></div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div><div>I think this would be a nice addition, to allow more fluid work between R and the various gdal utilities.</div><div><br></div><div>I'm cc'ing the developer of the format on this email, as an FYI.  Cheers!</div><div><br></div><div>--jonathan</div></div>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" rel="noreferrer" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a></blockquote></div></div>