New answers tagged

0

Hi Dave, this post is not an answer and must be deleted. None of the posts have been answers, but of all the posts you delete mine. Well done. Goodbye.


2

It looks like your problem is that the polygon data is multipart geometry. This means that you have multiple features (polygons) associated with single rows (attributes). Even if raster::extract works, this makes very little sense from a results standpoint. For your data to match, you need to explode your geometry into single part. Here is an example, you ...


0

I would like to suggest you to try mapnik node tile_server from this repo: https://github.com/pocketIlmatto/mapnik_node_tile_server


5

First off, a 100 byte shapefile indicates, that this file is just empty. The header of a shapefile is 100 byte. See https://www.esri.com/library/whitepapers/pdfs/shapefile.pdf on page 3. Secondly, the message "When you run QGIS and open a project file (.qgs), a new version of QGIS will read the .shp files and update them for format changes that go with the ...


1

Right-click the layer, click on "Filter..." and enter the expression LABEL % 10 != 5 into the text field labelled "Provider specific filter expression", then press OK.


6

Another way, using the same method as presented above by Joseph, with a slightly shorter code: if(LABEL % 10 = 5,"","LABEL") When the remainder of dividing by 10 is 5 (so 5, 15, 25...) it won't show the point.


3

We can use the modulus operator % to calculate the remainder values when dividing by 5. If there is no remainder then it is removed from the display. However, because multiples of 10 can be divided by 5 with no remainder, we will have to exclude this by adding some logic to the expression: if("LABEL" % 5 = 0 AND NOT "LABEL" % 10 = 0, "LABEL" % 5 AND NOT "...


0

Something might have happened to your data when converting to GeoJSON, because all the info you need is already in that dataset. I would try using QGIS to not only explore the data, but using it to export to GeoJSON and exploring the result of that export as well.


0

In general for cases like these refer to the GeoTools axis order FAQ - so you are correct that your axis order in the shapefile is lat/lon. Next, there are some issues with your code: CoordinateReferenceSystem geo = CRS.decode("EPSG:4326", true); creates a CRS with lon/lat if (and only if) the system hint FORCE_LONGITUDE_FIRST_AXIS_ORDER is set true, ...


5

GDAL supports two SQL dialects which can be used for all data sources. In addition to those two for some databases (Oracle, PostgeSQL etc.) the native SQL can be utilized. The OGR SQL dialect https://gdal.org/user/ogr_sql_dialect.html is made by the GDAL project and it is not at all a full featured SQL language. The other SQL dialect https://gdal.org/user/...


0

whyzar is right that Google Earth does not have a way to save to a shapefile, but can save to a KML which can then be converted to a shapefile with other software. Since it looks like you're using the new web version of Google Earth (not the Earth Pro desktop version), the UI is a bit different than described. You'll want to use the 3-dot menu at the top of ...


1

Edit: I looked around the Census TIGER repository and found these files which contain the census geographies and their polygons: https://www2.census.gov/geo/tiger/TIGER2018/ Old Answer: Try here: https://www2.census.gov/geo/tiger/TIGER_DP/2017ACS/ This is the entire ACS data, but the files do contain a listing of the census blocks as well as polygons. ...


1

Using the data from ?st_zm, here's a thing with Z and M coordinates: > a Simple feature collection with 2 features and 1 field geometry type: LINESTRING dimension: **XYZM** bbox: xmin: 1 ymin: 3 xmax: 8 ymax: 16 epsg (SRID): NA proj4string: NA a geom 1 1 LINESTRING ZM (1 9 17 25, 2... 2 2 LINESTRING ZM (1 ...


8

There is no export to shapefile via Google Earth. These steps when using the Google Earth web version: Once you are in the project section, you can click on the three vertical dots and then export as KML file as mentioned on this thread These steps when using Earth Pro desktop version: You will need to right-click the polygon layer, then click on Save As ...


3

While you country data seems to cover the whole planet (distinguishable are Cuba to the west, tilted at 90 °, Greenland overly large at the top middle, and something like Sri Lanka to the east, in the middle), your project, or more precisely, your projects coordinate reference system (CRS) covers only a small portion of the globe. While these "restricted" ...


1

Try re-creating the issue with your own polygon layer: Create a new polygon featureClass, Edit it in ArcMAP Draw all sorts of polygons : Hourglass, Islands, multi-polygon (user the cut/split editing functions). Draw polygons clockwise and anticlockwise. Save all edits Publish as a map service AND as a feature service Consume the map Service in a test ...


0

It seems like Overlay from Geopandas should work as well, through intersection http://geopandas.org/set_operations.html https://nbviewer.jupyter.org/github/geopandas/geopandas/blob/master/examples/overlays.ipynb But you need first: 1. Convert your netcdf into a dataframe, 2. convert latitude and longitud into a polygon like in this example, all the way to ...


1

In this case it looks like there's zero overlap between your habitats (red hashed) and the protected areas (orange solid). This map produced in QGIS: but in general once you've got both data sets as sf class objects and as polygons (there's some points in the protected area data) then you do st_intersection of the habitat with the protected area to get the ...


2

You can achieve this by setting rule based 3d rendering. Navigate to Properties > 3D View > Rule-based and create your own rules and styling. I did this really quick crude example:


0

You can achieve that using an SQL Query in the DB Manager of QGIS. You choose Database / Database Manager / Database Manager then Virtual Layers / Qgis Layers you can then try queries like the following : select zone, direction, count (*) as total from my_table group by zone having total>1 The result will show all zones having more than one direction (=...


1

To hide the labels outside feature - change to 'rule based labelling' and insert within($geometry, @atlas_geometry )=1 as filter


1

wfgeo's comment reminded me of a join I made between the shapefile and a text file/table I had added in QGIS. All of the column names in the joined shapefile that wouldn't import are: OBJECTID,GEOID,SUFFIX,NAME,ALAND,AWATER,INTPTLAT,INTPTLON,OR_pop_fie,OR_pop_sta,OR_pop_cou,OR_pop_tra,OR_pop_blo,OR_pop_P01,area_ha,popdens_ha,area_sqmi,popd_sqmi I undid the ...


0

You can download Road data without digitize but not in bulk you have to select which road you have to download. https://youtu.be/OsGBjHikAqs


-2

Use writeOGR() function to write the .shp file.


2

In d3 the correct way to do this is to use a composite projection. Here are plug-ins that specify composite projections for France: https://github.com/rveciana/d3-composite-projections (d3v4 and d3v3) https://github.com/seeschloss/francedom (unmaintained as of 20.12.2019)


5

Yes, you can physically move the polygons of the outer territories in geographic data, but this would degrade data unnecessarily. The much better solution hinted at by @ThingumaBob is to simply use multiple map frames when designing your export. Create a map frame for each outlying territory, then move those around on your sheet layout for printing/viewing ...


Top 50 recent answers are included