Re: [R-sig-Geo] Leaflet map nested in RShiny App - Improving speed & portability

Apologies- looks like that only supports raster tiles, as does the R package 'tiler'. Sorry for the multiple messages on this. Cheers, Mike Please pardon any typos, this message was sent from a mobile device. On Wed, Sep 5, 2018, 2:04 PM Michael Treglia wrote: > I'll just second Barry's i

Re: [R-sig-Geo] Leaflet map nested in RShiny App - Improving speed & portability

As a quick follow up, I just came across this blog post on serving up map tiles via a GitHub repository, generating tiles via a QGIS plugin. https://khufkens.com/2018/09/19/github-tile-server/ Not an R solution, but perhaps something that might help in a case like this. As explained, the

Re: [R-sig-Geo] Leaflet map nested in RShiny App - Improving speed & portability

I'll just second Barry's idea in particular, to set up as a standalone webpage. You could even use QGIS and the QGIS2Web Plugin to create that, and host via GitHub pages or similar. >From R, after creating a map via leaflet and similar packages, you can use htmlwidgets::saveWidget() to export as

Re: [R-sig-Geo] raster::clump not working?

Check out SDMTools::ConnCompLabel I think that will get you what you're after. Hope that helps! Mike Please pardon any typos, this message was sent from a mobile device. On Fri, May 25, 2018, 7:15 AM João Carreiras wrote: > Dear Ben, > > Thank you for your prompt

Re: [R-sig-Geo] Reading tables without geometry from gdb?

D > 1 24.42419 22.76134 77.23866101 > 2 20.64160 14.38314 85.61686102 > > > > > On 02/11/2018 05:22 AM, Michael Treglia wrote: > > Hi All, > > > > I have a file geodatabase with non-spatial tables that can be joined to > > other spatial objects.

[R-sig-Geo] Reading tables without geometry from gdb?

Hi All, I have a file geodatabase with non-spatial tables that can be joined to other spatial objects. Is there a best/easiest way to import the non-spatial tables in R these days? I've seen some solutions here:

Re: [R-sig-Geo] raster and oceancolor L2 netcdf data

Hi All, Just saw this thread, and wanted to share some code I've got, in case it's useful, though I'll give the disclaimer I haven't looked carefully at the focal dataset here, though it sounds somewhat similar. Basically, I'm working with modeled data for NY Harbor, across multiple time points,

Re: [R-sig-Geo] Maintain SRID with st_write to postgis db

ill apply. */ >> > | #ifdef __cplusplus >> > | extern "C" >> > | #endif >> > | char lwgeom_version (); >> > | int >> > | main () >> > | { >> > | return lwgeom_version (); >>

Re: [R-sig-Geo] Maintain SRID with st_write to postgis db

or: Command failed (1) On Wed, Mar 15, 2017 at 6:14 PM, Michael Treglia <mtreg...@gmail.com> wrote: > Wow - thanks so much for this quick fix, Edzer! I like your solution to > having syntatically different but semantically identical proj4string. Will > try this in a bit, and write bac

Re: [R-sig-Geo] Maintain SRID with st_write to postgis db

