add ref to geosphere::regularCoordinates

Author: edzer

R/sample.R

@@ -35,6 +35,8 @@ st_sample = function(x, size, ...) UseMethod("st_sample")
 #' 
 #' Fibonacci sampling see: Alvaro Gonzalez, 2010. Measurement of Areas on a Sphere Using Fibonacci and Latitude-Longitude Lattices. 
 #' Mathematical Geosciences 42(1), p. 49-64
+#' 
+#' For regular sampling on the sphere, see also \code{geosphere::regularCoordinates}.
 #'
 #' Sampling methods from package \code{spatstat} are interfaced (see examples), and need their own parameters to be set. 
 #' For instance, to use \code{spatstat.random::rThomas()}, set \code{type = "Thomas"}.

man/st_sample.Rd

@@ -366,6 +366,7 @@ List CPL_read_gdal(CharacterVector fname, CharacterVector options, CharacterVect
 	CharacterVector descriptions(bands.size());
 	NumericMatrix ranges(bands.size(), 4);
 	IntegerMatrix blocksizes(bands.size(), 2);
+	IntegerVector colorInterp(bands.size());
 	for (int i = 0; i < bands.size(); i++) {
 		if ((poBand = poDataset->GetRasterBand(bands(i))) == NULL)
 			stop("trying to read a band that is not present");
@@ -391,6 +392,7 @@ List CPL_read_gdal(CharacterVector fname, CharacterVector options, CharacterVect
 		poBand->GetBlockSize(&nBlockXSize, &nBlockYSize);
 		blocksizes(i, 0) = nBlockXSize;
 		blocksizes(i, 1) = nBlockYSize;
+		colorInterp(i) = (int) poBand->GetColorInterpretation();
 	}
 
 	// get metadata items:
@@ -462,7 +464,8 @@ List CPL_read_gdal(CharacterVector fname, CharacterVector options, CharacterVect
 		_["blocksizes"] = blocksizes,
 		_["descriptions"] = descriptions,
 		_["default_geotransform"] = default_geotransform,
-		_["proxy"] = LogicalVector::create(!read_data)
+		_["proxy"] = LogicalVector::create(!read_data),
+		_["colorInterp"] = colorInterp
 	);
 	if (read_data) {
 		ReturnList.attr("data") = read_gdal_data(poDataset, nodatavalue, nXOff, nYOff,