tmaptools | Tools for thematic maps | Map library

Lidar data is simply 3d coordinates, usually in las file format. Сontent example

I am trying to get a raster layer (5-metre Digital Elevation Model (DEM)) from the Australian government through the ArcGIS REST API

url: https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/

I can perform a query for the boundary and footprint layers with good results

boundary query:

"https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/1/query?where=&text=&objectIds=&time=&geometry=xmin:133.226100657682, ymin:-23.6165145646468, xmax:133.235899342318, ymax:-23.6074854353532&geometryType=esriGeometryEnvelope&inSR=4283&spatialRel=esriSpatialRelEnvelopeIntersects&relationParam=&outFields=&returnGeometry=true&maxAllowableOffset=&geometryPrecision=&outSR=&returnIdsOnly=false&returnCountOnly=false&orderByFields=&groupByFieldsForStatistics=&outStatistics=&returnZ=false&returnM=false&gdbVersion=&returnDistinctValues=false&f=json"

response (plotted in R):

R plot through jsonlite response

Footprint query:

"https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/2/query?where=&text=&objectIds=&time=&geometry=xmin:133.226100657682, ymin:-23.6165145646468, xmax:133.235899342318, ymax:-23.6074854353532&geometryType=esriGeometryEnvelope&inSR=4283&spatialRel=esriSpatialRelEnvelopeIntersects&relationParam=&outFields=&returnGeometry=true&maxAllowableOffset=&geometryPrecision=&outSR=&returnIdsOnly=false&returnCountOnly=false&orderByFields=&groupByFieldsForStatistics=&outStatistics=&returnZ=false&returnM=false&gdbVersion=&returnDistinctValues=false&f=json"

response:

jsonlite response in R

where I got the name of TIF files I am looking for, nevertheless all my queryes to the raster layer: https://services.ga.gov.au/site_9/rest/services/DEM_LiDAR_5m/MapServer/3 does not give me any result

I have tried the rcran package "arcpullr" as well, but no results

any help?

I have a map with an elevation raster that is in lat/long coordinates but I now need to display the map in easting/northing. How can I accomplish this?

I have tried reprojecting the raster from lat/long to UTM but this warps the map (I assume for reasons discussed in this SO post).

A minimal working example follows a description of the code that I am using to produce the map in lat/long coordinates. I used a manually defined bounding box, extent, to download the elevation raster of interest from FedData. I then plotted the raster using tmap.

I would like to record a series of maps using the sf and tmap packages in R. My data is an sf object with different attributes linked to square polygons. The breaks for the maps are the same for each attribute that I want to map, that's why I thought of coding a simple loop.

The code below will give you a sample of my data. Sorry, the dput is a bit long, it's because sf objects integrate geometries.

I'm trying to work out how to create boundary boxes from sf point geometry data on every row in an sf object. I'm trying to use the 'bb' function from the tmaptools package, along with dplyr and rowwise(). However, the output I get is just the same bounday box values copied onto every row, rather than a specific boundary box calculated from the point data on each row.

Here is a snip of the data frame sf object:

I am trying to add some pie charts inside a map in the following way:

I am building a shiny app and I want to have a static dashboardHeader() title with a logo to the right of it. I have looked through stackoverflow to figure out how to do this and it seems like this requires HTML tags. I am not a computer programmer - just an R user, so I don't really understand how these work. But following advice of others, it seems like I need tag$li() to set the the height of the header bar and the logo. When I run the app, the header appears as I want in the RStudio shiny viewer, but when I view it in a browser (Chrome), the title is sunken below the logo and cut off. Below is a reproducible example. You will first need to download the [R logo] https://www.r-project.org/logo/Rlogo.png and save it as "Rlogo.png" in a subdirectory called "www" in the same directory as the code below (which in turn must be saved as "app.R"):

I am building a Shiny app that provides canned reports and maps from a database of field-collected vegetation data. The app has a series of SQL scripts that draw in the data from appropriate database tables to generate the reports and to query the spatial data stored as blobs. The user interface allows the user to select a project and then a unit within that project on which to report. When my app runs, Shiny freezes when I try to change one of the projects. I suspect that this has to do either with how I have set up the reactivity to update the "Unit" drop down menu, or the interface between R and MS SQL. Unfortunately, I was unable to recreate the problem outside of my organization's enterprise data, so I will present two scripts which will hopefully be enough information. My apologies ahead of time for not being able to come up with a script that reproduces the error. Below is a working minimal reproducible example of how I have structured the app - this one works without error, but will give you a sense of how I have built the app and how it should work:

I am trying plot county level maps. I came across a tutorial that uses the sf package and the tmap package. Using ?st_read, it shows an example of how to get county level data using North Carolina as an example, using the following code -

nc = st_read(system.file("shape/nc.shp", package="sf"))

I used the code above to plot county level data for NC, however I am not sure how to get the data for other states. I have tried replacing the nc.shp with sc.shp or va.shp, etc, however I get the following error - Error: Fill argument neither colors nor valid variable name(s).

How do I extract the data for other states, or preferable "all" states. I would like to be able to plot the map for the entire US, regionally (north, south, east, west) and also multiple states within a region.

Here is some reproducible code for the map I created -