t_0=40.16 +lon_0=-74 +x_0=30.01 +y_0=0 > +datum=NAD83 +units=us-ft +no_defs <> > > > The change I made to sf ( > https://github.com/edzer/sfr/commit/d620f3b748d487bae910e3884c554f > a9e3ce7090) > now first asks gdal to resolve the epsg (in case it is not

Re: [R-sig-Geo] Maintain SRID with st_write to postgis db

in PostGIS to settle the table > parameters. > > HTH, > > Chris > > > > > > On Wed, Mar 15, 2017 at 3:50 PM, Michael Treglia <mtreg...@gmail.com> > wrote: > >> Hi All, >> >> Been working to import and manipulate a CSV file wit

[R-sig-Geo] Maintain SRID with st_write to postgis db

Hi All, Been working to import and manipulate a CSV file with point data (lat/long), and then export to a PostGIS db. Overall, successful, but one thing I'd like to fix - when I write out the layer to postgis, the SRID is not maintained. The final EPSG/SRID should be 2263, but when I check in

Re: [R-sig-Geo] Suggestions for irregular grid netCDF data?

e deal with scads of this > kind of data and it doesn't fit any good mould except the "code it > yourself" club. It's not a huge amount of code but the "many coordinate > systems" can be bewildering and it's really specifically bound to the > intricacies of raster and s

[R-sig-Geo] Suggestions for irregular grid netCDF data?

​Hi all, I've been trying to work with a somewhat complex netCDF dataset, and wanted to see if anybody has suggestions. The data are for marine variables in the form of a remote netCDF served through OPENDAP. One of the complexities is that it comes in an irregular grid, as fine as 25m and as

Re: [R-sig-Geo] Specifying Schema for PostGIS Layer with Simple Features

IS geometries in R, and not in Simple Features > themselves, you may want to give a look at rpostgis, which present a > standard and flexible solution (using rgeos and RPostgreSQL in the > background): > > https://cran.r-project.org/package=rpostgis > > Mathieu. > > > On 12

Re: [R-sig-Geo] Specifying Schema for PostGIS Layer with Simple Features

Oh, cool Lee! That's great to have around! Thanks, Mike On Fri, Dec 9, 2016 at 7:47 PM, Lee Hachadoorian < lee.hachadooria...@gmail.com> wrote: > > > On Fri, Dec 9, 2016 at 6:06 PM, Michael Treglia <mtreg...@gmail.com> > wrote: > >> > > >> > >

Re: [R-sig-Geo] Specifying Schema for PostGIS Layer with Simple Features

NE 1: select * from 'staging.ipis' ^ ) Error in CPL_hex_to_raw(sapply(x, skip0x, USE.NAMES = FALSE)) : not compatible with STRSXP On Fri, Dec 9, 2016 at 4:04 PM, Edzer Pebesma <edzer.pebe...@uni-muenster.de > wrote: > > > On 09/12/16 21:02, Michael

Re: [R-sig-Geo] Specifying Schema for PostGIS Layer with Simple Features

; You can check if the execute permission is granted to your user id on the > schema you are trying to access? > > Thanks & Regards, > Avipsa Roy > > Michael Treglia wrote on 2016-12-09: > > Hi All, > > > Getting going with Simple Features - awesome job Dev

[R-sig-Geo] Specifying Schema for PostGIS Layer with Simple Features

Hi All, Getting going with Simple Features - awesome job Developers, and thanks for your time! Quick question - I'm working on bringing in PostGIS layers to work with in R - I've got the db connection made using the RPostgreSQL connection and have successfully imported data from the public

Re: [R-sig-Geo] look for advice on handling 3D/4D trajectories in R

