https://community.wolfram.com RSS Feed for Wolfram Community showing questions tagged with Geographic Information System sorted by new. https://community.wolfram.com/groups/-/m/t/3645626 I used Locator in Manipulate[] to update mouse click point in GeoGraphics[]. The returned GeoPositionp[ from Locator has error in Latitude largely. See the following code. In this Poland map, Warsaw has {52.23,21.01} in latitude,longitude. Initialization of Locator to Warsaw did not work and there is 9.13 degree error in the latitude. It is huge error in location. Manipulate[ Column[ { "LocPos", LocPos, GeoGraphics[GeoRange -> {{48, 56}, {12, 27}}, GeoBackground -> "StreetMap", ImageSize -> Large] } ] , {{LocPos, GeoPosition[{52.23, 21.01}]}, Locator} ] Yeonjoon Park 2026-02-25T20:35:53Z https://community.wolfram.com/groups/-/m/t/3550646 In the release notes of WL 14.3, it is [said that][1] "An important feature of these maps is that they’re all produced with resolution-independent vector graphics." How can I export to svg or pdf a vector graphics of a geographic map of the Eiffel Tower, for example? GeoGraphics[Entity["Building", "EiffelTower::5h9w8"], GeoBackground -> "VectorClassic", GeoRange -> Quantity[400, "Meters"]] I right clicked on the output and clicked "Save Graphic as...", pdf format. I clicked the settings button and set "Use Highest Quality vector representation" as in the screenshot. Still, the pdf file is a raster graphics and not a vector graphics. Same when exporting to svg file. ![enter image description here][2] [1]: https://writings.stephenwolfram.com/2025/08/new-features-everywhere-launching-version-14-3-of-wolfram-language-mathematica/#maps-become-more-beautiful [2]: https://community.wolfram.com//c/portal/getImageAttachment?filename=ScreenShot2025-09-24at18.16.18.png&userId=1344988 Ehud Behar 2025-09-24T16:59:00Z https://community.wolfram.com/groups/-/m/t/3536075 I tried to show 3D terrain of United Arab Emirate (UAE) with GeoElevationGraphics3D function. I found there is significant shift error in Digital Terrain Elevation Data in GeoElevationGraphics3D [] function. The 3D mountain is drawn far to the right and does not match 2D "Reliefmap: or 2D "Streetmap". The 3D mountain is located on the ocean water to the east. The 3D mountain has to be on the peninsular. How can we manually fix the shift error of 3D Digital Terrain Elevation Data? Can we access the elevation data matrix? You can see the attached images. (* Default Map Range*) lattiLow = 23; lattiUp = 29; longiLow = 50.5; longiUp = 61; ResourceFunction["GeoElevationGraphics3D"][ Polygon[Entity["Country", "UnitedArabEmirates"]], GeoBackground -> GeoStyling["ReliefMap"], GeoRange -> {{lattiLow, lattiUp}, {longiLow, longiUp}}, Axes -> True, AxesLabel -> {"Longitude", "Latitude", "2.4xElevation"}, ImageSize -> {800, 800}, Ticks -> Automatic] Yeonjoon Park 2025-08-26T17:52:09Z https://community.wolfram.com/groups/-/m/t/3532941 In 14.3 at least, when I make a large map with any projection other than the default Equirectangular, the country borders that are overlaid by default are rough and fuzzy. Anyone else? E.g.: GeoGraphics["World", GeoProjection -> "Mollweide", ImageSize -> 1600] ![enter image description here][1] [1]: https://community.wolfram.com//c/portal/getImageAttachment?filename=Map.png&userId=1385683 Gareth Russell 2025-08-19T12:54:28Z https://community.wolfram.com/groups/-/m/t/3509967 I have a large number of locations (fiber end points from a large database) that I would like to calculate the distance between. In general, fibers tend to follow rail and highway right of ways so in the absence of actual fiber locations, I'm using the driving distance as the fiber length between two given points and the travel distance as the fiber length. For fibers where I have actual distances, the results are within 10% which is good enough for my application. I'm having trouble with TravelDirections/TravelDistance only intermittently working. Even using the documentation example of: TravelDistance[{Entity["City", {"Miami", "Florida", "UnitedStates"}], Entity["City", {"NewYork", "NewYork", "UnitedStates"}]}] Often returns $Failed without a reason. 2 minutes later, it will work for a while and then it will stop working again. The problem is that sometimes there is no good driving direction between the fiber endpoints even when TravelDirections IS working so there is no way for me to differentiate between the service not working and bad input points. My current approach is to first check if the example works and then if it does, try it on the fiber end points and hope that the service stays working for at least that long. Is anyone else experiencing this problem or have a clever solution? Thanks. Troy Rockwood 2025-07-12T18:58:19Z https://community.wolfram.com/groups/-/m/t/3438066 Given the following code: ```wls Import["https://www.gloriouslyawkward.com/wp-content/uploads/2025/03/2025-03-26-225828-Outdoor-Walking-Stevens-Apple-Watch.gpx", GeoBackground -> "SatelliteWithLabels"] ``` How do I make the line a different color instead of the default black? Steven Buehler 2025-04-02T16:52:38Z https://community.wolfram.com/groups/-/m/t/3413025 In the documentation of SpatialTrendFunction, I saw that variables can be introduced to implement the universal kriging method! <br />But many attempts have failed! <br />Can anyone with experience give me some guidance? &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/1f6d8218-4e1b-45cb-b4de-18b8801c8dbf Tsai Ming-Chou 2025-03-09T02:46:27Z https://community.wolfram.com/groups/-/m/t/3377021 Consider travel times between three cities in Romania : TravelTime[Entity["City", {"AlbaIulia", "Alba", "Romania"}], Entity["City", {"Harman", "Brasov", "Romania"}]] TravelTime[Entity["City", {"AlbaIulia", "Alba", "Romania"}], Entity["City", {"Prejmer", "Brasov", "Romania"}]] Three hours between "AlbaIulia" and "Harman", which seems correct, and a bit more between "AlbaIulia" and "Prejmer". But, for further computation, I need to convert these durations in minutes. In[200]:= Quantity[ TravelTime[Entity["City", {"AlbaIulia", "Alba", "Romania"}], Entity["City", {"Prejmer", "Brasov", "Romania"}]], "Minute"] Quantity[ TravelTime[Entity["City", {"AlbaIulia", "Alba", "Romania"}], Entity["City", {"Harman", "Brasov", "Romania"}]], "Minute"] Out[200]= Quantity[188, ("Minutes")^2] Out[201]= Quantity[1/20, ("Hours")^2] There are two problems. First, why square minutes? In addition, while the first "distance" is correct (188 minutes), the second one is a nonsense : 1/20 and the unit is square hour instead of minute! Is there a bug? Best regards, Claude Claude Mante 2025-02-03T17:48:26Z https://community.wolfram.com/groups/-/m/t/3293403 Does anyone have any kind of tutorial on how to work with GRIB/GRIB2 files (extracting and mapping the data with `GeoGraphics`), for example, https://tgftp.nws.noaa.gov/SL.us008001/ST.opnl/DF.gr2/DC.ndfd/AR.oceanic/VP.004-007/ds.wspd.bin ? The online documentation is insanely vague. Steven Buehler 2024-10-08T00:49:12Z https://community.wolfram.com/groups/-/m/t/3274221 How to specify that I want to use IGRF model instead of automatic when calling GeomagneticModelData ? The help is not clear and it is missing example. I tried various options like following but none of them is correct. geomagneticNorth = GeomagneticModelData["NorthGeomagneticPole", dates, Method->"IGRF"]; geomagneticNorth = GeomagneticModelData["NorthGeomagneticPole", dates, Method->IGRF]; Ján Kolár 2024-09-14T12:46:14Z https://community.wolfram.com/groups/-/m/t/3249891 me again. ;-) I having an issue creating a web form that I want to use for people to submit messages with a location whenever they find an item that I've left behind. I have the following code for the form, which you can try: ``` CloudDeploy@FormFunction[ { "Your Location" -> <|"Required" -> True, "Interpreter" -> "Location", "Input" :> Here["LatitudeLongitude"]|>, "Your Name" -> <| "Required" -> True, "Interpreter" -> "String"|>, "Phone Number" -> <|"Required" -> False, "Interpreter" -> "PhoneNumber"|>, "Email" -> <|"Required" -> True, "Interpreter" -> "EmailAddress"|>, "Message" -> <|"Required" -> False, "Interpreter" -> "String"|> } ``` The "Your Location" field automatically is filled in with the device location as `{lat,lng}` (with values for the `lat` and `lng`, of course). However, when submitting the form, it returns the error message, **"No geographic coordinates** [or **location**] **interpretation found. Try again."** on that field. What am I doing wrong? I've reviewed the documentation [here](https://reference.wolfram.com/language/workflow/ProvideInitialDataToAForm.html), but it's not working with the particular function I'm using for the prefill, which is `Here["LatitudeLongitude"]`. I have also tried `"Location"`, `"ComputedLocation"`, `"GeoCoordinates"`, and `"GeoPosition"` as values for `Interpreter`. Steven Buehler 2024-08-16T21:54:21Z https://community.wolfram.com/groups/-/m/t/3241392 I am brand new to the Wolfram ecosystem (as of today!) and I need to start with some basic research on selected towns or cities that are county seats in the contiguous 48 states. My selection criteria are: (1) Current county seat (2) population between 4,000 and 75,000 (3) not the state capital (4) not in Alaska or Hawaii or any US territories. I want to list them by region and state. I'd appreciate any steering you might offer to a total neophyte in this ecosystem. Philip Ullom 2024-08-05T17:24:19Z https://community.wolfram.com/groups/-/m/t/3237976 I am currently studying Professor André Dauphiné's book "Geographical Models with Mathematica." This book is full of concrete examples, and I find it very interesting and meaningful, as it allows me to combine the practice and theory of geography using Mathematica. However, some of the data is not included. Could someone please tell me the following data used in the examples in this book? (1) In the following 15 places, the syntax "data = Import[SystemDialogInput["FileOpen"]]" is used to read data from an external file, but since the external file is not attached, Mathematica cannot analyze it. p.38, 2.2.5.2. , Program 2.6. p.71, 2.5. , Program 2.16 p.170, 7.1.4.1. , Program 7.4. p.180, 7.2.1.1. , Program 7.9. p.184, 7.2.2.2. , Program 7.11. p.193, 1.2. , Program 7.14 p.193, 7.2.3.2. , Program 7.15. p.205, 8.2.1. , Program 8.1 p.209, 8.2.1. , Program 8.1. p.227, 9.1.1.2. , Program 9.2. p.233, 9.1.2.3. , Program 9.6. p.240, 9.3.1. , Program 9.7. p.243, 9.3.2. , Program 9.8. image p.247, 9.3.4. , Program 9.9. p.250, 10.1.1. , Program 10.1. (2) In p.34, 2.2.3.4, it says "line= LinearModelFit[d3,{1, v, x, y, z},{v, x, y, z}]", but Mathematica cannot perform the analysis because the data for the variables v, y, and z is not displayed. (3) There is currently no corresponding Wolfram Demonstrations Project for p.131 5.1., "Figure 5.1. The von Thünen–A A Alonso model" and p.132 "Figure 5.2. Formation of a polycentrric conurbation". (4) p.279, 10.5., Program 10.8., no image data is entered into the variable imageca in the first line, so analysis by Mathematica cannot be performed. (5) This may not be related to the data, but when I run p.260, 10.2.2.44. A realistic example: thee diffusion off an innovation in France, Program 10.44. Mesh of the map of France, I get the error "ToElementMesh::femimq: Element mesh quality 0.` is insufficient. If the quality estimate is less than 0., this may be due to the wrong order of element instances or self-intersections." and the output is "ElementMesh[{{-4.79537, 9.56036}, {45.5126, 59.623}}, {TriangleElement["<" 12032 ">"]}]", and I cannot proceed to the next Program 10.5. How should I deal with this? Any advice would be appreciated. kim hiro 2024-07-31T08:43:49Z https://community.wolfram.com/groups/-/m/t/3213254 The SpatialEstimatorFunction created using the dataset of 3D points is unable to estimate the 3D point. This is strange! &[Wolfram Notebook][1] [1]: https://www.wolframcloud.com/obj/ebbceb56-2725-44fd-a907-5379a70b0fd9 Tsai Ming-Chou 2024-07-11T07:23:44Z https://community.wolfram.com/groups/-/m/t/3185572 How do I get this: FlightData[Entity["Airport", "KMCO"] -> All, "FlightPath", Today] onto a _SINGLE_ `GeoGraphics` map? Steven Buehler 2024-05-29T20:24:43Z https://community.wolfram.com/groups/-/m/t/3101562 Hello, I have a Dataset with around 42k observation points from various meteorological stations. I run a Query[All,DateObject[{#year,#month,#day}]]@ds and it's around 0.25 seconds. Now I run Query[All,GeoPosition[{#latitude,#longitude}]]@ds and it's a stunning 20 seconds! What's wrong with that? Is it that complex to create a geoposition object? I don't see this mentioned in the docs. Best, Przemek Przemyslaw K. 2024-01-12T09:21:48Z https://community.wolfram.com/groups/-/m/t/3100673 I would like to use Laplacian, Div, and Grad in geophysical spherical coordinates rather than the usual spherical coordinates. The difference is that standard spherical coordinates uses the polar angle, θp, the angle from the North Pole with southwards being positive, while in geophysics it is standard to use latitude, θ, the angle from the Equator with northwards being positive. The transformation is relatively simple, θp = π/2 - θ. Does Mathematica have this as a defined coordinate system? If not, can I define a new coordinate system that will work with Laplacian, Div, and Grad? Jeffrey Weiss 2024-01-11T19:23:33Z https://community.wolfram.com/groups/-/m/t/3065876 If we have the following code: x = {-63., 45}; GeoGraphics[Locator[x]] This shows the locator in the correct place. But if we add a dynamic, then the locator ends up in the wrong place, and the graph boundaries are huge. GeoGraphics[Locator[Dynamic[x]]] Any suggestions? Thanks Vilis Nams 2023-11-16T13:08:30Z https://community.wolfram.com/groups/-/m/t/3063181 I have two points with a name and a latitude and longitude that I would like to plot on a map. The map would show the location as a red dot tagged by the name. &[Wolfram Notebook][1] bw = GeoPosition[{51.7091931533413, -1.25457498392513}]; lba = GeoPosition[{51.4976561183666, -0.127636445324134}]; When I use GeoListPlot I get a map with the points and tag, but it doesn’t have geographic information. points = GeoListPlot[{ Labeled[Point[bw], "BW", Left], Labeled[Point[lba], "LBA", Left] }] When I use GeoGraphics I get a map with the points and geographic information, but I can’t get the tags. When I put used Labeled I get the message, Labeled is not a Graphics primitive or directive. I’ve tried other approaches, but none work. GeoGraphics[{Red, PointSize[Large], Point[bw], Point[lba] }] How can I get a map with the points, tags, and geographic information? [1]: https://www.wolframcloud.com/obj/0e546223-8217-4b82-8a8c-02ce31d872a7 Dan O'Leary 2023-11-10T15:21:02Z https://community.wolfram.com/groups/-/m/t/3019459 Hello all dear, I have the satellite image of Canton dam reservoir and I want to get its volume and area. The actual specifications of the dam reservoir: volume=0.472km^3 & Surface area=32.01km^2 I did the following command to estimate the area, but the result was not acceptable and the unit is not clear. Please guide me to estimate the area and volume. Thank you ![enter image description here][1] First, I manually obtained the shape below from the Mask Tool command. ![enter image description here][2] imge =(image); img = Binarize@imge~Erosion~1; (m = MorphologicalComponents[img]) // Colorize ![enter image description here][3] linecount = 3 /. (Binarize@imge // ColorNegate // Thinning // ComponentMeasurements[#, "PerimeterCount"] &); areacount = 2 /. ComponentMeasurements[MorphologicalComponents[img], "Count"]; areacount/linecount^2 // N out put is : 0.111111 [1]: https://community.wolfram.com//c/portal/getImageAttachment?filename=cantondam.jpg&userId=2971053 [2]: https://community.wolfram.com//c/portal/getImageAttachment?filename=canton.jpg&userId=2971053 [3]: https://community.wolfram.com//c/portal/getImageAttachment?filename=fgg.jpg&userId=2971053 Erfan Abdi 2023-09-22T19:21:16Z