epsg.io | EPSG.io : Coordinate Systems Worldwide | Map library

I am trying to use kernel density smoothing to map the intensity of possible pest escape from vehicle traffic. Each route has been broken down into straight lines with each line having an integer attribute for the number of times the segment was travelled upon. However, when I use this attribute as the weight in kernel density smoothing, the weights don't seem to be used.

I've created a simplified reprex below with two abutting straight lines. Can anyone explain to me how I can make density.psp() account for the fact that one segment has an attribute 2x the magnitude of the other?

Many thanks for your help,

Josh

I am new to postgis, so please excuse what may be a very misguided question. I'm trying to reproject coordinates into a standard datum, but in some cases I don't think I'm getting a realistic result. The most obvious is from https://epsg.io/4813 to https://epsg.io/4326.

I have geotiff files load into xarray with a crs = EPSG:31467. I want to transform/reproject (don't know if there is a difference) these files into EPSG:4326. To do that, I use rasterio.warp.transform function which needs 1D arrays for x,y. To generate these i use numpy.meshgrid and flatten functions. Here is a small example with my data:

I have tried to follow the instructions here but I get wild results compared to this site.

Here is my code.

I am trying to calculate distance between two points using NetTopologySuite. Since I am referencing Microsoft documentation, I Came up with following GeometryExtension and GeometryHelper classes:

I want to check the area of geometry values.

The geometry values are POLYGON or POINT or MULTI POLYGON. The field has the type of geometry

I check the srid of the geom field:

I want to make grids (in the sense of data frames of x- and y-coordinates) over the US, or regions of the US, throwing out any points in the original grid that are beyond the borders of the US. I have some code that seems to work, but it's quite slow when I shrink the grid increment down to 1 km (1e3) or so. How can this be done faster? Perhaps there's a way to build the simple feature collection that I need without a lapply or loop, or perhaps this can be done with a MULTIPOINT instead of a simple feature collection of POINTs.

What I have:

• points in the arctic and antarctic
• raster data from various geophysical entities in arctic and antarctic

What I want:

A map in stereographic or any other polar projection with background map or coastlines, cropped to the extent of the points. In other words: A map like above with base map of my own choice.

What I did so far:

I loaded all the data (including land surface data from naturalearthdata; see MWE), projected them into stereographic and plotted that. The result including the polygon data looks then like this:

My MWE:

jsFiddle is given here: here

I am new to Open Layers 6 and I am trying to display Vector tile data on a map based, more or less, on the example given in the Open Layers workshop.

The URL for the vector tile source supplied in the above example code was not working so I am using the Vector Tile Source described in this page. There I read that the source is defined using using the RGF93 / Lambert-93 (EPSG:2154) coordinate system and then, using Google, I found that coordinate system's definition and bounds on this page.

In the code that follows I am using the projection's definition and the projected bounds from that last link.

Even though data do appear on the map, they appear only on the farthest left side of the screen and only at zoom level 2 as shown below:

If I change the zoom level, nothing is plotted on the screen.

The code is given below (see also link to JsFiddle above):