Hi Paedro, Might want to look into the mkde package ( https://cran.rstudio.com/web/packages/mkde/index.html) - it's designed to use 3d tracking data to estimate movement in 3 dimensions. Hope that helps, Mike On Fri, Nov 4, 2016 at 1:09 PM, Pedro C. Alvarez wrote: > Hi all,

[R-sig-Geo] Surprising performance difference across computers with raster operation

Hi All, Was doing some work across a couple of different computers with the raster package, and noticed a pretty considerable difference in processing time to simply merge two Bioclim tiles (reproducible code below, plus system.time() and sessionInfo() details). Both machines have plenty of RAM

Re: [R-sig-Geo] LiDAR LAS file import into R via RSAGA

In case it's of any use to others, I've been scripting some LiDAR processing from R using LASTools, focused on dealing with many individual las files. At this point we've just been extracting a bunch of metadata fields to get a sense of what we have, fixing some metadata errors that we're aware

Re: [R-sig-Geo] LiDAR LAS file import into R via RSAGA

Hi Clay, At some point I was looking into something similar, and don't think you can do that with RSAGA. rsaga.geoprocessor() is running things as a system commands (like running from the command line directly), so it is not bringing the Lidar data into R, if I recall correctly. One package you

Re: [R-sig-Geo] How to assess orientation of a topographic feature?

Roland, Do you have these features already defined and such? One thing that might help in some of this is calculating the Multiresolution Index of Valley Bottom Flatness (MRVBF) - it identify valleys and ridges in a quantitative way. I'm not aware of an implementation of this directly in R, but

Re: [R-sig-Geo] A simple correlogram using sp.correlogram with data(meuse)

Hi all, This is a bit a tangential to the question, but related to pgirmess::correlog used in this thread, so I thought it appropriate to follow-up here [if not, apologies - I can start a new thread if that's preferred]. In recent exploration with the pgirmess::correlog, I was unable to change

Re: [R-sig-Geo] Creating a raster with an hexagonal grid

Manuel, I might be mistaken, but as far as I understand, rasters in R use square pixels... and thats how the raster data are typically stored in raster file types and used by other GIS packages, with cell size referring to the square dimensions of pixels. Thus, all of the ideas are based on that

Re: [R-sig-Geo] Problem with using multicore resample in Raster package

Haven't tried to troubleshoot the problem, but your code works fine for me. Looks like I'm running a slightly older versions of packages (notably, for raster). My session info is below. Here's my session info in case it helps: R version 3.1.2 (2014-10-31) Platform: x86_64-w64-mingw32/x64

Re: [R-sig-Geo] Quadrat Variance Analyses in R?

analysis. It accepts command line input so I wrote a shell script that can be called from and return results to R. -Original Message- From: R-sig-Geo [mailto:r-sig-geo-boun...@r-project.org] On Behalf Of Michael Treglia Sent: Monday, February 02, 2015 10:31 AM To: r-sig-geo@r

[R-sig-Geo] Quadrat Variance Analyses in R?

Hi All, I was wondering if anybody has suggestions for packages with-which to do quadrat-based contiguous units analysis (e.g., paired quadrat variance, two-term local quadrat variance, etc). Or would these analyses involve writing my own functions, probably based on quadrat functions in

Re: [R-sig-Geo] Extract with Large Rasters

causes it. Could you perhaps email me the points? Thanks, Robert On Wed, Dec 17, 2014 at 2:27 AM, Michael Sumner mdsum...@gmail.com wrote: Hello, see response inline. On Wed Dec 17 2014 at 10:03:23 Michael Treglia mtreg...@gmail.com wrote: Hi All, I'm working with some pretty

[R-sig-Geo] Extract with Large Rasters

Hi All, I'm working with some pretty huge rasters and trying to extract data to SpatialPointsDataFrames, and running into some issues. I'm using the raster package for this. (sessionInfo at the end of this e-mail) Here are the specs of my raster: class : RasterLayer dimensions : 59901,

Re: [R-sig-Geo] zonal statistics as table

I've had the same experiences as described above... Depending on exactly what stats you need to calculate, you might want to check out SAGA GIS ( http://saga-gis.org/). I've generally found it to work great on giant rasters, though you might hit a limit depending on your hardware. The tool you'd

Re: [R-sig-Geo] raster plotting question: figured it out

Hey Erin, Could you please include the solution? I haven't had to do that, but I'm sure that I (and others) will probably have to in the future, so it may be useful to have the solution in this thread Thanks! Mike On Fri, Jul 25, 2014 at 2:25 PM, Hodgess, Erin hodge...@uhd.edu wrote: Sorry for

Re: [R-sig-Geo] a KMZ question, please

Hi Erin, I also use 7zip for most of my .zip-ing needs and such. I don't typically deal with .kmz files, but found this batch script using 7zip that might be useful [DISCLAIMER: I have not tried it... just a suggestion that might help]:

Re: [R-sig-Geo] RasterBrick: Change resolution, keep the values

Hi Michael, Looks like you're using the Raster package. In that package, check out the resample() function ( http://www.inside-r.org/packages/cran/raster/docs/resample) - looks like that does what you want. Make sure both layers are in the same projection - otherwise you can use projectRaster()

Re: [R-sig-Geo] Environmental Risk Surface

, to specific locations, to wetlands, etc) to assess the risk of a disease spread that is not yet in the country but could enter through migratory birds. Best, Manuel 2014-03-19 19:08 GMT-06:00 Michael Treglia mtreg...@gmail.com: Hi Manuel, There may be multiple types of environmental risk

Re: [R-sig-Geo] Environmental Risk Surface

Hi Manuel, There may be multiple types of environmental risk surfaces, based on different variables, and using different methods to effectively create indices of risk. Can you provide some more detail about what you are trying to do? If there is a paper that explains the general methodology, that

Re: [R-sig-Geo] How I make 2 rasters with equal extents

Manuel, Given that your layers also have different resolutions, I think you need to use resample(). Since you need to set the resolutions the same to stack the rasters, give it some thought to make sure you are doing what your really want - which resolution is more appropriate and such, and which

Re: [R-sig-Geo] How I make 2 rasters with equal extents

Also, regarding changing resolutions, you might want to check out this vignette: http://cran.r-project.org/web/packages/raster/vignettes/Raster.pdf In particular, see section 5.1, starting on page 12. Cheers, mike On Mon, Oct 21, 2013 at 11:25 PM, Michael Treglia mtreg...@gmail.comwrote

Re: [R-sig-Geo] Plotting trend of a pixel

Hi Eddie, There are a couple of ways you can do it using the raster package. For example, create a point at the desired location, and use extract() from a raster to a SpatialPoints data type. Or, you could have a data.frame or matrix with x/y coordinates for which you want the points. Or possibly

Re: [R-sig-Geo] Extract function - too much data for R?

Thanks Gareth, Next time I do this I'll give the spsample + extract approach a whirl. Also, I've learned about the beginCluster()/endCluster() commands - not sure if extract takes advantage of those at all, but I'll give a check. Best regards, mike On Sun, Sep 8, 2013 at 2:41 AM, Gareth Davies

Re: [R-sig-Geo] Extract function - too much data for R?

Thanks for the suggestion! I ended up using the centroids of my polygons, and doing the overlay in SAGA GIS. The whole process ran quickly. GD- Ill keep your suggestion in mind for the future. A quick question though - your method doesn't directly allow computation of weighted average of pixel

[R-sig-Geo] Extract function - too much data for R?

Hi All, I've been trying to use the extract function (in the raster packages) to get values of a raster entered as attributes of series of polygons, and have tested it successfully with a reduced dataset, but am concerned that my full dataset might be too big. My full datasets involve ~46000 and

Re: [R-sig-Geo] Extract function - too much data for R?

be much faster for specific cases.. With more information, folks on the list might be able to suggest specific solutions. Tom On Fri, Aug 30, 2013 at 10:13 AM, Michael Treglia mtreg...@gmail.comwrote: Hi All, I've been trying to use the extract function (in the raster packages) to get

Re: [R-sig-Geo] Extract function - too much data for R?

. On Fri, Aug 30, 2013 at 12:36 PM, Michael Treglia mtreg...@gmail.comwrote: Hi Tom, List members, Thanks for the reply. My polygons constitute a tiling of 100m or 200m square polygons, with no overlaps. I'll be trying to make predictions for the smaller scale, but some of my input variables

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

, it actually housed gdal commands or something of the sort, whereas now it does not? Thanks for your help! Mike On Fri, Jul 12, 2013 at 1:00 AM, Roger Bivand roger.biv...@nhh.no wrote: On Fri, 12 Jul 2013, Michael Treglia wrote: Thanks for trying Erin! I appreciate knowing its not just me

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

=+proj=utm +zone=11 +datum=WGS84,OutRes=800) On Fri, Jul 12, 2013 at 11:39 AM, Michael Treglia mtreg...@gmail.comwrote: Hi Roger, In R 2.13.0, it is Raster version 1.9-70; in R 3.0.0 it is Raster version 2.1-48 When running in R 3.0.0, upon running the 'raster' command, it loads the rgdal

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

if it is unknown (as with esri ascii files), but not to change it if the file supplies one. Robert On Fri, Jul 12, 2013 at 2:55 PM, Michael Treglia mtreg...@gmail.com wrote: Hi All, Thank you again for the suggestions and assistance. I used the work-around provided by Roger and Erin, and I've

[R-sig-Geo] Problem Setting Projection of Rasters upon Import

Hi All, I am using he raster package, and in particular the raster function. I am dealing with a number of layers, and it would be easiest to set the CRS upon import. Looking through the documentation, it seems that it is doable, and I have been using the following code to get started, though it

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

? That might do it. Thanks, Erin From: r-sig-geo-boun...@r-project.org [r-sig-geo-boun...@r-project.org] on behalf of Michael Treglia [mtreg...@gmail.com] Sent: Thursday, July 11, 2013 7:16 PM To: r-sig-geo@r-project.org Subject: [R-sig-Geo] Problem Setting

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

again! Mike On Thu, Jul 11, 2013 at 5:59 PM, Michael Treglia mtreg...@gmail.com wrote: Hi Erin, Thanks for the quick response - I did try it with the plus sign, and same result. (I also tried it with 'crs=+proj=longlat +ellps=GRS80' with the same result too.) test-raster

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

can check with library(rgdal) CRS(+ellps=GRS80) Error in CRS(+ellps=GRS80) : projection not named On Fri, Jul 12, 2013 at 11:10 AM, Michael Treglia mtreg...@gmail.com wrote: I just tried running it in an older version of R, 2.13.0 (64-bit) and it ran with no problem... the exact same

Re: [R-sig-Geo] Problem Setting Projection of Rasters upon Import

, but no luck. Weird. Thanks Erin -- *From:* Michael Treglia [mtreg...@gmail.com] *Sent:* Thursday, July 11, 2013 8:21 PM *To:* Michael Sumner *Cc:* Hodgess, Erin; r-sig-geo@r-project.org *Subject:* Re: [R-sig-Geo] Problem Setting Projection of Rasters upon